New implementation of molecular double point‐group symmetry in four‐component relativistic Gaussian‐type spinors

A new, practical implementation of double‐group symmetry to relativistic Gaussian spinors is presented for four‐component relativistic molecular calculations. We show that the systematic adaptability to irreducible representations under arbitrary point‐group symmetry, as well as Kramers (time‐reversal) symmetry, is inherent in the present basis spinors, which possess the analytic structure of Dirac atomic spinors. The implementation of double‐group symmetry entails significant computational efficiencies in the relativistic second‐order Møller–Plesset perturbation calculation on Au₂ and the density functional theory (DFT) calculation with the B3LYP functional on octahedral UF₆, in which the highest symmetries used are, respectively, C and D. The four‐component B3LYP equilibrium geometry of UF₆ is reported. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Nonadiabatic relaxation dynamics of water anion cluster and its isotope effects by ring‐polymer molecular dynamics simulation

We have applied a recently developed hybrid quantum ring‐polymer molecular dynamics method to the nonadiabatic p → s relaxation dynamics in water anion clusters to understand the isotope effects observed in previous experiments. The average relaxation times for (H₂O)₅₀^? and (D₂O)₅₀^? were calculated at 120 and 207 fs, respectively, and are comparable to the experimental results. Therefore, we conclude that nuclear quantum effects play an essential role in understanding the observed isotope effects for water anion cluster nonadiabatic dynamics. The nonadiabatic relaxation mechanisms are also discussed in detail. © 2014 Wiley Periodicals, Inc.     

Excited states of OsO4: A comprehensive time‐dependent relativistic density functional theory study

A large number of scalar as well as spinor excited states of OsO₄, in the experimentally accessible energy range of 3–11 eV, have been captured by time‐dependent relativistic density functional linear response theory based on an exact two‐component Hamiltonian resulting from the symmetrized elimination of the small component. The results are grossly in good agreement with those by the singles and doubles coupled‐cluster linear response theory in conjunction with relativistic effective core potentials. The simulated‐excitation spectrum is also in line with the available experiment. Furthermore, combined with detailed analysis of the excited states, the nature of the observed optical transitions is clearly elucidated. It is found that a few scalar states of ³T₁ and ³T₂ symmetries are split significantly by the spin‐orbit coupling. The possible source for the substantial spin‐orbit splittings of ligand molecular orbitals is carefully examined, leading to a new interpretation on the primary valence photoelectron ionization spectrum of OsO₄. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Four‐component relativistic DFT calculations of 77Se NMR chemical shifts: A gateway to a reliable computational scheme for the medium‐sized organoselenium molecules

A versatile high‐accuracy computational scheme for the ⁷⁷Se nuclear magnetic resonance (NMR) chemical shifts of the medium‐sized organoselenium compounds is suggested within a framework of a full four‐component relativistic density functional theory (DFT). The main accuracy factors (DFT functionals, relativistic geometry, vibrational corrections, and solvent effects) are addressed. The best result is achieved with NMR‐oriented KT2 functional of Keal–Tozer characterized by a fairly small error of only 30 ppm for the span of about 1700 ppm (<2%). © 2015 Wiley Periodicals, Inc.     

Local physical quantities for spin based on the relativistic quantum field theory in molecular systems

Local physical quantities for spin are investigated on the basis of the four‐ and two‐component relativistic quantum theory. In the quantum field theory, local physical quantities for spin such as the spin angular momentum density, spin torque density, zeta force density, and zeta potential play important roles in spin dynamics. We discuss how to calculate these local physical quantities based on the two‐component relativistic quantum theory. Some different types of relativistic numerical calculations of local physical quantities in Li atom and C₆H₆ are demonstrated and compared. Local physical quantities for each orbital are also discussed, and it is seen that a total local zeta potential is given as a result of some cancellation of large contributions from each orbital. © 2016 Wiley Periodicals, Inc.     

Two‐component calculations of spin–orbit effects for a van der Waals molecule Rn2

Electronic structures of the weakly bound Rn₂ were calculated by the two‐component Møller–Plesset second‐order perturbation and coupled‐cluster methods with relativistic effective core potentials including spin–orbit operators. The calculated spin–orbit effects are small, but depend strongly on the size of basis sets and the amount of electron correlations. Magnitudes of spin–orbit effects on D_e (0.7–3.0 meV) and R_e (−0.4∼−2.2 Å) of Rn₂ are comparable to previously reported values based on configuration interaction calculations. A two‐component approach seems to be a promising tool to investigate spin–orbit effects for the weak‐bonded systems containing heavy elements. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 139–143, 1999     

Implementation of the diagonalization‐free algorithm in the self‐consistent field procedure within the four‐component relativistic scheme

A recently developed Thouless‐expansion‐based diagonalization‐free approach for improving the efficiency of self‐consistent field (SCF) methods (Noga and ?imunek, J. Chem. Theory Comput. 2010, 6, 2706) has been adapted to the four‐component relativistic scheme and implemented within the program package ReSpect. In addition to the implementation, the method has been thoroughly analyzed, particularly with respect to cases for which it is difficult or computationally expensive to find a good initial guess. Based on this analysis, several modifications of the original algorithm, refining its stability and efficiency, are proposed. To demonstrate the robustness and efficiency of the improved algorithm, we present the results of four‐component diagonalization‐free SCF calculations on several heavy‐metal complexes, the largest of which contains more than 80 atoms (about 6000 4‐spinor basis functions). The diagonalization‐free procedure is about twice as fast as the corresponding diagonalization. © 2014 Wiley Periodicals, Inc.     

New implementation of high‐level correlated methods using a general block tensor library for high‐performance electronic structure calculations

This article presents an open‐source object‐oriented C++ library of classes and routines to perform tensor algebra. The primary purpose of the library is to enable post‐Hartree–Fock electronic structure methods; however, the code is general enough to be applicable in other areas of physical and computational sciences. The library supports tensors of arbitrary order (dimensionality), size, and symmetry. Implemented data structures and algorithms operate on large tensors by splitting them into smaller blocks, storing them both in core memory and in files on disk, and applying divide‐and‐conquer‐type parallel algorithms to perform tensor algebra. The library offers a set of general tensor symmetry algorithms and a full implementation of tensor symmetries typically found in electronic structure theory: permutational, spin, and molecular point group symmetry. The Q‐Chem electronic structure software uses this library to drive coupled‐cluster, equation‐of‐motion, and algebraic‐diagrammatic construction methods. © 2013 Wiley Periodicals, Inc.     

Single‐ and double‐electron reductions of CO2 by using superalkalis: An ab initio study

CO₂, a major contributor to global warming, can be balanced by converting it into fuels. The reduction of CO₂ has been difficult due to its extremely high stability. Recently, single‐electron reduction of CO₂ by superalkalis has been proposed using quantum chemical methods. Herein, we report a systematic study on the single‐reduction of CO₂ by using typical superalkalis. Superalkalis are hypervalent species possessing lower ionization energies than alkali atoms. We have studied the interaction of CO₂ with FLi₂, OLi₃, and NLi₄ superalkalis using ab initio MP2 calculations. We notice that this interaction leads to stable superalkali‐CO₂ complexes in which the structure of CO₂ is bent due to electron transfer from superalkalis. This clearly reveals that the CO₂ can successfully be reduced to the anion. It has been also noticed that the size of superalkalis plays a crucial in the single‐electron reduction of CO₂. For instance, the binding energy of superalkali‐CO₂ complex and charge transfer to CO₂ decreases monotonically with the increase in the size of superalkali. We have also proposed that CO₂ can be further reduced to in case of the anionic complex such as (FLi₂ CO₂)‾. Thus, FLi₂ superalkali is also capable of double‐electron reduction of CO₂. These findings should provide new insights into CO₂‐activation as well as motivate further research in this direction.     

Normal halogen dependence of 13C NMR chemical shifts of halogenomethanes revisited at the four‐component relativistic level

The ‘Normal Halogen Dependence’ of ¹³C NMR chemical shifts in the series of halogenomethanes is revisited at the four‐component relativistic level. Calculations of ¹³C NMR chemical shifts of 70 halogenomethanes have been carried out at the density functional theory (DFT) and MP2 levels with taking into account relativistic effects using the four‐component relativistic theory of Dirac‐Coulomb within the different computational methods (4RPA, 4OPW91) and hybrid computational schemes (MP2 + 4RPA, MP2 + 4OPW91). The most efficient computational protocols are derived for practical purposes. Relativistic shielding effect reaches as much as several hundreds of ppm for heavy halogenomethanes, and to account for this effect in comparison with experiment at the qualitative level, relativistic Dyall's basis sets of triple‐zeta quality or higher are to be used within the framework of the four‐component relativistic theory taking into account solvent effects. Relativistic geometrical optimization (as compared with the non‐relativistic level) is essential for the molecules containing at least two iodines at one carbon atom. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of exponent values in Gaussian‐type functions for development of protonic and deuteronic basis functions

We analyzed the exponent (α) values in Gaussian‐type functions (GTF) for protons and deuterons in BH₃, CH₄, NH₃, H₂O, HF, and their deuterated molecules for the development of nuclear basis functions, which are used for molecular orbital (MO) calculations that directly include nuclear quantum effects. The optimized α (α_opt) value in the single s‐type ([1s]) GTF for protons is changed due to the difference in flexibility of the electronic basis sets. The difference between the energy obtained by using the α_opt value for each molecule and that obtained by using the average α (α_ave) value for these exponents with the 6‐31G(d,p) electronic basis function is only 2 × 10^?5 a.u. The α_ave values of protonic and deuteronic [1s] GTFs by the present calculation are 24.1825 and 35.6214, respectively. We found that the α_ave values enable the evaluation of the total energy and the geometrical changes in hydrogen bonding, such as O…H? O, O…H? N, and O…H? C, while the α_opt value became small by forming a hydrogen bond. The result using only the [1s] GTF for the protonic and deuteronic basis functions is sufficient to explain the differences of energy and geometry induced by the H/D isotope effect, although the total energy of ～5 × 10^?4 a.u. was improved by using the s‐, p‐, and d‐type ([1s1p1d]) GTFs for protons and deuterons. We clearly demonstrate that the protonic and deuteronic basis functions based on the α_ave value enable us to apply the method to other sample molecules (glycine, malonaldehyde, and formic acid dimer). The protonic and deuteronic basis functions we developed treat the quantum effects of protons and deuterons effectively and extend the application range of the MO calculation to include nuclear quantum effects. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Ab initio path‐integral calculations of kinetic and equilibrium isotope effects on base‐catalyzed RNA transphosphorylation models

Detailed understandings of the reaction mechanisms of RNA catalysis in various environments can have profound importance for many applications, ranging from the design of new biotechnologies to the unraveling of the evolutionary origin of life. An integral step in the nucleolytic RNA catalysis is self‐cleavage of RNA strands by 2′‐O‐transphosphorylation. Key to elucidating a reaction mechanism is determining the molecular structure and bonding characteristics of transition state. A direct and powerful probe of transition state is measuring isotope effects on biochemical reactions, particularly if we can reproduce isotope effect values from quantum calculations. This article significantly extends the scope of our previous joint experimental and theoretical work in examining isotope effects on enzymatic and nonenzymatic 2′‐O‐transphosphorylation reaction models that mimic reactions catalyzed by RNA enzymes (ribozymes), and protein enzymes such as ribonuclease A (RNase A). Native reactions are studied, as well as reactions with thio substitutions representing chemical modifications often used in experiments to probe mechanism. Here, we report and compare results from eight levels of electronic‐structure calculations for constructing the potential energy surfaces in kinetic and equilibrium isotope effects (KIE and EIE) computations, including a “gold‐standard” coupled‐cluster level of theory [CCSD(T)]. In addition to the widely used Bigeleisen equation for estimating KIE and EIE values, internuclear anharmonicity and quantum tunneling effects were also computed using our recently developed ab initio path‐integral method, that is, automated integration‐free path‐integral method. The results of this work establish an important set of benchmarks that serve to guide calculations of KIE and EIE for RNA catalysis. © 2014 Wiley Periodicals, Inc.     

Indirect relativistic bridge and substituent effects from the ‘heavy’ environment on the one‐bond and two‐bond 13C1H spin–spin coupling constants

Indirect relativistic bridge effect (IRBE) and indirect relativistic substituent effect (IRSE) induced by the ‘heavy’ environment of the IV‐th, V‐th and VI‐th main group elements on the one‐bond and geminal ¹³C? ¹H spin–spin coupling constants are observed, and spin‐orbit parts of these two effects were interpreted in terms of the third‐order Rayleigh–Schrödinger perturbation theory. Both effects, IRBE and IRSE, rapidly increase with the total atomic charge of the substituents at the coupled carbon. The accumulation of IRSE for geminal coupling constants is not linear with respect to the number of substituents in contrast to the one‐bond couplings where IRSE is an essentially additive quantity. Copyright © 2015 John Wiley & Sons, Ltd.     

Temperature‐dependent conformational analysis of cyclononane: An ab initio study

Quantum chemical methods were used for the theoretical determination of the conformational population for the relevant conformers of cyclononane, i.e., TBC, TCB, TCC, and M4 (or C1), which have been previously investigated experimentally through detailed examination of the nuclear magnetic resonance (NMR) spectrum. Our best Gibbs free energy result, evaluated with MP4(SDTQ)/6‐31G(d,p)//MP2/6‐31G(d,p) energy differences and MP2/6‐31G(d,p) thermal corrections, lead to a temperature‐dependent population in excellent agreement with the experimental results based on the analysis of the low temperature ¹³C NMR spectrum. The nice agreement with experiment is achieved using MP2 harmonic frequencies for the evaluation of vibration partition functions within the standard statistic thermodynamics formalism. Theoretical temperature‐dependent infrared (IR) and ¹³C NMR spectra were simulated and compared with experimental data, which confirmed the ab initio conformational population reported here. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Development of spin‐dependent relativistic open‐shell Hartree–Fock theory with time‐reversal symmetry (I): The unrestricted approach

An open‐shell Hartree–Fock (HF) theory for spin‐dependent, two‐component relativistic calculations, termed the Kramers‐unrestricted HF (KUHF) method, is developed. The present KUHF method, which is formulated as a relativistic counterpart of nonrelativistic UHF, is based on quaternion algebra and partly uses time‐reversal symmetry. The fundamental characteristics of KUHF are discussed in this study. From numerical assessments, it was revealed that KUHF gives a corresponding solution to nonrelativistic UHF; furthermore, KUHF properly describes spin‐orbit interactions. In addition, KUHF can improve the self‐consistent field convergence behavior in spin‐dependent calculations, for example, for f‐block elements.     

Development of spin‐dependent relativistic open‐shell Hartree–Fock theory with time‐reversal symmetry (II): The restricted open‐shell approach

An open‐shell Hartree–Fock (HF) theory for spin‐dependent two‐component relativistic calculations, termed the Kramers‐restricted open‐shell HF (KROHF) method, is developed. The present KROHF method is defined as a relativistic analogue of ROHF using time‐reversal symmetry and quaternion algebra, based on the Kramers‐unrestricted HF (KUHF) theory reported in our previous study (Int. J. Quantum Chem., doi: 10.1002/qua.25356 ). As seen in the nonrelativistic ROHF theory, the ambiguity of the KROHF Fock operator gives physically meaningless spinor energies. To avoid this problem, the canonical parametrization of KROHF to satisfy Koopmans' theorem is also discussed based on the procedure proposed by Plakhutin et al. (J. Chem. Phys. 2006 , 125, 204110). Numerical assessments confirmed that KROHF using Plakhutin's canonicalization procedure correctly gives physical spinor energies within the frozen‐orbital approximation under spin–orbit interactions.     

Evaluation of two‐center overlap and nuclear attraction integrals over Slater‐type orbitals with integer and noninteger principal quantum numbers

A general formula has been established for the expansion of the product of two normalized associated Legendre functions centered on the nuclei a and b. This formula has been utilized for the evaluation of two‐center overlap and nuclear attraction integrals over Slater‐type orbitals (STOs) with integer and noninteger principal quantum numbers. The formulas given in this study for the evaluation of two‐center overlap and nuclear attraction integrals show good rate of convergence and great numerical stability under wide range of quantum numbers, orbital exponents, and internuclear distances. © 2001 Wiley Periodicals, Inc. Int J Quantum Chem, 2001     

Comparison of basis set effects and the performance of ab initio and DFT methods for probing equilibrium fluctuations

The electronic absorption and emission spectra of large molecules reflect the extent and timescale of electron-vibration coupling and therefore the extent and timescale of relaxation/reorganization in response to a perturbation. In this paper, we present a comparison of the calculated absorption and emission spectra of NADH in liver alcohol dehydrogenase (LADH), using quantum mechanical/molecular mechanical methods, in which we vary the QM component. Specifically, we have looked at the influence of basis set (STO-3G, 3-21G*, 6-31G*, CC-pVDZ, and 6-311G**), as well as the influence of applying the DFT TD-B3LYP and ab initio TD-HF and CIS methods to the calculation of absorption/emission spectra and the reorganization energy (Stokes shift). The ab initio TD-HF and CIS methods reproduce the experimentally determined Stokes shift and spectral profiles to a high level of agreement, while the TD-B3LYP method significantly underestimates the Stokes shift, by 45%. We comment on the origin of this problem and suggest that ab initio methods may be naturally more suited to predicting molecular behavior away from equilibrium geometries.     

ChemInform Abstract: An ab initio Theoretical Study of the Electronic Structure of UO2+ and [UO2(CO3)3]5‐.

ChemInform is a weekly Abstracting Service, delivering concise information at a glance that was extracted from about 200 leading journals. To access a ChemInform Abstract of an article which was published elsewhere, please select a “Full Text” option. The original article is trackable via the “References” option.