Property calculations using perturbed orbitals via state-specific multireference coupled-cluster and perturbation theories

Negative Ion Resonance States: The Fock-Space Coupled-Cluster Way

The negative ion resonance states, which are electron-molecule metastable compound states, play the most important role in free-electron controlled molecular reactions and low-energy free-electron-induced DNA damage. Their electronic structure is often only poorly described but crucial to an understanding of their reaction dynamics. One of the most important challenges to current electronic structure theory is the computation of negative ion resonance states. As a major step forward, coupled-cluster theories, which are well-known for their ability to produce the best approximate bound state electronic eigen solutions, are upgraded to offer the most accurate and effective approximations for negative ion resonance states. The existing Fock-space coupled-cluster (FSCC) and the equation-of-motion coupled-cluster (EOM-CC) approaches that compute bound states are redesigned for the direct and simultaneous determination of both the kinetic energy of the free electron at which the electron-molecule compound states are resonantly formed and the corresponding autodetachment decay rate of the electron from the metastable compound state. This Feature Article reviews the computation of negative ion resonances using the FSCC approach and, in passing, provides the highlights of the equivalent EOM-CC approach.

Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer

In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate formal properties of the used method, such as energy size consistency and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of what we have defined as a complete active space-external space gap describing the basis splitting between the complete active space and the external part generalizing the concept of a HOMO–LUMO gap. Furthermore, the behavior of the energy error for an optimal basis splitting, i.e., an active space choice minimizing the density matrix renormalization group-tailored coupled-cluster singles doubles error, is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose a complete active space. Moreover, the dependence of the ground-state energy error on the complete active space has been analyzed numerically in order to find an optimal split between the complete active space and external space by minimizing the density matrix renormalization group-tailored coupled-cluster error.

Extending spin-symmetry projected coupled-cluster to large model spaces using an iterative null-space projection technique

Recently, we introduced an orbital-invariant approximate coupled-cluster (CC) method in the spin-projection manifold. The multi-determinantal property of spin-projection means that the parametrization in the spin-extended CC (ECC) ansatz is nonorthogonal and overcomplete. Therefore, the linear dependencies must be removed by an orthogonalization procedure to obtain meaningful solutions. Multi-reference methods often achieve this by diagonalizing a metric of the equation system, but this is not feasible with ECC because of the enormous size of the metric, a consequence of the incomplete active space of the spin-projected Hartree-Fock reference. As a result, the applicability of ECC has been limited to small benchmark systems, for which the ansatz was shown to be superior to the configuration interaction and linearized approximations. In this article, we provide a solution to this problem that completely avoids the metric diagonalization by iteratively projecting out its null-space from the working equations. As the additional computational cost required for this iterative projection is only marginal, it greatly expands the application range of ECC. We demonstrate the potential of approximate ECC by studying the complete basis set limit of F₂ and transition metal complexes such as NiO, Mn₂ , and [Cu₂ O₂ ]²⁺ , which have all been hindered by the prohibitively large metric size. We also identify the potential inadequacy of the molecular orbitals given by spin-projected Hartree-Fock in some cases, and propose possible solutions. © 2018 Wiley Periodicals, Inc.

Orbital-invariant spin-extended approximate coupled-cluster for multi-reference systems

We present an approximate treatment of spin-extended coupled-cluster (ECC) based on the spin-projection of the broken-symmetry coupled-cluster (CC) ansatz. ECC completely eliminates the spin-contamination of unrestricted CC and is therefore expected to provide better descriptions of dynamical and static correlation effects, but introduces two distinct problems. The first issue is the emergence of non-terminating amplitude equations, which are caused by the de-excitation effects inherent in symmetry projection operators. In this study, we take a minimalist approach and truncate the Taylor series of the exponential ansatz at a certain order such that the approximation safely recovers the traditional CC without spin-projection. The second issue is that the nonlinear equations of ECC become underdetermined, although consistent, yielding an infinitude of solutions. This problem arises because of the redundancies in the excitation manifold, as is common in other multi-reference approaches. We remove the linear dependencies in ECC by employing an orthogonal projection manifold. We also propose an efficient solver for our method, in which the components are usually sparse but not diagonal-dominant. It is shown that our approach is rigorously orbital-invariant and provides more accurate results than its configuration interaction and linearized CC analogues for chemical systems.

Structure, stability, and properties of the trans peroxo nitrate radical: the importance of nondynamic correlation

We report a comparative single-reference and multireference coupled-cluster investigation on the structure, potential energy surface, and IR spectroscopic properties of the trans peroxo nitrate radical, one of the key intermediates in stratospheric NOX chemistry. The previous single-reference ab initio studies predicted an unbound structure for the trans peroxo nitrate radical. However, our Fock space multireference coupled-cluster calculation confirms a bound structure for the trans peroxo nitrate radical, in accordance with the experimental results reported earlier. Further, the analysis of the potential energy surface in FSMRCC method indicates a well-behaved minima, contrary to the shallow minima predicted by the single-reference coupled-cluster method. The harmonic force field analysis, of various possible isomers of peroxo nitrate also reveals that only the trans structure leads to the experimentally observed IR peak at 1840 cm(-1). The present study highlights the critical importance of nondynamic correlation in predicting the structure and properties of high-energy stratospheric NOx radicals.

Effect of Triples to Dipole Moments in Fock-Space Multireference Coupled Cluster Method

In this paper, we present the new implementation of partial triples for the dipole moment of doublet radicals in Lagrangian formulation of Fock-space multireference coupled cluster (Λ-FSMRCC) response method. We have implemented a specific scheme of noniterative triples, in addition to singles and doubles schemes, which accounts for the effects appearing at least at the third order in dipole moments. The method is applied to the ground states of OH, OOH, HCOO, CN, CH, and PO radicals.

Coupling term derivation and general implementation of state-specific multireference coupled cluster theories

Simple closed-form expressions are derived for the "same vacuum" renormalization terms that arise in state-specific multireference coupled cluster (MRCC) theories. Explicit equations are provided for these coupling terms through the triple excitation level of MRCC theory, and a general expression is included for arbitrary-order excitations. The first production-level code (PSIMRCC) for state-specific and rigorously size-extensive Mukherjee multireference coupled cluster singles and doubles (MkCCSD) computations has been written. This code is also capable of evaluating analogous Brillouin-Wigner multireference energies (BWCCSD), including a posteriori size-extensivity corrections. Using correlation-consistent basis sets (cc-pVXZ, X=D,T,Q), MkCCSD and BWCCSD were tested and compared on two classic multireference problems: (1) the dissociation potential curve of molecular fluorine (F(2)) and (2) the structure and vibrational frequencies of ozone. Comparison with experimental data shows that the Mukherjee method is generally superior to the Brillouin-Wigner theory in predicting energies, structures, and vibrational frequencies. Particularly accurate results for F(2) are obtained by applying the MkCCSD method with localized molecular orbitals. Although the MkCCSD theory greatly improves upon single-reference CCSD for the geometric parameters and a(1) vibrational frequencies of ozone, the antisymmetric stretching frequency omega(3)(b(2)) remains pathological and cannot be properly treated without the inclusion of connected triple excitations. Finally, preliminary multireference MkCCSD results are reported for the singlet-triplet splittings in ortho-, meta-, and para-benzyne, coming within 1.5 kcal mol(-1) of experiment in all cases.

Multiconfiguration perturbation theory: Size consistency at second order

A modified version of a previously elaborated multiconfiguration perturbation theory (MCPT) [Rolik et al. J. Chem. Phys. 119, 1922 (2003)] is presented. In the modified formulation size consistency is ensured at second order in energy, by omitting projectors from the zero order Hamiltonian operator. This MCPT formulation is abbreviated as SC2-MCPT (size consistent at second order). To ensure proper separability, we also require that energy denominators are constructed as differences of some one-particle energies. A similar choice for energy denominators also renders the well-known multireference Moller-Plesset (MRMP) energy size consistent at second order. The same thing applies to the related multireference perturbation theory by Witek, Nakano, and Hirao.

Multireference State-Specific Coupled-Cluster Theory and Multiconfigurationality Index. BH Dissociation

Multiconfigurationality index calculated for the coupled-cluster wave function based on an algorithm developed using a computer-aided generation approach is applied to analyze the multireference state-specific coupled-cluster method with the CAS reference (i.e. the so called the CAS(n,m)CCSD approach). The numerical results concern dissociation of the BH molecule where at larger displacement from the equilibrium significant quasi-degeneracy arises. The analysis shows that the CAS(n,m)CCSD approach performs very well in such a situation.

A state-specific approach to multireference coupled electron-pair approximation like methods: Development and applications

The traditional multireference (MR) coupled-cluster (CC) methods based on the effective Hamiltonian are often beset by the problem of intruder states, and are not suitable for studying potential energy surface (PES) involving real or avoided curve crossing. State-specific MR-based approaches obviate this limitation. The state-specific MRCC (SS-MRCC) method developed some years ago can handle quasidegeneracy of varying degrees over a wide range of PES, including regions of real or avoided curve-crossing. Motivated by its success, we have suggested and explored in this paper a suite of physically motivated coupled electron-pair approximations (SS-MRCEPA) like methods, which are designed to capture the essential strength of the parent SS-MRCC method without significant sacrificing its accuracy. These SS-MRCEPA theories, like their CC counterparts, are based on complete active space, treat all the reference functions on the same footing and provide a description of potentially uniform precision of PES of states with varying MR character. The combining coefficients of the reference functions are self-consistently determined along with the cluster amplitudes themselves. The newly developed SS-MRCEPA methods are size-extensive, and are also size-consistent with localized orbitals. Among the various versions, there are two which are invariant with respect to the restricted rotations among doubly occupied and active orbitals separately. Similarity of performance of this latter and the noninvariant versions at the crossing points of the degenerate orbitals imply that the all the methods presented are rather robust with respect to the rotations among degenerate orbitals. Illustrative numerical applications are presented for PES of the ground state of a number of difficult test cases such as the model H4, H(8) problems, the insertion of Be into H(2), and Li(2), where intruders exist and for a state of a molecule such as CH(2), with pronounced MR character. Results obtained with SS-MRCEPA methods are found to be comparable in accuracy to the parent SS-MRCC and FCI/large scale CI results throughout the PES, which indicates the efficacy of our SS-MRCEPA methods over a wide range of geometries, despite their neglect of a host of complicated nonlinear terms, even when the traditional MR-based methods based on effective Hamiltonians fail due to intruders.