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.     

Fock space multi-reference coupled cluster study of transition moment and oscillator strength

Summary A method of calculating transition moment and oscillator strength within the framework of the Fock space multi-reference coupled cluster method is described. Diagrammatic technique is used to obtain coupled cluster equations. The general form of equations for the transition moment betweenN-electron ground and excited states is obtained. MBPT analysis of the final equations is done. The excitation energies, dipole transition moments and oscillator strengths for theCH ⁺ molecule are calculated.     

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     

Fock space multireference coupled cluster theory: noniterative inclusion of triples for excitation energies

In this paper we propose and numerically implement a specific scheme for calculating the excitation energies (EEs) within the Fock space multireference coupled cluster framework, which includes the contributions from noniterative triples cluster amplitudes. These contribute to the EEs at the third order. We present results for CH⁺ and N2, and study the effects of these noniterative triples on EEs. Received: 28 July 1997 / Accepted: 8 December 1997     

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.     

On the calculation of expectation values and transition matrix elements by coupled cluster method

Summary Finite order expressions are derived for expectation values and transition matrix elements within the framework of the coupled cluster method.     

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     

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.     

A variational coupled cluster theory for closed shells using a propagator modification procedure

A general partial summation method for including arbitrary classes of diagrams to all orders in the coupled cluster based size consistent energy functional for closed shell states is developed. Since the various reduced density matrices which appear in the energy functional are essentially the time-independent analogues of the corresponding many body Green functions, it is possible to derive Dyson-like equations for these quantities. By expanding the associated proper self energy parts in terms of the T-amplitudes, one can carry out partial summations in the reduced density matrices and thus in energy. At a higher level, higher order terms in a proper self energy can also be generated by renormalizing the internal propagators in it, and considering only the irreducible self-energy terms.     

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     

On expectation value calculations of one-electron properties using the coupled cluster wave functions

The ability of various approximate coupled cluster (CC) methods to provide accurate first-order one-electron properties calculated as expectation values is theoretically analysed and computationally examined for BH and CO. For actual calculations the infinite number of terms of the expectation value expansion (O=¦exp (T ⁺)O exp (T)¦_c) was truncated so that T ₁ T ₂, T ₃, and (1/2) T ₂T₂ clusters were retained on both sides of O. The role of individual clusters is carefully discussed. Inclusion of T ₁, is unavoidable, but if triples are essential in the energy evaluation, they may play an even more important role in the property expansion, as shown in the case of CO. It is shown that the CC wave function, which is exact to second order, effectively satisfies the Hellmann-Feynman theorem.     

Direct evaluation of one-electron properties in coupled cluster methods

Summary An analysis of a method for approximate calculations of expectation values for one-electron operators from available coupled cluster amplitudes is presented and illustrated numerically for the polarizability of the Be atom. The one-particle density matrix resulting from the present approach is accurate through the fourth order in the electron correlation perturbation. It has been found that, in order to obtain quantitative agreement between the energy derivative results and the approximate expectation value formalism, the third orderT ₁ T ₂⁽⁰⁾ wave function term must be included into the calculation of the one-particle density matrix. The present method is also considered as a promising tool for calculations of higher-order atomic and molecular properties from high level correlated wave functions.     

Valence bond corrected single reference coupled cluster approach

Summary An extension of the single reference coupled cluster method truncated to 1- and 2-body cluster components (CCSD) to quasidegenerate systems, where 3-and 4-body connected cluster components play an important role, is proposed. The basic idea is to extract the information concerning the 3- and 4-body clusters from some independent source, similarly as was implicitly done in the so-called ACPQ or ACC(S)D methods, and correct accordingly the absolute term in the CCSD equations. As a source of these approximate 3- and 4-body clusters, simple valence bond (VB) type wave functions are employed, since they are capable of describing electronic structure of various molecular systems for a wide range of nuclear conformations including their dissociation. The cluster analysis of these VB wave functions, that provides the desired information concerning the connected 3- and 4-body cluster components, is outlined and the explicit form of required correction terms to the CCSD equations is given.     

The saddle point of the nucleophilic substitution reaction Cl− + CH3Cl: results of large-scale coupled cluster calculations

On the basis of large-scale coupled cluster calculations including connectedz triple substitutions in a perturbative way, the geometrical parameters of the D _{3 h} saddle point of the Walden inversion reaction Cl⁻ + CH₃Cl′→ ClCH₃ + Cl′⁻ are predicted to be R _s (C—Cl) = 2.301 ? and r _s (C—H) = 1.069 ?. The barrier height with respect to the reactants is recommended to be 11.5 ± 1.0 kJ mol⁻¹. Connected triple substitutions lower the barrier height by almost a factor of 2, but have very little influence on the geometric structure of the saddle point. Received: 26 June 1998 / Accepted: 15 July 1998 / Published online: 28 September 1998     

A holomorphic representation approach to the regularization of model field theories in coupled cluster form

Summary We explain in detail the so-called Bargmann or holomorphic representation, and apply it to the general class of single-mode bosonic field theories. Since these model field theories have no attribute of separability and are, in some sense, maximally nonlocal, they are an especially severe test of the capability of coupled cluster methods to parametrize them satisfactorily. They include the cases of anharmonic oscillators of order 2K (K=2, 3,...), for which ordinary perturbation theory is known to diverge, and we therefore make a special study of such systems. We demonstrate for the first time for any quantum-mechanical problem with infinite Hilbert space that both the normal and extended coupled cluster methods (NCCM and ECCM) have phase spaces which rigorously exist. We analyze completely the asymptotic properties of the complete sets of the NCCM and ECCM amplitudes, either of which fully characterizes the system. It is thereby shown how the holomorphic representation can be used to regularize completely all otherwise formally divergent series that appear. We demonstrate in detail how the entire NCCM and ECCM programmes can be carried through for these systems, including the diagonalization of the classically mapped Hamilitonians in the respective classical NCCM and ECCM phase spaces.     

Adiabatic multi-step separation method and its application to coupled oscillators

Extension of the adiabatic approach to a multi-step separation method is presented. This method step by step reduces the multi-dimensional Schrödinger equation to the effective equations of lower dimensions. The reduction procedure allows to take advantage of multi-level hierarchy of various physical systems. The multi-step separation method is applied in the calculation of vibrational energies of coupled oscillators. The new method is found to be very effective and accurate.     

Multireference Fock space coupled cluster method in the effective and intermediate Hamiltonian formulation for the (2,0) sector

The effective and intermediate Hamiltonian (IH) multireference coupled cluster (CC) method with singles (S) and doubles (D) within the double electron attached (2,0) sector of the Fock space (FS) is formulated and implemented. The intermediate Hamiltonian realization of the (2,0) FS problem allows to replace the iterative scheme of the FS-CC equations based on the effective Hamiltonian with the diagonalization of the properly constructed matrix. The proposed method, IH-FS-CCSD (2,0), is rigorously size-extensive, easy to code, and numerically very efficient with the results comparable or slightly better than equation-of-motion ones at the CCSDT (T--triples) level. The performance of the method is discussed on the basis of test calculations for potential energy curves of the systems for which double positive ions dissociate into closed shell fragments (e.g., Na(2) dimer). The double electron attachment (DEA) scheme can be also useful in determination of the excitation spectra for difficult cases. The example is a carbon atom which has two electrons out of the closed shell structure. The newly implemented method is also analyzed by plotting potential energy curve for twisted ethylene case as a function of a dihedral angle between two methylene groups. Using DEA scheme one obtains a smooth, cusp free curve.     

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