A relativistic 4-component general-order multi-reference coupled cluster method: initial implementation and application to HBr

We present the initial implementation of a determinant-based general-order coupled cluster method which fully accounts for relativistic effects within the four-component framework. The method opens the way for the treatment of multi-reference problems through a state-selective expansion of the model space. The evaluation of the coupled cluster vector function is carried out via relativistic configuration interaction expansions. The implementation is based on a large-scale configuration interaction technique, which may efficiently treat long determinant expansions of more than 10⁸ terms. We demonstrate the capabilities of the new method in calculations of complete potential energy curves of the HBr molecule. The inclusion of spin–orbit interaction and higher excitations than coupled cluster double excitations, either by multi-reference model spaces or the inclusion of full iterative triple excitations, lead to highly accurate results for spectral constants of HBr. An erratum to this article can be found at     

Ab initio computed diabatic potential energy surfaces of OH-HCl

The two four-dimensional diabatic potential energy surfaces (DPESs) for OH-HCl are computed that correlate with the twofold degenerate (2)Pi ground state of the free OH radical. About 20 000 points on the surface are obtained by the ab initio coupled-cluster and multi-reference configuration interaction methods. Analytic forms for the diabatic potential energy surfaces are derived as expansions in complete sets of orthogonal functions depending on the three intermolecular angles. The numeric computation of the angular expansion coefficients is discussed. The distance-dependence of the angular coefficients is represented by the reproducing kernel Hilbert space method. It is checked that both diabatic potentials converge for large intermolecular separations to the values computed directly from the electrostatic multipole expansion. The final DPESs are discussed and illustrated by some physically meaningful one- and two-dimensional cuts through them.     

Automated generation of coupled-cluster diagrams: implementation in the multireference state-specific coupled-cluster approach with the complete-active-space reference

An algorithm for generation of the spin-orbital diagrammatic representation, the corresponding algebraical formulas, and the computer code of the coupled-cluster (CC) method with an arbitrary level of the electronic excitations has been developed. The method was implemented in the general case as well as for specific application in the state-specific multireference coupled-cluster theory (SSMRCC) based on the concept of a "formal reference state." The algorithm was tested in SSMRCC calculations describing dissociation of a single bond and in calculations describing simultaneous dissociation of two single bonds--the problem requiring up to six-particle excitations in the CC operator.     

Calculation of excited-state properties using general coupled-cluster and configuration-interaction models

Using string-based algorithms excitation energies and analytic first derivatives for excited states have been implemented for general coupled-cluster (CC) models within CC linear-response (LR) theory which is equivalent to the equation-of-motion (EOM) CC approach for these quantities. Transition moments between the ground and excited states are also considered in the framework of linear-response theory. The presented procedures are applicable to both single-reference-type and multireference-type CC wave functions independently of the excitation manifold constituting the cluster operator and the space in which the effective Hamiltonian is diagonalized. The performance of different LR-CC/EOM-CC and configuration-interaction approaches for excited states is compared. The effect of higher excitations on excited-state properties is demonstrated in benchmark calculations for NH(2) and NH(3). As a first application, the stationary points of the S(1) surface of acetylene are characterized by high-accuracy calculations.     

Analytic second derivatives for general coupled-cluster and configuration-interaction models

Analytic second derivatives of energy for general coupled-cluster (CC) and configuration-interaction (CI) methods have been implemented using string-based many-body algorithms. Wave functions truncated at an arbitrary excitation level are considered. The presented method is applied to the calculation of CC and CI harmonic frequencies and nuclear magnetic resonance chemical shifts up to the full CI level for some selected systems. The present benchmarks underline the importance of higher excitations in high-accuracy calculations.     

A Hilbert space multi-reference coupled-cluster study of the H4 model system

Summary Employing the Hilbert space ansatz, a fully quadratic coupled-cluster method with a multidimensional reference space is applied to a DZP basis study of the model system, H₄. The reference space is described by two to four configurations at the level of single and double excitations, and single and double excitation operators are included in the expansions for the cluster and wave operator through quadratic terms. The performance of quadratic MRCCSD is investigated for the ground and three excited states of the H₄ system consisting of two stretched hydrogen molecules in a trapezoidal configuration where the degree of quasidegeneracy is varied from a nondegenerate situation to a completely degenerate one. Compared to full CI, in the highly degenerate region, the MRCCSD works quite well. In less degenerate regions, the accuracy is less satisfactory.     

Superior performance of Mukherjee’s state-specific multi-reference coupled-cluster theory at the singles and doubles truncation scheme with localized active orbitals

The state-specific multi-reference coupled-cluster (SS-MRCC) theory of Mukherjee et al., in its singles and doubles truncation scheme (SS-MRCCSD), misses important couplings between the virtual functions reached by single and double excitations from different model functions. Since the SS-MRCC theory is not invariant with respect to the transformations among the active orbitals, the results are dependent on the active orbitals chosen. We demonstrate in this paper with results for potential energy curves for several example molecules involving single and multiple bond dissociation that the performance of SS-MRCCSD is significantly improved if localized active orbitals are used. The improvement is remarkable both in terms of the non-parallelity error and the magnitude of correlation energy recovered vis-a-vis the full configuration interaction results with the same basis set. The results bolster our claim that SS-MRCCSD with localized orbitals is an accurate general theory for potential energy surfaces.     

Approximate treatment of higher excitations in coupled-cluster theory

The possibilities for the approximate treatment of higher excitations in coupled-cluster (CC) theory are discussed. Potential routes for the generalization of corresponding approximations to lower-level CC methods are analyzed for higher excitations. A general string-based algorithm is presented for the evaluation of the special contractions appearing in the equations specific to those approximate CC models. It is demonstrated that several iterative and noniterative approximations to higher excitations can be efficiently implemented with the aid of our algorithm and that the coding effort is mostly reduced to the generation of the corresponding formulas. The performance of the proposed and implemented methods for total energies is assessed with special regard to quadruple and pentuple excitations. The applicability of our approach is illustrated by benchmark calculations for the butadiene molecule. Our results demonstrate that the proposed algorithm enables us to consider the effect of quadruple excitations for molecular systems consisting of up to 10-12 atoms.     

Coupled-cluster methods including noniterative corrections for quadruple excitations

A new method is presented for treating the effects of quadruple excitations in coupled-cluster theory. In the approach, quadruple excitation contributions are computed from a formula based on a non-Hermitian perturbation theory analogous to that used previously to justify the usual noniterative triples correction used in the coupled cluster singles and doubles method with a perturbative treatment of the triple excitations (CCSD(T)). The method discussed in this paper plays a parallel role in improving energies obtained with the full coupled-cluster singles, doubles, and triples method (CCSDT) by adding a perturbative treatment of the quadruple excitations (CCSDT(Q)). The method is tested for an extensive set of examples, and is shown to provide total energies that compare favorably with those obtained with the full singles, doubles, triples, and quadruples (CCSDTQ) method.     

General implementation of the relativistic coupled-cluster method

We report the development of a general order relativistic coupled-cluster (CC) code. Our implementation is based on Kramers-paired molecular spinors, utilizes double group symmetry, and is applicable with the full Dirac-Coulomb and several approximate relativistic Hamiltonians. The available methods include iterative and perturbative single-reference CC approaches with arbitrary excitations as well as a state-selective multi-reference CC ansatz. To illustrate the performance of the new code, benchmark calculations have been performed for the total energies, bond lengths, and vibrational frequencies of the monoxides of Group IVa elements. The trends due to the simultaneous inclusion of relativity as well as higher-order electron correlation effects are analyzed. The newly developed code significantly widens the scope of the ab initio relativistic calculations, for both molecules and atoms alike, surpassing the accuracy and reliability of the currently available implementations in the literature.     

The coupled-cluster method in high sectors of the Fock space

The coupled-cluster method is applied to high sectors of the Fock space (up to five electrons outside a closed shell). Ionization potentials and excitation energies are calculated for atomic N, O, F, P, and S and their ions. The effect of virtual triple excitations is included, exactly or at several levels of appoximation. These terms are important for accurate results in the high sectors. © 1995 John Wiley & Sons, Inc.     

Universal state-selective corrections to multi-reference coupled-cluster theories with single and double excitations

The recently proposed universal state-selective (USS) corrections [K. Kowalski, J. Chem. Phys. 134, 194107 (2011)] to approximate multi-reference coupled-cluster (MRCC) energies can be commonly applied to any type of MRCC theory based on the Jeziorski-Monkhorst [B. Jeziorski and H. J. Monkhorst, Phys. Rev. A 24, 1668 (1981)] exponential ansatz. In this paper we report on the performance of a simple USS correction to the Brillouin-Wigner and Mukherjee's MRCC approaches employing single and double excitations (USS-BW-MRCCSD and USS-Mk-MRCCSD). It is shown that the USS-BW-MRCCSD correction, which employs the manifold of single and double excitations, can be related to a posteriori corrections utilized in routine BW-MRCCSD calculations. In several benchmark calculations we compare the USS-BW-MRCCSD and USS-Mk-MRCCSD results with the results obtained with the full configuration interaction method.     

Spin-flip configuration interaction: an electronic structure model that is both variational and size-consistent

A new formulation of the configuration interaction (CI) method is presented. It is based on the recently introduced spin-flip (SF) approach. SF-CI target states are described as spin-flipping excitations from the reference Hartree–Fock high-spin, e.g., M_s=1 (|), determinant. The resulting model is both variational and size-consistent. Moreover, the SF-CI model can describe within a single-reference formalism some inherently multi-reference situations, such as single bond-breaking and diradicals. Initial benchmarks for the SF-CI model with single and double substitutions (SF-CISD) are presented.     

pCCSD: parameterized coupled-cluster theory with single and double excitations

The primary characteristics of single reference coupled-cluster (CC) theory are size-extensivity and size-consistency, invariance under orbital rotations of the occupied or virtual space, the exactness of CC theory for N electron systems when the cluster operator is truncated to N-tuple excitations, and the relative insensitivity of CC theory to the choice of the reference determinant. In this work, we propose a continuous class of methods which display the desirable features of the coupled-cluster approach with single and double excitations (CCSD). These methods are closely related to the CCSD method itself and are inspired by the coupled electron pair approximation (CEPA). It is demonstrated that one can systematically improve upon CCSD and obtain geometries, harmonic vibrational frequencies, and total energies from a parameterized version of CCSD or pCCSD(α,β) by selecting a specific member from this continuous family of approaches. In particular, one finds that one such approach, the pCCSD(-1,1) method, is a significant improvement over CCSD for the calculation of equilibrium structures and harmonic frequencies. Moreover, this method behaves surprisingly well in the calculation of potential energy surfaces for single bond dissociation. It appears that this methodology has significant promise for chemical applications and may be particularly useful in applications to larger molecules within the framework of a high accuracy local correlation approach.     

Relativistic calculations on thallium hydride

The spectroscopic parameters of the ground state of thallium hydride are obtained using the four-component relativistic Dirac–Coulomb–Gaunt–coupled-cluster wave function with single and double excitations and an estimated triples correction method. Core correlation effects make the bond weaker but have little effect on the bond length. Inclusion of the distance dependence of the Gaunt part of the two-electron Breit interaction has an opposite but smaller effect on these properties. Received: 8 September 2000 / Revised version: 5 October 2000 / Published online: 21 December 2000     

Calculation of current densities using gauge-including atomic orbitals

A method for calculating the various components of the magnetically induced current-density tensor using gauge-including atomic orbitals is described. The method is formulated in the framework of analytical derivative theory, thus enabling implementation at the Hartree-Fock self-consistent-field (HF-SCF) as well as at electron-correlated levels. First-order induced current densities have been computed up to the coupled-cluster singles and doubles level (CCSD) augmented by a perturbative treatment of triple excitations [CCSD(T)] for carbon dioxide and benzene and up to the full coupled-cluster singles, doubles, and triples (CCSDT) level in the case of ozone. The applicability of the gauge including magnetically induced current method to larger molecules is demonstrated by computing first-order current densities for porphin and hexabenzocoronene at the HF-SCF and density-functional theory level. Furthermore, a scheme for obtaining quantitative values for the induced currents in a molecule via numerical integration over the current flow is presented. For benzene, a perpendicular magnetic field induces a (field dependent) ring current of 12.8 nA T(-1) at the HF-SCF level using a triple-zeta basis set augmented with polarization functions (TZP). At the CCSD(T)/TZP level the induced current was found to be 11.4 nA T(-1). Gauge invariance and its relation to charge-current conservation is discussed.     

Coupled-cluster with active space selected higher amplitudes: performance of seminatural orbitals for ground and excited state calculations

The active space approach for coupled-cluster models is generalized using the general active space concept and implemented in a string-based general coupled-cluster code. Particular attention is devoted to the choice of orbitals on which the subspace division is based. Seminatural orbitals are proposed for that purpose. These orbitals are obtained by diagonalizing only the hole-hole and particle-particle block of the one-electron density of a lower-order method. The seminatural orbitals are shown to be a good replacement for complete active space self-consistent field orbitals and avoid the ambiguities with respect to the reference determinant introduced by the latter orbitals. The seminatural orbitals also perform well in excited state calculations, including excited states with strong double excitation contributions, which usually are difficult to describe with standard coupled-cluster methods. A set of vertical excitation energies is obtained and benchmarked against full configuration interaction calculations, and alternative hierarchies of active space coupled-cluster models are proposed. As a simple application the spectroscopic constants of the C(2) B (1)Delta(g) and B(') (1)Sigma(g) (+) states are calculated using active space coupled-cluster methods and basis sets up to quadruple-zeta quality in connection with extrapolation and additivity schemes.     

Application of the MC SCF method to the π → π* excitation energies of ethylene

The MC SCF method is employed to calculate the N → T and N → V π → π^* vertical excitation energies of ethylene. To obtain accurate excitation energies it is found to be necessary to utilize an expanded valence space containing two π and two π^* orbitals. Relatively small MC SCF calculations, allowing at most one-electron excitations from the sigma space, are found to yield excitation energies and spatial extents of the excited states in excellent agreement with the predictions of large multi-reference or iterative-natural-orbital CI calculations. These results show that within an MC SCF framework σ-σ correlation is unimportant for describing the π → π^* processes. We also conclude that the neglect of the effects of unlinked cluster terms in some of the CI calculations may have introduced small, but important, errors in the excitation energies and predictions of the spatial extent of the V state.