High-order electron-correlation methods with scalar relativistic and spin-orbit corrections

An assortment of computer-generated, parallel-executable programs of ab initio electron-correlation methods has been fitted with the ability to use relativistic reference wave functions. This has been done on the basis of scalar relativistic and spin-orbit effective potentials and by allowing the computer-generated programs to handle complex-valued, spinless orbitals determined by these potentials. The electron-correlation methods that benefit from this extension are high-order coupled-cluster methods (up to quadruple excitation operators) for closed- and open-shell species, coupled-cluster methods for excited and ionized states (up to quadruples), second-order perturbation corrections to coupled-cluster methods (up to triples), high-order perturbation corrections to configuration-interaction singles, and active-space (multireference) coupled-cluster methods for the ground, excited, and ionized states (up to active-space quadruples). A subset of these methods is used jointly such that the dynamical correlation energies and scalar relativistic effects are computed by a lower-order electron-correlation method with more extensive basis sets and all-electron relativistic treatment, whereas the nondynamical correlation energies and spin-orbit effects are treated by a higher-order electron-correlation method with smaller basis sets and relativistic effective potentials. The authors demonstrate the utility and efficiency of this composite scheme in chemical simulation wherein the consideration of spin-orbit effects is essential: ionization energies of rare gases, spectroscopic constants of protonated rare gases, and photoelectron spectra of hydrogen halides.     

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.     

Correction to constrained coupled cluster doubles models based on the second coupled cluster central moment

We develop a correction for the coupled cluster version of the perfect pairing (PP) model. The correction is based on finding modified values of the PP amplitudes such that the second coupled cluster central moment defined in the space of all valence single and double substitutions vanishes and, subject to this constraint, minimizing the deviation between the modified and unmodified PP amplitudes with respect to a chosen metric. We discuss how this correction can be generalized to other constrained doubles models, such as local correlation and active-space models. While the correction is not strictly size consistent and retains some of the deficiencies of the PP model, numerical results indicate that much of the missing active-space coupled cluster singles and doubles correlation energy is recovered.     

Coupled‐cluster theory for the polarizable continuum model. III. A response theory for molecules in solution

A coupled‐cluster (CC) response functions theory for molecular solutes described with the framework of the polarizable continuum model (PCM) is presented. The theory is an extension to the dynamical molecular properties of the PCM‐CC analytic derivatives recently proposed for the calculation of static molecular properties (Cammi, Jr Chem Phys 2009, 131, 164104). The theory is presented for linear and quadratic response functions, and the operative expressions of these response functions can accurately account for the nonequilibrium solvation effects. The excitation energies and transition moments of the solvated chromophores have been determined from the linear response functions. Accurate expressions for gradients of excitation energies for the evaluation of the excited state properties have been also discussed. © 2012 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Can large active-space CASSCF calculation make sense to the reaction analysis of iron complex? A benchmark study of methane oxidation reaction by FeO+

A methane oxidation reaction by FeO⁺ cation was theoretically investigated based on the density functional theory (DFT) and the complete active-space self-consistent field (CASSCF) method as well as the coupled-cluster singles, doubles, and perturbative triples (CCSD(T)) to explore the active-space dependency to computational analyses in such strongly correlated reaction systems. A small active-space CASSCF(5e in 5o) calculation, which only includes five 3d orbitals of the Fe atom in the active-space, showed remarkable difference both in energy and geometry compared to those computed by the DFT and CCSD(T) methods. Interestingly, a large active-space CASSCF(17e in 17o) calculation, which includes almost all the valence orbitals gives a qualitative agreement with either the DFT or the CCSD(T) results in the first half part of the reaction, although it varies from them in the latter half part. Therefore, it is indicated that the active-space dependency is serious in some part of the reaction and the small active-space CASSCF might lead a wrong discussion. We further investigated the optimized geometry of the intermediate complex with the small and the large active-space CASSCF methods as well as the CCSD(T) method, and found that the CASSCF(5e in 5o)-optimized geometry is considerably different from the others. In consequence, a small active-space CASSCF/CASPT2 calculation does not really work for such a strongly correlated reaction system even qualitatively, and a sophisticated assessment using the large active-space CASSCF/CASPT2 method will be indispensable. © 2018 Wiley Periodicals, Inc.     

Thermodynamics of titanium and vanadium reduction in non-aqueous environment calculated at various levels of theory

Reduction of titanium and vanadium compounds is a process accompanying the activation of coordinative olefin polymerization catalysts. Four density functional theory (DFT) functionals, coupled cluster with single, double, and perturbative triple excitations method CCSD(T) as well as complete active-space second-order perturbation theory method CASPT2 with a complete active-space self-consistent field CASSCF reference wave function were applied to investigate the thermodynamics of titanium and vanadium reduction. The performance of these theoretical methods was assessed and compared with experimental values. The calculations indicate that vanadium(IV) chloride is more easily reduced by trimethylaluminum than the corresponding titanium compound; the energies of reaction calculated at the CCSD(T) level are equal -57.21 and -33.10 kcal/mol, respectively. The calculations deal with the redox reactions of metal chlorides in the gas phase, rather than solvated ions in the aqueous solution. This approach may be more appropriate for olefin polymerization, usually carried out in nonpolar solvents.     

Convergence of triples energy in CCSD(T) and CC3 calculations with correlation-consistent basis sets

The convergence behavior of connected triples correlation energy in CCSD(T) and CC3 calculations with (aug-)cc-pVXZ basis sets has been accurately described in terms of a power law of the type E(X)=E(infinity)+AX(-4). Calculations ranging from double-Z through septuple-Z attest the validity of this X(-4) convergence model. Extrapolations generated from (X-1,X)-Z calculations yield energies of nearly (X+1)-Z quality. Typically, the fraction of triples correlation energy recovered is 0.92+/-0.05 in (D,T) extrapolations; 0.98+/-0.01 in (T,Q) extrapolations; 1.0002+/-0.0012 in (Q,5) extrapolations; and 0.9995+/-0.0005 in (5,6) extrapolations.     

Orbital-optimized coupled-cluster theory does not reproduce the full configuration-interaction limit

It is shown that due to the mixing of the usual projection approach of coupled cluster with variational orbital optimization, orbital-optimized coupled cluster (OCC) fails to reproduce the full configuration-interaction (full CI) limit when the cluster operator becomes complete. It is pointed out that the fulfillment of the projected singles equations, which define the orbital gradient in Brueckner coupled cluster (BCC), is mandatory for a correct behavior. As numerical examples we present general OCC and BCC calculations up to the full CI limit on CH(2) and an active-space model of ozone. The observed deviations of OCC from full CI are of the order of the correlation error obtained in calculations with up to quadruples excitations. Thus the failure of OCC may be considered tolerable in more approximate calculations but clearly prohibitive for any benchmark application. For applications to active-space models a hybrid approach for OCC is suggested in which for active particle-hole rotations the Brueckner orbital gradient is employed, whereas for the remaining orbital rotations the variational orbital gradient is retained.     

Application of renormalized coupled-cluster methods to potential function of water

The goal of this paper is to examine the performance of the conventional and renormalized single-reference coupled-cluster (CC) methods in calculations of the potential energy surface of the water molecule. A comparison with the results of the internally contracted multi-reference configuration interaction calculations including the quasi-degenerate Davidson correction (MRCI(Q)) and the spectroscopically accurate potential energy surface of water resulting from the use of the energy switching (ES) approach indicates that the relatively inexpensive completely renormalized (CR) CC methods with singles (S), doubles (D), and a non-iterative treatment of triples (T) or triples and quadruples (TQ), such as CR-CCSD(T), CR-CCSD(TQ), and the recently developed rigorously size extensive extension of CR-CCSD(T), termed CR-CC(2,3), provide substantial improvements in the results of conventional CCSD(T) and CCSD(TQ) calculations at larger internuclear separations. It is shown that the CR-CC(2,3) results corrected for the effect of quadruply excited clusters through the CR-CC(2,3)+Q approach can compete with the highly accurate MRCI(Q) data. The excellent agreement between the CR-CC(2,3)+Q and MRCI(Q) results suggests ways of improving the global potential energy surface of water resulting from the use of the ES approach in the regions of intermediate bond stretches and intermediate energies connecting the region of the global minimum with the asymptotic regions. Contribution to the Mark S. Gordon 65th Birthday Festschrift Issue.     

An overview of coupled cluster theory and its applications in physics

Summary What has since become known as the normal coupled cluster method (NCCM) was invented about thirty years ago to calculate ground-state energies of closed-shell atomic nuclei. Coupled cluster (CC) techniques have since been developed to calculate excited states, energies of open-shell systems, density matrices and hence other properties, sum rules, and the sub-sum-rules that follow from imbedding linear response theory within the NCCM. Further extensions deal both with systems at nonzero temperature and with general dynamical behaviour. More recently, a new version of CC theory, the so-called extended coupled cluster method (ECCM) has been introduced. It has the potential to describe such global phenomena as phase transitions, spontaneous symmetry breaking, states of topological excitation, and nonequilibrium behaviour. CC techniques are now widely recognized as providing one of the most universally applicable, most powerful, and most accurate of all microscopicab initio methods in quantum many-body theory. The number of successful applications within physics is now impressively large. In most such cases the numerical results are either the best or among the best available. A typical case is the electron gas, where the CC results for the correlation energy agree over the entire metallic density range to within less than 1 millihartree (or <1%) with the essentially exact Green's function Monte Carlo results. The role of CC theory within modern quantum many-body theory is first surveyed, by a comparison with other techniques. Its full range of applications in physics is then reviewed. These include problems in nuclear physics, both for finite nuclei and infinite nuclear matter; the electron gas; various integrable and nonintegrable models; various relativistic quantum field theories; and quantum spin chain and lattice models. Particular applications of the ECCM include the quantum hydrodynamics of a zero-temperature, strongly-interacting condensed Bose fluid; a charged impurity in a polarizable medium (e.g., positron annihilation in metals); and various anharmonic oscillator and spin systems.     

Assessing the performance of density functional theory for the electronic structure of metal-salens: the M06 suite of functionals and the d⁴-metals

We have systematically investigated the electronic structure of the d? metal-salen complexes including the Cr(II)-, Mn(III)-, Fe(IV)-, Mo(II)-, Tc(III)-, and Ru(IV)-salen complexes. Density functional theory (DFT) has been employed, using the BP86 and B3LYP functionals, and the entire M05 and M06 suites of meta-generalized gradient functionals. These results are compared to robust complete active-space self-consistent field (CASSCF) optimized geometries and complete active-space third-order perturbation theory (CASPT3) energies for the lowest singlet, triplet, and quintet states. Although the M06 and M06-L DFT functionals have been generally recommended for computations on complexes that contain main group and transition metals, none of the M0-functionals appear statistically better than the B3LYP functional for the computation of spin-state energy gaps. DFT- and CASSCF-optimized geometries normally agree to within 0.3 ? least root mean squared deviation (LRMSD) for the singlet and triplet structures and less than 0.1 ? LRMSD for the quintet structures. It can be concluded that no DFT functional tested here can be considered reliable over all metal-salen complexes and it is highly recommended that the accuracy of any given DFT functional should be assessed on a case-by-case basis.     

High-level ab initio studies of the electronic excited states of the hydroxyl radical and water-hydroxyl complex

The lowest-energy electronic transitions in the hydroxyl radical and the hydrogen bound complex H(2)O.HO have been studied using ab initio methods. We have used the complete active-space self-consistent field and multireference configuration interaction (MRCI) methods to calculate vertical excitation energies and oscillator strengths. At the MRCI level the lowest-lying (2)Sigma(+)<--(2)Pi electronic transition is redshifted by about 2500 cm(-1) upon formation of the H(2)O.HO complex. We propose that this transition could be used to identify the complex in the gas phase, which in turn could be used to examine the role of H(2)O.HO in atmospheric reactions.     

Energetics, transition states, and intrinsic reaction coordinates for reactions associated with O(3P) processing of hydrocarbon materials

Electronic structure calculations based on multiconfiguration wave functions are used to investigate a set of archetypal reactions relevant to O(3P) processing of hydrocarbon molecules and surfaces. These include O(3P) reactions with methane and ethane to give OH plus methyl or ethyl radicals, O(3P) + ethane to give CH3O + CH3, and secondary reactions of the OH product radical with ethane and the ethyl radical. Geometry optimization is carried out with CASSCF/cc-pVTZ for all reactions, and with CASPT2/cc-pVTZ for O(3P) + methane/ethane. Single-point energy corrections are applied with CASPT2, CASPT3, and MRCI + Q with the cc-pVTZ and cc-pVQZ basis sets, and the energies extrapolated to the complete basis set limit (CBL). Where comparison of computed barriers and energies of reaction with experiment is possible, the agreement is good to excellent. The best agreement (within experimental error) is found for MRCI + Q/CBL applied to O(3P) + methane. For the other reactions, CASPT2/CBL and MRCI + Q/CBL predictions differ from experiment by 1-5 kcal/mol for 0 K enthalpies of reaction, and are within 1 kcal/mol of the best-estimate experimental range of 0 K barriers for O(3P) + ethane and OH + ethane. The accuracy of MRCI + Q/CBL is limited mainly by the quality of the active space. CASPT2/CBL barriers are consistently lower than MRCI + Q/CBL barriers with identical reference spaces.     

Method of moments of coupled-cluster equations: a new formalism for designing accurate electronic structure methods for ground and excited states

The method of moments of coupled-cluster equations (MMCC), which provides a systematic way of improving the results of the standard coupled-cluster (CC) and equation-of-motion CC (EOMCC) calculations for the ground- and excited-state energies of atomic and molecular systems, is described. The MMCC theory and its generalized MMCC (GMMCC) extension that enables one to use the cluster operators resulting from the standard as well as nonstandard CC calculations, including those obtained with the extended CC (ECC) approaches, are based on rigorous mathematical relationships that define the many-body structure of the differences between the full configuration interaction (CI) and CC or EOMCC energies. These relationships can be used to design the noniterative corrections to the CC/EOMCC energies that work for chemical bond breaking and potential energy surfaces of excited electronic states, including excited states dominated by double excitations, where the standard single-reference CC/EOMCC methods fail. Several MMCC and GMMCC approximations are discussed, including the renormalized and completely renormalized CC/EOMCC methods for closed- and open-shell states, the quadratic MMCC approaches, the CI-corrected MMCC methods, and the GMMCC approaches for multiple bond breaking based on the ECC cluster amplitudes.     

C2 in a Box: Determining Its Intrinsic Bond Strength for the X1Σg+ Ground State

The intrinsic bond strength of C₂ in its ¹Σ_g⁺ ground state is determined from its stretching force constant utilizing MR‐CISD+Q(8,8), MR‐AQCC(8,8), and single‐determinant coupled cluster calculations with triple and quadruple excitations. By referencing the CC stretching force constant to its local counterparts of ethane, ethylene, and acetylene, an intrinsic bond strength half way between that of a double bond and a triple bond is obtained. Diabatic MR‐CISD+Q results do not change this. Confinement of C₂ and suitable reference molecules in a noble gas cage leads to compression, polarization, and charge transfer effects, which are quantified by the local CC stretching force constants and differences of correlated electron densities. These results are in line with two π bonds and a partial σ bond. Bond orders and bond dissociation energies of small hydrocarbons do not support quadruple bonding in C₂.     

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.     

Multireference self-consistent-field energies without the many-electron wave function through a variational low-rank two-electron reduced-density-matrix method

The variational two-electron reduced-density-matrix (2-RDM) method allows for the computation of accurate ground-state energies and 2-RDMs of atoms and molecules without the explicit construction of an N-electron wave function. While previous work on variational 2-RDM theory has focused on calculating full configuration-interaction energies, this work presents the first application toward approximating multiconfiguration self-consistent-field (MCSCF) energies via low-rank restrictions on the 1- and 2-RDMs. The 2-RDM method with two- or three-particle N-representability conditions reduces the exponential active-space scaling of MCSCF methods to a polynomial scaling. Because the first-order algorithm [Mazziotti, Phys. Rev. Lett. 93, 213001 (2004)] represents each form of the 1- and 2-RDMs by a matrix factorization, the RDMs are readily defined to have a low rank rather than a full rank by setting the matrix factors to be rectangular rather than square. Results for the potential energy surfaces of hydrogen fluoride, water, and the nitrogen molecule show that the low-rank 2-RDM method yields accurate approximations to the MCSCF energies. We also compute the energies along the symmetric stretch of a 20-atom hydrogen chain where traditional MCSCF calculations, requiring more than 17x10(9) determinants in the active space, could not be performed.