Hyperpolarizabilities vibrational wavefunction

The difference equation or numerical integration method for vibrational wavefunctions usually referred to as the Numerov-Cooley method [111] has been extended by Dykstra and Malik [116] to an open-ended method for the analytical differentiation of the vibrational Schrodinger equation of a diatomic. This is particularly important for high-order derivatives (i.e., hyperpolarizabilities) where numerical difficulties may limit the use of finite-field treatments. As in Numerov-Cooley, this is a procedure that invokes the Born-Oppenheimer approximation. The accuracy of the results are limited only by the quality of the electronic wavefunction s description of the stretching potential and of the electrical property functions and by the adequacy of the Born-Oppenheimer approximation. 

The CCSD model gives for static and frequency-dependent hyperpolarizabilities usually results close to the experimental values, provided that the effects of vibrational averaging and the pure vibrational contributions have been accounted for. Zero point vibrational corrections for the static and the electric field induced second harmonic generation (ESHG) hyperpolarizability of methane have recently been calculated by Bishop and Sauer using SCF and MCSCF wavefunctions [51]. 

For nonlinear (magneto-) optical properties, calculations of an accuracy close to that of modern gas phase experiments require - similar to what has also been found for other properties like structures [79, 109], reaction enthalpies [79, 110, 111], vibrational frequencies [112, 113], NMR chemical shifts [114], etc. - at least an approximate inclusion of connected triple excitations in the wavefunction. This has been known for years now from calculations of static hyperpolarizabilities with the CCSD(T) approximation [9-13]. CCSD(T) accounts rather efficiently for connected triples through a perturbative correction on top of CCSD. For the reasons pointed out in Section 2.1 CCSD(T) is, as a two-step approach, not suitable for the calculation of frequency-dependent properties. Therefore, the CC3 model has been proposed [56, 58] as an alternative to CCSD(T) especially designed for use in connection with response theory. CC3 is an approximation to CCSDT - alike CCSDT-la and related methods - where the triples equations are truncated such that the scaling of the computational efforts with system size is reduced to as for CCSD(T),... 

To obtain hyperpolarizabilities of calibrational quality, a number of standards must be met. The wavefunctions used must be of the highest quality and include electronic correlation. The frequency dependence of the property must be taken into account from the start and not be simply treated as an ad hoc add-on quantity. Zero-point vibrational averaging coupled with consideration of the Maxwell-Boltzmann distribution of populations amongst the rotational states must also be included. The effects of the electric fields (static and dynamic) on nuclear motion must likewise be brought into play (the results given in this section include these effects, but exactly how will be left until Section 3.2.). All this is obviously a tall order and can (and has) only been achieved for the simplest of species He, H2, and D2. Comparison with dilute gas-phase dc-SHG experiments on H2 and D2 (with the helium theoretical values as the standard) shows the challenge to have been met. 

---