Hybrid molecular dynamics‐quantum mechanics simulations of solute spectral properties in the condensed phase: Evaluation of simulation parameters

We have explored the impact of a number of basic simulation parameters on the results of a recently developed hybrid molecular dynamics‐quantum mechanics (MD‐QM) method (Mercer et al., J Phys Chem B 1999, 103, 7720). The method utilizes MD simulations to explore the ground‐state configuration space of the system and QM evaluation of those structures to yield the time‐dependent electronic transition energy, which is transformed into the optical line‐broadening function using the second‐order cumulant expansion. Both linear and nonlinear optical spectra can then be generated for comparison to experiment. The dependence of the resulting spectra on the length of the MD trajectory, the QM sampling rate, and the QM model chemistry have all been examined. In particular, for the system of oxazine‐4 in methanol studied here, at least 20 ps of MD trajectory are needed for qualitative convergence of linear spectral properties, and > 100 ps is needed for quantitative convergence. Surprisingly, little difference is found between the 3‐21G and 6‐31G(d) basis sets, and the CIS and TD‐B3LYP methods yield remarkably similar spectra. The semiempirical INDO/s method yields the most accurate results, reproducing the experimental Stokes shift to within 5% and the FWHM to within 20%. Nonlinear 3‐pulse photon echo peak shift (3PEPS) decays have also been simulated. Decays are generally poorly reproduced, though the initial peak shift which depends on the overall coupling of motions to the solute transition energy is within 15% of experiment for all model chemistries other than those using the STO‐3G basis. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007

[1]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

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

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

[4]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. I. Use of Gaussian Expansions of Slater‐Type Atomic Orbitals , 1969 .

[5]  M. Maroncelli.  Computer simulations of solvation dynamics in acetonitrile , 1991 .

[6]  Michael C. Zerner,et al.  An intermediate neglect of differential overlap technique for spectroscopy: Pyrrole and the azines , 1973 .

[7]  Klaus Schulten,et al.  Excitons in a photosynthetic light-harvesting system: a combined molecular dynamics, quantum chemistry, and polaron model study. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  G. Scuseria,et al.  An efficient implementation of time-dependent density-functional theory for the calculation of excitation energies of large molecules , 1998 .

[9]  T. Okada,et al.  Solute Dependence of Three Pulse Photon Echo Peak Shift Measurements in Methanol Solution , 2003 .

[10]  D. Wiersma,et al.  Ultrafast electronic fluctuation and solvation in liquids , 1994 .

[11]  P. Kollman,et al.  Structure and Properties of Neat Liquids Using Nonadditive Molecular Dynamics: Water, Methanol, and N-Methylacetamide , 1995 .

[12]  Mark A. Berg,et al.  Ultrafast Raman Echo Measurements of Vibrational Dephasing and the Nature of Solvent−Solute Interactions , 1997 .

[13]  Dennis R. Salahub,et al.  Molecular excitation energies to high-lying bound states from time-dependent density-functional response theory: Characterization and correction of the time-dependent local density approximation ionization threshold , 1998 .

[14]  S. Bradforth,et al.  Excited States of Iodide Anions in Water: A Comparison of the Electronic Structure in Clusters and in Bulk Solution , 2002 .

[15]  K. Mikkelsen,et al.  A quantum mechanical method for calculating nonlinear optical properties of condensed phase molecules coupled to a molecular mechanics field: A quadratic multiconfigurational self-consistent-field/molecular mechanics response method , 2001 .

[16]  S. Mukamel,et al.  Photon echoes and related four-wave-mixing spectroscopies using phase-locked pulses , 1992 .

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

[18]  Mark S. Gordon,et al.  Self-consistent molecular-orbital methods. 22. Small split-valence basis sets for second-row elements , 1980 .

[19]  Graham R. Fleming,et al.  Two-dimensional spectroscopy of electronic couplings in photosynthesis , 2005, Nature.

[20]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[21]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[22]  Jacopo Tomasi,et al.  An Integrated Effective Fragment—Polarizable Continuum Approach to Solvation: Theory and Application to Glycine , 2002 .

[23]  Wei Zhao,et al.  Interference, dephasing, and vibrational coupling effects between coherence pathways in doubly vibrationally enhanced nonlinear spectroscopies , 2001 .

[24]  David M. Jonas,et al.  Two-dimensional Fourier transform electronic spectroscopy , 2001 .

[25]  E. Carter,et al.  Solvation dynamics for an ion pair in a polar solvent: Time‐dependent fluorescence and photochemical charge transfer , 1991 .

[26]  B. Ladanyi,et al.  Breakdown of linear response for solvation dynamics in methanol , 1991 .

[27]  Michael J. Frisch,et al.  Toward a systematic molecular orbital theory for excited states , 1992 .

[28]  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 .

[29]  Graham R. Fleming,et al.  Third-order nonlinear time domain probes of solvation dynamics , 1996 .

[30]  S. Mukamel.  Principles of Nonlinear Optical Spectroscopy , 1995 .

[31]  David R. Klug,et al.  Comparison of basis set effects and the performance of ab initio and DFT methods for probing equilibrium fluctuations , 2007, J. Comput. Chem..

[32]  P. C. Hariharan,et al.  The effect of d-functions on molecular orbital energies for hydrocarbons , 1972 .

[33]  J. Fourkas.  Higher-order optical correlation spectroscopy in liquids. , 2002, Annual review of physical chemistry.

[34]  Junmei Wang,et al.  Automatic parameterization of force field by systematic search and genetic algorithms , 2001, J. Comput. Chem..

[35]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[36]  Wiersma,et al.  Chirped four-wave mixing. , 1993, Physical review. A, Atomic, molecular, and optical physics.

[37]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[38]  R. M. Stratt,et al.  Rotational Coherence and a Sudden Breakdown in Linear Response Seen in Room-Temperature Liquids , 2006, Science.

[39]  O. Golonzka,et al.  Coupling and orientation between anharmonic vibrations characterized with two-dimensional infrared vibrational echo spectroscopy , 2001 .

[40]  David Chandler,et al.  Importance sampling and theory of nonequilibrium solvation dynamics in water , 2000 .

[41]  J. Pople,et al.  Self-consistent molecular orbital methods. 21. Small split-valence basis sets for first-row elements , 2002 .

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

[43]  M. Maroncelli,et al.  Computer simulation of the dynamics of aqueous solvation , 1988 .

[44]  R. F. Loring,et al.  Calculation of the photon echo with mixed-state propagation , 1998 .

[45]  G. Fleming,et al.  The Integrated Photon Echo and Solvation Dynamics , 1996 .

[46]  K. Nelson,et al.  Time-resolved vibrational spectroscopy in the impulsive limit , 1994 .