Cubic response function, dispersion

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. 

In the normal dispersion region below the first pole, response functions can be expanded in power series in their frequency arguments. The four frequencies, associated with the operator arguments of the cubic response function are related by the matching condition a -fwfl +wc -t-U ) — 0. Thus second hyperpoiarizabiiities or in general cubic response properties are functions of only three independent frequency variables, which may be chosen as u>b, ljc and U > ... 

---