Coupled-cluster method tailored by configuration interaction

A method is presented which combines coupled cluster (CC) and configuration interaction (CI) to describe accurately potential-energy surfaces (PESs). We use the cluster amplitudes extracted from the complete active space CI calculation to manipulate nondynamic correlation to tailor a single reference CC theory (TCC). The dynamic correlation is then incorporated through the framework of the CC method. We illustrate the method by describing the PESs for HF, H2O, and N2 molecules which involve single, double, and triple bond-breaking processes. To the dissociation limit, this approach yields far more accurate PESs than those obtained from the conventional CC method and the additional computational cost is negligible compared with the CC calculation steps. We anticipate that TCC offers an effective and generally applicable approach for many problems.     

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.     

A Simple Approximate Perturbation Approach to Quasi-degenerate Systems

In this paper a simple approximate approach for the study of quasi-degenerate systems is presented in the frame of multireference perturbation theory. The formulation can be considered as an approximation of the quasi-degenerate perturbation theory (QDPT) with the simplification that only the state specific (diagonal) perturbation corrections to the energy have to be computed. The new approach is discussed and compared with previous QDPT formulations using the weakly avoided crossing model (for which new properties are here presented) and applied to the case of the neutral/ionic energy crossing in the LiF molecule.     

Reduced‐cost sparsity‐exploiting algorithm for solving coupled‐cluster equations

We present an algorithm for reducing the computational work involved in coupled‐cluster (CC) calculations by sparsifying the amplitude correction within a CC amplitude update procedure. We provide a theoretical justification for this approach, which is based on the convergence theory of inexact Newton iterations. We demonstrate by numerical examples that, in the simplest case of the CCD equations, we can sparsify the amplitude correction by setting, on average, roughly 90% nonzero elements to zeros without a major effect on the convergence of the inexact Newton iterations.     

Coupled-cluster theory in a projected atomic orbital basis

We present a biorthogonal formulation of coupled-cluster (CC) theory using a redundant projected atomic orbital (PAO) basis. The biorthogonal formulation provides simple equations, where the projectors involved in the definition of the PAO basis are absorbed in the integrals. Explicit expressions for the coupled-cluster singles and doubles equations are derived in the PAO basis. The PAO CC equations can be written in a form identical to the standard molecular orbital CC equations, only with integrals that are related to the atomic orbital integrals through different transformation matrices. The dependence of cluster amplitudes, integrals, and correlation energy contributions on the distance between the participating atomic centers and on the number of involved atomic centers is illustrated in numerical case studies. It is also discussed how the present reformulation of the CC equations opens new possibilities for reducing the number of involved parameters and thereby the computational cost.     

Making More Extensive Use of the Coupled-cluster Wave Function: from the Standard Energy Expression to the Energy Expectation Value

The standard coupled-cluster (CC) scheme with single and double excitations in the cluster operator (CCSD) includes only up to quadruple excitations in the equations. The CCSD exponential expansion generates, however, all possible excitations out of the reference function through products of the cluster operators. Clearly, in all standard approximate CC approaches only a part of the CC wave function is used in the equations. If the standard CCSD wave function is inserted into the energy expectation value expression then the complete CCSD wave function contributes to the energy. Such an energy expectation value expression can be presented as a sum of the standard CCSD energy formula plus correction terms. The correction terms provide an information about the quality of the total CC function. Contributions associated with the presence of higher than double excitations in the bra CCSD wave function supplement the CCSD energy obtained within the standard scheme. These contributions can be generated in a sequential way by considering intermediate excitation levels for the bra CCSD wave function in the expectation value expression before reaching the highest excitation level. In this way the importance of particular components differing in the standard and expectation value CCSD energies can be investigated. Some of the contributions can be recognized as close to or identical with the so-called renormalized noniterative corrections to the CC methods. We try to see to what an extent the nonstandard energy expressions, like the energy expectation value or the asymmetric energy formula, can be used to extend the applicability of the CCSD method illustrating our considerations with some numerical examples.Dedicated to Professor Jean-Paul Malrieu to honor his contribution to quantum chemistry and physics     

Model study of the impact of orbital choice on the accuracy of coupled-cluster energies. I. Single-reference-state formulation

The impact of the choice of molecular orbital sets on the results of single-reference-state coupled-cluster (CC) methods was studied for the H4 model. This model offers a straightforward way of taking into account all possible symmetry-adapted orbitals. Moreover, the degree of quasi-degeneracy of its ground state can be varied over a wide range by changing its geometry. The CCD, CCSD, and CCSDT approaches are considered. Surfaces representing the dependence of the energy on the parameters defining the orbitals are obtained. It is documented that for every method there exist alternative orbital sets which allow one to obtain more accurate energies than the standard (HF, BO, and NO) ones. However, for many of the former orbital sets, one obtains relatively large one-body amplitudes or one may encounter problems with solving the CC equations by conventional methods. An interesting variety of orbitals which might be useful for studies of quasi-degenerate states by the CCD method was found. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 205–219, 1998     

Quasispin symmetry for the derivation of coupled cluster equations for the Hubbard model of benzene

The concept of quasispin is applied to a special case of the Pariser–Parr–Pople (PPP ) model of the benzene molecule, namely, the Hubbard Hamiltonian. Added to the spin, space, and alternancy symmetries already taken into account in the PPP Hamiltonian, this new symmetry, called quasispin symmetry, has the effect of reducing the size of the CI matrix. Coupled cluster (CC ) equations are then obtained after applying the CC approach with doubles as well as its extension that accounts for triexcited clusters (CCSDT -1). The derivation of these equations following the use of quasispin to the Hubbard model of benzene constitutes the most simple nontrivial example of CC results. In addition, the CC equations can be written in explicit algebraic form using the symbolic computation language MAPLE.     

CCSD-PCM: improving upon the reference reaction field approximation at no cost

The combination of the coupled cluster (CC) method with the polarizable continuum model (PCM) of solvation requires a much larger computational effort than gas phase CC calculations, since the PCM contribution depends nonlinearly on the CC reduced density: perturbation theory energy and density (PTED) scheme. An approximation can be introduced that neglects the "correlation" PCM contribution and only considers the "reference" PCM contribution to the free energy: PTE scheme. The PTE scheme is a computationally efficient strategy, since the cost is comparable to gas phase CC, but the difference in the free energy with respect to the PTED scheme can be significant. In this work, two intermediate approximations are presented, PTE(S) and PTES (where S stands for singles), which retain the computational efficiency of the PTE scheme while reducing the energy gap with the PTED scheme. PTE(S) only introduces an energy correction to the PTE free energy, while PTES introduces explicit PCM terms in the iterative solution of the CC equations. PTE(S) improves the PTE free energy, although such correction is small. PTES recovers 50%-80% of the PTE-PTED difference and represents a promising approach to perform calculations in solution of CC quality at a cost comparable to gas phase CC. The expressions for the CC-PTE(S) and -PTES wave functions, free energy, and free energy analytical gradients are presented, and the methods are tested with numerical examples.     

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.     

Orthogonally spin-adapted multi-reference Hilbert space coupled-cluster formalism: diagrammatic formulation

Summary The problem of spin-adaptation of the multi-reference (MR) coupled-cluster (CC) formalism, employing Jeziorski-Monkhorst ansatz, is addressed. The diagrammatic technique based on graphical methods of spin algebras is generalized to the MR case, so that both direct and coupling terms can be determined. Usefulness of this fully diagrammatic spin-adaptation approach is illustrated on a derivation of explicit expressions for the linear and bilinear coupling terms that are required in the special two-reference MR-CC theory involving singly and doubly excited states (MR-CCSD formalism). Results obtained with the diagrammatic approach are compared with those derived earlier using the algebraic technique and relative advantages of both procedures are compared.To Prof. Klaus Ruedenberg at the occasion of his 70th anniversary     

Improved embedded molecular cluster model

We demonstrate that boundary effects (i.e., displacements of the cluster boundary atoms from their lattice sites and the difference between effective charges of the perfect crystal atoms and those of the cluster atoms in the case of a cluster with no point defect in it) in an embedded molecular cluster (EMC) model can be radically reduced. A new embedding scheme is proposed. It includes search for the structural elements (SE) of which perfect crystal is composed, use of corresponding to these SE expression for the total energy, and application of the degree of localization of equations consistent with the wave functions of the cluster. To get equations for the cluster wave functions, the problem of varied subsystem in the field of the frozen remaining part of the whole electron system” is investigated in the framework of a one-electron approximation. The consideration is general for every task of this type. Homogeneous equations resulting directly from variation of the total energy expression are obtained and transformed to the eigenvalue problem equations. Orthogonality constraints are not imposed during variation. A particular case of the equations describing mutually orthogonal one-electron wave functions of the cluster staying nonorthogonal to those of the remaining crystal is found. A proposed embedding scheme is realized in the CLUSTER code based on the calculation scheme of the semiempirical INDO method. Boundary effects both in the standard (cluster in the field of the infinite lattice of nonpoint spherical charges) and new embedding scheme are investigated, calculating the clusters of LiF, MgO, NaCl, KCl, and AgCl crystals. Significant reduction of the boundary effects in the new embedding scheme is achieved. Reasons for the boundary effects are discussed. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

A comparison of singlet and triplet states for one- and two-dimensional graphene nanoribbons using multireference theory

This study examines the radical nature and spin symmetry of the ground state of the quasi-linear acene and two-dimensional periacene series. For this purpose, high-level ab initio calculations have been performed using the multireference averaged quadratic coupled cluster theory and the COLUMBUS program package. A reference space consisting of restricted and complete active spaces is taken for the π-conjugated space, correlating 16 electrons with 16 orbitals with the most pronounced open-shell character for the acenes and a complete active-space reference approach with eight electrons in eight orbitals for the periacenes. This reference space is used to construct the total configuration space by means of single and double excitations. By comparison with more extended calculations, it is shown that a focus on the π space with a 6-31G basis set is sufficient to describe the major features of the electronic character of these compounds. The present findings suggest that the ground state is a singlet for the smaller members of these series, but that for the larger ones, singlet and triplet states are quasi-degenerate. Both the acenes and periacenes exhibit significant polyradical character beyond the traditional diradical.     

Applicability of multi-reference many-body perturbation theory to the determination of potential energy surfaces: A model study

The recently written CI -based multi-reference many-body perturbation theory (MR-MBPT ) program package is exploited to study a simple ab initio minimum basis set model involving four hydrogen atoms in a rectangular configuration. This model was examined earlier by several authors using both coupled cluster (CC ) and finite-order MBPT approaches. Here we present the MR-MBPT results up to the 50th order and examine the effect of various shifting techniques on the convergence behavior of this approach. It is shown that in contrast with CC methods, both single and MR finite-order MBPT potential energy calculations are plagued with convergency and intruder state problems, which can be particularly severe when the latter approach is employed for non-degenerate situations.     

A contemporary linear representation theory for ordinary differential equations: probabilistic evolutions and related approximants for unidimensional autonomous systems

In this paper, we build on our previous research on probabilistic foundations of dynamical systems and introduce a theory of linear representation for ordinary differential equations. The theory is developed for explicit ODEs and can be further extended to cover implicit cases. In this report, we investigate the case of a canonical single unknown autonomous system. First we construct a linear representation to get an infinite linear ODE set with a constant coefficient matrix which can be transformed into an upper triangular form. Then we find its approximate truncated solutions. We describe a number of properties of the theory using this framework. The companion of this paper expands this canonical approach to cover multidimensional cases using the theory of folded arrays which is another line of research established by our research group.