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.     

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.     

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.     

One-electron contributions to the g-tensor for second-order Douglas-Kroll-Hess theory

The electric g-tensor is a central quantity for the interpretation of electron paramagnetic resonance spectra. In this paper, a detailed derivation of the 1-electron contributions to the g-tensor is presented in the framework of linear response theory and the second-order Douglas-Kroll-Hess (DKH) transformation. Importantly, the DKH transformation in the presence of a magnetic field is not unique. Whether or not the magnetic field is included in the required Foldy-Wouthuysen transformation, different transformation matrices and, consequently, Hamiltonians result. In this paper, a detailed comparison of both approaches is presented, paying particular attention to the mathematical properties of the resulting Hamiltonians. In contrast to previous studies that address the g-tensor in the framework of DKH theory, the resulting terms are compared to those of the conventional Pauli theory and are given a physical interpretation. Based on these mathematical and physical arguments, we establish that the proper DKH transformation for systems with constant magnetic fields is based on a gauge-invariant Foldy-Wouthuysen transformation, i.e., a Foldy-Wouthuysen transformation including the magnetic field. Calculations using density functional theory (DFT) are carried out on a set of heavy, diatomic molecules, and a set of transition-metal complexes. Based on these calculations, the performance of the relativistic calculation with and without inclusion of picture-change effects is compared. Additionally, the g-tensor is calculated for the Lanthanide dihydrides. Together with the results from the other two molecular test sets, these calculations serve to quantify the magnitude of picture-change effects and elucidate trends across the periodic table.     

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. IV. Automated evaluation of molecular properties within the Douglas-Kroll-Hess theory up to arbitrary order

In Part III [J. Chem. Phys. 124, 064102 (2005)] of this series of papers on exact decoupling of the Dirac Hamiltonian within transformation theory, we developed the most general account on how to treat magnetic and electric properties in a unitary transformation theory on the same footing. In this paper we present an implementation of a general algorithm for the calculation of magnetic as well as electric properties within the framework of Douglas-Kroll-Hess theory. The formal and practical principles of this algorithm are described. We present the first high-order Douglas-Kroll-Hess results for property operators. As for model properties we propose to use the well-defined radial moments, i.e., expectation values of r(k), which can be understood as terms of the Taylor-series expansion of any property operator. Such moments facilitate a rigorous comparison of methods free of uncertainties which may arise in a direct comparison with experiment. This is important in view of the fact that various approaches to two-component molecular properties may yield numerically very small terms whose approximate or inaccurate treatment would not be visible in a direct comparison to experimental data or to another approximate computational reference. Results are presented for various degrees of decoupling of the model properties within the Douglas-Kroll-Hess scheme.     

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.     

Atomic approximation to the projection on electronic states in the Douglas-Kroll-Hess approach to the relativistic Kohn-Sham method

We suggest an approximate relativistic model for economical all-electron calculations on molecular systems that exploits an atomic ansatz for the relativistic projection transformation. With such a choice, the projection transformation matrix is by definition both transferable and independent of the geometry. The formulation is flexible with regard to the level at which the projection transformation is approximated; we employ the free-particle Foldy-Wouthuysen and the second-order Douglas-Kroll-Hess variants. The (atomic) infinite-order decoupling scheme shows little effect on structural parameters in scalar-relativistic calculations; also, the use of a screened nuclear potential in the definition of the projection transformation shows hardly any effect in the context of the present work. Applications to structural and energetic parameters of various systems (diatomics AuH, AuCl, and Au(2), two structural isomers of Ir(4), and uranyl dication UO(2) (2+) solvated by 3-6 water ligands) show that the atomic approximation to the conventional second-order Douglas-Kroll-Hess projection (ADKH) transformation yields highly accurate results at substantial computational savings, in particular, when calculating energy derivatives of larger systems. The size-dependence of the intrinsic error of the ADKH method in extended systems of heavy elements is analyzed for the atomization energies of Pd(n) clusters (n相似文献   

Four-component relativistic theory for nuclear magnetic shielding constants: critical assessments of different approaches

Both formal and numerical analyses have been carried out on various exact and approximate variants of the four-component relativistic theory for nuclear magnetic shielding constants. These include the standard linear response theory (LRT), the full or external field-dependent unitary transformations of the Dirac operator, as well as the orbital decomposition approach. In contrast with LRT, the latter schemes take explicitly into account both the kinetic and magnetic balances between the large and small components of the Dirac spinors, and are therefore much less demanding on the basis sets. In addition, the diamagnetic contributions, which are otherwise "missing" in LRT, appear naturally in the latter schemes. Nevertheless, the definitions of paramagnetic and diamagnetic terms are not the same in the different schemes, but the difference is only of O(c(-2)) and thus vanishes in the nonrelativistic limit. It is shown that, as an operator theory, the full field-dependent unitary transformation approach cannot be applied to singular magnetic fields such as that due to the magnetic point dipole moment of a nucleus. However, the inherent singularities can be avoided by the corresponding matrix formulation (with a partial closed summation). All the schemes are combined with the Dirac-Kohn-Sham ansatz for ground state calculations, and by using virtually complete basis sets a new and more accurate set of absolute nuclear magnetic resonance shielding scales for the rare gases He-Rn have been established.     

Quasirelativistic theory. II. Theory at matrix level

The Dirac operator in a matrix representation in a kinetically balanced basis is transformed to the matrix representation of a quasirelativistic Hamiltonian that has the same electronic eigenstates as the original Dirac matrix (but no positronic eigenstates). This transformation involves a matrix X, for which an exact identity is derived and which can be constructed either in a noniterative way or by various iteration schemes, not requiring an expansion parameter. Both linearly convergent and quadratically convergent iteration schemes are discussed and compared numerically. The authors present three rather different schemes, for each of which even in unfavorable cases convergence is reached within three or four iterations, for all electronic eigenstates of the Dirac operator. The authors present the theory both in terms of a non-Hermitian and a Hermitian quasirelativistic Hamiltonian. Quasirelativistic approaches at the matrix level known from the literature are critically analyzed in the frame of the general theory.     

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.     

Flying onto global minima on potential energy surfaces: A swarm intelligence guided route to molecular electronic structure

A novel quantum‐classical recipe for locating the global minimum on the potential energy surface of a large molecule and simultaneously predicting the associated electronic charge distribution is developed by interfacing the classical particle swarm optimization with a near optimal unitary evolution scheme for the trial one electron density matrix. The unitary transformation is generated by an antihermitian matrix linked to the molecular electronic Hamiltonian at the instantaneous nuclear configurations discovered by the swarm as it flies. The algorithm is used to predict the extensive reorganization of electronic charge distribution and bond lengths in polythiophene oligomers on doping at various levels.     

Relation between different variants of the generalized Douglas-Kroll transformation through sixth order

Wolf et al. have recently investigated a generalized Douglas-Kroll transformation. From a general class of unitary transformations that can be used in the Douglas-Kroll transformation, they pick one which is supposed to give, at a given order, an optimal transformed Dirac Hamiltonian. Results were presented through the fifth order. However, no data were given to demonstrate to which extent the so-called "optimal" Douglas-Kroll transformation is superior to other choices. In this work, the Douglas-Kroll transformation is extended to the sixth order for the first time, using computer algebra algorithms to obtain the working equations. It is shown how, at a given order, different variants of the Douglas-Kroll Hamiltonians are related. Various choices of the generalized transformation are examined numerically for the ground states of the one-electron atomic ions with nuclear charges Z=20, 40, 60, 80, 100, and 120. It is shown that compared to the improvement obtained by including the next order, the differences between various choices for the generalized Douglas-Kroll transformation are almost negligible. Results closest to the Dirac eigenvalues are not obtained with the optimal Douglas-Kroll transformation given by Wolf et al., but with the parametrization originally suggested by Douglas and Kroll.     

Analytic second derivatives for the spin-free exact two-component theory

The formulation and implementation of the spin-free (SF) exact two-component (X2c) theory at the one-electron level (SFX2c-1e) is extended in the present work to the analytic evaluation of second derivatives of the energy. In the X2c-1e scheme, the four-component one-electron Dirac Hamiltonian is block diagonalized in its matrix representation and the resulting "electrons-only" two-component Hamiltonian is then used together with untransformed two-electron interactions. The derivatives of the two-component Hamiltonian can thus be obtained by means of simple manipulations of the parent four-component Hamiltonian integrals and derivative integrals. The SF version of X2c-1e can furthermore exploit available nonrelativistic quantum-chemical codes in a straightforward manner. As a first application of analytic SFX2c-1e second derivatives, we report a systematic study of the equilibrium geometry and vibrational frequencies for the bent ground state of the copper hydroxide (CuOH) molecule. Scalar-relativistic, electron-correlation, and basis-set effects on these properties are carefully assessed.     

Density functional based structure optimization for molecules containing heavy elements: analytical energy gradients for the Douglas-Kroll-Hess scalar relativistic approach to the LCGTO-DF method

The self-consistent scalar-relativistic linear combination of Gaussian-type orbitals density functional (LCGTO-DF) method has been extended to calculate analytical energy gradients. The method is based on the use of a unitary second order Douglas-Kroll-Hess (DKH) transformation for decoupling large and small components of the full four-component Dirac-Kohn-Sham equation. The approximate DKH transformation most common in molecular calculations has been implemented; this variant employs nuclear potential based projectors and it leaves the electron-electron interaction untransformed. Examples are provided for the geometry optimization of a series of heavy metal systems which feature a variety of metal-ligand bonds, like Au₂, AuCl, AuH, Mo(CO)₆ and W(CO)₆ as well as the d¹⁰ complexes [Pd(PH₃)₂O₂] and [Pt(PH₃)₂O₂]. The calculated results, obtained with several gradient-corrected exchange-correlation potentials, compare very well with experimental data and they are of similar or even better accuracy than those of other high quality relativistic calculations reported so far.     

A paradox of grid-based representation techniques: accurate eigenvalues from inaccurate matrix elements

Several approximately variational grid-based representation techniques devised to solve the time-independent nuclear-motion Schrödinger equation share a similar behavior: while the computed eigenpairs, the only results which are of genuine interest, are accurate, many of the underlying Hamiltonian matrix elements are inaccurate, deviating substantially from their values in a variational basis representation. Examples are presented for the discrete variable representation and the Lagrange-mesh approaches, demonstrating that highly accurate eigenvalues and eigenfunctions can be obtained even if some or even all of the Hamiltonian matrix elements in these grid-based representations are inaccurate. It is shown how the apparent contradiction of obtaining accurate eigenpairs with far less accurate individual matrix elements can be resolved by considering the unitary transformation between the representations. Furthermore, the relations connecting orthonormal bases and the corresponding Lagrange bases are generalized to relations connecting nonorthogonal, regularized bases and the corresponding nonorthogonal, regularized Lagrange bases.     

Use of a unitary wavefunction in the calculation of static electronic properties

A unitary coupled cluster method is advocated in this paper for the calculation of static properties. Corresponding to the perturbed Hamiltonian H() including the relevant static property, a suitable unitary wavefunction is envisaged. It is shown that a specific nonvariational model of calculating various order static properties utilising this unitary ansatz results in simplifications compared to the previous Coupled Cluster Theories using only hole-particle excitation parameters formulated for this purpose.NCL Communications No. 3533