1
|
Lin S, Peng D, Yang W, Gu FL, Lan Z. Theoretical studies on triplet-state driven dissociation of formaldehyde by quasi-classical molecular dynamics simulation on machine-learning potential energy surface. J Chem Phys 2021; 155:214105. [PMID: 34879677 PMCID: PMC8654486 DOI: 10.1063/5.0067176] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2021] [Accepted: 11/09/2021] [Indexed: 11/15/2022] Open
Abstract
The H-atom dissociation of formaldehyde on the lowest triplet state (T1) is studied by quasi-classical molecular dynamic simulations on the high-dimensional machine-learning potential energy surface (PES) model. An atomic-energy based deep-learning neural network (NN) is used to represent the PES function, and the weighted atom-centered symmetry functions are employed as inputs of the NN model to satisfy the translational, rotational, and permutational symmetries, and to capture the geometry features of each atom and its individual chemical environment. Several standard technical tricks are used in the construction of NN-PES, which includes the application of clustering algorithm in the formation of the training dataset, the examination of the reliability of the NN-PES model by different fitted NN models, and the detection of the out-of-confidence region by the confidence interval of the training dataset. The accuracy of the full-dimensional NN-PES model is examined by two benchmark calculations with respect to ab initio data. Both the NN and electronic-structure calculations give a similar H-atom dissociation reaction pathway on the T1 state in the intrinsic reaction coordinate analysis. The small-scaled trial dynamics simulations based on NN-PES and ab initio PES give highly consistent results. After confirming the accuracy of the NN-PES, a large number of trajectories are calculated in the quasi-classical dynamics, which allows us to get a better understanding of the T1-driven H-atom dissociation dynamics efficiently. Particularly, the dynamics simulations from different initial conditions can be easily simulated with a rather low computational cost. The influence of the mode-specific vibrational excitations on the H-atom dissociation dynamics driven by the T1 state is explored. The results show that the vibrational excitations on symmetric C-H stretching, asymmetric C-H stretching, and C=O stretching motions always enhance the H-atom dissociation probability obviously.
Collapse
Affiliation(s)
| | | | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Feng Long Gu
- Authors to whom correspondence should be addressed: and
| | - Zhenggang Lan
- Authors to whom correspondence should be addressed: and
| |
Collapse
|
2
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 154] [Impact Index Per Article: 51.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
3
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
4
|
Abstract
Theoretical simulations of electronic excitations and associated processes in molecules are indispensable for fundamental research and technological innovations. However, such simulations are notoriously challenging to perform with quantum mechanical methods. Advances in machine learning open many new avenues for assisting molecular excited-state simulations. In this Review, we track such progress, assess the current state of the art and highlight the critical issues to solve in the future. We overview a broad range of machine learning applications in excited-state research, which include the prediction of molecular properties, improvements of quantum mechanical methods for the calculations of excited-state properties and the search for new materials. Machine learning approaches can help us understand hidden factors that influence photo-processes, leading to a better control of such processes and new rules for the design of materials for optoelectronic applications.
Collapse
|
5
|
Alborzpour JP, Tew DP, Habershon S. Efficient and accurate evaluation of potential energy matrix elements for quantum dynamics using Gaussian process regression. J Chem Phys 2016; 145:174112. [DOI: 10.1063/1.4964902] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jonathan P. Alborzpour
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - David P. Tew
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
6
|
Kim CW, Rhee YM. Constructing an Interpolated Potential Energy Surface of a Large Molecule: A Case Study with Bacteriochlorophyll a Model in the Fenna-Matthews-Olson Complex. J Chem Theory Comput 2016; 12:5235-5246. [PMID: 27760297 DOI: 10.1021/acs.jctc.6b00647] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Constructing a reliable potential energy surface (PES) is a key step toward computationally studying the chemical dynamics of any molecular system. The interpolation scheme is a useful tool that can closely follow the accuracy of quantum chemical means at a dramatically reduced computational cost. However, applying interpolation to building a PES of a large molecule is not a straightforward black-box approach, as it frequently encounters practical difficulties associated with its large dimensionality. Here, we present detailed courses of applying interpolation toward building a PES of a large chromophore molecule. We take the example of S0 and S1 electronic states of bacteriochlorophyll a (BChla) molecules in the Fenna-Matthews-Olson light harvesting complex. With a reduced model molecule that bears BChla's main π-conjugated ring, various practical approaches are designed for improving the PES quality in a stable manner and for fine-tuning the final surface such that the surface can be adopted for long time molecular dynamics simulations. Combined with parallel implementation, we show that interpolated mechanics/molecular mechanics (IM/MM) simulations of the entire complex in the nanosecond time scale can be conducted readily without any practical issues. With 1500 interpolation data points for each chromophore unit, the PES error relative to the reference quantum chemical calculation is found to be ∼0.15 eV in the thermally accessible region of the conformational space, together with ∼0.01 eV error in S0 - S1 transition energies. The performance issue related to the use of a large interpolation database within the framework of our parallel routines is also discussed.
Collapse
Affiliation(s)
- Chang Woo Kim
- Center for Self-assembly and Complexity, Institute for Basic Science (IBS) , Pohang 37673, Korea.,Department of Chemistry, Pohang University of Science and Technology (POSTECH) , Pohang 37673, Korea
| | - Young Min Rhee
- Center for Self-assembly and Complexity, Institute for Basic Science (IBS) , Pohang 37673, Korea.,Department of Chemistry, Pohang University of Science and Technology (POSTECH) , Pohang 37673, Korea
| |
Collapse
|
7
|
An overview of nonadiabatic dynamics simulations methods, with focus on the direct approach versus the fitting of potential energy surfaces. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1526-1] [Citation(s) in RCA: 102] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
8
|
Lin SY, Zhang P, Zhang JZ. Hybrid many-body-expansion/Shepard-interpolation method for constructing ab initio potential energy surfaces for quantum dynamics calculations. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2012.11.089] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
9
|
Li Q, Chen XY, Qiu YX, Wang SG. Investigation of spin-flip reactions of Zr + CH3CN by relativistic density functional theory. J Phys Chem A 2012; 116:5019-25. [PMID: 22578173 DOI: 10.1021/jp211975e] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To explore the details of the reaction mechanisms of Zr atoms with acetonitrile molecules, the triplet and singlet spin-state potential energy surfaces have been investigated. Density functional theory (DFT) with the relativistic zero-order regular approximation at the PW91/TZ2P level has been applied. The complicated minimum energy reaction path involves four transition states (TS), stationary states 1-5 and one spin inversion (indicated by ⇒): (3)Zr + NCCH(3) → (3)Zr-η(1)-NCCH(3) ((3)1) → (3)TS(1/2) → (3)Zr-η(2)-(NC)CH(3) ((3)2) → (3)TS(2/3) → (3)ZrH-η(3)-(NCCH(2)) ((3)3) → (3)TS(3/4) → CNZrCH(3) ((3)4) ⇒ (1)TS(4/5) → CN(ZrH)CH(2) ((1)5). The minimum energy crossing point was determined with the help of the DFT fractional-occupation-number approach. The spin inversion leading from the triplet to the singlet state facilitates the activation of a C-H bond, lowering the rearrangement-barrier by 78 kJ/mol. The overall reaction is calculated to be exothermic by about 296 kJ/mol. All intermediate and product species were frequency and NBO analyzed. The species can be rationalized with the help of Lewis type formulas.
Collapse
Affiliation(s)
- Qiang Li
- School of Chemistry and Chemical Engineering, Shanghai Jiao Tong University, 200240 Shanghai, China
| | | | | | | |
Collapse
|
10
|
Kalinowski J, Räsänen M, Gerber RB. Mechanism and electronic transition in the reaction: On the fly dynamics simulations with multi-reference potentials. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.03.068] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
11
|
DOLTSINIS NIKOSL, MARX DOMINIK. FIRST PRINCIPLES MOLECULAR DYNAMICS INVOLVING EXCITED STATES AND NONADIABATIC TRANSITIONS. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2012. [DOI: 10.1142/s0219633602000257] [Citation(s) in RCA: 126] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Extensions of traditional molecular dynamics to excited electronic states and non-Born–Oppenheimer dynamics are reviewed focusing on applicability to chemical reactions of large molecules, possibly in condensed phases. The latter imposes restrictions on both the level of accuracy of the underlying electronic structure theory and the treatment of nonadiabaticity. This review, therefore, exclusively deals with ab initio "on the fly" molecular dynamics methods. For the same reason, mainly mixed quantum-classical approaches to nonadiabatic dynamics are considered.
Collapse
Affiliation(s)
- NIKOS L. DOLTSINIS
- Lehrstuhl für Theoretische Chemie, Ruhr–Universität Bochum, 44780 Bochum, Germany
| | - DOMINIK MARX
- Lehrstuhl für Theoretische Chemie, Ruhr–Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
12
|
Evenhuis C, Martínez TJ. A scheme to interpolate potential energy surfaces and derivative coupling vectors without performing a global diabatization. J Chem Phys 2011; 135:224110. [DOI: 10.1063/1.3660686] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
13
|
PARK SEUNGC, BRAAMS BASTIAANJ, BOWMAN JOELM. CONSTRUCTION OF A GLOBAL POTENTIAL ENERGY SURFACE FROM NOVEL AB INITIO MOLECULAR DYNAMICS FOR THE O(3P) + C3H3 REACTION. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2011. [DOI: 10.1142/s0219633605001350] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
We present a global potential energy surface (PES) for the 2 A state of the O (3 P ) + C 3 H 3 radical reaction. The global PES is constructed mainly using direct ab initio molecular dynamics and further sampling is done using the Diffusion Monte Carlo method. The potential is fully invariant with respect to permutational symmetry of like atoms. Special techniques, based on invariant theory of finite groups, have been used to develop basis functions for fitting that display this symmetry. The resultant potential energy surface shows multiple reaction paths with six different product channels. The products of the reactions are CO + C 2 H 3 radicals H + C 3 H 2 O radicals (with two isomers, propynal and propa-1,2-dien-1-one) and OH + C 3 H 2 radicals (with three isomers, vinylidenecarbene, propargylene and cyclopropenylidene). Energies of the PES are in excellent agreement with ab initio energies for each stationary point, the reactants and the products. Most stationary points are fitted at the sub Kcal/mol level. The global potential surface represents all the stationary points and six different product channels correctly. Preliminary dynamics calculations show abstraction and insertion mechanisms for the OH + vinylidenecarbene channel and the H + propynal channel, respectively.
Collapse
Affiliation(s)
- SEUNG C. PARK
- Department of Chemistry, Sungkyunkwan University, Suwon 440-746 KOREA and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, GA 30322, USA
| | - BASTIAAN J. BRAAMS
- Department of Mathematics and Computer Science and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, GA 30322, USA
| | - JOEL M. BOWMAN
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, GA 30322, USA
| |
Collapse
|
14
|
Isotope Effects in Photodissociation: Chemical Reaction Dynamics and Implications for Atmospheres. ADVANCES IN QUANTUM CHEMISTRY 2008. [DOI: 10.1016/s0065-3276(07)00207-9] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
15
|
Chen X, Batista VS. The MP/SOFT methodology for simulations of quantum dynamics: Model study of the photoisomerization of the retinyl chromophore in visual rhodopsin. J Photochem Photobiol A Chem 2007. [DOI: 10.1016/j.jphotochem.2007.05.015] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
16
|
Fernandez-Ramos A, Miller JA, Klippenstein SJ, Truhlar DG. Modeling the kinetics of bimolecular reactions. Chem Rev 2007; 106:4518-84. [PMID: 17091928 DOI: 10.1021/cr050205w] [Citation(s) in RCA: 474] [Impact Index Per Article: 27.9] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Antonio Fernandez-Ramos
- Departamento de Quimica Fisica, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | | | | | | |
Collapse
|
17
|
Chen X, Batista VS. Matching-pursuit/split-operator-Fourier-transform simulations of excited-state nonadiabatic quantum dynamics in pyrazine. J Chem Phys 2006; 125:124313. [PMID: 17014180 DOI: 10.1063/1.2356477] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A simple approach for numerically exact simulations of nonadiabatic quantum dynamics in multidimensional systems is introduced and applied to the description of the photoabsorption spectroscopy of pyrazine. The propagation scheme generalizes the recently developed matching-pursuit/split-operator-Fourier-transform (MP/SOFT) method [Y. Wu and V. S. Batista, J. Chem. Phys. 121, 1676 (2004)] to simulations of nonadiabatic quantum dynamics. The time-evolution operator is applied, as defined by the Trotter expansion to second order accuracy, in dynamically adaptive coherent-state expansions. These representations are obtained by combining the matching-pursuit algorithm with a gradient-based optimization method. The accuracy and efficiency of the resulting computational approach are demonstrated in calculations of time-dependent survival amplitudes and photoabsorption cross sections, using a model Hamiltonian that allows for direct comparisons with benchmark calculations. Simulations in full-dimensional potential energy surfaces involve the propagation of a 24-dimensional wave packet to describe the S(1)S(2) interconversion of pyrazine after S(0)-->S(2) photoexcitation. The reported results show that the generalized MP/SOFT method is a practical and accurate approach to model nonadiabatic reaction dynamics in polyatomic systems.
Collapse
Affiliation(s)
- Xin Chen
- Department of Chemistry, Yale University, New Haven, CT 06520-8107, USA
| | | |
Collapse
|
18
|
Wu Y, Herman MF, Batista VS. Matching-pursuit∕split-operator Fourier-transform simulations of nonadiabatic quantum dynamics. J Chem Phys 2005; 122:114114. [PMID: 15836208 DOI: 10.1063/1.1881132] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A rigorous and practical approach for simulations of nonadiabatic quantum dynamics is introduced. The algorithm involves a natural extension of the matching-pursuitsplit-operator Fourier-transform (MPSOFT) method [Y. Wu and V. S. Batista, J. Chem. Phys. 121, 1676 (2004)] recently developed for simulations of adiabatic quantum dynamics in multidimensional systems. The MPSOFT propagation scheme, extended to nonadiabatic dynamics, recursively applies the time-evolution operator as defined by the standard perturbation expansion to first-, or second-order, accuracy. The expansion is implemented in dynamically adaptive coherent-state representations, generated by an approach that combines the matching-pursuit algorithm with a gradient-based optimization method. The accuracy and efficiency of the resulting propagation method are demonstrated as applied to the canonical model systems introduced by Tully for testing simulations of dual curve-crossing nonadiabatic dynamics.
Collapse
Affiliation(s)
- Yinghua Wu
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118, USA
| | | | | |
Collapse
|
19
|
Chen X, Wu Y, Batista VS. Matching-pursuit/split-operator-Fourier-transform computations of thermal correlation functions. J Chem Phys 2005; 122:064102. [PMID: 15740362 DOI: 10.1063/1.1848513] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A rigorous and practical methodology for evaluating thermal-equilibrium density matrices, finite-temperature time-dependent expectation values, and time-correlation functions is described. The method involves an extension of the matching-pursuit/split-operator-Fourier-transform method to the solution of the Bloch equation via imaginary-time propagation of the density matrix and the evaluation of Heisenberg time-evolution operators through real-time propagation in dynamically adaptive coherent-state representations.
Collapse
Affiliation(s)
- Xin Chen
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, USA
| | | | | |
Collapse
|
20
|
Evenhuis CR, Collins MA. Interpolation of diabatic potential energy surfaces. J Chem Phys 2004; 121:2515-27. [PMID: 15281848 DOI: 10.1063/1.1770756] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
A method is presented for constructing diabatic potential energy matrices from ab initio quantum chemistry data. The method is similar to that reported previously for single adiabatic potential energy surfaces, but correctly accounts for the nuclear permutation symmetry of diabatic potential energy matrices and other complications that arise from the derivative coupling of electronic states. The method is tested by comparison with an analytic model for the two lowest energy states of H(3).
Collapse
Affiliation(s)
- Christian R Evenhuis
- Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia
| | | |
Collapse
|
21
|
|
22
|
|
23
|
|
24
|
Quenneville J, Ben-Nun M, Martı́nez TJ. Photochemistry from first principles — advances and future prospects. J Photochem Photobiol A Chem 2001. [DOI: 10.1016/s1010-6030(01)00452-x] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
25
|
Hack MD, Wensmann AM, Truhlar DG, Ben-Nun M, Martı́nez TJ. Comparison of full multiple spawning, trajectory surface hopping, and converged quantum mechanics for electronically nonadiabatic dynamics. J Chem Phys 2001. [DOI: 10.1063/1.1377030] [Citation(s) in RCA: 126] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
26
|
Soo Kim M, Tae Park S, June Sung B, Hee Moon J. Scaling of classical rate constants on scaled potential-energy surfaces. J Chem Phys 2001. [DOI: 10.1063/1.1374938] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
27
|
Hayes RL, Fattal E, Govind N, Carter EA. Long live vinylidene! A new view of the H(2)C=C: --> HC triple bond CH rearrangement from ab initio molecular dynamics. J Am Chem Soc 2001; 123:641-57. [PMID: 11456576 DOI: 10.1021/ja000907x] [Citation(s) in RCA: 73] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present complete active space self-consistent field (CASSCF) ab initio molecular dynamics (AIMD) simulations of the preparation of the metastable species vinylidene, and its subsequent, highly exothermic isomerization to acetylene, via electron removal from vinylidene anion (D(2)C=C(-) --> D(2)C=C: --> DC triple bond CD). After equilibrating vinylidene anion-d(2) at either 600 +/- 300 K (slightly below the isomerization barrier) or 1440 K +/- 720 K (just above the isomerization barrier), we remove an electron to form a vibrationally excited singlet vinylidene-d(2) and follow its dynamical evolution for 1.0 ps. Remarkably, we find that none of the vinylidenes equilibrated at 600 K and only 20% of the vinylidenes equilibrated at 1440 K isomerized, suggesting average lifetimes >1 ps for vibrationally excited vinylidene-d(2). Since the anion and neutral vinylidene are structurally similar, and yet extremely different geometrically from the isomerization transition state (TS), neutral vinylidene is not formed near the TS so that it must live until it has sufficient instantaneous kinetic energy in the correct vibrational mode(s). The origin of the delay is explained via both orbital rearrangement and intramolecular vibrational energy redistribution (IVR) effects. Unique signatures of the isomerization dynamics are revealed in the anharmonic vibrational frequencies extracted from the AIMD, which should be observable by ultrafast vibrational spectroscopy and in fact are consistent with currently available experimental spectra. Most interestingly, of those trajectories that did isomerize, every one of them violated conventional transition-state theory by recrossing back to vinylidene multiple times, against conventional notions that expect highly exothermic reactions to be irreversible. The dynamical motion responsible for the multiple barrier recrossings involves strong mode-coupling between the vinylidene CD(2) rock and a local acetylene DCC bend mode that has been recently observed experimentally. The multiple barrier recrossings can be used, via a generalized definition of lifetime, to reconcile extremely disparate experimental estimates of vinylidene's lifetime (differing by at least 6 orders of magnitude). Last, a caveat: These results are constrained by the approximations inherent in the simulation (classical nuclear motion, neglect of rotation-vibration coupling, and restriction to C(s) symmetry); refinement of these predictions may be necessary when more exact simulations someday become feasible.
Collapse
Affiliation(s)
- R L Hayes
- Contribution from the Department of Chemistry and Biochemistry, Box 951569, University of California, Los Angeles, California 90095-1569, USA
| | | | | | | |
Collapse
|
28
|
Yu HG, Nyman G. Interpolatedab initioquantum scattering for the reaction of OH with HCl. J Chem Phys 2000. [DOI: 10.1063/1.1319999] [Citation(s) in RCA: 38] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
29
|
Sung BJ, Kim MS. Accuracy of the energy partitioning data obtained by classical trajectory calculations on potential energy surfaces constructed by interpolation: H2CO→H2 + CO as an example. J Chem Phys 2000. [DOI: 10.1063/1.1287326] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
|
30
|
Volobuev YL, Hack MD, Topaler MS, Truhlar DG. Continuous surface switching: An improved time-dependent self-consistent-field method for nonadiabatic dynamics. J Chem Phys 2000. [DOI: 10.1063/1.481609] [Citation(s) in RCA: 101] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
31
|
Ben-Nun M, Quenneville J, Martínez TJ. Ab Initio Multiple Spawning: Photochemistry from First Principles Quantum Molecular Dynamics. J Phys Chem A 2000. [DOI: 10.1021/jp994174i] [Citation(s) in RCA: 644] [Impact Index Per Article: 26.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- M. Ben-Nun
- Department of Chemistry and The Beckman Institute, University of Illinois, Urbana, Illinois 61801
| | - Jason Quenneville
- Department of Chemistry and The Beckman Institute, University of Illinois, Urbana, Illinois 61801
| | - Todd J. Martínez
- Department of Chemistry and The Beckman Institute, University of Illinois, Urbana, Illinois 61801
| |
Collapse
|
32
|
Kobrak MN, Bittner ER. A quantum molecular dynamics study of exciton self-trapping in conjugated polymers: Temperature dependence and spectroscopy. J Chem Phys 2000. [DOI: 10.1063/1.481379] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
|
33
|
A generalized discrete variable representation approach to interpolating or fitting potential energy surfaces. Chem Phys Lett 2000. [DOI: 10.1016/s0009-2614(00)00359-6] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
34
|
Kobrak MN, Bittner ER. A dynamic model for exciton self-trapping in conjugated polymers. I. Theory. J Chem Phys 2000. [DOI: 10.1063/1.481109] [Citation(s) in RCA: 17] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
35
|
Kobrak MN, Bittner ER. A dynamic model for exciton self-trapping in conjugated polymers. II. Implementation. J Chem Phys 2000. [DOI: 10.1063/1.481126] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
|
36
|
Ben-Nun M, Martínez TJ. Electronic Absorption and Resonance Raman Spectroscopy from Ab Initio Quantum Molecular Dynamics. J Phys Chem A 1999. [DOI: 10.1021/jp992197r] [Citation(s) in RCA: 54] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- M. Ben-Nun
- Department of Chemistry and The Beckman Institute, University of Illinois, Urbana, Illinois 61801
| | - Todd J. Martínez
- Department of Chemistry and The Beckman Institute, University of Illinois, Urbana, Illinois 61801
| |
Collapse
|
37
|
Micha DA. Time-Dependent Many-Electron Treatment of Electronic Energy and Charge Transfer in Atomic Collisions. J Phys Chem A 1999. [DOI: 10.1021/jp9906839] [Citation(s) in RCA: 62] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- David A. Micha
- Quantum Theory Project, Departments of Chemistry and of Physics, University of Florida, Gainesville, Florida 32611-8435
| |
Collapse
|
38
|
Lin SY, Park SC, Kim MS. Construction of an accurate potential energy surface by interpolation for quantum dynamics studies of a three-body system. J Chem Phys 1999. [DOI: 10.1063/1.479681] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
39
|
Ben-Nun M, Martı́nez TJ. Exploiting temporal nonlocality to remove scaling bottlenecks in nonadiabatic quantum dynamics. J Chem Phys 1999. [DOI: 10.1063/1.478297] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
|