Non-equilibrium reaction and relaxation dynamics in a strongly interacting explicit solvent: F + CD3CN treated with a parallel multi-state EVB model.

We describe a parallelized linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM and TINKER. Forces are obtained using the Hellmann-Feynman relationship, giving continuous gradients, and good energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to explicitly correlated coupled cluster theory, we built a 64-state MS-EVB model designed to study the F + CD3CN → DF + CD2CN reaction in CD3CN solvent (recently reported in Dunning et al. [Science 347(6221), 530 (2015)]). This approach allows us to build a reactive potential energy surface whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We ran molecular dynamics simulations to examine a range of observables which follow in the wake of the reactive event: energy deposition in the nascent reaction products, vibrational relaxation rates of excited DF in CD3CN solvent, equilibrium power spectra of DF in CD3CN, and time dependent spectral shifts associated with relaxation of the nascent DF. Many of our results are in good agreement with time-resolved experimental observations, providing evidence for the accuracy of our MS-EVB framework in treating both the solute and solute/solvent interactions. The simulations provide additional insight into the dynamics at sub-picosecond time scales that are difficult to resolve experimentally. In particular, the simulations show that (immediately following deuterium abstraction) the nascent DF finds itself in a non-equilibrium regime in two different respects: (1) it is highly vibrationally excited, with ∼23 kcal mol(-1) localized in the stretch and (2) its post-reaction solvation environment, in which it is not yet hydrogen-bonded to CD3CN solvent molecules, is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of the post-reaction solvation environment results in a red shift. These two competing effects mean that the post-reaction relaxation profile is distinct from what is observed when Franck-Condon vibrational excitation of DF occurs within a microsolvation environment initially at equilibrium. Our conclusions, along with the theoretical and parallel software framework presented in this paper, should be more broadly applicable to a range of complex reactive systems.

[1]  Michael A. Bellucci,et al.  Molecular dynamics of excited state intramolecular proton transfer: 3-hydroxyflavone in solution. , 2012, The Journal of chemical physics.

[2]  G. Stock.  Classical simulation of quantum energy flow in biomolecules. , 2009, Physical review letters.

[3]  Takehiko Shimanouchi,et al.  Tables of molecular vibrational frequencies. Consolidated volume II , 1972 .

[4]  David R. Glowacki,et al.  Quantum dynamics simulations of energy redistribution in HO-SO2 , 2010 .

[5]  R. Hochstrasser,et al.  Vibrational relaxation dynamics in solutions. , 1994, Annual review of physical chemistry.

[6]  H. Bakker,et al.  Vibrational excitation induced proton transfer in hydrated Nafion membranes. , 2015, The journal of physical chemistry. B.

[7]  Arieh Warshel,et al.  An empirical valence bond approach for comparing reactions in solutions and in enzymes , 1980 .

[8]  P. H. Berens,et al.  Molecular Dynamics and Spectra. II. Diatomic Raman. , 1981 .

[9]  C. Elles,et al.  Connecting chemical dynamics in gases and liquids. , 2006, Annual review of physical chemistry.

[10]  P. Morse.  Diatomic Molecules According to the Wave Mechanics. II. Vibrational Levels , 1929 .

[11]  J. Polanyi.  Some Concepts in Reaction Dynamics , 1987, Science.

[12]  Romelia Salomón-Ferrer,et al.  Dynamics and dissipation in enzyme catalysis , 2011, Proceedings of the National Academy of Sciences.

[13]  N. Slater.  Classical Motion under a Morse Potential , 1957, Nature.

[14]  D. Dlott.  Vibrational energy redistribution in polyatomic liquids: 3D infrared–Raman spectroscopy , 2001 .

[15]  V. Sundström,et al.  Femtochemistry and femtobiology : ultrafast reaction dynamics at atomic-scale resolution : Nobel Symposium 101, Björkborn, Sweden, September 9-12, 1996 , 1997 .

[16]  M. Towrie,et al.  Vibrational relaxation and microsolvation of DF after F-atom reactions in polar solvents , 2015, Science.

[17]  William H. Miller,et al.  An empirical valence bond model for constructing global potential energy surfaces for chemical reactions of polyatomic molecular systems , 1990 .

[18]  G. M. Greetham,et al.  Comparing molecular photofragmentation dynamics in the gas and liquid phases. , 2013, Physical chemistry chemical physics : PCCP.

[19]  Todd J. Martinez,et al.  Graphical Processing Units for Quantum Chemistry , 2008, Computing in Science & Engineering.

[20]  G. M. Greetham,et al.  Vibrationally resolved dynamics of the reaction of Cl atoms with 2,3-dimethylbut-2-ene in chlorinated solvents , 2013 .

[21]  G. Schatz,et al.  Theoretical studies of polyatomic bimolecular reaction dynamics. , 1995, Annual review of physical chemistry.

[22]  G. Fleming,et al.  COMMUNICATIONS Femtosecond solvation dynamics in acetonitrile: Observation of the inertial contribution to the solvent response , 1991 .

[23]  A. Nitzan,et al.  Computing vibrational energy relaxation for high-frequency modes in condensed environments , 1997 .

[24]  David R Glowacki,et al.  Product energy deposition of CN + alkane H abstraction reactions in gas and solution phases. , 2011, The Journal of chemical physics.

[25]  Frederick R. Manby,et al.  Explicitly correlated coupled cluster methods with pair-specific geminals , 2011 .

[26]  Arieh Warshel,et al.  The empirical valence bond model: theory and applications , 2011 .

[27]  Gregory A. Voth,et al.  The computer simulation of proton transport in water , 1999 .

[28]  Gregory A. Voth,et al.  A second generation multistate empirical valence bond model for proton transport in aqueous systems , 2002 .

[29]  J. Hynes,et al.  VB resonance theory in solution. I. Multistate formulation , 1995 .

[30]  Michael D. Fayer,et al.  Ultrafast Infrared And Raman Spectroscopy , 2001 .

[31]  Robin M. Hochstrasser,et al.  Transition State Dynamics and Relaxation Processes in Solutions: A Frontier of Physical Chemistry , 1996 .

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

[33]  Stewart K. Reed,et al.  Classical, quantum and statistical simulations of vibrationally excited HOSO(2): IVR, dissociation, and implications for OH + SO(2) kinetics at high pressures. , 2009, Physical chemistry chemical physics : PCCP.

[34]  A. Tokmakoff,et al.  Local and collective reaction coordinates in the transport of the aqueous hydroxide ion. , 2014, The journal of physical chemistry. B.

[35]  Emanuele Paci,et al.  Boxed Molecular Dynamics: Decorrelation Time Scales and the Kinetic Master Equation. , 2011, Journal of chemical theory and computation.

[36]  D. Rapp,et al.  A REVIEW OF THE THEORY OF VIBRATIONAL ENERGY TRANSFER BETWEEN SIMPLE MOLECULES IN NONREACTIVE COLLISIONS. , 1968 .

[37]  R. Hochstrasser,et al.  Vibrational Relaxation of the CO Stretch Vibration in Hemoglobin-CO, Myoglobin-CO, and Protoheme-CO , 1995 .

[38]  D. Pines,et al.  Solvent-induced red-shifts for the proton stretch vibrational frequency in a hydrogen-bonded complex. 1. A valence bond-based theoretical approach. , 2014, The journal of physical chemistry. B.

[39]  J. Skinner,et al.  Vibrational energy relaxation of small molecules and ions in liquids , 2011 .

[40]  Martin Schütz,et al.  Molpro: a general‐purpose quantum chemistry program package , 2012 .

[41]  J. Harvey,et al.  Non-equilibrium phenomena and molecular reaction dynamics: mode space, energy space and conformer space , 2013 .

[42]  L. Andrews,et al.  Infrared spectra of hydrogen-bonded complexes between methyl cyanide and hydrogen fluoride in solid argon at 12 K , 1983 .

[43]  P. H. Berens,et al.  Molecular dynamics and spectra. I. Diatomic rotation and vibration , 1981 .

[44]  D. Case,et al.  Proton Transfer Studied Using a Combined Ab Initio Reactive Potential Energy Surface with Quantum Path Integral Methodology. , 2010, Journal of chemical theory and computation.

[45]  Yatsandra Oyola,et al.  Dynamics and the failure of transition state theory in alkene hydroboration. , 2009, Journal of the American Chemical Society.

[46]  Markus Meuwly,et al.  Multisurface Adiabatic Reactive Molecular Dynamics. , 2014, Journal of chemical theory and computation.

[47]  Aaron Sisto,et al.  Ab initio nonadiabatic dynamics of multichromophore complexes: a scalable graphical-processing-unit-accelerated exciton framework. , 2014, Accounts of chemical research.

[48]  T. Elsaesser,et al.  Ultrafast hydrogen bonding dynamics and proton transfer processes in the condensed phase , 2002 .

[49]  Kisam Park,et al.  Calculating Vibrational Energy Relaxation Rates from Classical Molecular Dynamics Simulations: Quantum Correction Factors for Processes Involving Vibration−Vibration Energy Transfer† , 2001 .

[50]  B. Carpenter,et al.  Do we fully understand what controls chemical selectivity? , 2011, Physical chemistry chemical physics : PCCP.

[51]  G. Voth,et al.  Proton conduction in exchange membranes across multiple length scales. , 2012, Accounts of chemical research.

[52]  V. Batista,et al.  Ultrafast branching of reaction pathways in 2-(2'-hydroxyphenyl)benzothiazole in polar acetonitrile solution. , 2011, The journal of physical chemistry. A.

[53]  Barbara Kirchner,et al.  Computing vibrational spectra from ab initio molecular dynamics. , 2013, Physical chemistry chemical physics : PCCP.

[54]  P. Wolynes,et al.  Vibrational energy flow and chemical reactions. , 2004, Accounts of chemical research.

[55]  M. Stei,et al.  Single solvent molecules can affect the dynamics of substitution reactions. , 2012, Nature chemistry.

[56]  S. Bradforth,et al.  Contrasting the excited state reaction pathways of phenol and para-methylthiophenol in the gas and liquid phases. , 2012, Faraday discussions.

[57]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[58]  T. J. Preston,et al.  Probing the photoisomerization of CHBr3 and CHI3 in solution with transient vibrational and electronic spectroscopy. , 2013, The journal of physical chemistry. A.

[59]  D. Truhlar,et al.  Combined Electrostatically Embedded Multiconfiguration Molecular Mechanics and Molecular Mechanical Method: Application to Molecular Dynamics Simulation of a Chemical Reaction in Aqueous Solution with Hybrid Density Functional Theory. , 2008, Journal of chemical theory and computation.

[60]  L. Raff.  Projection methods for obtaining intramolecular energy transfer rates from classical trajectory results: Application to 1,2‐difluoroethane , 1988 .

[61]  G. M. Greetham,et al.  Reaction dynamics of CN radicals with tetrahydrofuran in liquid solutions. , 2012, Physical chemistry chemical physics : PCCP.

[62]  Michael Towrie,et al.  Vibrationally Quantum-State–Specific Reaction Dynamics of H Atom Abstraction by CN Radical in Solution , 2011, Science.

[63]  James T. Yardley,et al.  Introduction to Molecular Energy Transfer , 1980 .

[64]  D. Glowacki,et al.  Chemical Reaction Dynamics in Liquid Solutions. , 2011, The journal of physical chemistry letters.

[65]  A. Orr-Ewing.  Perspective: Bimolecular chemical reaction dynamics in liquids. , 2014, The Journal of chemical physics.

[66]  S. Marsden,et al.  Alkene hydroboration: hot intermediates that react while they are cooling. , 2010, Journal of the American Chemical Society.

[67]  S. Bradforth,et al.  Linking photochemistry in the gas and solution phase: S-H bond fission in p-methylthiophenol following UV photoexcitation. , 2011, Faraday discussions.

[68]  Emanuele Paci,et al.  Boxed molecular dynamics: a simple and general technique for accelerating rare event kinetics and mapping free energy in large molecular systems. , 2009, The journal of physical chemistry. B.

[69]  A. Tokmakoff,et al.  TEMPERATURE-DEPENDENT VIBRATIONAL RELAXATION IN POLYATOMIC LIQUIDS : PICOSECOND INFRARED PUMP-PROBE EXPERIMENTS , 1994 .

[70]  D. W. Noid,et al.  A spectral analysis method of obtaining molecular spectra from classical trajectories , 1977 .

[71]  S. Marsden,et al.  Significance of nonstatistical dynamics in organic reaction mechanisms: time-dependent stereoselectivity in cyclopentyne-alkene cycloadditions. , 2009, Journal of the American Chemical Society.

[72]  S. Okazaki,et al.  A surface hopping method for chemical reaction dynamics in solution described by diabatic representation: an analysis of tunneling and thermal activation. , 2006, The Journal of chemical physics.

[73]  A. Hagler,et al.  Derivation of force fields for molecular mechanics and dynamics from ab initio energy surfaces. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[74]  J. Harvey,et al.  Ultrafast energy flow in the wake of solution-phase bimolecular reactions. , 2011, Nature chemistry.

[75]  M. Fayer.  Vibrational Echo Chemical Exchange Spectroscopy , 2013 .

[76]  J. Hynes.  Chemical Reaction Dynamics in Solution , 1985 .

[77]  David R Glowacki,et al.  Prediction of enhanced solvent-induced enantioselectivity for a ring opening with a bifurcating reaction path. , 2015, Physical chemistry chemical physics : PCCP.

[78]  J. Hynes,et al.  Activation to the transition state: reactant and solvent energy flow for a model SN2 reaction in water. , 1991, Journal of the American Chemical Society.

[79]  H. Schlegel,et al.  Empirical Valence-Bond Models for Reactive Potential Energy Surfaces Using Distributed Gaussians. , 2006, Journal of chemical theory and computation.

[80]  Lawrence M. Goldman,et al.  Nonstatistical dynamics in unlikely places: [1,5] hydrogen migration in chemically activated cyclopentadiene. , 2011, Journal of the American Chemical Society.