FREQUENCY-DEPENDENT HYPERPOLARIZABILITIES IN THE MOLLER-PLESSET PERTURBATION-THEORY

A formulation for calculating frequency-dependent hyperpolarizabilities in the Moller-Plesset perturbation theory is presented as the correlation correction to the TDHF approximation. Our quasienergy derivative (QED) method is applied, and the difference between the QED method and the pseudoenergy derivative (PED) method by Rice and Handy is discussed. The Lagrangian technique is utilized to obtain simple and practical expressions for response properties in which the TDHF orbital rotation parameters satisfy the 2n + 1 rule and the Lagrange multipliers satisfy the 2n + 2 rule. Explicit expressions for response properties up to third order [mu, alpha(omega1;omega1), beta(omega(sigma);omega1,omega2)] are derived in the second-order Moller-Plesset perturbation theory.