Resonant-convergent second-order nonlinear response functions at the levels of Hartree-Fock and Kohn-Sham density functional theory.

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.

[1]  Patrick Norman,et al.  Non-linear electric and magnetic properties obtained from cubic response functions in the random phase approximation , 1996 .

[2]  D. M. Bishop.  General dispersion formulas for molecular third‐order nonlinear optical properties , 1989 .

[3]  N. Handy,et al.  A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP) , 2004 .

[4]  M Gu,et al.  Use of two-photon excitation for erasable-rewritable three-dimensional bit optical data storage in a photorefractive polymer. , 1999, Optics letters.

[5]  H. Ågren,et al.  Cubic response functions in the random phase approximation , 1995 .

[6]  A. Dvornikov,et al.  Two-photon three-dimensional optical storage memory. , 2009, The journal of physical chemistry. A.

[7]  S. Karna,et al.  Frequency dependent nonlinear optical properties of molecules: Formulation and implementation in the HONDO program , 1991 .

[8]  Trygve Helgaker,et al.  Recent advances in wave function-based methods of molecular-property calculations. , 2012, Chemical reviews.

[9]  J. Olsen,et al.  Multiconfigurational Hartree–Fock response functions , 1983 .

[10]  J. Olsen,et al.  Linear and nonlinear response functions for an exact state and for an MCSCF state , 1985 .

[11]  Poul Jørgensen,et al.  Quasienergy formulation of damped response theory. , 2009, The Journal of chemical physics.

[12]  D A Parthenopoulos,et al.  Three-Dimensional Optical Storage Memory , 1989, Science.

[13]  D. M. Bishop.  General dispersion formulae for atomic third-order non-linear optical properties , 1988 .

[14]  D. M. Bishop,et al.  Near-resonant absorption in the time-dependent self-consistent field and multiconfigurational self-consistent field approximations , 2001 .

[15]  D. M. Bishop.  Dispersion formula for the average first hyperpolarizability β , 1991 .

[16]  Poul Jørgensen,et al.  Response functions from Fourier component variational perturbation theory applied to a time-averaged quasienergy , 1998 .

[17]  Daniel C. Flynn,et al.  Synthesis and Optical Properties of Two-Photon Absorbing GFP-type Probes , 2011 .

[18]  J. Bolger,et al.  Linear and nonlinear optical properties of platinum-ethynyl , 1998 .

[19]  Patrick Norman,et al.  A Polarization Propagator for Nonlinear X-ray Spectroscopies. , 2016, The journal of physical chemistry letters.

[20]  P. Jørgensen,et al.  Damped response theory description of two-photon absorption. , 2011, The Journal of chemical physics.

[21]  P. Jørgensen,et al.  Frequency-dependent second hyperpolarizabilities using coupled cluster cubic response theory , 1998 .

[22]  Patrick Norman,et al.  A perspective on nonresonant and resonant electronic response theory for time-dependent molecular properties. , 2011, Physical chemistry chemical physics : PCCP.

[23]  L. Jensen,et al.  Simulating Third-Order Nonlinear Optical Properties Using Damped Cubic Response Theory within Time-Dependent Density Functional Theory. , 2016, Journal of chemical theory and computation.

[24]  W. M. McClain,et al.  Polarization Dependence of the Two-Photon Absorption of Tumbling Molecules with Application to Liquid 1-Chloronaphthalene and Benzene , 1970 .

[25]  W. M. McClain.  Excited State Symmetry Assignment Through Polarized Two‐Photon Absorption Studies of Fluids , 1971 .

[26]  Bishop.  Dispersion formulas for certain nonlinear optical processes. , 1988, Physical review letters.

[27]  Dan Jonsson,et al.  A general, recursive, and open‐ended response code , 2014, J. Comput. Chem..

[28]  Patrick Norman,et al.  CUBIC RESPONSE FUNCTIONS IN THE MULTICONFIGURATION SELF-CONSISTENT FIELD APPROXIMATION , 1996 .

[29]  K. Kneipp,et al.  Two-photon vibrational spectroscopy for biosciences based on surface-enhanced hyper-Raman scattering , 2006, Proceedings of the National Academy of Sciences.

[30]  Fumihiko Aiga,et al.  Frequency‐dependent hyperpolarizabilities in the Mo/ller–Plesset perturbation theory , 1993 .

[31]  Michael P. Gaj,et al.  Third-harmonic generation and its applications in optical image processing , 2009 .

[32]  Luca Frediani,et al.  The Dalton quantum chemistry program system , 2013, Wiley interdisciplinary reviews. Computational molecular science.

[33]  Michel Dupuis,et al.  Ab initio analytic polarizability, first and second hyperpolarizabilities of large conjugated organic molecules: Applications to polyenes C4H6 to C22H24 , 1988 .

[34]  C. Hättig,et al.  Derivation of coupled cluster excited states response functions and multiphoton transition moments between two excited states as derivatives of variational functionals , 1998 .

[35]  Trygve Helgaker,et al.  Efficient parallel implementation of response theory: Calculations of the second hyperpolarizability of polyacenes , 1996 .

[36]  M. Jungen,et al.  Universal Gaussian basis sets for an optimum representation of Rydberg and continuum wavefunctions , 1989 .

[37]  J. Strickler,et al.  Three-dimensional optical data storage in refractive media by two-photon point excitation. , 1991, Optics letters.

[38]  Hideo Sekino,et al.  Frequency dependent nonlinear optical properties of molecules , 1986 .

[39]  W. Webb,et al.  Exchange of protein molecules through connections between higher plant plastids. , 1997, Science.

[40]  Lee W. Tutt,et al.  Optical limiting performance of C60 and C70 solutions , 1992, Nature.

[41]  Evert Jan Baerends,et al.  Accurate density functional calculations on frequency-dependent hyperpolarizabilities of small molecules , 1998 .

[42]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[43]  O. Christiansen,et al.  Static and frequency-dependent polarizabilities of excited singlet states using coupled cluster response theory , 1998 .

[44]  P. Norman,et al.  Theoretical study of linear and nonlinear absorption in platinum-organic compounds , 2002 .

[45]  D. M. Bishop,et al.  The frequency dependence of nonlinear optical processes , 1996 .

[46]  R. Bartlett,et al.  Hyperpolarizabilities of molecules with frequency dependence and electron correlation , 1991 .

[47]  Branislav Jansik,et al.  Calculations of static and dynamic polarizabilities of excited states by means of density functional theory. , 2004, The Journal of chemical physics.

[48]  Richard L. Sutherland,et al.  Understanding the one-photon photophysical properties of a two-photon absorbing chromophore , 2004 .

[49]  E. F. Hayes,et al.  Time‐Dependent Hellmann‐Feynman Theorems , 1965 .

[50]  Paweł Sałek,et al.  Cubic response functions in time-dependent density functional theory. , 2005, The Journal of chemical physics.

[51]  D. M. Bishop,et al.  Nonlinear response theory with relaxation: the first-order hyperpolarizability. , 2005, The Journal of chemical physics.

[52]  P. W. Langhoff,et al.  Aspects of Time-Dependent Perturbation Theory , 1972 .

[53]  Kenneth Ruud,et al.  Benchmarking two-photon absorption cross sections: performance of CC2 and CAM-B3LYP. , 2015, Physical chemistry chemical physics : PCCP.

[54]  R. Bartlett,et al.  New algorithm for high‐order time‐dependent hartree–fock theory for nonlinear optical properties , 1992 .

[55]  Kenneth Ruud,et al.  A density matrix-based quasienergy formulation of the Kohn-Sham density functional response theory using perturbation- and time-dependent basis sets. , 2008, The Journal of chemical physics.

[56]  Evert Jan Baerends,et al.  Calculating frequency-dependent hyperpolarizabilities using time-dependent density functional theory , 1998 .

[57]  K. Ruud,et al.  The ab initio calculation of molecular electric, magnetic and geometric properties. , 2011, Physical chemistry chemical physics : PCCP.

[58]  F. Aiga,et al.  Higher‐order response theory based on the quasienergy derivatives: The derivation of the frequency‐dependent polarizabilities and hyperpolarizabilities , 1993 .

[59]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[60]  D. M. Bishop,et al.  The frequency dependence of hyperpolarizabilities for noncentrosymmetric molecules , 1996 .