Response function coupled cluster

In a recent publication [22] we reported the implementation of dispersion coefficients for first hyperpolarizabiiities based on the coupled cluster quadratic response approach. In the present publication we extend the work of Refs. [22-24] to the analytic calculation of dispersion coefficients for cubic response properties, i.e. second hyperpolarizabiiities. We define the dispersion coefficients by a Taylor expansion of the cubic response function in its frequency arguments. Hence, this approach is... 

In the next section we derive the Taylor expansion of the coupled cluster cubic response function in its frequency arguments and the equations for the required expansions of the cluster amplitude and Lagrangian multiplier responses. For the experimentally important isotropic averages 7, 7i and yx we give explicit expressions for the A and higher-order coefficients in terms of the coefficients of the Taylor series. In Sec. 4 we present an application of the developed approach to the second hyperpolarizability of the methane molecule. We test the convergence of the hyperpolarizabilities with respect to the order of the expansion and investigate the sensitivity of the coefficients to basis sets and correlation treatment. The results are compared with dispersion coefficients derived by least square fits to experimental hyperpolarizability data or to pointwise calculated hyperpolarizabilities of other ab inito studies. 

An implementation of the cubic response function Eq. (30) for the coupled cluster model hierachy CCS, CC2 and CCSD was reported in Ref. [24]. 

Koch H, Jensen HJA, Jorgensen P, Helgaker T (1990) Excitation-energies from the coupled cluster singles and doubles linear response function (CCSDLR) - applications to be, CH+, CO, and H2O. J Chem Phys 93 3345... 

Dunning s correlation consistent basis sets cc-pVAZ [27] augmented with diffuse functions [28] were used in the calculations. We considered cardinal numbers X—D, T, Q, 5, 6 and single (s), double (d), triple (t), and quadruple (q) augmentations. The orbitals were not allowed to relax in the coupled cluster response calculations. 

Static charge-density susceptibilities have been computed ab initio by Li et al (38). The frequency-dependent susceptibility x(r, r cd) can be calculated within density functional theory, using methods developed by Ando (39 Zang-will and Soven (40 Gross and Kohn (4I and van Gisbergen, Snijders, and Baerends (42). In ab initio work, x(r, r co) can be determined by use of time-dependent perturbation techniques, pseudo-state methods (43-49), quantum Monte Carlo calculations (50-52), or by explicit construction of the linear response function in coupled cluster theory (53). Then the imaginary-frequency susceptibility can be obtained by analytic continuation from the susceptibility at real frequencies, or by a direct replacement co ico, where possible (for example, in pseudo-state expressions). 

Also in response theory the summation over excited states is effectively replaced by solving a system of linear equations. Spin-orbit matrix elements are obtained from linear response functions, whereas quadratic response functions can most elegantly be utilized to compute spin-forbidden radiative transition probabilities. We refrain from going into details here, because an excellent review on this subject has been published by Agren et al.118 While these authors focus on response theory and its application in the framework of Cl and multiconfiguration self-consistent field (MCSCF) procedures, an analogous scheme using coupled-cluster electronic structure methods was presented lately by Christiansen et al.124... 

LINEAR RESPONSE THEORY IN CONNECTION TO DENSITY FUNCTIONAL THEORY/MOLECULAR DYNAMICS AND COUPLED CLUSTER/MOLECULAR DYNAMICS METHODS... 

Our present focus is on density functional theory and coupled cluster methods for describing molecular systems interacting with a structured environment, and we focus on the derivation of linear response properties and compare the expressions that we obtain for the two different electronic structure methods. Based on linear response... 

In this section we outline the coupled cluster-molecular mechanics response method, the CC/MM response method. Ref. [51] considers CC response functions for molecular systems in vacuum and for further details we refer to this article. The identification of response functions is closely connected to time-dependent perturbation theory [51,65,66,67,68,69,70], Our starting point is the quasienergy and we identify the response functions as simple derivatives of the quasienergy. For a molecular system in vacuum where Hqm is the vacuum Hamiltonian for the unperturbed molecule and V" is a time-dependent perturbation we have the following time-dependent Hamiltonian, H,... 

Dalgaard E, Monkhorst HJ (1983) Some aspects of the time-dependent coupled-cluster approach to dynamic response functions. Phys Rev A 28 1217—1222. 

Koch H, Jprgensen P (1990) Coupled cluster response functions. J Chem Phys 93 3333—3344. 

---