scholarly journals Quasi-Relativistic Calculation of EPR g-Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance

Author(s):  
Yannick J. Franzke ◽  
Jason M. Yu

We present an exact two-component (X2C) ansatz for the EPR g-tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and "fully" relativistic ansätze for the g-tensor, this implementation results in a gauge-origin invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are considered to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU) and the (multipole-accelerated) resolution of the identity approximation. The X2C ansatz is compared to Douglas-Kroll-Hess theory and the zeroth-order regular approximation for 11 diatomic molecules. The impact of the relativistic Hamiltonian, the basis set, and the density functional approximation is subsequently assessed for a set of 17 transition-metal complexes to complement our previous work on the hyperfine coupling constant [DOI: 10.33774/chemrxiv-2021-wnz1v-v2]. In total, 24 basis sets and 22 density functional approximations are considered. The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent "fully" relativistic four-component theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. Generally, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations are needed to match the experimental findings. Here, hybrids based on the meta -generalized gradient approximation are not an a priori improvement. Compared to the other computational parameters, the impact of the GIAOs and the magnetic balance on the actual results in standard calculations is less pronounced. Routine calculations of large molecules are possible with widely available and comparably low- cost hardware as demonstrated for [Pt(C6Cl5)4]− with 3360 basis functions and three spin-(1/2) La(II) and Lu(II) compounds. Both approaches based on a common gauge origin and GIAOs using triple- ζ basis sets lead to a good agreement with the experimental findings. The best agreement is found with hybrid functionals such as PBE0 and ωB97X-D.

2010 ◽  
Vol 7 (2) ◽  
pp. 449-455
Author(s):  
S. D. S. Chauhan ◽  
A.K. Sharma ◽  
R. Kumar ◽  
D. Kulshreshtha ◽  
R. Gupta ◽  
...  

Vibrational frequencies of aniline in gas phase have been calculated and each of their modes of vibration assigned properly at RHF and DFT with 6-31G(d) basis set. In the present study, it has been observed that the 6-31G(d) basis set at both RHF and DFT levels of calculations provides better agreement to the experimental findings as compared to other basis sets. Simultaneously, Density functional theory is found to be superior to its counterpart Hartree Fock method.


2019 ◽  
Author(s):  
Kamal Batra ◽  
Stefan Zahn ◽  
Thomas Heine

<p>We thoroughly benchmark time-dependent density- functional theory for the predictive calculation of UV/Vis spectra of porphyrin derivatives. With the aim to provide an approach that is computationally feasible for large-scale applications such as biological systems or molecular framework materials, albeit performing with high accuracy for the Q-bands, we compare the results given by various computational protocols, including basis sets, density-functionals (including gradient corrected local functionals, hybrids, double hybrids and range-separated functionals), and various variants of time-dependent density-functional theory, including the simplified Tamm-Dancoff approximation. An excellent choice for these calculations is the range-separated functional CAM-B3LYP in combination with the simplified Tamm-Dancoff approximation and a basis set of double-ζ quality def2-SVP (mean absolute error [MAE] of ~0.05 eV). This is not surpassed by more expensive approaches, not even by double hybrid functionals, and solely systematic excitation energy scaling slightly improves the results (MAE ~0.04 eV). </p>


2007 ◽  
Vol 62 (12) ◽  
pp. 711-715 ◽  
Author(s):  
Ahmad Seif ◽  
Mahmoud Mirzaei ◽  
Mehran Aghaie ◽  
Asadollah Boshra

Density functional theory (DFT) calculations were performed to calculate the electric field gradient (EFG) tensors at the sites of aliminium (27Al) and nitrogen (14N) nuclei in an 1 nm of length (6,0) single-walled aliminium nitride nanotube (AlNNT) in three forms of the tubes, i. e. hydrogencapped, aliminium-terminated and nitrogen-terminated as representatives of zigzag AlNNTs. At first, each form was optimized at the level of the Becke3,Lee-Yang-Parr (B3LYP) method, 6-311G∗∗ basis set. After, the EFG tensors were calculated at the level of the B3LYP method, 6-311++G∗∗ and individual gauge for localized orbitals (IGLO-II and IGLO-III) types of basis sets in each of the three optimized forms and were converted to experimentally measurable nuclear quadrupole resonance (NQR) parameters, i. e. quadrupole coupling constant (qcc) and asymmetry parameter (ηQ). The evaluated NQR parameters revealed that the considered model of AlNNT can be divided into four equivalent layers with similar electrostatic properties.With the exception of Al-1, all of the three other Al layers have almost the same properties, however, N layers show significant differences in the magnitudes of the NQR parameters in the length of the nanotube. Furthermore, the evaluated NQR parameters of Al-1 in the Al-terminated form and N-1 in the N-terminated form revealed the different roles of Al (base agent) and of N (acid agent) in AlNNT. All the calculations were carried out using the GAUSSIAN 98 package program.


Author(s):  
Bole Chen ◽  
Gennady L. Gutsev ◽  
Weiguo Sun ◽  
Xiao-Yu Kuang ◽  
Cheng Lu ◽  
...  

The coalescence of two Fe8N as well as the structure of the Fe16N2 cluster were studied using density functional theory with the generalized gradient approximation and a basis set of...


2021 ◽  
Author(s):  
Marco Caricato

<div> <div> <div> <p>We present an origin-invariant approach to compute the full optical rotation tensor (Buckingham/Dunn tensor) in the length dipole gauge without recourse to London atomic orbitals, called LG(OI). The LG(OI) approach is simpler and less computationally demanding than the more common LG-London and modified velocity gauge (MVG) approaches and it can be used with any approximate wave function or density functional method. We report an implementation at coupled cluster with single and double excitations level (CCSD), for which we present the first simulations of the origin-invariant Buckingham/Dunn tensor in the length gauge. With this method, we attempt to decouple the effects of electron correlation and basis set incompleteness on the choice of gauge for optical rotation calculations on simple test systems. The simulations show a smooth convergence of the LG(OI) and MVG results with the basis set size towards the complete basis set limit. However, these preliminary results indicate that CCSD may not be close to a complete description of the electron correlation effects on this property even for small molecules, and that basis set incompleteness may be a less important cause of discrepancy between choices of gauge than electron correlation incompleteness. </p> </div> </div> </div>


Author(s):  
Rabiu Nuhu Muhammad ◽  
N. M. Mahraz ◽  
A. S Gidado ◽  
A. Musa

Tetrathiafulvalene () is an organosulfur compound used in the production of molecular devices such as switches, sensors, nonlinear optical devices and rectifiers. In this work, a theoretical study on the effects of solvent on TTF molecule was investigated and reported based on Density Functional Theory (DFT) as implemented in Gaussian 03 package using B3LYP/6-31++G(d,p) basis set. Different solvents were introduced as a bridge to investigate their effects on the electronic structure. The HUMO, LUMO, energy gap, global chemical index, thermodynamic properties, NLO and DOS analysis of the TTF molecule in order to determine the reactivity and stability of the molecule were obtained. The results obtained showed that the solvents have effects on the electronic and non-linear-optical properties of the molecule. The optimized bond length revealed that the molecule has strong bond in gas phase with smallest bond length of about 1.0834Å than in the rest of the solvents. It was observed that the molecule is more stable in acetonitrile with HOMO-LUMO gap and chemical hardness of 3.6373eV and 1.8187eV respectively. This indicates that the energy gap and chemical hardness of TTF molecule increases with the increase in polarity and dielectric constant of the solvents. The computed results agreed with the results in the literature. The thermodynamics and NLO properties calculation also indicated that TTF molecule has highest value of specific heat capacity (Cv), total dipole moment () and first order hyperpolarizability () in acetonitrile, while acetone has the highest value of entropy and toluene has a slightly higher value of zero point vibrational energy (ZPVE) than the rest of the solvents. The results show that careful selection of the solvents and basis sets can tune the frontier molecular orbital energy gap of the molecule and can be used for molecular device applications.


2021 ◽  
Vol 2090 (1) ◽  
pp. 012078
Author(s):  
A Paolone ◽  
S Brutti

Abstract In this paper we investigated the calculation of the anodic limit of two anions of ionic liquids, largely used as electrolyte of lithium batteries. Starting from a model based on calculations performed on single ions at the MP2 level of theory, we showed that the matching between calculation and experiments decreases while using more expanded basis set with respect to 6-31G**, possibly because of the destabilization of the neutral species when larger basis sets are considered. Additionally, in order to decrease the computational time, the performances for the calculation of the anodic limit obtained by means of a series of DFT functionals with increasing level of complexity (from the Generalized Gradient Approximation to the Range Separated Hybrid meta-Generalized Gradient Approximation) were compared. Overall, the best performing functionals are BMK, ωB97M-V and MN12-SX, while acceptable results can be obtained by M06-2X, M11, M08-HX and M11-L. Some less computationally expensive functionals, like CAM-B3LYP and ωB97X-D, also provide reasonable values of the anodic limit.


Author(s):  
Qin Yang ◽  
Marco Mendolicchio ◽  
Vincenzo Barone ◽  
Julien Bloino

Vibrational spectroscopy represents an active frontier for the identification and characterization of molecular species in the context of astrochemistry and astrobiology. As new missions will provide more data over broader ranges and at higher resolution, especially in the infrared region, which could be complemented with new spectrometers in the future, support from laboratory experiments and theory is crucial. In particular, computational spectroscopy is playing an increasing role in deepening our understanding of the origin and nature of the observed bands in extreme conditions characterizing the interstellar medium or some planetary atmospheres, not easily reproducible on Earth. In this connection, the best compromise between reliability, feasibility and ease of interpretation is still a matter of concern due to the interplay of several factors in determining the final spectral outcome, with larger molecular systems and non-covalent complexes further exacerbating the dichotomy between accuracy and computational cost. In this context, second-order vibrational perturbation theory (VPT2) together with density functional theory (DFT) has become particularly appealing. The well-known problem of the reliability of exchange-correlation functionals, coupled with the treatment of resonances in VPT2, represents a challenge for the determination of standardized or “black-box” protocols, despite successful examples in the literature. With the aim of getting a clear picture of the achievable accuracy and reliability of DFT-based VPT2 calculations, a multi-step study will be carried out here. Beyond the definition of the functional, the impact of the basis set and the influence of the resonance treatment in VPT2 will be analyzed. For a better understanding of the computational aspects and the results, a short summary of vibrational perturbation theory and the overall treatment of resonances for both energies and intensities will be given. The first part of the benchmark will focus on small molecules, for which very accurate experimental and theoretical data are available, to investigate electronic structure calculation methods. Beyond the reliability of energies, widely used for such systems, the issue of intensities will also be investigated in detail. The best performing electronic structure methods will then be used to treat larger molecular systems, with more complex topologies and resonance patterns.


2014 ◽  
Vol 13 (04) ◽  
pp. 1450023 ◽  
Author(s):  
Reza Ghiasi ◽  
Morteza Zaman Fashami ◽  
Amir Hossein Hakimioun

In this work, the interaction of C 20 with N 2 X 2 ( X = H , F , Cl , Br , Me ) molecules has been explored using the B3LYP, M062x methods and 6-311G(d,p) and 6-311+G(d,p) basis sets. The interaction energies (IEs) obtained with standard method were corrected by basis set superposition error (BSSE) during the geometry optimization for all molecules at the same levels of theory. It was found C 20… N 2 H 2 interaction is stronger than the interaction of other N 2 X 2 ( X = F , Cl , Br , Me ) with C 20. Highest occupied and lowest unoccupied molecular orbitals (HOMO and LUMO, respectively) levels are illustrated by density of states spectra (DOS). The nucleus-independent chemical shifts (NICSs) confirm that C 20… N 2 X 2 molecules exhibit aromatic characteristics. Geometries obtained from DFT calculations were used to perform NBO analysis. Also, 14 N NQR parameters of the C 20… N 2 X 2 molecules are predicted.


Molecules ◽  
2020 ◽  
Vol 25 (21) ◽  
pp. 5174
Author(s):  
Frederick Stein ◽  
Jürg Hutter ◽  
Vladimir V. Rybkin

Intermolecular interactions play an important role for the understanding of catalysis, biochemistry and pharmacy. Double-hybrid density functionals (DHDFs) combine the proper treatment of short-range interactions of common density functionals with the correct description of long-range interactions of wave-function correlation methods. Up to now, there are only a few benchmark studies available examining the performance of DHDFs in condensed phase. We studied the performance of a small but diverse selection of DHDFs implemented within Gaussian and plane waves formalism on cohesive energies of four representative dispersion interaction dominated crystal structures. We found that the PWRB95 and ωB97X-2 functionals provide an excellent description of long-ranged interactions in solids. In addition, we identified numerical issues due to the extreme grid dependence of the underlying density functional for PWRB95. The basis set superposition error (BSSE) and convergence with respect to the super cell size are discussed for two different large basis sets.


Sign in / Sign up

Export Citation Format

Share Document