Vibrational excitation energies from vibrational coupled cluster response theory

Response theory in the context of vibrational coupled cluster (VCC) theory is introduced and used to obtain vibrational excitation energies. The relation to the vibrational configuration interaction (VCI) approach is described, and the increase in accuracy of VCC response energies relative to VCI energies is discussed theoretically in terms of a perturbational order expansion and demonstrated numerically. To illustrate the theory, a pilot implementation is used to obtain anharmonic vibrational frequencies for fundamental, first overtone and combination excitations of formaldehyde as well as for the fundamental transitions of ethylene.     

Vibrational coupled cluster response theory: a general implementation

The calculation of vibrational contributions to molecular properties using vibrational coupled cluster (VCC) response theory is discussed. General expressions are given for expectation values, linear response functions, and transition moments. It is shown how these expressions can be evaluated for arbitrary levels of excitation in the wave function parameterization as well as for arbitrary coupling levels in the potential and property surfaces. The convergence of the method is assessed by benchmark calculations on formaldehyde. Furthermore, excitation energies and infrared intensities are calculated for the fundamental vibrations of furan using VCC limited to up to two-mode and up to three-mode excitations, VCC[2] and VCC[3], as well as VCC with full two-mode and approximate three-mode couplings, VCC[2pt3].     

Quasi-variational coupled cluster theory

We extend our previous work on the construction of new approximations of the variational coupled cluster method. By combining several linked pair functional transformations in such a way as to give appropriately balanced infinite-order contributions, in order to approximate (L) well at all orders, we formulate a new quantum chemical method, which we name quasi-variational coupled cluster. We demonstrate this method to be particularly robust in the regime of strong static electron correlation, improving significantly on our earlier approximate variational coupled cluster approach.     

无自旋混杂的耦簇理论

本文用APCCSD(T)，自旋投影方法消除CCSD中的自旋混杂，自旋投影算符^P用自旋湮灭算符^As 1和^As 2的乘积近似表出，提高了计算精度。     

Approximate variational coupled cluster theory

We show that it is possible to construct an accurate approximation to the variational coupled cluster method, limited to double substitutions, from the minimization of a functional that is rigorously extensive, exact for isolated two-electron subsystems and invariant to transformations of the underlying orbital basis. This approximate variational coupled cluster theory is a modification and enhancement of our earlier linked pair functional theory. It is first motivated by the constraint that the inverse square root of the matrix that transforms the cluster amplitudes must exist. Low-order corrections are then included to enhance the accuracy of the approximation of variational coupled cluster, while ensuring that the computational complexity of the method never exceeds that of the standard traditional coupled cluster method. The effects of single excitations are included by energy minimization with respect to the orbitals defining the reference wavefunction. The resulting quantum chemical method is demonstrated to be a robust approach to the calculation of molecular electronic structure and performs well when static correlation effects are strong.     

Multireference coupled cluster theory in Fock space

Summary A Fock space multireference coupled cluster method based on incomplete model spaces is described. Some of the essential computational aspects of the theory are discussed with the aid of the diagrammatic representation of the equations. An application to the calculation of ionization potentials and excitation energies ofs-tetrazine is presented along with comparisons with conventionalab initio calculations and experimental results.     

Method of moments approach and coupled cluster theory

Summary The single reference coupled cluster (CC) approach to the many-electron correlation problem is examined from the viewpoint of the method of moments (MM). This yields generally an inconsistent (overcomplete) set of equations for cluster amplitudes, which can be solved either in the least squares sense or by selective projection process restricting the number of equations to that of the unknowns. These resulting generalized MM-CC equations always contain the standard CC equations as a special case. Since, in the MM-CC formalism, the Schrödinger equation will be approximately satisfied on a subspace spanned by non-canonical configurations, this procedure may be helpful in extending the standard single reference CC theory to quasi-degenerate situations. To examine the potential usefulness of this idea, we explore the linear version of the CC approach for systems with a quasi-degenerate reference, in which case the standard linear theory is plagued with singularities due to the intruder states. Implications of this analysis for the structure of the wavefunction are also briefly discussed.Killam Research Fellow 1987–89     

Two-particle density matrix cumulant of coupled cluster theory

A two-particle density matrix obtained from the expectation value expression for the coupled cluster wave function is separated into the antisymmetrized product of the one-particle density matrices and the remaining cumulant part. It is demonstrated that the proposed formula for the coupled cluster cumulant is a valid cumulant expression, since it is a connected, and therefore size-extensive quantity. It is also shown that the density matrices from coupled cluster gradient theory, when used to define a cumulant, result in the expression containing disconnected terms. The proposed formulation of the coupled cluster cumulant makes it easy to develop size-extensive truncation schemes. As an example, explicit equations for the cumulant at the coupled cluster single, double and triple excitation level are presented.     

Unitary group approach to spin-adapted open-shell coupled cluster theory

We show that the irreducible tensor operators of the unitary group provide a natural operator basis for the exponential Ansatz which preserves the spin symmetry of the reference state, requires a minimal number of independent cluster amplitudes for each substitution order, and guarantees the invariance of the correlation energy under unitary transformations of core, open-shell, and virtual orbitals. When acting on the closed-shell reference state with n_c doubly occupied and n_v unoccupied (virtual) orbitals, the irreducible tensor operators of the group U(n_c) ? U(n_V) generate all Gelfand-Tsetlin (GT) states corresponding to appropriate irreducible representation of U(n_c + n_v). The tensor operators generating the M-tuply excited states are easily constructed by symmetrizing products of M unitary group generators with the Wigner operators of the symmetric group S_M. This provides an alternative to the Nagel-Moshinsky construction of the GT basis. Since the corresponding cluster amplitudes, which are also U(n_c) ? U(n_s) tensors, can be shown to be connected, the irreducible tensor operators of U(n_c) ? U(n_v) represent a convenient basis for a spin-adapted full coupled cluster calculation for closed-shell systems. For a high-spin reference determinant with n, singly occupied open-shell orbitals, the corresponding representation of U(n), n=n_c + n_v + n_s is not simply reducible under the group U(n_c) ? U(n_s) ? U(n_v). The multiplicity problem is resolved using the group chain U(n) ? U(n_c + n_v) ? U(n_s) ? U(n_c) ?U(n_s)? U(n_v) ? U(n_v). The labeling of the resulting configuration-state functions (which, in general, are not GT states when n_c > 1) by the irreducible representations of the intermediate group U(n_c + n_v) ?U(n_s) turns out to be equivalent to the classification based on the order of interaction with the reference state. The irreducible tensor operators defined by the above chain and corresponding to single, double, and triple substitutions from the first-, second-, and third-order interacting spaces are explicitly constructed from the U(n) generators. The connectedness of the corresponding cluster amplitudes and, consequently, the size extensivity of the resulting spin-adapted open-shell coupled cluster theory are proved using group theoretical arguments. The perturbation expansion of the resulting coupled cluster equations leads to an explicitly connected form of the spin-restricted open-shell many-body perturbation theory. Approximation schemes leading to manageable computational procedures are proposed and their relation to perturbation theory is discussed. © 1995 John Wiley & Sons, Inc.     

The Fock space coupled cluster method: theory and application

Summary The Fock space coupled cluster method and its application to atomic and molecular systems are described. The importance of conserving size extensivity is demonstrated by the electron affinities of the alkali atoms. Two types of intruder states are discussed, one attributable to the orbital energy spectrum and the other caused by two-electron interactions. They are illustrated by the excited states of Li₂ and by¹ S states of Be, respectively. It is shown how both problems may be solved using incomplete model spaces. The selection of the model space in a moderately dense spectrum is discussed in connection with N₂ excited states.Supported in part by the U.S.-Israel Binational Science Foundation     

Dispersion coefficients for first hyperpolarizabilities using coupled cluster quadratic response theory

The frequency dependence of third-order properties can in the normal dispersion region be expanded in a Taylor series in the frequency arguments. The dispersion coefficients thus obtained provide an efficient way of expressing the dispersion of frequency-dependent properties and are transferable between different optical processes. We derive analytic expressions for the dispersion coefficients of third-order properties in coupled cluster quadratic response theory and report an implementation for the three coupled cluster models CCS, CC2, and CCSD. Calculations are performed for the first hyperpolarizability of the NH₃ molecule. The convergence of the dispersion expansion with the order of the coefficients is examined and we find good convergence up to about half the frequency at which the first pole in the hyperpolarizability occurs. Padé approximants improve the convergence dramatically and extend the application range of the dispersion expansion to frequencies close to the first pole. The sensitivity of the dispersion coefficients on the dynamic correlation treatment and on the choice of the one-electron basis set is investigated. The results demonstrate that, contrary to presumptions in the literature, the dispersion coefficients are sensitive to basis set effects and correlation treatment similar to the static hyperpolarizabilities. Received: 26 March 1998 / Accepted: 21 July 1998 / Published online: 19 October 1998     

Fock space multireference coupled cluster theory: Study of shape resonance

The complex absorbing potential along with correlated independent particle potential (CIP) Fock space multireference coupled cluster method is used for the study of resonances. We have studied shape resonance of e^?‐ F₂, e^?‐ N₂O and e^?‐CO molecules. In particular, we have studied e^?‐ F₂ scattering at different bond lengths to know whether is bound at the equilibrium bond length of F₂. © 2013 Wiley Periodicals, Inc.     

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.     

Application of open-shell coupled cluster theory to the ground state of GaAs

Summary Theoretical calculations at the coupled cluster level of theory including all single, double and perturbative triple excitations, CCSD(T), are carried out for the^3– ground state of GaAs. Employing a (7s5p3d1f) basis set, the theoretical predictions forr _e (2.560 Å), _e (217 cm^–1),D _e (1.84 eV), and IP (7.80 eV), are in good agreement with recent experimental results. The importance of includingf-type polarization functions in the basis set and the effect of correlating 3d electrons are discussed in detail.     

Localized optimized orbitals, coupled cluster theory, and chiroptical response properties

The impact of orbital localization on the efficiency and accuracy of the optimized-orbital coupled cluster model is examined for the prediction of chiroptical properties, in particular optical rotation. The specific rotations of several test cases-(P)-[4]triangulane, (S)-1-phenylethanol, and chiral conformers of 1-fluoropentane, heptane, and nonane-were computed using an approach in which localization is enforced throughout the orbital optimization and subsequent linear response computation. This method provides a robust local-correlation scheme for future production-level implementation. Although the cross-over point between the canonical and localized coupled cluster approach lies at larger molecules than for ground-state energies, the scheme presented should still provide reduced scaling sufficient to investigate much larger molecules than are presently accessible.     

Explicitly correlated coupled cluster F12 theory with single and double excitations

Full explicitly correlated F12 coupled cluster theory with single and double excitations and with Slater-type geminal as a correlation factor is introduced and implemented within the standard approximation. The variant "C" that does not require integrals over the commutator between the kinetic operator and the correlation factor has been used. All the necessary integrals are analytically calculated. With variant C also, first results are reported for the correlation factor being the interelectronic distance coordinate, i.e., for original R12 method. Calculations have been performed for a set of eight molecules including CH(2)((1)A(1)), CH(4), NH(3), H(2)O, HF, CO, N(2), and F(2), as well as for the constituting atoms. Atomization energies are reported too.     

Perturbative treatment of triple excitations in internally contracted multireference coupled cluster theory

Internally contracted multireference coupled cluster (ic-MRCC) methods with perturbative treatment of triple excitations are formulated based on Dyall's definition of a zeroth-order Hamiltonian. The iterative models ic-MRCCSDT-1, ic-MRCC3, and their variants ic-MRCCSD(T), ic-MRCC(3) which determine the energy correction from triples by a non-iterative step are consistent in the single-reference limit with CCSDT-1a, CC3, CCSD(T), and CC(3), respectively. Numerical tests on the potential energy surfaces of BeH(2), H(2)O, and N(2) as well as on the structure and harmonic vibrational frequencies of the ozone molecule show that these methods account very well for higher order correlation effects. The ic-MRCCSD(T) method is further applied to the geometry optimization and harmonic frequencies of the symmetric vibrational modes of the binuclear transition metal oxide Ni(2)O(2), to the singlet-triplet splittings of o-, m-, and p-benzyne and to a ring-opening reaction of an azirine compound with the molecular formula C(6)H(7)NO. The size of the active spaces used in this study ranges from CAS(2,2) to CAS(8,8). Comparisons of results based on differently sized active spaces indicate that the ic-MRCCSD(T) method provides a highly accurate and efficient treatment of both static and dynamic electron correlation in connection with minimal active spaces.     

Vibrational absorption spectra from vibrational coupled cluster damped linear response functions calculated using an asymmetric Lanczos algorithm

We report the theory and implementation of vibrational coupled cluster (VCC) damped response functions. From the imaginary part of the damped VCC response function the absorption as function of frequency can be obtained, requiring formally the solution of the now complex VCC response equations. The absorption spectrum can in this formulation be seen as a matrix function of the characteristic VCC Jacobian response matrix. The asymmetric matrix version of the Lanczos method is used to generate a tridiagonal representation of the VCC response Jacobian. Solving the complex response equations in the relevant Lanczos space provides a method for calculating the VCC damped response functions and thereby subsequently the absorption spectra. The convergence behaviour of the algorithm is discussed theoretically and tested for different levels of completeness of the VCC expansion. Comparison is made with results from the recently reported [P. Seidler, M. B. Hansen, W. Gyo?rffy, D. Toffoli, and O. Christiansen, J. Chem. Phys. 132, 164105 (2010)] vibrational configuration interaction damped response function calculated using a symmetric Lanczos algorithm. Calculations of IR spectra of oxazole, cyclopropene, and uracil illustrate the usefulness of the new VCC based method.