Density functional self-consistent quantum mechanics/molecular mechanics theory for linear and nonlinear molecular properties: Applications to solvated water and formaldehyde.

A combined quantum mechanics/molecular mechanics (QM/MM) method is described, where the polarization between the solvent and solute is accounted for using a self-consistent scheme linear in the solvent polarization. The QM/MM method is implemented for calculation of energies and molecular response properties including the calculation of linear and quadratic response functions using the density-functional theory (DFT) and the Hartree-Fock (HF) theory. Sample calculations presented for ground-state energies, first-order ground-state properties, excitation energies, first-order excited state properties, polarizabilities, first-hyperpolarizabilities, and two-photon absorptions strengths of formaldehyde suggests that DFT may in some cases be a sufficiently reliable alternative to high-level theory, such as coupled-cluster (CC) theory, in modeling solvent shifts, whereas results obtained with the HF wave function deviate significantly from the CC results. Calculations carried out on water gives results that also are comparable with CC calculations in accuracy for ground-state and first-order properties. However, to obtain such accuracy an exchange-correlation functional capable of describing the diffuse Rydberg states must be chosen.

[1]  F.J.Olivares del Valle,et al.  A multiconfiguration self-consistent field/molecular dynamics study of the (n→π*)1 transition of carbonyl compounds in liquid water , 2000 .

[2]  David R. Klug,et al.  Large and Fast Relaxations inside a Protein: Calculation and Measurement of Reorganization Energies in Alcohol Dehydrogenase , 2002 .

[3]  P. Prasad,et al.  Observation of stimulated emission by direct three-photon excitation , 2002, Nature.

[4]  Jiali Gao,et al.  HYBRID AB INITIO QM/MM SIMULATION OF N-METHYLACETAMIDE IN AQUEOUS SOLUTION , 1997 .

[5]  D. Klug,et al.  A Quantum Mechanical/Molecular Mechanical Approach to Relaxation Dynamics: Calculation of the Optical Properties of Solvated Bacteriochlorophyll-a , 1999 .

[6]  J. G. Snijders,et al.  A discrete solvent reaction field model within density functional theory , 2003 .

[7]  Gregory K. Schenter,et al.  Excited States of the Bacteriochlorophyll b Dimer of Rhodopseudomonas viridis: A QM/MM Study of the Photosynthetic Reaction Center That Includes MM Polarization , 1995 .

[8]  Arieh Warshel,et al.  Nature of the Surface Crossing Process in Bacteriorhodopsin: Computer Simulations of the Quantum Dynamics of the Primary Photochemical Event , 2001 .

[9]  Arieh Warshel,et al.  Calculations of chemical processes in solutions , 1979 .

[10]  Y. Liu,et al.  Two-Photon Photochromism of an Organic Material for Holographic Recording , 2002 .

[11]  U. Narang,et al.  Characterization of a New Solvent-Sensitive Two-Photon-Induced Fluorescent (Aminostyryl)pyridinium Salt Dye , 1996 .

[12]  Evert Jan Baerends,et al.  Molecular calculations of excitation energies and (hyper)polarizabilities with a statistical average of orbital model exchange-correlation potentials , 2000 .

[13]  W. Denk,et al.  Two-photon laser scanning fluorescence microscopy. , 1990, Science.

[14]  D. D. Yue,et al.  Theory of Electric Polarization , 1974 .

[15]  I. Thormählen,et al.  Refractive Index of Water and Its Dependence on Wavelength, Temperature, and Density , 1985 .

[16]  Piet Th. van Duijnen,et al.  A discrete solvent reaction field model for calculating molecular linear response properties in solution , 2003 .

[17]  A. Warshel,et al.  Microscopic Calculations of Solvent Effects on Absorption Spectra of Conjugated Molecules , 1991 .

[18]  K. Mikkelsen,et al.  Second harmonic generation second hyperpolarizability of water calculated using the combined coupled cluster dielectric continuum or different molecular mechanics methods. , 2004, The Journal of chemical physics.

[19]  K. Mikkelsen,et al.  A CC2 dielectric continuum model and a CC2 molecular mechanics model , 2003 .

[20]  Michele Parrinello,et al.  Water Molecule Dipole in the Gas and in the Liquid Phase , 1999 .

[21]  Kajzar,et al.  Third-harmonic generation in liquids. , 1985, Physical review. A, General physics.

[22]  Kurt V. Mikkelsen,et al.  Linear response functions for coupled cluster/molecular mechanics including polarization interactions , 2003 .

[23]  M. Ratner.  Molecular electronic-structure theory , 2000 .

[24]  W. Webb,et al.  Two‐photon molecular excitation provides intrinsic 3‐dimensional resolution for laser‐based microscopy and microphotochemistry , 1994, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[25]  Jiali Gao,et al.  Solvent effects on the nπ* transition of pyrimidine in aqueous solution , 1997 .

[26]  M. Grüning,et al.  On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies , 2002 .

[27]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

[28]  J. Gao,et al.  A priori evaluation of aqueous polarization effects through Monte Carlo QM-MM simulations. , 1992, Science.

[29]  Jacob Kongsted,et al.  Statistical mechanically averaged molecular properties of liquid water calculated using the combined coupled cluster/molecular dynamics method. , 2006, The Journal of chemical physics.

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

[31]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[32]  Martin Karplus,et al.  Molecular properties from combined QM/MM methods. I. Analytical second derivative and vibrational calculations , 2000 .

[33]  K. Mikkelsen,et al.  Nonlinear optical response properties of molecules in condensed phases using the coupled cluster/dielectric continuum or molecular mechanics methods , 2003 .

[34]  Kurt V. Mikkelsen,et al.  The combined multiconfigurational self-consistent-field/molecular mechanics wave function approach , 2001 .

[35]  Kurt V. Mikkelsen,et al.  The QM/MM approach for wavefunctions, energies and response functions within self-consistent field and coupled cluster theories , 2002 .

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

[37]  M. Thompson,et al.  QM/MMpol: A Consistent Model for Solute/Solvent Polarization. Application to the Aqueous Solvation and Spectroscopy of Formaldehyde, Acetaldehyde, and Acetone , 1996 .

[38]  S. Canuto,et al.  Solvent effects in emission spectroscopy: A Monte Carlo quantum mechanics study of the n←π* shift of formaldehyde in water , 2000 .

[39]  Seth R. Marder,et al.  Photoresponsive Hydrogel Microstructure Fabricated by Two‐Photon Initiated Polymerization , 2002 .

[40]  Seth R. Marder,et al.  Two-photon polymerization initiators for three-dimensional optical data storage and microfabrication , 1999, Nature.

[41]  K. Mikkelsen,et al.  Two-photon absorption cross sections: an investigation of solvent effects. Theoretical studies on formaldehyde and water. , 2006, The Journal of chemical physics.

[42]  C. Bethea,et al.  Effects on hyperpolarizabilities of molecular interactions in associating liquid mixtures , 1976 .

[43]  Jacob Kongsted,et al.  Automatic generation of force fields and property surfaces for use in variational vibrational calculations of anharmonic vibrational energies and zero-point vibrational averaged properties. , 2006, The Journal of chemical physics.

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

[45]  Paras N. Prasad,et al.  Highly Active Two-Photon Dyes: Design, Synthesis, and Characterization toward Application , 1998 .

[46]  Sylvio Canuto,et al.  A Monte Carlo–quantum mechanics study of the spectroscopic properties of molecules in solution , 2001 .

[47]  S. Keiding,et al.  Two-photon dissociation and ionization of liquid water studied by femtosecond transient absorption spectroscopy , 1999 .

[48]  K. Mikkelsen,et al.  Dipole and quadrupole moments of liquid water calculated within the coupled cluster/molecular mechanics method , 2002 .

[49]  Kristian O. Sylvester-Hvid,et al.  Refractive index of liquid water in different solvent models. , 2005, The journal of physical chemistry. A.

[50]  P. Prasad,et al.  Phase-conjugation properties of two-photon-pumped backward-stimulated emission. I. Experimental studies , 1998 .

[51]  Nicholas C. Handy,et al.  Density functional calculations of the hyperpolarisabilities of small molecules , 1999 .

[52]  C. Coulson,et al.  Interactions of H2O molecules in ice I. The dipole moment of an H2O molecule in ice , 1966, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[53]  I. Gould,et al.  Optical properties of solvated molecules calculated by a QMMM method Chlorophyll a and bacteriochlorophyll a , 1997 .

[54]  K. Mikkelsen,et al.  Solvent effects on the n-->pi* electronic transition in formaldehyde: a combined coupled cluster/molecular dynamics study. , 2004, The Journal of chemical physics.

[55]  Melvin B. Robin,et al.  Higher excited states of polyatomic molecules , 1974 .

[56]  R. Leeuwen,et al.  Exchange-correlation potential with correct asymptotic behavior. , 1994, Physical review. A, Atomic, molecular, and optical physics.

[57]  Paweł Sałek,et al.  Benchmarking two-photon absorption with CC3 quadratic response theory, and comparison with density-functional response theory. , 2006, The Journal of chemical physics.

[58]  Yi Luo,et al.  Density functional response theory calculations of three-photon absorption. , 2004, The Journal of chemical physics.

[59]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[60]  M. Karplus,et al.  Molecular Properties from Combined QM/MM Methods. 2. Chemical Shifts in Large Molecules , 2000 .

[61]  Kurt V. Mikkelsen,et al.  Coupled Cluster/Molecular Mechanics Method: Implementation and Application to Liquid Water , 2003 .

[62]  Kurt V. Mikkelsen,et al.  A theoretical study of the electronic spectrum of water , 2000 .

[63]  R. N. Hamm,et al.  Collective oscillation in liquid water , 1974 .

[64]  J. G. Snijders,et al.  A discrete solvent reaction field model for calculating frequency-dependent hyperpolarizabilities of molecules in solution , 2003 .

[65]  S. Canuto,et al.  From hydrogen bond to bulk: Solvation analysis of then-?* transition of formaldehyde in water , 2000 .

[66]  Paweł Sałek,et al.  Density-functional theory of linear and nonlinear time-dependent molecular properties , 2002 .

[67]  K. Mikkelsen,et al.  Coupled cluster calculation of the n --> pi* electronic transition of acetone in aqueous solution. , 2005, The journal of physical chemistry. A.

[68]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[69]  K. Morokuma,et al.  A simple model of solvation within the molecular orbital theory , 1975 .