Calculation of electric-field gradients based on higher-order generalized Douglas-Kroll transformations

In this paper, the calculation of electric-field-like properties based on higher-order Douglas-Kroll-Hess (DKH) transformations is discussed. The electric-field gradient calculated within the Hartree-Fock self-consistent field framework is used as a representative property. The properties are expressed as an analytic first derivative of the four-component Dirac energy and the nth-order DKH energy, respectively. The differences between a "forward" transformation of the relativistic energy or the "back transformation" of the wave function is discussed in some detail. Detailed test calculations were carried out on the electric-field gradient at the halogen nucleus in the series HX (X=F,Cl,Br,I,At) for which extensive reference data are available. The DKH method is shown to reproduce (spin-free) four-component Dirac-Fock results to an accuracy of better than 99% which is significantly closer than previous DKH studies. The calculations of both the Hamiltonian and the property operator are shown to be essentially converged after the second-order transformation, even for elements as heavy as At. In addition, we have obtained results within the density-functional framework using the DKHZ and zeroth-order regular approximation (ZORA) methods. The latter results included picture-change effects at the scalar relativistic variant of the ZORA-4 level and were shown to be in quantitative agreement with earlier results obtained by van Lenthe and Baerends. The picture-change effects are somewhat smaller for the ZORA method compared to DKH. For heavier elements significant differences in the field gradients predicted by the two methods were found. Based on comparison with four-component Dirac-Kohn-Sham calculations, the DKH results are more accurate. Compared to the spin-free Dirac-Kohn-Sham reference values, the ZORA-4 formalism did not improve the results of the ZORA calculations.     

Analytic high-order Douglas-Kroll-Hess electric field gradients

In this work we present a comprehensive study of analytical electric field gradients in hydrogen halides calculated within the high-order Douglas-Kroll-Hess (DKH) scalar-relativistic approach taking picture-change effects analytically into account. We demonstrate the technical feasibility and reliability of a high-order DKH unitary transformation for the property integrals. The convergence behavior of the DKH property expansion is discussed close to the basis set limit and conditions ensuring picture-change-corrected results are determined. Numerical results are presented, which show that the DKH property expansion converges rapidly toward the reference values provided by four-component methods. This shows that in closed-shell cases, the scalar-relativistic DKH(2,2) approach which is of second order in the external potential for both orbitals and property operator yields a remarkable accuracy. As a parameter-dependence-free high-order DKH model, we recommend DKH(4,3). Moreover, the effect of a finite-nucleus model, different parametrization schemes for the unitary matrices, and the reliability of standard basis sets are investigated.     

Exact decoupling of the Dirac Hamiltonian. I. General theory

Exact decoupling of positive- and negative-energy states in relativistic quantum chemistry is discussed in the framework of unitary transformation techniques. The obscure situation that each scheme of decoupling transformations relies on different, but very special parametrizations of the employed unitary matrices is critically analyzed. By applying the most general power series ansatz for the parametrization of the unitary matrices it is shown that all transformation protocols for decoupling the Dirac Hamiltonian have necessarily to start with an initial free-particle Foldy-Wouthuysen step. The purely numerical iteration scheme applying X-operator techniques to the Barysz-Sadlej-Snijders (BSS) Hamiltonian is compared to the analytical schemes of the Foldy-Wouthuysen (FW) and Douglas-Kroll-Hess (DKH) approaches. Relying on an illegal 1/c expansion of the Dirac Hamiltonian around the nonrelativistic limit, any higher-order FW transformation is in principle ill defined and doomed to fail, irrespective of the specific features of the external potential. It is shown that the DKH method is the only valid analytic unitary transformation scheme for the Dirac Hamiltonian. Its exact infinite-order version can be realized purely numerically by the BSS scheme, which is only able to yield matrix representations of the decoupled Hamiltonian but no analytic expressions for this operator. It is explained why a straightforward numerical iterative extension of the DKH procedure to arbitrary order employing matrix representations is not feasible within standard one-component electronic structure programs. A more sophisticated ansatz based on a symbolical evaluation of the DKH operators via a suitable parser routine is needed instead and introduced in Part II of this work.     

Local relativistic exact decoupling

We present a systematic hierarchy of approximations for local exact decoupling of four-component quantum chemical Hamiltonians based on the Dirac equation. Our ansatz reaches beyond the trivial local approximation that is based on a unitary transformation of only the atomic block-diagonal part of the Hamiltonian. Systematically, off-diagonal Hamiltonian matrix blocks can be subjected to a unitary transformation to yield relativistically corrected matrix elements. The full hierarchy is investigated with respect to the accuracy reached for the electronic energy and for selected molecular properties on a balanced test molecule set that comprises molecules with heavy elements in different bonding situations. Our atomic (local) assembly of the unitary exact-decoupling transformation--called local approximation to the unitary decoupling transformation (DLU)--provides an excellent local approximation for any relativistic exact-decoupling approach. Its order-N(2) scaling can be further reduced to linear scaling by employing a neighboring-atomic-blocks approximation. Therefore, DLU is an efficient relativistic method well suited for relativistic calculations on large molecules. If a large molecule contains many light atoms (typically hydrogen atoms), the computational costs can be further reduced by employing a well-defined nonrelativistic approximation for these light atoms without significant loss of accuracy. We also demonstrate that the standard and straightforward transformation of only the atomic block-diagonal entries in the Hamiltonian--denoted diagonal local approximation to the Hamiltonian (DLH) in this paper--introduces an error that is on the order of the error of second-order Douglas-Kroll-Hess (i.e., DKH2) when compared with exact-decoupling results. Hence, the local DLH approximation would be pointless in an exact-decoupling framework, but can be efficiently employed in combination with the fast to evaluate DKH2 Hamiltonian in order to speed up calculations for which ultimate accuracy is not the major concern.     

Douglas–Kroll–Hess Theory: a relativistic electrons-only theory for chemistry

A unitary transformation allows to separate (block-diagonalize) the Dirac Hamiltonian into two parts one part: solely describes electrons, while the other gives rise to negative-energy states, which are the so-called positronic states. The block-diagonal form of the Hamiltonian no longer accounts for the coupling of both kinds of states. The positive-energy (‘electrons-only’) part can serve as a ‘fully’ relativistic electrons-only theory, which can be understood as a rigorous basis for chemistry. Recent developments of the Douglas–Kroll–Hess (DKH) method allowed to derive a sequence of expressions, which approximate this electrons-only Hamiltonian up to arbitrary-order. While all previous work focused on the numerical stability and accuracy of these arbitrary-order DKH Hamiltonians, conceptual issues and paradoxa of the method were mostly left aside. In this work, the conceptual side of DKH theory is revisited in order to identify essential aspects of the theory to be distinguished from purely computational consideration.     

A second-quantization framework for the unified treatment of relativistic and nonrelativistic molecular perturbations by response theory

A formalism is presented for the calculation of relativistic corrections to molecular electronic energies and properties. After a discussion of the Dirac and Breit equations and their first-order Foldy-Wouthuysen [Phys. Rev. 78, 29 (1950)] transformation, we construct a second-quantization electronic Hamiltonian, valid for all values of the fine-structure constant alpha. The resulting alpha-dependent Hamiltonian is then used to set up a perturbation theory in orders of alpha(2), using the general framework of time-independent response theory, in the same manner as for geometrical and magnetic perturbations. Explicit expressions are given to second order in alpha(2) for the Hartree-Fock model. However, since all relativistic considerations are contained in the alpha-dependent Hamiltonian operator rather than in the wave function, the same approach may be used for other wave-function models, following the general procedure of response theory. In particular, by constructing a variational Lagrangian using the alpha-dependent electronic Hamiltonian, relativistic corrections can be calculated for nonvariational methods as well.     

Analytic energy gradients for the spin-free exact two-component theory using an exact block diagonalization for the one-electron Dirac Hamiltonian

We report the implementation of analytic energy gradients for the evaluation of first-order electrical properties and nuclear forces within the framework of the spin-free (SF) exact two-component (X2c) theory. In the scheme presented here, referred to in the following as SFX2c-1e, the decoupling of electronic and positronic solutions is performed for the one-electron Dirac Hamiltonian in its matrix representation using a single unitary transformation. The resulting two-component one-electron matrix Hamiltonian is combined with untransformed two-electron interactions for subsequent self-consistent-field and electron-correlated calculations. The "picture-change" effect in the calculation of properties is taken into account by considering the full derivative of the two-component Hamiltonian matrix with respect to the external perturbation. The applicability of the analytic-gradient scheme presented here is demonstrated in benchmark calculations. SFX2c-1e results for the dipole moments and electric-field gradients of the hydrogen halides are compared with those obtained from nonrelativistic, SF high-order Douglas-Kroll-Hess, and SF Dirac-Coulomb calculations. It is shown that the use of untransformed two-electron interactions introduces rather small errors for these properties. As a first application of the analytic geometrical gradient, we report the equilibrium geometry of methylcopper (CuCH(3)) determined at various levels of theory.     

Exact decoupling of the Dirac Hamiltonian. II. The generalized Douglas-Kroll-Hess transformation up to arbitrary order

In order to achieve exact decoupling of the Dirac Hamiltonian within a unitary transformation scheme, we have discussed in part I of this series that either a purely numerical iterative technique (the Barysz-Sadlej-Snijders method) or a stepwise analytic approach (the Douglas-Kroll-Hess method) are possible. For the evaluation of Douglas-Kroll-Hess Hamiltonians up to a pre-defined order it was shown that a symbolic scheme has to be employed. In this work, an algorithm for this analytic derivation of Douglas-Kroll-Hess Hamiltonians up to any arbitrary order in the external potential is presented. We discuss how an estimate for the necessary order for exact decoupling (within machine precision) for a given system can be determined from the convergence behavior of the Douglas-Kroll-Hess expansion prior to a quantum chemical calculation. Once this maximum order has been accomplished, the spectrum of the positive-energy part of the decoupled Hamiltonian, e.g., for electronic bound states, cannot be distinguished from the corresponding part of the spectrum of the Dirac operator. An efficient scalar-relativistic implementation of the symbolic operations for the evaluation of the positive-energy part of the block-diagonal Hamiltonian is presented, and its accuracy is tested for ground-state energies of one-electron ions over the whole periodic table. Furthermore, the first many-electron calculations employing sixth up to fourteenth order DKH Hamiltonians are presented.     

Perturbational relativistic theory of electron spin resonance g-tensor

We carry out a complete treatment of the leading-order relativistic one-electron contributions, arising from the Breit-Pauli Hamiltonian, to the g-tensor of electron spin resonance spectroscopy. We classify the different terms and discuss their interpretation as well as give numerical ab initio estimates for the F2(-), Cl2(-), Br2(-), and I2(-) series, using analytical response theory calculations with a multiconfigurational self-consistent field reference state. The results are compared to available experimental data. (c) 2004 American Institute of Physics     

Understanding the EPR parameters of glycine-derived radicals: the case of N-acetylglycyl in the N-acetylglycine single-crystal environment

As a step toward an in-depth understanding of the electron paramagnetic resonance parameters of glycyl radicals in proteins, the hyperfine tensors and, particularly, the g-tensor of N-acetylglcyl in the environment of a single crystal of N-acetylglycine have been studied by systematic state-of-the-art quantum chemical calculations on various suitable model systems. The quantitative computation of the g-tensors for such glycyl-derived radicals is a veritable challenge, mainly because of the very small g-anisotropy combined with a nonsymmetrical, delocalized spin-density distribution and several atoms with comparable spin-orbit contributions to the g-tensors. The choice of gauge origin of the magnetic vector potential, and of approximate spin-orbit operators, both turn out to be more critical than found in previous studies of g-tensors for organic radicals. Environmental effects, included by supermolecular hydrogen-bonded models, were found to be moderate, because of a partial compensation between the influences from intramolecular and intermolecular hydrogen bonds. The largest effects on the g-tensor are caused by the conformation of the radical. The density functional theory methods employed systematically overestimate both the Delta gx and Delta gy components of the g-tensor. This is important for parallel investigations on the protein-glycyl radicals. The 1H alpha and 13C alpha hyperfine couplings depend only slightly on the supermolecular model chosen and appear less sensitive probes of detailed structure and environment.     

A fully relativistic method for calculation of nuclear magnetic shielding tensors with a restricted magnetically balanced basis in the framework of the matrix Dirac-Kohn-Sham equation

A new relativistic four-component density functional approach for calculations of NMR shielding tensors has been developed and implemented. It is founded on the matrix formulation of the Dirac-Kohn-Sham (DKS) method. Initially, unperturbed equations are solved with the use of a restricted kinetically balanced basis set for the small component. The second-order coupled perturbed DKS method is then based on the use of restricted magnetically balanced basis sets for the small component. Benchmark relativistic calculations have been carried out for the (1)H and heavy-atom nuclear shielding tensors of the HX series (X=F,Cl,Br,I), where spin-orbit effects are known to be very pronounced. The restricted magnetically balanced basis set allows us to avoid additional approximations and/or strong basis set dependence which arises in some related approaches. The method provides an attractive alternative to existing approximate two-component methods with transformed Hamiltonians for relativistic calculations of chemical shifts and spin-spin coupling constants of heavy-atom systems. In particular, no picture-change effects arise in property calculations.     

Harnessing the Combined Power of Theoretical and Experimental Data through Multiscale Informatics

Monumental, recent and rapidly continuing, improvements in the capabilities of ab initio theoretical kinetics calculations provides reason to believe that progress in the field of chemical kinetics can be accelerated through a corresponding evolution of the role of theory in kinetic modeling and its relationship with experiment. The present article reviews and provides additional demonstrations of the unique advantages that arise when theoretical and experimental data across multiple scales are considered on equal footing, including the relevant uncertainties of both, within a single mathematical framework. Namely, the multiscale informatics framework simultaneously integrates information from a wide variety of sources and scales: ab initio electronic structure calculations of molecular properties, rate constant determinations for individual reactions, and measured global observables of multireaction systems. The resulting model representation consists of a set of theoretical kinetics parameters (with constrained uncertainties) that are related through elementary kinetics models to rate constants (with propagated uncertainties) that in turn are related through physical models to global observables (with propagated uncertainties). An overview of the approach and typical implementation is provided along with a brief discussion of the major uncertainties (parametric and structural) in theoretical kinetics calculations, kinetic models for complex chemical mechanisms, and physical models for experiments. Higher levels of automation in all aspects, including closed‐loop autonomous mixed‐experimental‐and‐computational model improvement, are advocated for facilitating scalability of the approach to larger systems with reasonable human effort and computational cost. The unique advantages of combining theoretical and experimental data across multiple scales are illustrated through a series of examples. Previous results demonstrating the utility of simultaneous interpretation of theoretical and experimental data for assessing consistency in complex systems and for reliable, physics‐based extrapolation of limited data are briefly summarized. New results are presented to demonstrate the high predictive accuracy of multiscale informed models for both small (molecular properties) and large (global observables) scales. These new results provide examples where the optimization yields physically realistic parameter adjustments and where physical model uncertainties in experiments are larger than kinetic model uncertainties. New results are also presented to demonstrate the utility of the multiscale informatics approach for design of experiments and theoretical calculations, accounting for both theoretical and experimental existing knowledge as well as relevant parametric and structural uncertainties in interpreting potential new data. These new results provide examples where neglecting structural uncertainties in design of experiments results in failure to identify the most worthwhile experiment. Further progress in the chemical kinetics field (particularly at the intersection of theory, kinetic modeling, and experiment) would benefit from increased attention to understanding parametric and structural uncertainties for all three—the uncertainty magnitude and cross‐correlations among model parameters as well as limitations of the model structures themselves.     

Gauge invariance of the spin-other-orbit contribution to the g-tensors of electron paramagnetic resonance

The spin-other-orbit (SOO) contribution to the g-tensor (DeltagSOO) of electron paramagnetic resonance arises due to the interaction of electron-spin magnetic moment with the magnetic field produced by the orbital motion of other electrons. A similar mechanism is responsible for the leading term in nuclear magnetic-shielding tensors sigma. We demonstrate that analogous to sigma, paramagnetic DeltagSOO contribution exhibits a pronounced dependence on the choice of the magnetic-field gauge. The gauge corrections to DeltagSOO are similar in magnitude, and opposite in sign, to the paramagnetic SOO term. We calculate gauge-invariant DeltagSOO values using gauge-including atomic orbitals and density-functional theory. For organic radicals, complete gauge-invariant DeltagSOO values typically amount to less than 500 parts per million (ppm), and are small compared to other g-tensor contributions. For the first-row transition-metal compounds, DeltagSOO may contribute several thousand ppm to the g-tensor, but are negligible compared to the remaining deviations from experiment. With popular choices for the magnetic-field gauge, the individual gauge-variant contributions may be an order of magnitude higher, and do not provide a reliable estimation of DeltagSOO.     

Connections between perturbation theory and the unitary transformation methods of deriving effective Hamiltonians

The connections between open shell Brillouin–Wigner perturbation theory and the Van Vleck unitary transformation formalisms for generating effective Hamiltonians are explored. An explicit expression is obtained relating the generator ? of the unitary transformation e^i? with the amplitudes to be found from perturbation theory. The “renormalization effects” needed to produce the explicit “orthogonal-Hermitian” form of the effective Hamiltonian in perturbation theory are related directly to the generator of the unitary transformation. The conclusions reached previously by Jørgensen and Brandow regarding the identity of the effective Hamiltonians of the formalisms are explicitly verified for the case that the generator ? satisfied the Kemble condition. The procedure suggests how the powerful techniques of perturbation theory can be used within the unitary transformation framework to guarantee properly renormalized wave functions.     

Resolution of identity Dirac-Kohn-Sham method using the large component only: Calculations of g-tensor and hyperfine tensor

A new relativistic two-component density functional approach, based on the Dirac-Kohn-Sham method and an extensive use of the technique of resolution of identity (RI), has been developed and is termed the DKS2-RI method. It has been applied to relativistic calculations of g and hyperfine tensors of coinage-metal atoms and some mercury complexes. The DKS2-RI method solves the Dirac-Kohn-Sham equations in a two-component framework using explicitly a basis for the large component only, but it retains all contributions coming from the small component. The DKS2-RI results converge to those of the four-component Dirac-Kohn-Sham with an increasing basis set since the error associated with the use of RI will approach zero. The RI approximation provides a basis for a very efficient implementation by avoiding problems associated with complicated integrals otherwise arising from the elimination of the small component. The approach has been implemented in an unrestricted noncollinear two-component density functional framework. DKS2-RI is related to Dyall's [J. Chem. Phys. 106, 9618 (1997)] unnormalized elimination of the small component method (which was formulated at the Hartree-Fock level and applied to one-electron systems only), but it takes advantage of the local Kohn-Sham exchange-correlation operators (as, e.g., arising from local or gradient-corrected functionals). The DKS2-RI method provides an attractive alternative to existing approximate two-component methods with transformed Hamiltonians (such as Douglas-Kroll-Hess [Ann. Phys. 82, 89 (1974); Phys. Rev. A 33, 3742 (1986)] method, zero-order regular approximation, or related approaches) for relativistic calculations of the structure and properties of heavy-atom systems. In particular, no picture-change effects arise in the property calculations.     

Perturbation theory of relativistic corrections

Methods for perturbation theory of relativistic corrections for an electron in a Coulomb field are divided into three categories: (1) in terms of 4-component spinors; (2) in terms of the ‘large components’ of the Dirac spinor; (3) involving a Foldy-Wouthuysen type transformation, where one attempts to obtain a two-component spinor different from the ‘large component’. In methods of category 1 (the ‘direct perturbation theory’ of paper I of this series, the related approaches by Rutkowski as well as by Gesteszy, Grosse, and Thaller and a somewhat different one by Moore) the wave function, the energy and the Hamiltonian are analytic inc ^?2. No divergent terms arise. In methods of category 2 (that of the elemination of the small component as well as a similarity transformation in intermediate normalization) wave function and energy are still analytic inc ^?2, but the effective Hamiltonian no longer is. Regularized results can be obtained by controlled cancellation of divergent terms. In category 3 both the effective Hamiltonian and the wave function are highly singular and non-analytic inc ^?1. A controlled cancellation of divergent terms is at least very difficult. These pathologic feature survive in the non-relativistic limit and have hence little to do with relativistic effects. They are related to the fact that forr → 0 the sign of the quantum number κ rather than that of the energy determines which component of the Dirac spinor is large and which is small. In the limitr → 0 andc → ∞ the Foldy-Wouthuysen wave function of a 2p _1/2 state is a 1p wave function. Hierarchies of transformations of the Dirac equation and its non-relativistic limit are presented and discussed. Finally the problem of the regularization of effective Hamiltonians on 2-component level ‘for electrons only’ is addressed.     

Electronic g-tensors of solvated molecules using the polarizable continuum model

We present the implementation of density functional response theory combined with the polarizable continuum model (PCM), enabling first principles calculations of molecular g-tensors of solvated molecules. The calculated g-tensor shifts are compared with experimental g-tensor shifts obtained from electron paramagnetic resonance spectra for a few solvated species. The results indicate qualitative agreement between the calculations and the experimental data for aprotic solvents, whereas PCM fails to reproduce the electronic g-tensor behavior for protic solvents. This failure of PCM for protic solvents can be resolved by including into the model those solvent molecules which are involved in hydrogen bonding with the solute. The results for the protic solvents show that the explicit inclusion of the solvent molecules of the first solvation sphere is not sufficient in order to reproduce the behavior of the electronic g-tensor in protic solvents, and that better agreement with experimental data can be obtained by including the long-range electrostatic effects accounted for by the PCM approach on top of the explicit hydrogen-bonded complexes.