Electronic absorption spectra and solvatochromic shifts by the vertical excitation model: solvated clusters and molecular dynamics sampling.

A physically realistic treatment of solvatochromic shifts in liquid-phase electronic absorption spectra requires a proper account for various short- and long-range equilibrium and nonequilibrium solute-solvent interactions. The present article demonstrates that such a treatment can be accomplished using a mixed discrete-continuum approach based on the two-time-scale self-consistent state-specific vertical excitation model (called VEM) for electronic excitation in solution. We apply this mixed approach in combination with time-dependent density functional theory to compute UV/vis absorption spectra in solution for the n → π* ((1)A2) transition for acetone in methanol and in water, the π → π* ((1)A1) transition for para-nitroaniline (PNA) in methanol and in water, the n → π* ((1)B1) transition for pyridine in water, and the n → π* ((1)B1) transition for pyrimidine in water. Hydrogen bonding and first-solvation-shell-specific complexation are included by means of explicit solvent molecules, and solute-solvent dispersion is included by using the solvation model with state-specific polarizability (SMSSP). Geometries of microsolvated clusters were treated in two different ways, (i) using single liquid-phase global-minimum solute-solvent clusters containing up to two explicit solvent molecules and (ii) using solute-solvent cluster snapshots derived from molecular dynamics (MD) trajectories. The calculations in water involve using VEM/TDDFT excitation energies and oscillator strengths computed over 200 MD-derived solute-solvent clusters and convoluted with Gaussian functions. We also calculate ground- and excited-state dipole moments for interpretation. We find that inclusion of explicit solvent molecules generally improves the agreement with experiment and can be recommended as a way to include the effect of hydrogen bonding in solvatochromic shifts.

[1]  Rudolph A. Marcus,et al.  Electrostatic Free Energy and Other Properties of States Having Nonequilibrium Polarization. I , 1956 .

[2]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[3]  J. Tomasi,et al.  Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects , 1981 .

[4]  Ruifang Li,et al.  PAPER www.rsc.org/pccp | Physical Chemistry Chemical Physics Density functional approximations for charge transfer excitations with intermediate spatial overlapw , 2010 .

[5]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[6]  Walter Thiel,et al.  Benchmarks for electronically excited states: CASPT2, CC2, CCSD, and CC3. , 2008, The Journal of chemical physics.

[7]  Donald G. Truhlar,et al.  Effectiveness of Diffuse Basis Functions for Calculating Relative Energies by Density Functional Theory , 2003 .

[8]  A. Millefiori,et al.  Electronic spectra and structure of nitroanilines , 1977 .

[9]  Kurt V. Mikkelsen,et al.  Failures of TDDFT in describing the lowest intramolecular charge-transfer excitation in para-nitroaniline , 2013 .

[10]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

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

[12]  S. Canuto,et al.  A Monte Carlo–quantum mechanical study of the solvatochromism of pyrimidine in water and in carbon tetrachloride , 2001 .

[13]  Donald G Truhlar,et al.  Density functionals with broad applicability in chemistry. , 2008, Accounts of chemical research.

[14]  V. Barone,et al.  Effective method to compute Franck-Condon integrals for optical spectra of large molecules in solution. , 2007, The Journal of chemical physics.

[15]  K. Hahn,et al.  Density Functional Vertical Self-Consistent Reaction Field Theory for Solvatochromism Studies of Solvent-Sensitive Dyes , 2004 .

[16]  Jacopo Tomasi,et al.  A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics , 1997 .

[17]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[18]  D. Sherrington,et al.  Correlations of electrons in small molecules , 1964 .

[19]  G. Chang,et al.  Macromodel—an integrated software system for modeling organic and bioorganic molecules using molecular mechanics , 1990 .

[20]  S. Hoffmann,et al.  The electronic states of pyrimidine studied by VUV photoabsorption and electron energy-loss spectroscopy. , 2010, Physical chemistry chemical physics : PCCP.

[21]  Richard A. Friesner,et al.  Integrated Modeling Program, Applied Chemical Theory (IMPACT) , 2005, J. Comput. Chem..

[22]  C. Cramer,et al.  Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. , 2009, The journal of physical chemistry. B.

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

[24]  Henry F. Schaefer,et al.  On the evaluation of analytic energy derivatives for correlated wave functions , 1984 .

[25]  Mark E. Casida,et al.  Time-dependent density-functional theory for molecules and molecular solids , 2009 .

[26]  K. Hahn,et al.  A theoretical study of the UV/visible absorption and emission solvatochromic properties of solvent-sensitive dyes. , 2003, Chemphyschem : a European journal of chemical physics and physical chemistry.

[27]  F. Lipparini,et al.  Vertical Electronic Excitations in Solution with the EOM-CCSD Method Combined with a Polarizable Explicit/Implicit Solvent Model. , 2013, Journal of chemical theory and computation.

[28]  I. Renge.  Solvent dependence of n-pi* absorption in acetone. , 2009, The journal of physical chemistry. A.

[29]  H. Berendsen,et al.  A LEAP-FROG ALGORITHM FOR STOCHASTIC DYNAMICS , 1988 .

[30]  Giovanni Scalmani,et al.  Practical computation of electronic excitation in solution: vertical excitation model , 2011 .

[31]  M. Born.  Volumen und Hydratationswärme der Ionen , 1920 .

[32]  Jacopo Tomasi,et al.  Nonequilibrium solvation: An ab initio quantum‐mechanical method in the continuum cavity model approximation , 1993 .

[33]  S. Dancoff.  Non-Adiabatic Meson Theory of Nuclear Forces , 1950 .

[34]  M. E. Casida,et al.  Progress in time-dependent density-functional theory. , 2011, Annual review of physical chemistry.

[35]  Jacopo Tomasi,et al.  Geometries and properties of excited states in the gas phase and in solution: theory and application of a time-dependent density functional theory polarizable continuum model. , 2006, The Journal of chemical physics.

[36]  Jeffrey R. Reimers,et al.  The Low-Lying Excited States of Pyridine , 2000 .

[37]  L. M. L. Daku,et al.  Ab initio static and molecular dynamics study of the absorption spectra of the 4-styrylpyridine photoswitch in its cis and trans forms. , 2010, Physical chemistry chemical physics : PCCP.

[38]  Donald G. Truhlar,et al.  Generalized born fragment charge model for solvation effects as a function of reaction coordinate , 1989 .

[39]  S. Kovalenko,et al.  Femtosecond relaxation of photoexcited para-nitroaniline: solvation, charge transfer, internal conversion and cooling , 2000 .

[40]  E. Gross,et al.  Density-Functional Theory for Time-Dependent Systems , 1984 .

[41]  W. P. Weijland,et al.  Reduction potentials of various aromatic hydrocarbons and their univalent anions , 1956 .

[42]  M. B. Pinto,et al.  Optimized δ expansion for relativistic nuclear models , 1997, nucl-th/9709049.

[43]  C. Cramer,et al.  Uniform Treatment of Solute-Solvent Dispersion in the Ground and Excited Electronic States of the Solute Based on a Solvation Model with State-Specific Polarizability. , 2013, Journal of chemical theory and computation.

[44]  Donald G. Truhlar,et al.  Two-response-time model based on CM2/INDO/S2 electrostatic potentials for the dielectric polarization component of solvatochromic shifts on vertical excitation energies , 2000 .

[45]  M. Caricato.  A comparison between state-specific and linear-response formalisms for the calculation of vertical electronic transition energy in solution with the CCSD-PCM method. , 2013, The Journal of chemical physics.

[46]  L. Slipchenko,et al.  Solvent effects on the electronic transitions of p-nitroaniline: a QM/EFP study. , 2011, The journal of physical chemistry. A.

[47]  Donald G Truhlar,et al.  Sorting Out the Relative Contributions of Electrostatic Polarization, Dispersion, and Hydrogen Bonding to Solvatochromic Shifts on Vertical Electronic Excitation Energies. , 2010, Journal of chemical theory and computation.

[48]  S. Kunte,et al.  Solvent effects on electronic absorption spectra of nitrochlorobenzenes, nitrophenols and nitroanilines. I: Studies in nonpolar solvents , 1986 .

[49]  S. Mason.  245. The electronic spectra of N-heteroaromatic systems. Part I. The n?? transitions of monocyclic azines , 1959 .

[50]  Jacopo Tomasi,et al.  Approximate evaluations of the electrostatic free energy and internal energy changes in solution processes , 1982 .

[51]  E. Condon,et al.  Nuclear Motions Associated with Electron Transitions in Diatomic Molecules , 1928 .