Second- and third-order triples and quadruples corrections to coupled-cluster singles and doubles in the ground and excited states

Second- and third-order perturbation corrections to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) incorporating excited configurations in the space of triples [EOM-CCSD(2)T and (3)T] or in the space of triples and quadruples [EOM-CCSD(2)TQ] have been implemented. Their ground-state counterparts--third-order corrections to coupled-cluster singles and doubles (CCSD) in the space of triples [CCSD(3)T] or in the space of triples and quadruples [CCSD(3)TQ]--have also been implemented and assessed. It has been shown that a straightforward application of the Rayleigh-Schrodinger perturbation theory leads to perturbation corrections to total energies of excited states that lack the correct size dependence. Approximations have been introduced to the perturbation corrections to arrive at EOM-CCSD(2)T, (3)T, and (2)TQ that provide size-intensive excitation energies at a noniterative O(n(7)), O(n(8)), and O(n(9)) cost (n is the number of orbitals) and CCSD(3)T and (3)TQ size-extensive total energies at a noniterative O(n(8)) and O(n(10)) cost. All the implementations are parallel executable, applicable to open and closed shells, and take into account spin and real Abelian point-group symmetries. For excited states, they form a systematically more accurate series, CCSD1 eV) and the ground-state wave function has single-determinant character. In other cases, however, the corrections tend to overestimate the triples and quadruples effects, the origin of which is discussed. For ground states, the third-order corrections lead to a rather small improvement over the highly effective second-order corrections [CCSD(2)T and (2)TQ], which is a manifestation of the staircase convergence of perturbation series.     

Cost reduction of high-order coupled-cluster methods via active-space and orbital transformation techniques

We discuss several techniques which have the potential to decrease the computational expenses of high-order coupled-cluster (CC) methods with a reasonable loss in accuracy. In particular, the CC singles, doubles, and triples (CCSDT) as well as the CC singles, doubles, triples, and perturbative quadruples [CCSDT(Q)] methods are considered, which are frequently used in high-precision model chemistries for the calculation of iterative triples and quadruples corrections. First, we study the possibilities for using active spaces to decrease the computational costs. In this case, an active space is defined and some indices of cluster amplitudes are restricted to be in the space. Second, the application of transformed virtual orbitals is investigated. In this framework, to reduce the computation time the dimension of the properly transformed virtual one-particle space is truncated. We have found that the orbital transformation techniques outperform the active-space approaches. Using the transformation techniques, the computational time can be reduced in average by an order of magnitude without significant loss in accuracy. It is demonstrated that high-order CC calculations are possible for considerably larger systems than before using the implemented techniques.     

HEAT: High accuracy extrapolated ab initio thermochemistry

A theoretical model chemistry designed to achieve high accuracy for enthalpies of formation of atoms and small molecules is described. This approach is entirely independent of experimental data and contains no empirical scaling factors, and includes a treatment of electron correlation up to the full coupled-cluster singles, doubles, triples and quadruples approach. Energies are further augmented by anharmonic zero-point vibrational energies, a scalar relativistic correction, first-order spin-orbit coupling, and the diagonal Born-Oppenheimer correction. The accuracy of the approach is assessed by several means. Enthalpies of formation (at 0 K) calculated for a test suite of 31 atoms and molecules via direct calculation of the corresponding elemental formation reactions are within 1 kJ mol(-1) to experiment in all cases. Given the quite different bonding environments in the product and reactant sides of these reactions, the results strongly indicate that even greater accuracy may be expected in reactions that preserve (either exactly or approximately) the number and types of chemical bonds.     

A theoretical study of the excited states of AmO2n+, n=1,2,3

The ground and excited states of the AmO(2) (+), AmO(2) (2+), and AmO(2) (3+) ions have been studied using the four-component configuration interaction singles doubles, spin-orbit complete active space self-consistent field, and spin-orbit complete active space-order perturbation theory methods. The roles of scalar relativistic effects and spin-orbit coupling are analyzed; results with different methods are carefully compared by a precise analysis of the wave functions. A molecular spinor diagram is used in relation to the four-component calculations while a ligand field model is used for the two-step method. States with the same number of electrons in the four nonbonding orbitals are in very good agreement with the two methods while ligand field and charge transfer states do not have the same excitation energies.     

Ab initio study of a Bi3+ impurity in Cs2NaYCl6 and Y2O3: comparison of perturbative and variational electron correlation methods

Ab initio study of excitation energies and oscillator strengths for absorption towards the (3)P(1) and (1)P(1) states of the Bi(3+) ion has been performed for the Bi(3+) ion in gas phase and as a dopant of the cubic elpasolite Cs(2)NaYCl(6) and the yttria Y(2)O(3) crystal using the ab initio embedded-cluster method. The ground and excited states were computed with a relativistic spin-orbit configuration interaction approach suited for heavy elements. Electron correlation was treated in the scalar relativistic scheme with perturbative, variational, and coupled-cluster methods. Intermediate coupling is included via an effective-Hamiltonian based spin-orbit configuration interaction approach. Small-core (60 electrons) and large-core (78 electrons) relativistic effective core potentials (ECPs) have been used to describe the bismuth ion. The best match with experiment was obtained with the small-core ECP. The accuracy of excitation energies strongly depends on the electron correlation method used. The agreement between experimental data and the results obtained using second-order multiconfigurational perturbation theory is greatly improved with the shifted zeroth-order Hamiltonian proposed by Ghido et al. [Chem. Phys. Lett. 396, 142 (2004)]. Although quite time consuming, coupled-cluster and variational methods yield good agreement with experimental data. The first absorption band recorded for the doped elpasolite crystal is positioned with an excellent accuracy while the computed energy of the second absorbing manifold is in poorer agreement with experimental data. This suggests that interactions with neglected close-lying excited states with a ligand-to-metal charge transfer character may be significant. Calculations of the spectrum of Bi(3+) doping yttria in both the S(6) and C(2) site symmetries indicate that the absorbing manifold arises from electronic excitations localized on the Bi(3+) doping ion with main triplet 6s6p character. Our results predict the first absorbing peak to lie about 0.5 eV lower for the S(6) sites than for the C(2) site, thus attributing the violet and the green emission wavelengths to the S(6) and C(2) sites, respectively. A subsequent study of Stokes shift and emission wavelength should hopefully lead to a final assignment of the measured excitation spectra.     

The barrier height of the F+H2 reaction revisited: coupled-cluster and multireference configuration-interaction benchmark calculations

Large scale coupled-cluster benchmark calculations have been carried out to determine the barrier height of the F+H2 reaction as accurately as possible. The best estimates for the barrier height of the linear and bent transition states amount to 2.16 and 1.63 kcal/mol, respectively. These values include corrections for core correlation, scalar-relativistic effects, spin-orbit effects, as well as the diagonal Born-Oppenheimer correction. The CCSD(T) basis-set limits are estimated using extrapolation techniques with augmented quintuple and sextuple-zeta basis sets, and remaining N-electron errors are determined using coupled-cluster singles, doubles, triples, quadruples calculations with up to augmented quintuple-zeta basis sets. The remaining uncertainty is estimated to be less than 0.1 kcal/mol. The coupled-cluster results are used to calibrate multireference configuration-interaction calculations with empirical scaling of the correlation energy.     

Accurate relativistic energy-consistent pseudopotentials for the superheavy elements 111 to 118 including quantum electrodynamic effects

Energy-consistent two-component semi-local pseudopotentials for the superheavy elements with atomic numbers 111-118 have been adjusted to fully relativistic multi-configuration Dirac-Hartree-Fock calculations based on the Dirac-Coulomb Hamiltonian, including perturbative corrections for the frequency-dependent Breit interaction in the Coulomb gauge and lowest-order quantum electrodynamic effects. The pseudopotential core includes 92 electrons corresponding to the configuration [Xe]4f(14)5d(10)5f(14). The parameters for the elements 111-118 were fitted by two-component multi-configuration Hartree-Fock calculations in the intermediate coupling scheme to the total energies of 267 up to 797 J levels arising from 31 up to 62 nonrelativistic configurations, including also anionic and highly ionized states, with mean absolute errors clearly below 0.02 eV for averages corresponding to nonrelativistic configurations. Primitive basis sets for one- and two-component pseudopotential calculations have been optimized for the ground and excited states and exhibit finite basis set errors with respect to the finite-difference Hartree-Fock limit below 0.01 and 0.02 eV, respectively. General contraction schemes have been applied to obtain valence basis sets of polarized valence double- to quadruple-zeta quality. Results of atomic test calculations in the intermediate coupling scheme at the Fock-space coupled-cluster level are in good agreement with those of corresponding fully relativistic all-electron calculations based on the Dirac-Coulomb-Breit Hamiltonian. The results demonstrate besides the well-known need of a relativistic treatment at the Dirac-Coulomb level also the necessity to include higher-order corrections for the superheavy elements.     

Two- and four-component relativistic generalized-active-space coupled cluster method: implementation and application to BiH

A string-based coupled-cluster method of general excitation rank and with optimal scaling which accounts for special relativity within the four-component framework is presented. The method opens the way for the treatment of multi-reference problems through an active-space inspired single-reference based state-selective expansion of the model space. The evaluation of the coupled-cluster vector function is implemented by considering contractions of elementary second-quantized operators without setting up the amplitude equations explicitly. The capabilities of the new method are demonstrated in application to the electronic ground state of the bismuth monohydride molecule. In these calculations simulated multi-reference expansions with both doubles and triples excitations into the external space as well as the regular coupled-cluster hierarchy up to full quadruples excitations are compared. The importance of atomic outer core-correlation for obtaining accurate results is shown. Comparison to the non-relativistic framework is performed throughout to illustrate the additional work of the transition to the four-component relativistic framework both in implementation and application. Furthermore, an evaluation of the highest order scaling for general-order expansions is presented.     

Relativistic spin-orbit effects on hyperfine coupling tensors by density-functional theory

A second-order perturbation theory treatment of spin-orbit corrections to hyperfine coupling tensors has been implemented within a density-functional framework. The method uses the all-electron atomic mean-field approximation and/or spin-orbit pseudopotentials in incorporating one- and two-electron spin-orbit interaction within a first-principles framework. Validation of the approach on a set of main-group radicals and transition metal complexes indicates good agreement between all-electron and pseudopotential results for hyperfine coupling constants of the lighter nuclei in the system, except for cases in which scalar relativistic effects become important. The nonrelativistic Fermi contact part of the isotropic hyperfine coupling constants is not always accurately reproduced by the exchange-correlation functionals employed, particularly for the triplet and pi-type doublet radicals in the present work. For this reason, ab initio coupled-cluster singles and doubles with perturbative triples results for the first-order contributions have been combined in the validation calculations with the density-functional results for the second-order spin-orbit contributions. In the cases where spin-orbit corrections are of significant magnitude relative to the nonrelativistic first-order terms, they improve the agreement with experiment. Antisymmetric contributions to the hyperfine tensor arise from the spin-orbit contributions and are discussed for the IO2 radical, whereas rovibrational effects have been evaluated for RhC, NBr, and NI.     

Ab initio potential energy surface and spectrum of the B(3Pi) state of the HeI2 complex

The three-dimensional interaction potential for I2(B 3Pi0u+)+He is computed using accurate ab initio methods and a large basis set. Scalar relativistic effects are accounted for by large-core relativistic pseudopotentials for the iodine atoms. Using multireference configuration interaction calculations with subsequent treatment of spin-orbit coupling, it is shown for linear and perpendicular structures of the complex that the interaction potential for I2(B 3Pi0u+)+He is very well approximated by the average of the 3A' and 3A" interaction potentials obtained without spin-orbit coupling. The three-dimensional 3A' and 3A" interaction potentials are computed at the unrestricted open-shell coupled-cluster level of theory using large basis sets. Bound state calculations based on the averaged surface are carried out and binding energies, vibrationally averaged structures, and frequencies are determined. These results are found to be in excellent accord with recent experimental measurements from laser-induced fluorescence and action spectra of HeI2. Furthermore, in combination with a recent X-state potential, the spectral blueshift is obtained and compared with available experimental values.     

Third-order Douglas-Kroll relativistic coupled-cluster theory through connected single, double, triple, and quadruple substitutions: applications to diatomic and triatomic hydrides

Coupled-cluster methods including through and up to the connected single, double, triple, and quadruple substitutions have been derived and implemented automatically for sequential and parallel executions by an algebraic and symbolic manipulation program TCE (TENSOR CONTRACTION ENGINE) for use in conjunction with a one-component third-order Douglas-Kroll approximation for relativistic corrections. A combination of the converging electron-correlation methods, the accurate relativistic reference wave functions, and the use of systematic basis sets tailored to the relativistic approximation has been shown to predict the experimental singlet-triplet separations within 0.02 eV (0.5 kcal/mol) for five triatomic hydrides (CH2, NH2+, SiH2, PH2+, and AsH2+), the experimental bond lengths (re or r0) within 0.002 angstroms, rotational constants (Be or B0) within 0.02 cm(-1), vibration-rotation constants (alphae) within 0.01 cm(-1), centrifugal distortion constants (De) within 2%, harmonic vibration frequencies (omegae) within 8 cm(-1) (0.4%), anharmonic vibrational constants (xomegae) within 2 cm(-1), and dissociation energies (D0(0)) within 0.02 eV (0.4 kcal/mol) for twenty diatomic hydrides (BH, CH, NH, OH, FH, AlH, SiH, PH, SH, ClH, GaH, GeH, AsH, SeH, BrH, InH, SnH, SbH, TeH, and IH) containing main-group elements across the second through fifth rows of the periodic table. In these calculations, spin-orbit effects on dissociation energies, which were assumed to be additive, were estimated from the measured spin-orbit coupling constants of atoms and diatomic molecules, and an electronic energy in the complete-basis-set, complete-electron-correlation limit has been extrapolated in two ways to verify the robustness of the results: One assuming Gaussian-exponential dependence of total energies on double through quadruple zeta basis sets and the other assuming n(-3) dependence of correlation energies on double through quintuple zeta basis sets.     

Heats of formation of MHxCly (M = Si, P, As, Sb) compounds and main group fluorides from high level electronic structure calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for the MH(x)Cl(y) compounds (M = Si, P, As, and Sb) and for a number of trivalent, tetravalent, and pentavalent fluorides (SbF(3), BiF(3), GeF(4), SnF(4), PbF(4), AsF(5), SbF(5)) from coupled cluster theory (CCSD(T)) calculations using correlation consistent basis sets and extrapolation to the complete basis set limit. Small-core, relativistic effective core potentials were used for the heavier elements (Ge, As, Sn, Sb, Pb, and Bi), including correlation of the outer core electrons. Additional scalar relativistic (for the lighter elements) and atomic spin-orbit corrections are included in order to achieve near chemical accuracy of ±1.5 kcal/mol. Vibrational zero point energies were computed from scaled harmonic frequencies at the second order M?ller-Plesset perturbation theory (MP2) level where possible. Agreement between theory and the available experimental data is excellent. We present a revised heat of formation of the antimony atom in the gas phase. The calculated values will be of use in predicting the behavior of chemical vapor deposition systems.     

Efficient formulation and computer implementation of the active-space electron-attached and ionized equation-of-motion coupled-cluster methods

The efficient, general-purpose implementations of the active-space electron-attached (EA) and ionized (IP) equation-of-motion coupled-cluster (EOMCC) methods including up to 3p-2h and 3h-2p excitations, called EA-EOMCCSDt and IP-EOMCCSDt, respectively, are discussed. The details of the algorithm that enables one to achieve a high degree of code vectorization for the active-space methods and the factorized forms of the EA- and IP-EOMCCSDt equations that maximize the benefits of using active orbitals in the process of selecting the dominant 3p-2h and 3h-2p excitations are presented. The results of benchmark calculations for the low-lying doublet and quartet states of the CH and SH radicals reveal that the active-space EA-EOMCCSDt and IP-EOMCCSDt methods are capable of producing results for the electronic excitations in open-shell systems that match the high accuracy of EA- and IP-EOMCC calculations with a full treatment of 3p-2h and 3h-2p excitations, even when the excited states of interest display a manifestly multideterminantal nature, with the costs that can be on the same order of those characterizing the basic EOMCC singles and doubles approach.     

On the spectroscopic and thermochemical properties of ClO, BrO, IO, and their anions

A coupled cluster composite approach has been used to accurately determine the spectroscopic constants, bond dissociation energies, and heats of formation for the X1(2)II(3/2) states of the halogen oxides ClO, BrO, and IO, as well as their negative ions ClO-, BrO-, and IO-. After determining the frozen core, complete basis set (CBS) limit CCSD(T) values, corrections were added for core-valence correlation, relativistic effects (scalar and spin-orbit), the pseudopotential approximation (BrO and IO), iterative connected triple excitations (CCSDT), and iterative quadruples (CCSDTQ). The final ab initio equilibrium bond lengths and harmonic frequencies for ClO and BrO differ from their accurate experimental values by an average of just 0.0005 A and 0.8 cm-1, respectively. The bond length of IO is overestimated by 0.0047 A, presumably due to an underestimation of molecular spin-orbit coupling effects. Spectroscopic constants for the spin-orbit excited X2(2)III(1/2) states are also reported for each species. The predicted bond lengths and harmonic frequencies for the closed-shell anions are expected to be accurate to within about 0.001 A and 2 cm-1, respectively. The dissociation energies of the radicals have been determined by both direct calculation and through use of negative ion thermochemical cycles, which made use of a small amount of accurate experimental data. The resulting values of D0, 63.5, 55.8, and 54.2 kcal/mol for ClO, BrO, and IO, respectively, are the most accurate ab initio values to date, and those for ClO and BrO differ from their experimental values by just 0.1 kcal/mol. These dissociation energies lead to heats of formation, DeltaH(f) (298 K), of 24.2 +/- 0.3, 29.6 +/- 0.4, and 29.9 +/- 0.6 kcal/mol for ClO, BrO, and IO, respectively. Also, the final calculated electron affinities are all within 0.2 kcal/mol of their experimental values. Improved pseudopotential parameters for the iodine atom are also reported, together with revised correlation consistent basis sets for this atom.     

Perturbative treatment of the electron-correlation contribution to the diagonal Born-Oppenheimer correction

A perturbative scheme for the treatment of electron-correlation effects on the diagonal Born-Oppenheimer correction (DBOC) is suggested. Utilizing the usual Moller-Plesset partitioning of the Hamiltonian formulas for first and second orders (termed as MP1 and MP2) are obtained by expanding the wave function in the corresponding coupled-cluster expressions for the DBOC[J. Gauss et al., J. Chem. Phys. 125, 144111 (2006)]. The obtained expressions are recast in terms of one- and two-particle density matrices in order to take advantage of existing analytic second-derivative implementations for many-body methods. Test calculations show that both MP1 and MP2 recover large fractions (on average 90% and 95%, respectively) of the coupled-cluster singles and doubles (CCSD) electron-correlation corrections to the DBOC and thus render the suggested MP treatments cost-effective (though still accurate) alternatives to high-level coupled cluster (CC) treatments. The applicability of the MP1 and MP2 schemes for treating DBOC is demonstrated in calculations for the atomization energies of benzene, naphthalene, anthracene, and tetracene. The corresponding corrections are surprisingly large (about 0.6 kJmol for benzene, 1.1 kJmol for naphthalene, 1.5 kJmol for anthracene, and 1.8 kJmol for tetracene) with the electron-correlation corrections reducing the corresponding Hartree-Fock self-consistent field values by 25%-30%.     

A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect

In the present work we have proposed an approximate time-dependent density-functional theory (TDDFT) formalism to deal with the influence of spin-orbit coupling effect on the excitation energies for closed-shell systems. In this formalism scalar relativistic TDDFT calculations are first performed to determine the lowest single-group excited states and the spin-orbit coupling operator is applied to these single-group excited states to obtain the excitation energies with spin-orbit coupling effects included. The computational effort of the present method is much smaller than that of the two-component TDDFT formalism and this method can be applied to medium-size systems containing heavy elements. The compositions of the double-group excited states in terms of single-group singlet and triplet excited states are obtained automatically from the calculations. The calculated excitation energies based on the present formalism show that this formalism affords reasonable excitation energies for transitions not involving 5p and 6p orbitals. For transitions involving 5p orbitals, one can still obtain acceptable results for excitations with a small truncation error, while the formalism will fail for transitions involving 6p orbitals, especially 6p1/2 spinors.     

Heats of formation of krypton fluorides and stability predictions for KrF4 and KrF6 from high level electronic structure calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for KrF+, KrF-, KrF2, KrF3+, KrF4, KrF5+, and KrF6 from coupled-cluster theory (CCSD(T)) calculations with effective core potential correlation-consistent basis sets for krypton. To achieve near chemical accuracy (+/-1 kcal/mol), three corrections were added to the complete basis set binding energies based on frozen core coupled-cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, and a correction for first-order atomic spin-orbit effects. Vibrational zero point energies were computed at the coupled-cluster level of theory. The calculated value for the heat of formation of KrF2 is in excellent agreement with the experimental value. Contrary to the analogous xenon fluorides, KrF2, KrF4, and KrF6 are predicted to be thermodynamically unstable with respect to loss of F2. An analysis of the energetics of KrF4 and KrF6 with respect to fluorine atom loss together with calculations of the transition states for the intramolecular loss of F2 show that fluorine atom loss is the limiting factor determining the kinetic stabilities of these molecules. Whereas KrF4 possesses a marginal energy barrier of 10 kcal/mol toward fluorine atom loss and might be stable at moderately low temperatures, the corresponding barrier in KrF6 is only 0.9 kcal/mol, suggesting that it could exist only at very low temperatures. Although the simultaneous reactions of either two or four fluorine atoms with KrF2 to give KrF4 or KrF6, respectively, are exothermic, they do not represent feasible synthetic approaches because the attack of the fluorine ligands of KrF2 by the fluorine atoms, resulting in F2 abstraction, is thermodynamically favored over oxidative fluorination of the krypton central atom. Therefore, KrF6 could exist only at very low temperatures, and even the preparation of KrF4 will be extremely difficult.