Abstract

The second-order nonlinear (or cubic) response function is derived from the Ehrenfest theorem with inclusion made of the finite lifetimes of the excited states, representing the extension of the derivation of the quadratic response function in the same framework [P. Norman et al., J. Chem. Phys. 123, 194103 (2005)]. The resulting damped response functions are physically sound and converging also in near-resonance and resonance regions of the spectrum. Being an accurate approximation for small complex frequencies (defined as the sum of an optical frequency and an imaginary damping parameter), the polynomial expansion of the complex cubic response function in terms of the said frequencies is presented and used to validate the program implementation. In terms of approximate state theory, the computationally tractable expressions of the damped cubic response function are derived and implemented at the levels of Hartree-Fock and Kohn-Sham density functional theory. Numerical examples are provided in terms of studies of the intensity-dependent refractive index of para-nitroaniline and the two-photon absorption cross section of neon. For the latter property, a numerical comparison is made against calculations of the square of two-photon matrix elements that are identified from a residue analysis of the resonance-divergent quadratic response function.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.