Automated calculation of anharmonic vibrational contributions to first hyperpolarizabilities: quadratic response functions from vibrational configuration interaction wave functions.

Quadratic response functions are derived and implemented for a vibrational configuration interaction state. Combined electronic and vibrational quadratic response functions are derived using Born-Oppenheimer vibronic product wave functions. Computational tractable expressions are derived for determining the total quadratic response contribution as a sum of contributions involving both electronic and vibrational linear and quadratic response functions. In the general frequency-dependent case this includes a new and more troublesome type of electronic linear response function. Pilot calculations for the FH, H(2)O, CH(2)O, and pyrrole molecules demonstrate the importance of vibrational contributions for accurate comparison to experiment and that the vibrational contributions in some cases can be very large. The calculation of transition properties between vibrational states is combined with sum-over-states expressions for analysis purposes. On the basis of this some simple analysis methods are suggested. Also, a preliminary study of the effect of finite lifetimes on quadratic response functions is presented.