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.     

Infinite-order quasirelativistic density functional method based on the exact matrix quasirelativistic theory

The exact one-electron matrix quasirelativistic theory [Kutzelnigg and Liu, J. Chem. Phys. 123, 241102 (2005)] is extended to the effective one-particle Kohn-Sham scheme of density functional theory. Several variants of the resultant theory are discussed. Although they are in principle equivalent, consideration of computational efficiency strongly favors the one (F(+)) in which the effective potential remains untransformed. Further combined with the atomic approximation for the matrix X relating the small and large components of the Dirac spinors as well as a simple ansatz for correcting the two-electron picture change errors, a very elegant, accurate, and efficient infinite-order quasirelativistic approach is obtained, which is far simpler than all existing quasirelativistic theories and must hence be regarded as a breakthrough in relativistic quantum chemistry. In passing, it is also shown that the Dirac-Kohn-Sham scheme can be made as efficient as two-component approaches without compromising the accuracy. To demonstrate the performance of the new methods, atomic calculations on Hg and E117 are first carried out. The spectroscopic constants (bond length, vibrational frequency, and dissociation energy) of E117(2) are then reported. All the results are in excellent agreement with those of the Dirac-Kohn-Sham calculations.     

Accurate and efficient treatment of two-electron contributions in quasirelativistic high-order Douglas-Kroll density-functional calculations

Two-component quasirelativistic approaches are in principle capable of reproducing results from fully relativistic calculations based on the four-component Dirac equation (with fixed particle number). For one-electron systems, this also holds in practice, but in many-electron systems one has to transform the two-electron interaction, which is necessary because a picture change occurs when going from the Dirac equation to a two-component method. For one-electron properties, one can take full account of picture change in a manageable way, but for the electron interaction, this would spoil the computational advantages which are the main reason to perform quasirelativistic calculations. Exploiting those picture change effects are largest in the atomic cores, which in molecular applications do not differ too much from the cores of isolated neutral atoms, we propose an elegant, efficient, and accurate approximation to the two-electron picture change problem. The new approach, called the "model potential" approach because it makes use of atomic (four- and two-component) data to estimate picture change effects in molecules, shares with the nuclear-only approach that the Douglas-Kroll operator needs to be constructed only once (not in each self-consistent-field iteration) and that no time-consuming multicenter relativistic two-electron integrals need to be calculated. The new approach correctly describes the screening of both the nearest nucleus and distant nuclei, for the scalar-relativistic as well as the spin-orbit parts of the Hamiltonian. The approach is tested on atomic and molecular-orbital energies as well as spectroscopic constants of the lead dimer.     

Gradients for two-component quasirelativistic methods. Application to dihalogenides of element 116

The authors report the implementation of geometry gradients for quasirelativistic two-component Hartree-Fock and density functional methods using either the zero-order regular approximation Hamiltonian or spin-dependent effective core potentials. The computational effort of the resulting program is comparable to that of corresponding nonrelativistic calculations, as it is dominated by the evaluation of derivative two-electron integrals, which is the same for both types of calculations. Besides the implementation of derivatives of matrix elements of the one-particle Hamiltonian with respect to nuclear displacements, the calculation of the derivative exchange-correlation energy for the open shell case involves complicated expressions because of the noncollinear approach chosen to define the spin density. A pilot application to dihalogenides of element 116 shows how spin-orbit coupling strongly affects the chemistry of the superheavy p-block elements. While these molecules are bent at a scalar-relativistic level, spin-orbit coupling is so strong that only the 7p3/2 atomic orbitals of element 116 are involved in bonding, which favors linear molecular geometries for dihalogenides with heavy terminal halogen atoms.     

A gauge-independent zeroth-order regular approximation to the exact relativistic Hamiltonian--formulation and applications

A simple modification of the zeroth-order regular approximation (ZORA) in relativistic theory is suggested to suppress its erroneous gauge dependence to a high level of approximation. The method, coined gauge-independent ZORA (ZORA-GI), can be easily installed in any existing nonrelativistic quantum chemical package by programming simple one-electron matrix elements for the quasirelativistic Hamiltonian. Results of benchmark calculations obtained with ZORA-GI at the Hartree-Fock (HF) and second-order Moller-Plesset perturbation theory (MP2) level for dihalogens X(2) (X=F,Cl,Br,I,At) are in good agreement with the results of four-component relativistic calculations (HF level) and experimental data (MP2 level). ZORA-GI calculations based on MP2 or coupled-cluster theory with single and double perturbations and a perturbative inclusion of triple excitations [CCSD(T)] lead to accurate atomization energies and molecular geometries for the tetroxides of group VIII elements. With ZORA-GI/CCSD(T), an improved estimate for the atomization energy of hassium (Z=108) tetroxide is obtained.     

A new parametrizable model of molecular electronic structure

A new electronic structure model is developed in which the ground state energy of a molecular system is given by a Hartree-Fock-like expression with parametrized one- and two-electron integrals over an extended (minimal + polarization) set of orthogonalized atom-centered basis functions, the variational equations being solved formally within the minimal basis but the effect of polarization functions being included in the spirit of second-order perturbation theory. It is designed to yield good dipole polarizabilities and improved intermolecular potentials with dispersion terms. The molecular integrals include up to three-center one-electron and two-center two-electron terms, all in simple analytical forms. A method to extract the effective one-electron Hamiltonian of nonlocal-exchange Kohn-Sham theory from the coupled-cluster one-electron density matrix is designed and used to get its matrix representation in a molecule-intrinsic minimal basis as an input to the parametrization procedure--making a direct link to the correlated wavefunction theory. The model has been trained for 15 elements (H, Li-F, Na-Cl, 720 parameters) on a set of 5581 molecules (including ions, transition states, and weakly bound complexes) whose first- and second-order properties were computed by the coupled-cluster theory as a reference, and a good agreement is seen. The model looks promising for the study of large molecular systems, it is believed to be an important step forward from the traditional semiempirical models towards higher accuracy at nearly as low a computational cost.     

Analytical evaluation of first-order electrical properties based on the spin-free Dirac-Coulomb Hamiltonian

We report an analytical scheme for the calculation of first-order electrical properties using the spin-free Dirac-Coulomb (SFDC) Hamiltonian, thereby exploiting the well-developed density-matrix formulations in nonrelativistic coupled-cluster (CC) derivative theory. Orbital relaxation effects are fully accounted for by including the relaxation of the correlated orbitals with respect to orbitals of all types, viz., frozen-core, occupied, virtual, and negative energy state orbitals. To demonstrate the applicability of the presented scheme, we report benchmark calculations for first-order electrical properties of the hydrogen halides, HX with X = F, Cl, Br, I, At, and a first application to the iodo(fluoro)methanes, CH(n)F(3 - n)I, n = 0-3. The results obtained from the SFDC calculations are compared to those from nonrelativistic calculations, those obtained via leading-order direct perturbation theory as well as those from full Dirac-Coulomb calculations. It is shown that the full inclusion of spin-free (SF) relativistic effects is necessary to obtain accurate first-order electrical properties in the presence of fifth-row elements. The SFDC scheme is also recommended for applications to systems containing lighter elements because it introduces no extra cost in the rate-determining steps of a CC calculation in comparison to the nonrelativistic case. On the other hand, spin-orbit contributions are generally small for first-order electrical properties of closed-shell molecules and may be handled efficiently by means of perturbation theory.     

Nonvariational time-dependent multiconfiguration self-consistent field equations for electronic dynamics in laser-driven molecules

A time-dependent multiconfiguration self-consistent field (TDMCSCF) scheme is developed to describe the time-resolved electron dynamics of a laser-driven many-electron atomic or molecular system, starting directly from the time-dependent Schrodinger equation for the system. This nonvariational formulation aims at the full exploitations of concepts, tools, and facilities of existing, well-developed quantum chemical MCSCF codes. The theory uses, in particular, a unitary representation of time-dependent configuration mixings and orbital transformations. Within a short-time, or adiabatic approximation, the TDMCSCF scheme amounts to a second-order split-operator algorithm involving generically the two noncommuting one-electron and two-electron parts of the time-dependent electronic Hamiltonian. We implement the scheme to calculate the laser-induced dynamics of the two-electron H2 molecule described within a minimal basis, and show how electron correlation is affected by the interaction of the molecule with a strong laser field.     

An infinite-order two-component relativistic Hamiltonian by a simple one-step transformation

The authors report the implementation of a simple one-step method for obtaining an infinite-order two-component (IOTC) relativistic Hamiltonian using matrix algebra. They apply the IOTC Hamiltonian to calculations of excitation and ionization energies as well as electric and magnetic properties of the radon atom. The results are compared to corresponding calculations using identical basis sets and based on the four-component Dirac-Coulomb Hamiltonian as well as Douglas-Kroll-Hess and zeroth-order regular approximation Hamiltonians, all implemented in the DIRAC program package, thus allowing a comprehensive comparison of relativistic Hamiltonians within the finite basis approximation.     

Perturbative treatment of scalar-relativistic effects in coupled-cluster calculations of equilibrium geometries and harmonic vibrational frequencies using analytic second-derivative techniques

An analytic scheme for the computation of scalar-relativistic corrections to nuclear forces is presented. Relativistic corrections are included via a perturbative treatment involving the mass-velocity and the one-electron and two-electron Darwin terms. Such a scheme requires mixed second derivatives of the nonrelativistic energy with respect to the relativistic perturbation and the nuclear coordinates and can be implemented using available second-derivative techniques. Our implementation for Hartree-Fock self-consistent field, second-order Moller-Plesset perturbation theory, as well as the coupled-cluster level is used to investigate the relativistic effects on the geometrical parameters and harmonic vibrational frequencies for a set of molecules containing light elements (HX, X=F, Cl, Br; H2X, X=O, S; HXY, X=O, S and Y=F, Cl, Br). The focus of our calculations is the basis-set dependence of the corresponding relativistic effects, additivity of electron correlation and relativistic effects, and the importance of core correlation on relativistic effects.     

Computation scheme for optimizing multiconfigurational wave-functions

A computation scheme is proposed to determine the wave-functions of molecular systems within the framework of the CMC SCF theory and the APSG SCF approach. The orbital optimization is carried out by the refined first-order one-electron Hamiltonian method. Explicit expressions of the first and second energy derivatives are obtained. In the suggested scheme all the calculations are based on using the matrices of the “partial” Coulomb and the exchange operators constructed over the orbitals at the current iteration cycle.     

Auxiliary basis sets for density-fitted correlated wavefunction calculations: weighted core-valence and ECP basis sets for post-d elements

We report optimised auxiliary basis sets for the resolution-of-the-identity (or density-fitting) approximation of two-electron integrals in second-order M?ller-Plesset perturbation theory (MP2) and similar electronic structure calculations with correlation-consistent basis sets for the post-d elements Ga-Kr, In-Xe, and Tl-Rn. The auxiliary basis sets are optimised such that the density-fitting error is negligible compared to the one-electron basis set error. To check to which extent this criterion is fulfilled we estimated for a test set of 80 molecules the basis set limit of the correlation energy at the MP2 level and evaluated the remaining density-fitting and the one-electron basis set errors. The resulting auxiliary basis sets are only 2-6 times larger than the corresponding one-electron basis sets and lead in MP2 calculations to speed-ups of the integral evaluation by one to three orders of magnitude. The density-fitting errors in the correlation energy are at least hundred times smaller than the one-electron basis set error, i.e. in the order of only 1-100 μH per atom.     

Efficient treatment of the Hartree interaction in the relativistic Kohn-Sham problem

We elaborate the two-component Douglas-Kroll reduction of the Dirac-Kohn-Sham problem of relativistic density-functional theory as introduced by Matveev and Rosch [J. Chem. Phys. 118, 3997 (2003)]. That method retains corrections to the Coulomb self-interaction (or Hartree) term of the energy functional that are due to the picture change. Using analytic expressions for the matrix elements, one is able to abandon the resolution of the identity approach for a crucial step of the relativistic transformation. Thus, a major source of uncertainties of the method is eliminated because basis sets no longer have to be extended by functions of higher angular momentum, previously required to ensure kinetic balance. This approach also relies on the electron charge-density fitting scheme via an auxiliary basis set. An efficient approximate implementation results if one restricts the relativistic transformation to the spherically symmetric atom-centered auxiliary functions. It provides accurate results while simplifying greatly the expressions for the matrix elements of the relativistically transformed operators and significantly reducing the computational effort. We demonstrate the performance of the method for the fine structure of one-electron levels of the Hg atom, the g-tensor shifts of NO2, and the properties of the diatomic molecules Bi2, Pb2, PbO, and TlH.     

The MC SCF theory: Method of one-electron Hamiltonian

The one-electron Hamiltonian method is developed to solve the variational equations of the MC SCF theory. The many-parameter family of the one-electron Hamiltonian is derived and conditions for parameters to provide convergence of the SCF procedure to the energy minimum are obtained. A computation scheme based on the use of the one-electron Hamiltonian is described.     

A closed-shell coupled-cluster treatment of the Breit-Pauli first-order relativistic energy correction

First-order relativistic corrections to the energy of closed-shell molecular systems are calculated, using all terms in the two-component Breit-Pauli Hamiltonian. In particular, we present the first implementation of the two-electron Breit orbit-orbit integrals, thus completing the first-order relativistic corrections within the two-component Pauli approximation. Calculations of these corrections are presented for a series of small and light molecules, at the Hartree-Fock and coupled-cluster levels of theory. Comparisons with four-component Dirac-Coulomb-Breit calculations demonstrate that the full Breit-Pauli energy corrections represent an accurate approximation to a fully relativistic treatment of such systems. The Breit interaction is dominated by the spin-spin interaction, the orbit-orbit interaction contributing only about 10% to the total two-electron relativistic correction in molecules consisting of light atoms. However, the relative importance of the orbit-orbit interaction increases with increasing nuclear charge, contributing more than 20% in H(2)S.     

Convergence of approximate two-component Hamiltonians: how far is the Dirac limit

A systematic elimination of the off-diagonal parts of the Dirac Hamiltonian is carried out in the spirit of the Douglas-Kroll [Ann. Phys. 82, 87 1974] approach and the recently proposed infinite-order two-component method. The present approach leads to a series of approximate two-component Hamiltonians which are exact through a certain order in the external potential. These Hamiltonians are used to study the convergence pattern of approximate two-component theories. It is shown that to achieve an acceptably high accuracy for low-lying one-electron levels in heavy and superheavy systems one needs to use approximate Hamiltonians of prohibitively high order in the external potential. One can conclude that the finite-order two-component Hamiltonians are of limited usefulness in accurate relativistic calculations for heavy and superheavy systems.     

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.     

Comments on adiabatic calculations of the electrical properties of molecules

The Hamiltonian used in recent calculations (by Silverman, Bishop and others) of the electrical properties of one-electron diatomic molecules is compared, in the adiabatic approximation, with the Stark Hamiltonian.