Optimizing conical intersections of solvated molecules: the combined spin-flip density functional theory/effective fragment potential method.

Solvent effects on a potential energy surface crossing are investigated by optimizing a conical intersection (CI) in solution. To this end, the analytic energy gradient has been derived and implemented for the collinear spin-flip density functional theory (SFDFT) combined with the effective fragment potential (EFP) solvent model. The new method is applied to the azomethane-water cluster and the chromophore of green fluorescent protein in aqueous solution. These applications illustrate not only dramatic changes in the CI geometries but also strong stabilization of the CI in a polar solvent. Furthermore, the CI geometries obtained by the hybrid SFDFT/EFP scheme reproduce those by the full SFDFT, indicating that the SFDFT/EFP method is an efficient and promising approach for understanding nonadiabatic processes in solution.

[1]  Mark S Gordon,et al.  Optimizing conical intersections by spin-flip density functional theory: application to ethylene. , 2009, The journal of physical chemistry. A.

[2]  M. Gordon,et al.  Gradients of the polarization energy in the effective fragment potential method. , 2006, The Journal of chemical physics.

[3]  David Casanova,et al.  The spin-flip extended single excitation configuration interaction method. , 2008, The Journal of chemical physics.

[4]  Yu Kay Law,et al.  DNA excited-state dynamics: from single bases to the double helix. , 2009, Annual review of physical chemistry.

[5]  Kevin M. Dunn,et al.  Ab initio molecular orbital study of the mechanism of photodissociation of trans‐azomethane , 1996 .

[6]  V. R. Saunders,et al.  On methods for converging open-shell Hartree-Fock wave-functions , 1974 .

[7]  R. Mitrić,et al.  Nonadiabatic dynamics within the time dependent density functional theory: Ultrafast photodynamics in pyrazine , 2008 .

[8]  A. Krylov,et al.  Effect of microhydration on the electronic structure of the chromophores of the photoactive yellow and green fluorescent proteins. , 2011, The Journal of chemical physics.

[9]  M. E. Martín,et al.  Location of conical intersections in solution using a sequential quantum mechanics/molecular dynamics method , 2007 .

[10]  Mark S. Gordon,et al.  Solvent-induced frequency shifts: configuration interaction singles combined with the effective fragment potential method. , 2010, The journal of physical chemistry. A.

[11]  Wenjian Liu,et al.  Spin-adapted open-shell time-dependent density functional theory. III. An even better and simpler formulation. , 2011, The Journal of chemical physics.

[12]  Filipp Furche,et al.  On the density matrix based approach to time-dependent density functional response theory , 2001 .

[13]  S. Meech,et al.  Excited state reactions in fluorescent proteins. , 2009, Chemical Society reviews.

[14]  Thierry Deutsch,et al.  Assessment of noncollinear spin-flip Tamm-Dancoff approximation time-dependent density-functional theory for the photochemical ring-opening of oxirane. , 2010, Physical chemistry chemical physics : PCCP.

[15]  A. Kahan,et al.  Photoreactivity of a push-pull merocyanine in static electric fields: a three-state model of isomerization reactions involving conical intersections. , 2009, The journal of physical chemistry. A.

[16]  M. E. Casida,et al.  Mixed time-dependent density-functional theory/classical trajectory surface hopping study of oxirane photochemistry. , 2008, The Journal of chemical physics.

[17]  B. A. Lindquist,et al.  Photodynamics in complex environments: ab initio multiple spawning quantum mechanical/molecular mechanical dynamics. , 2009, The journal of physical chemistry. B.

[18]  S. Kato,et al.  Locating the lowest free-energy point on conical intersection in polar solvent: reference interaction site model self-consistent field study of ethylene and CH2NH2(+). , 2005, The Journal of chemical physics.

[19]  Walter R. Duncan,et al.  Trajectory surface hopping in the time-dependent Kohn-Sham approach for electron-nuclear dynamics. , 2005, Physical review letters.

[20]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[21]  Evgeny Epifanovsky,et al.  Quantum Chemical Benchmark Studies of the Electronic Properties of the Green Fluorescent Protein Chromophore: 2. Cis-Trans Isomerization in Water. , 2009, Journal of chemical theory and computation.

[22]  Qin Wu,et al.  Configuration interaction based on constrained density functional theory: a multireference method. , 2007, The Journal of chemical physics.

[23]  G. Mazur,et al.  Time-dependent density functional theory gradients in the Amsterdam density functional package: geometry optimizations of spin-flip excitations , 2011 .

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

[25]  Yihan Shao,et al.  General formulation of spin-flip time-dependent density functional theory using non-collinear kernels: theory, implementation, and benchmarks. , 2012, The Journal of chemical physics.

[26]  F. Neese,et al.  Efficient time-dependent density functional theory approximations for hybrid density functionals: analytical gradients and parallelization. , 2011, The Journal of chemical physics.

[27]  Filipp Furche,et al.  Adiabatic time-dependent density functional methods for excited state properties , 2002 .

[28]  M. Gordon,et al.  Implementation of the analytic energy gradient for the combined time-dependent density functional theory/effective fragment potential method: application to excited-state molecular dynamics simulations. , 2011, The Journal of chemical physics.

[29]  Claudia Filippi,et al.  Absorption Spectrum of the Green Fluorescent Protein Chromophore: A Difficult Case for ab Initio Methods? , 2009, Journal of chemical theory and computation.

[30]  Massimo Olivucci,et al.  Probing the rhodopsin cavity with reduced retinal models at the CASPT2//CASSCF/AMBER level of theory. , 2003, Journal of the American Chemical Society.

[31]  Jacob Kongsted,et al.  Scrutinizing the effects of polarization in QM/MM excited state calculations. , 2011, Physical chemistry chemical physics : PCCP.

[32]  Mark S. Gordon,et al.  The Effective Fragment Potential Method: A QM-Based MM Approach to Modeling Environmental Effects in Chemistry , 2001 .

[33]  Kieron Burke,et al.  Time-dependent density functional theory: past, present, and future. , 2005, The Journal of chemical physics.

[34]  Anna I. Krylov,et al.  The spin–flip approach within time-dependent density functional theory: Theory and applications to diradicals , 2003 .

[35]  G. Seifert,et al.  Nonadiabatic dynamics within time-dependent density functional tight binding method. , 2009, The journal of physical chemistry. A.

[36]  M. Persico,et al.  Semiclassical simulations of azomethane photochemistry in the gas phase and in solution. , 2001, Journal of the American Chemical Society.

[37]  Shahnawaz Rafiq,et al.  Excited state relaxation dynamics of model green fluorescent protein chromophore analogs: evidence for cis-trans isomerism. , 2011, The journal of physical chemistry. A.

[38]  Toshifumi Mori,et al.  Conical intersections of free energy surfaces in solution: effect of electron correlation on a protonated Schiff base in methanol solution. , 2010, The Journal of chemical physics.

[39]  A. Zewail,et al.  Femtochemistry of trans-azomethane: a combined experimental and theoretical study. , 2003, Chemphyschem : a European journal of chemical physics and physical chemistry.

[40]  Tom Ziegler,et al.  Time-dependent density functional theory based on a noncollinear formulation of the exchange-correlation potential. , 2004, The Journal of chemical physics.

[41]  L. Cederbaum,et al.  Environmental effects on a conical intersection: a model study. , 2004, Faraday discussions.

[42]  Todd J Martínez,et al.  Optimizing conical intersections without derivative coupling vectors: application to multistate multireference second-order perturbation theory (MS-CASPT2). , 2008, The journal of physical chemistry. B.

[43]  Anna I. Krylov,et al.  A spin-complete version of the spin-flip approach to bond breaking: What is the impact of obtaining spin eigenfunctions? , 2003 .

[44]  S. Ruhman,et al.  Solvent tuning of a conical intersection: direct experimental verification of a theoretical prediction. , 2011, The journal of physical chemistry. A.

[45]  H. Lischka,et al.  The effect of hydration on the photo-deactivation pathways of 4-aminopyrimidine , 2010 .

[46]  O. Christiansen,et al.  Photoabsorption studies of neutral green fluorescent protein model chromophores in vacuo. , 2009, Physical chemistry chemical physics : PCCP.

[47]  G. Nienhaus.  The green fluorescent protein: a key tool to study chemical processes in living cells. , 2008, Angewandte Chemie.

[48]  Benjamin T. Miller,et al.  A parallel implementation of the analytic nuclear gradient for time-dependent density functional theory within the Tamm–Dancoff approximation , 1999 .

[49]  Michael J. Bearpark,et al.  A direct method for the location of the lowest energy point on a potential surface crossing , 1994 .

[50]  H. Ågren,et al.  Spin-flip time dependent density functional theory applied to excited states with single, double, or mixed electron excitation character. , 2010, The Journal of chemical physics.

[51]  I. Polyakov,et al.  Potential Energy Landscape of the Electronic States of the GFP Chromophore in Different Protonation Forms: Electronic Transition Energies and Conical Intersections. , 2010, Journal of chemical theory and computation.

[52]  Walter Thiel,et al.  Comparison of algorithms for conical intersection optimisation using semiempirical methods , 2007 .

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

[54]  Lars V. Schäfer,et al.  ADVANCES IN QUANTUM CHEMISTRY, VOL 59 , 2010 .

[55]  F. Buda,et al.  Bathochromic Shift in Green Fluorescent Protein: A Puzzle for QM/MM Approaches. , 2012, Journal of chemical theory and computation.

[56]  Mark S. Gordon,et al.  An effective fragment method for modeling solvent effects in quantum mechanical calculations , 1996 .

[57]  R. Mitrić,et al.  Time-dependent density functional theory excited state nonadiabatic dynamics combined with quantum mechanical/molecular mechanical approach: photodynamics of indole in water. , 2011, The Journal of chemical physics.

[58]  M. E. Casida.  Time-Dependent Density Functional Response Theory for Molecules , 1995 .

[59]  Mark S Gordon,et al.  Solvent effects on the S(N)2 reaction: Application of the density functional theory-based effective fragment potential method. , 2005, The journal of physical chemistry. A.

[60]  T. Ziegler,et al.  Introduction of a New Theory for the Calculation of Magnetic Coupling Based on Spin-Flip Constricted Variational Density Functional Theory. Application to Trinuclear Copper Complexes which Model the Native Intermediate in Multicopper Oxidases. , 2011, Journal of chemical theory and computation.

[61]  Wenjian Liu,et al.  Theoretical and numerical assessments of spin-flip time-dependent density functional theory. , 2012, The Journal of chemical physics.

[62]  Sean C. Smith,et al.  Bond selection in the photoisomerization reaction of anionic green fluorescent protein and kindling fluorescent protein chromophore models. , 2008, Journal of the American Chemical Society.

[63]  A. Krylov,et al.  Effect of protein environment on electronically excited and ionized states of the green fluorescent protein chromophore. , 2011, The journal of physical chemistry. B.

[64]  Todd J. Martínez,et al.  Conical intersections and double excitations in time-dependent density functional theory , 2006 .

[65]  R. Mitrić,et al.  Simulation of time resolved photoelectron spectra with Stieltjes imaging illustrated on ultrafast internal conversion in pyrazine. , 2010, The Journal of chemical physics.

[66]  Marco Garavelli,et al.  Structure of the intersection space associated with ZIE photoisomerization of retinal in rhodopsin proteins. , 2004, Faraday discussions.

[67]  Michael W. Schmidt,et al.  Solvent-induced shifts in electronic spectra of uracil. , 2011, The journal of physical chemistry. A.

[68]  S. Kato,et al.  Solvent effect on conical intersections in excited-state 9H-adenine: radiationless decay mechanism in polar solvent. , 2007, Journal of the American Chemical Society.

[69]  I. Tavernelli,et al.  Trajectory surface hopping within linear response time-dependent density-functional theory. , 2007, Physical review letters.

[70]  Pascal Krause,et al.  Nonadiabatic events and conical intersections. , 2011, Annual review of physical chemistry.

[71]  Mark S Gordon,et al.  Solvent-induced shift of the lowest singlet π → π* charge-transfer excited state of p-nitroaniline in water: an application of the TDDFT/EFP1 method. , 2011, The journal of physical chemistry. A.

[72]  M. Barbatti,et al.  Azomethane: nonadiabatic photodynamical simulations in solution. , 2010, The journal of physical chemistry. A.

[73]  Basile F. E. Curchod,et al.  Nonadiabatic molecular dynamics with solvent effects: A LR-TDDFT QM/MM study of ruthenium (II) tris (bipyridine) in water , 2011 .

[74]  D. Truhlar,et al.  Density functional study of multiplicity-changing valence and Rydberg excitations of p-block elements: delta self-consistent field, collinear spin-flip time-dependent density functional theory (DFT), and conventional time-dependent DFT. , 2011, The Journal of chemical physics.

[75]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[76]  E. Tapavicza,et al.  Unravelling the details of vitamin D photosynthesis by non-adiabatic molecular dynamics simulations. , 2011, Physical chemistry chemical physics : PCCP.

[77]  T. Van Voorhis,et al.  Communication: Conical intersections using constrained density functional theory-configuration interaction. , 2010, The Journal of chemical physics.

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

[79]  M. Barbatti,et al.  Nonadiabatic excited-state dynamics with hybrid ab initio quantum-mechanical/molecular-mechanical methods: solvation of the pentadieniminium cation in apolar media. , 2010, The journal of physical chemistry. A.

[80]  O. Sugino,et al.  A time-dependent density-functional approach to nonadiabatic electron-nucleus dynamics: formulation and photochemical application. , 2009, Physical chemistry chemical physics : PCCP.

[81]  Jingsong Wang,et al.  Electrostatic control of the regioselectivity in the photoisomerization of trans,trans-1-fluoro-2,4-hexadiene: evidence for competing conical intersections. , 2004, Journal of the American Chemical Society.

[82]  Hans Lischka,et al.  Excited-state properties and environmental effects for protonated schiff bases: a theoretical study. , 2006, Chemphyschem : a European journal of chemical physics and physical chemistry.

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

[84]  Potential energy surface crossings in organic photochemistry , 1996 .

[85]  Satoshi Maeda,et al.  Updated Branching Plane for Finding Conical Intersections without Coupling Derivative Vectors. , 2010, Journal of chemical theory and computation.

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

[87]  H. Ågren,et al.  Spin-flip time dependent density functional theory for singlet–triplet splittings in σ,σ-biradicals , 2010 .

[88]  M. Zimmer,et al.  Green fluorescent protein (GFP): applications, structure, and related photophysical behavior. , 2002, Chemical reviews.

[89]  Lars V. Schäfer,et al.  Computer simulations of photobiological processes: The effect of the protein environment , 2010 .

[90]  M. Persico,et al.  Diabatic and adiabatic potential-energy surfaces for azomethane photochemistry , 2000 .

[91]  Weitao Yang,et al.  Conical intersections in solution: formulation, algorithm, and implementation with combined quantum mechanics/molecular mechanics method. , 2011, The Journal of chemical physics.

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

[93]  T. Martínez,et al.  Conical Intersections in Solution: A QM/MM Study Using Floating Occupation Semiempirical Configuration Interaction Wave Functions , 2003 .

[94]  R. Mitrić,et al.  Nonadiabatic dynamics and simulation of time resolved photoelectron spectra within time-dependent density functional theory: Ultrafast photoswitching in benzylideneaniline. , 2008, The Journal of chemical physics.

[95]  C. Crespo-Hernández,et al.  Ultrafast excited-state dynamics in nucleic acids. , 2004, Chemical reviews.

[96]  Bingbing Suo,et al.  Spin-adapted open-shell time-dependent density functional theory. II. Theory and pilot application. , 2011, The Journal of chemical physics.

[97]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[98]  T. Martínez,et al.  Conical intersection dynamics in solution: the chromophore of Green Fluorescent Protein. , 2004, Faraday discussions.

[99]  2-Aminopurine non-radiative decay and emission in aqueous solution: A theoretical study , 2008 .

[100]  M. Gordon,et al.  Modeling Solvent Effects on Electronic Excited States , 2011 .

[101]  Hui Li,et al.  Analytic energy gradient in combined time-dependent density functional theory and polarizable force field calculation. , 2010, The Journal of chemical physics.

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

[103]  Angel Rubio,et al.  Time-dependent density-functional approach for biological chromophores: the case of the green fluorescent protein. , 2003, Physical review letters.

[104]  Tom Ziegler,et al.  The performance of time-dependent density functional theory based on a noncollinear exchange-correlation potential in the calculations of excitation energies. , 2005, The Journal of chemical physics.

[105]  M. Gordon,et al.  Solvent effects on optical properties of molecules: a combined time-dependent density functional theory/effective fragment potential approach. , 2008, The Journal of chemical physics.

[106]  T. Martínez,et al.  Protonic gating of excited-state twisting and charge localization in GFP chromophores: a mechanistic hypothesis for reversible photoswitching. , 2010, Journal of the American Chemical Society.

[107]  Mark S Gordon,et al.  Photoisomerization of stilbene: a spin-flip density functional theory approach. , 2011, The journal of physical chemistry. A.

[108]  P. Engel.  Mechanism of the thermal and photochemical decomposition of azoalkanes , 1980 .

[109]  Raymond Kapral,et al.  Quantum-classical description of environmental effects on electronic dynamics at conical intersections. , 2010, The Journal of chemical physics.

[110]  T. Vreven,et al.  Searching for conical intersections of potential energy surfaces with the ONIOM method: application to previtamin D. , 2008, The journal of physical chemistry. A.

[111]  M. Barbatti,et al.  Photodynamics of azomethane: a nonadiabatic surface-hopping study. , 2010, The journal of physical chemistry. A.

[112]  Wenjian Liu,et al.  Spin-adapted open-shell random phase approximation and time-dependent density functional theory. I. Theory. , 2010, The Journal of chemical physics.

[113]  Benjamin G. Levine,et al.  Isomerization through conical intersections. , 2007, Annual review of physical chemistry.