1
|
Toldo JM, Mattos RS, Pinheiro M, Mukherjee S, Barbatti M. Recommendations for Velocity Adjustment in Surface Hopping. J Chem Theory Comput 2024; 20:614-624. [PMID: 38207213 DOI: 10.1021/acs.jctc.3c01159] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2024]
Abstract
This study investigates velocity adjustment directions after hopping in surface hopping dynamics. Using fulvene and a protonated Schiff base (PSB4) as case studies, we investigate the population decay and reaction yields of different sets of dynamics with the velocity adjusted in either the nonadiabatic coupling, gradient difference, or momentum directions. For the latter, in addition to the conventional algorithm, we investigated the performance of a reduced kinetic energy reservoir approach recently proposed. Our evaluation also considered velocity adjustment in the directions of approximate nonadiabatic coupling vectors. While results for fulvene are susceptible to the adjustment approach, PSB4 is not. We correlated this dependence to the topography near the conical intersections. When nonadiabatic coupling vectors are unavailable, the gradient difference direction is the best adjustment option. If the gradient difference is also unavailable, a semiempirical vector direction or the momentum direction with a reduced kinetic energy reservoir becomes an excellent option to prevent an artificial excess of back hoppings. The precise velocity adjustment direction is less crucial for describing the nonadiabatic dynamics than the kinetic energy reservoir's size.
Collapse
Affiliation(s)
- Josene M Toldo
- Aix-Marseille University, CNRS, ICR, Marseille 13397, France
| | - Rafael S Mattos
- Aix-Marseille University, CNRS, ICR, Marseille 13397, France
| | - Max Pinheiro
- Aix-Marseille University, CNRS, ICR, Marseille 13397, France
| | | | - Mario Barbatti
- Aix-Marseille University, CNRS, ICR, Marseille 13397, France
- Institut Universitaire de France, Paris 75231, France
| |
Collapse
|
2
|
Avagliano D, Lorini E, González L. Sampling effects in quantum mechanical/molecular mechanics trajectory surface hopping non-adiabatic dynamics. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200381. [PMID: 35341304 PMCID: PMC8958275 DOI: 10.1098/rsta.2020.0381] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Accepted: 09/01/2021] [Indexed: 05/29/2023]
Abstract
The impact of different initial conditions in non-adiabatic trajectory surface hopping dynamics within a hybrid quantum mechanical/molecular mechanics scheme is investigated. The influence of a quantum sampling, based on a Wigner distribution, a fully thermal sampling, based on classical molecular dynamics, and a quantum sampled system, but thermally equilibrated with the environment, is investigated on the relaxation dynamics of solvated fulvene after light irradiation. We find that the decay from the first singlet excited state to the ground state shows high dependency on the initial condition and simulation parameters. The three sampling methods lead to different distributions of initial geometries and momenta, which then affect the fate of the excited state dynamics. We evaluated both the effect of sampling geometries and momenta, analysing how the ultrafast decay of fulvene changes accordingly. The results are expected to be of interest to decide how to initialize non-adiabatic dynamics in the presence of the environment. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
- Davide Avagliano
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria
| | - Emilio Lorini
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria
| | - Leticia González
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Straße 17, A-1180 Vienna, Austria
| |
Collapse
|
3
|
T. do Casal M, Toldo JM, Pinheiro Jr M, Barbatti M. Fewest switches surface hopping with Baeck-An couplings. OPEN RESEARCH EUROPE 2022; 1:49. [PMID: 37645211 PMCID: PMC10446015 DOI: 10.12688/openreseurope.13624.2] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/01/2022] [Indexed: 08/31/2023]
Abstract
In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.
Collapse
|
4
|
T. do Casal M, Toldo JM, Pinheiro Jr M, Barbatti M. Fewest switches surface hopping with Baeck-An couplings. OPEN RESEARCH EUROPE 2022; 1:49. [PMID: 37645211 PMCID: PMC10446015 DOI: 10.12688/openreseurope.13624.1] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/01/2022] [Indexed: 03/12/2024]
Abstract
In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.
Collapse
|
5
|
Lassmann Y, Curchod BFE. AIMSWISS-Ab initio multiple spawning with informed stochastic selections. J Chem Phys 2021; 154:211106. [PMID: 34240975 DOI: 10.1063/5.0052118] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
Ab initio multiple spawning (AIMS) offers a reliable strategy to describe the excited-state dynamics and nonadiabatic processes of molecular systems. AIMS represents nuclear wavefunctions as linear combinations of traveling, coupled Gaussians called trajectory basis functions (TBFs) and uses a spawning algorithm to increase as needed the size of this basis set during nonadiabatic transitions. While the success of AIMS resides in this spawning algorithm, the dramatic increase in TBFs generated by multiple crossings between electronic states can rapidly lead to intractable dynamics. In this Communication, we introduce a new flavor of AIMS, coined ab initio multiple spawning with informed stochastic selections (AIMSWISS), which proposes a parameter-free strategy to beat the growing number of TBFs in an AIMS dynamics while preserving its accurate description of nonadiabatic transitions. The performance of AIMSWISS is validated against the photodynamics of ethylene, cyclopropanone, and fulvene. This technique, built upon the recently developed stochastic-selection AIMS, is intended to serve as a computationally affordable starting point for multiple spawning simulations.
Collapse
Affiliation(s)
- Yorick Lassmann
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
6
|
Ibele LM, Lassmann Y, Martínez TJ, Curchod BFE. Comparing (stochastic-selection) ab initio multiple spawning with trajectory surface hopping for the photodynamics of cyclopropanone, fulvene, and dithiane. J Chem Phys 2021; 154:104110. [DOI: 10.1063/5.0045572] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Yorick Lassmann
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry, Stanford University, Stanford, California 94305, USA and PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
7
|
Tran T, Jenkins AJ, Worth GA, Robb MA. The quantum-Ehrenfest method with the inclusion of an IR pulse: Application to electron dynamics of the allene radical cation. J Chem Phys 2020; 153:031102. [PMID: 32716173 DOI: 10.1063/5.0015937] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We describe the implementation of a laser control pulse in the quantum-Ehrenfest method, a molecular quantum dynamics method that solves the time-dependent Schrödinger equation for both electrons and nuclei. The oscillating electric field-dipole interaction is incorporated directly in the one-electron Hamiltonian of the electronic structure part of the algorithm. We then use the coupled electron-nuclear dynamics of the π-system in the allene radical cation (•CH2=C=CH2)+ as a simple model of a pump-control experiment. We start (pump) with a two-state superposition of two cationic states. The resulting electron dynamics corresponds to the rapid oscillation of the unpaired electron between the two terminal methylenes. This electron dynamics is, in turn, coupled to the torsional motion of the terminal methylenes. There is a conical intersection at 90° twist, where the electron dynamics collapses because the adiabatic states become degenerate. After passing the conical intersection, the electron dynamics revives. The IR pulse (control) in our simulations is timed to have its maximum at the conical intersection. Our simulations show that the effect of the (control) pulse is to change the electron dynamics at the conical intersection and, as a consequence, the concomitant nuclear dynamics, which is dominated by the change in the torsional angle.
Collapse
Affiliation(s)
- Thierry Tran
- Department of Chemistry, University College London, 20, Gordon St., WC1H 0AJ London, United Kingdom
| | - Andrew J Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Graham A Worth
- Department of Chemistry, University College London, 20, Gordon St., WC1H 0AJ London, United Kingdom
| | - Michael A Robb
- Department of Chemistry, Molecular Sciences Research Hub, Imperial College London, White City Campus, 80 Wood Lane, W12 0BZ London, United Kingdom
| |
Collapse
|
8
|
Ibele LM, Curchod BFE. A molecular perspective on Tully models for nonadiabatic dynamics. Phys Chem Chem Phys 2020; 22:15183-15196. [DOI: 10.1039/d0cp01353f] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
We present a series of standardized molecular tests for nonadiabatic dynamics, reminiscent of the one-dimensional Tully models proposed in 1990.
Collapse
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry
- Durham University
- Durham DH1 3LE
- UK
| | | |
Collapse
|
9
|
Ibele LM, Nicolson A, Curchod BFE. Excited-state dynamics of molecules with classically driven trajectories and Gaussians. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1665199] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry, Durham University, Durham, UK
| | | | | |
Collapse
|
10
|
Xie W, Sapunar M, Došlić N, Sala M, Domcke W. Assessing the performance of trajectory surface hopping methods: Ultrafast internal conversion in pyrazine. J Chem Phys 2019; 150:154119. [DOI: 10.1063/1.5084961] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Weiwei Xie
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
- Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Marin Sapunar
- Department of Physical Chemistry, Ruđer Bošković Institute, HR-10000 Zagreb, Croatia
| | - Nađa Došlić
- Department of Physical Chemistry, Ruđer Bošković Institute, HR-10000 Zagreb, Croatia
| | - Matthieu Sala
- Laboratoire Interdisciplinaire Carnot de Bourgogne UMR 6303 CNRS, Université de Bourgogne, BP 47870, F-21078 Dijon, France and Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstr. 40, D-24098 Kiel, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Lichtenbergstr. 4, 85747 Garching, Germany
| |
Collapse
|
11
|
Mignolet B, Curchod BFE. Excited-State Molecular Dynamics Triggered by Light Pulses—Ab Initio Multiple Spawning vs Trajectory Surface Hopping. J Phys Chem A 2019; 123:3582-3591. [DOI: 10.1021/acs.jpca.9b00940] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Benoit Mignolet
- Theoretical Physical Chemistry, UR MolSYS, B6c, University of Liège, B4000 Liège, Belgium
| | - Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
12
|
Posenitskiy E, Rapacioli M, Lepetit B, Lemoine D, Spiegelman F. Non-adiabatic molecular dynamics investigation of the size dependence of the electronic relaxation in polyacenes. Phys Chem Chem Phys 2019; 21:12139-12149. [DOI: 10.1039/c9cp00603f] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Electronic relaxation from the brightest excited state has been investigated for neutral polyacenes ranging in size from naphthalene to heptacene.
Collapse
Affiliation(s)
- Evgeny Posenitskiy
- Laboratoire Collisions Agrégats et Réactivité (LCAR)
- IRSAMC UMR5589
- Université de Toulouse (UPS) and CNRS
- F-31062 Toulouse
- France
| | - Mathias Rapacioli
- Laboratoire de Chimie et Physique Quantiques (LCPQ)
- IRSAMC UMR5626
- Université de Toulouse (UPS) and CNRS
- F-31062 Toulouse
- France
| | - Bruno Lepetit
- Laboratoire Collisions Agrégats et Réactivité (LCAR)
- IRSAMC UMR5589
- Université de Toulouse (UPS) and CNRS
- F-31062 Toulouse
- France
| | - Didier Lemoine
- Laboratoire Collisions Agrégats et Réactivité (LCAR)
- IRSAMC UMR5589
- Université de Toulouse (UPS) and CNRS
- F-31062 Toulouse
- France
| | - Fernand Spiegelman
- Laboratoire de Chimie et Physique Quantiques (LCPQ)
- IRSAMC UMR5626
- Université de Toulouse (UPS) and CNRS
- F-31062 Toulouse
- France
| |
Collapse
|
13
|
Ma T, Bonfanti M, Eisenbrandt P, Martinazzo R, Burghardt I. Multi-configurational Ehrenfest simulations of ultrafast nonadiabatic dynamics in a charge-transfer complex. J Chem Phys 2018; 149:244107. [DOI: 10.1063/1.5062608] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tianji Ma
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt/Main, Germany
| | - Matteo Bonfanti
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt/Main, Germany
| | - Pierre Eisenbrandt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt/Main, Germany
| | - Rocco Martinazzo
- Department of Chemistry, Università degli Studi di Milano, Via Golgi 19, 20133 Milano, Italy
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry, Goethe University Frankfurt, Max-von-Laue-Str. 7, D-60438 Frankfurt/Main, Germany
| |
Collapse
|
14
|
Richings GW, Robertson C, Habershon S. Improved on-the-Fly MCTDH Simulations with Many-Body-Potential Tensor Decomposition and Projection Diabatization. J Chem Theory Comput 2018; 15:857-870. [DOI: 10.1021/acs.jctc.8b00819] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry, CV4 7AL, U.K
| | - Christopher Robertson
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry, CV4 7AL, U.K
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry, CV4 7AL, U.K
| |
Collapse
|
15
|
Han S, Xie D, Guo H. Modified Gaussian Wave Packet Method for Calculating Initial State Wave Functions in Photodissociation. J Chem Theory Comput 2018; 14:5527-5534. [PMID: 30234984 DOI: 10.1021/acs.jctc.8b00771] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A modified Gaussian wave packet relaxation method is proposed to calculate the ground state wave function using an expansion of frozen Gaussian wave packets. This new procedure consists of two steps. In the first step, a multidimensional Gaussian product placed at the ground state equilibrium geometry is propagated in imaginary time. The relaxation optimizes the widths of the one-dimensional Gaussians. In the second step, additional Gaussian wave packets with the same widths are placed near the equilibrium geometry, and the corresponding expansion coefficients are optimized using the same relaxation method. This new algorithm is tested in photodissociation of NOCl and NH3, and the results show good agreement with the exact results in the energy, wave function, and absorption spectrum. In particular, the highly structured absorption spectrum of NH3 is reproduced, underscoring the accuracy of both the initial wave packet and the excited state propagation.
Collapse
Affiliation(s)
- Shanyu Han
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering , Nanjing University , Nanjing 210093 , China.,Department of Chemistry and Chemical Biology , University of New Mexico , Albuquerque , New Mexico 87131 , United States
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering , Nanjing University , Nanjing 210093 , China
| | - Hua Guo
- Department of Chemistry and Chemical Biology , University of New Mexico , Albuquerque , New Mexico 87131 , United States
| |
Collapse
|
16
|
Bonfanti M, Petersen J, Eisenbrandt P, Burghardt I, Pollak E. Computation of the S 1 ← S 0 Vibronic Absorption Spectrum of Formaldehyde by Variational Gaussian Wavepacket and Semiclassical IVR Methods. J Chem Theory Comput 2018; 14:5310-5323. [PMID: 30141930 DOI: 10.1021/acs.jctc.8b00355] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The vibronic absorption spectrum of the electric dipole forbidden and vibronically allowed S1(1 A2) ← S0(1 A1) transition of formaldehyde is calculated by Gaussian wavepacket and semiclassical methods, along with numerically exact reference calculations, using the potential energy surface of Fu, Shepler, and Bowman ( J. Am. Chem. Soc. 2011, 133, 7957). Specifically, the variational multiconfigurational Gaussian (vMCG) approach and the Herman-Kluk semiclassical initial value representation (HK-SCIVR) are compared to assess the accuracy and convergence of these methods, benchmarked against numerically exact time-dependent wavepacket propagation (TDWP) on the reference potential energy surface. The vMCG calculation is shown to converge quite well with about 100 variationally evolving Gaussian functions and using a local cubic expansion instead of the conventional local harmonic approximation. By contrast, the HK-SCIVR approach with ∼105 trajectories reproduces the vibrationally structured spectral envelope correctly but yields a strongly broadened spectrum. The comparison of the computed absorption spectrum with experiment shows that the relevant vibronic progressions are reasonably reproduced by all computations, but deviations of the order of 10-100 cm-1 occur, underscoring that both electronic structure calculations and dynamical approaches remain challenging in the calculation of typical small-molecule excited-state spectra by trajectory-based methods.
Collapse
Affiliation(s)
- Matteo Bonfanti
- Institute of Physical and Theoretical Chemistry , Goethe University Frankfurt , Max-von-Laue-Str. 7 , D-60438 Frankfurt/Main , Germany
| | - Jakob Petersen
- Chemical and Biological Physics Department , Weizmann Institute of Science , 76000 Rehovot , Israel
| | - Pierre Eisenbrandt
- Institute of Physical and Theoretical Chemistry , Goethe University Frankfurt , Max-von-Laue-Str. 7 , D-60438 Frankfurt/Main , Germany
| | - Irene Burghardt
- Institute of Physical and Theoretical Chemistry , Goethe University Frankfurt , Max-von-Laue-Str. 7 , D-60438 Frankfurt/Main , Germany
| | - Eli Pollak
- Chemical and Biological Physics Department , Weizmann Institute of Science , 76000 Rehovot , Israel
| |
Collapse
|
17
|
Affiliation(s)
- Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
18
|
Chen L, Borrelli R, Zhao Y. Dynamics of Coupled Electron–Boson Systems with the Multiple Davydov D1 Ansatz and the Generalized Coherent State. J Phys Chem A 2017; 121:8757-8770. [DOI: 10.1021/acs.jpca.7b07069] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lipeng Chen
- Division
of Materials Science, Nanyang Technological University, 639798, Singapore
| | - Raffaele Borrelli
- Department
of Agricultural, Forestry and Food Science, Universitá di Torino, I-10095 Grugliasco, Turin, Italy
| | - Yang Zhao
- Division
of Materials Science, Nanyang Technological University, 639798, Singapore
| |
Collapse
|
19
|
Richings GW, Habershon S. Direct grid-based quantum dynamics on propagated diabatic potential energy surfaces. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.01.063] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
20
|
|
21
|
Richings GW, Habershon S. Direct Quantum Dynamics Using Grid-Based Wave Function Propagation and Machine-Learned Potential Energy Surfaces. J Chem Theory Comput 2017; 13:4012-4024. [DOI: 10.1021/acs.jctc.7b00507] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry and
Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry and
Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
22
|
Hirshberg B, Sagiv L, Gerber RB. Approximate Quantum Dynamics using Ab Initio Classical Separable Potentials: Spectroscopic Applications. J Chem Theory Comput 2017; 13:982-991. [DOI: 10.1021/acs.jctc.6b01129] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Barak Hirshberg
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
| | - Lior Sagiv
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
| | - R. Benny Gerber
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
23
|
Borrelli R, Gelin MF. The Generalized Coherent State ansatz: Application to quantum electron-vibrational dynamics. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.05.013] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
24
|
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
|
25
|
Vacher M, Bearpark MJ, Robb MA. Direct methods for non-adiabatic dynamics: connecting the single-set variational multi-configuration Gaussian (vMCG) and Ehrenfest perspectives. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1937-2] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
26
|
Borrelli R, Peluso A. Quantum dynamics of electronic transitions with Gauss-Hermite wave packets. J Chem Phys 2016; 144:114102. [DOI: 10.1063/1.4943538] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Raffaele Borrelli
- Department of Agricultural, Forestry and Food Science, Università di Torino, I-10095 Grugliasco, TO, Italy
| | - Andrea Peluso
- Department of Chemistry and Biology, Università di Salerno, I-84081 Fisciano, SA, Italy
| |
Collapse
|
27
|
Kong X, Markmann A, Batista VS. Time-Sliced Thawed Gaussian Propagation Method for Simulations of Quantum Dynamics. J Phys Chem A 2016; 120:3260-9. [DOI: 10.1021/acs.jpca.5b12192] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Xiangmeng Kong
- Department of Chemistry, Yale University, 225 Prospect
Street, New Haven, Connecticut 06520-8107, United States
| | - Andreas Markmann
- Department of Chemistry, Yale University, 225 Prospect
Street, New Haven, Connecticut 06520-8107, United States
| | - Victor S. Batista
- Department of Chemistry, Yale University, 225 Prospect
Street, New Haven, Connecticut 06520-8107, United States
| |
Collapse
|
28
|
Mendive-Tapia D, Kortekaas L, Steen JD, Perrier A, Lasorne B, Browne WR, Jacquemin D. Accidental degeneracy in the spiropyran radical cation: charge transfer between two orthogonal rings inducing ultra-efficient reactivity. Phys Chem Chem Phys 2016; 18:31244-31253. [DOI: 10.1039/c6cp06907j] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We unravel an original photoswitching mechanism in spiropyran radical cation.
Collapse
Affiliation(s)
- David Mendive-Tapia
- Institut Charles Gerhardt Montpellier
- UMR 5253
- CNRS-UM-ENSCM
- CTMM
- Université Montpellier
| | - Luuk Kortekaas
- Molecular Inorganic Chemistry
- Stratingh Institute for Chemistry
- Faculty of Mathematics and Natural Sciences
- University of Groningen
- 9747AG Groningen
| | - Jorn D. Steen
- Molecular Inorganic Chemistry
- Stratingh Institute for Chemistry
- Faculty of Mathematics and Natural Sciences
- University of Groningen
- 9747AG Groningen
| | - Aurélie Perrier
- Université Paris Diderot Sorbonne Paris Cité
- F-75205 Paris Cedex 13
- France
- Institut de Recherche de Chimie Paris
- PSL Research University
| | - Benjamin Lasorne
- Institut Charles Gerhardt Montpellier
- UMR 5253
- CNRS-UM-ENSCM
- CTMM
- Université Montpellier
| | - Wesley R. Browne
- Molecular Inorganic Chemistry
- Stratingh Institute for Chemistry
- Faculty of Mathematics and Natural Sciences
- University of Groningen
- 9747AG Groningen
| | - Denis Jacquemin
- Chimie Et Interdisciplinarité
- Synthèse, Analyse, Modélisation (CEISAM)
- UMR CNRS no. 6230
- 44322 Nantes Cedex 3
- France
| |
Collapse
|
29
|
Richings GW, Worth GA. A Practical Diabatisation Scheme for Use with the Direct-Dynamics Variational Multi-Configuration Gaussian Method. J Phys Chem A 2015; 119:12457-70. [DOI: 10.1021/acs.jpca.5b07921] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Gareth W. Richings
- School
of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, U.K
| | - Graham A. Worth
- School
of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, U.K
| |
Collapse
|
30
|
Skouteris D, Barone V. Nonadiabatic photodynamics of phenol on a realistic potential energy surface by a novel multilayer Gaussian MCTDH program. Chem Phys Lett 2015; 636:15-21. [PMID: 29551833 DOI: 10.1016/j.cplett.2015.06.065] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We report the main features of a new implementation of the Gaussian Multi-Configuration Time-Dependent Hartree (G-MCTDH) model. The code allows effective computations of time-dependent phenomena, including calculation of vibronic spectra (in one or more electronic states), relative state populations etc., with the possibility of a multilayer formulation. We have validated the code on the diabatic surfaces recently published by Truhlar and coworkers to study the nonadiabatic photodynamics of phenol. Using an Ehrenfest-like, single-nuclear-configuration (but in a fully quantum formalism) model we calculate the optical spectrum and relative state populations of the system as a function of time.
Collapse
Affiliation(s)
- D Skouteris
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - V Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
31
|
Polyak I, Allan CSM, Worth GA. A complete description of tunnelling using direct quantum dynamics simulation: Salicylaldimine proton transfer. J Chem Phys 2015; 143:084121. [DOI: 10.1063/1.4929478] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Iakov Polyak
- School of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
| | - Charlotte S. M. Allan
- School of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
| | - Graham A. Worth
- School of Chemistry, University of Birmingham, Edgbaston, Birmingham B15 2TT, United Kingdom
| |
Collapse
|
32
|
Richings G, Polyak I, Spinlove K, Worth G, Burghardt I, Lasorne B. Quantum dynamics simulations using Gaussian wavepackets: the vMCG method. INT REV PHYS CHEM 2015. [DOI: 10.1080/0144235x.2015.1051354] [Citation(s) in RCA: 174] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
33
|
Tavernelli I. Nonadiabatic molecular dynamics simulations: synergies between theory and experiments. Acc Chem Res 2015; 48:792-800. [PMID: 25647401 DOI: 10.1021/ar500357y] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Recent developments in nonadiabatic dynamics enabled ab inito simulations of complex ultrafast processes in the condensed phase. These advances have opened new avenues in the study of many photophysical and photochemical reactions triggered by the absorption of electromagnetic radiation. In particular, theoretical investigations can be combined with the most sophisticated femtosecond experimental techniques to guide the interpretation of measured time-resolved observables. At the same time, the availability of experimental data at high (spatial and time) resolution offers a unique opportunity for the benchmarking and the improvement of those theoretical models used to describe complex molecular systems in their natural environment. The established synergy between theory and experiments can produce a better understanding of new ultrafast physical and chemical processes at atomistic scale resolution. Furthermore, reliable ab inito molecular dynamics simulations can already be successfully employed as predictive tools to guide new experiments as well as the design of novel and better performing materials. In this paper, I will give a concise account on the state of the art of molecular dynamics simulations of complex molecular systems in their excited states. The principal aim of this approach is the description of a given system of interest under the most realistic ambient conditions including all environmental effects that influence experiments, for instance, the interaction with the solvent and with external time-dependent electric fields, temperature, and pressure. To this end, time-dependent density functional theory (TDDFT) is among the most efficient and accurate methods for the representation of the electronic dynamics, while trajectory surface hopping gives a valuable representation of the nuclear quantum dynamics in the excited states (including nonadiabatic effects). Concerning the environment and its effects on the dynamics, the quantum mechanics/molecular mechanics (QM/MM) approach has the advantage of providing an atomistic (even though approximated) description of the solvent molecules, which is crucial for the characterization of all ultrafast relaxation phenomena that depend on the geometrical arrangement at the interface between a molecule and the solvent, for example, the hydrogen bond network. After a short description of the TDDFT-based implementation of Ehrenfest and trajectory surface hopping dynamics, I will present applications in different domains of molecular chemistry and physics: the analysis and the understanding of (time-resolved) X-ray absorption spectra, the interpretation of the ultrafast relaxation dynamics of photoexcited dyes in solution, and the design of specific laser pulses (capable of inducing desired chemical reactions) using local control theory.
Collapse
Affiliation(s)
- Ivano Tavernelli
- Laboratory
of Computational
Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, Laussane 1015, Switzerland
| |
Collapse
|
34
|
Saller MAC, Habershon S. Basis Set Generation for Quantum Dynamics Simulations Using Simple Trajectory-Based Methods. J Chem Theory Comput 2014; 11:8-16. [DOI: 10.1021/ct500657f] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Maximilian A. C. Saller
- Department
of Chemistry and
Centre for Scientific Computing, University of Warwick, Coventry, CV4 7AL, United Kingdom
| | - Scott Habershon
- Department
of Chemistry and
Centre for Scientific Computing, University of Warwick, Coventry, CV4 7AL, United Kingdom
| |
Collapse
|
35
|
Malhado JP, Bearpark MJ, Hynes JT. Non-adiabatic dynamics close to conical intersections and the surface hopping perspective. Front Chem 2014; 2:97. [PMID: 25485263 PMCID: PMC4240175 DOI: 10.3389/fchem.2014.00097] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2014] [Accepted: 10/16/2014] [Indexed: 11/13/2022] Open
Abstract
Conical intersections play a major role in the current understanding of electronic de-excitation in polyatomic molecules, and thus in the description of photochemistry and photophysics of molecular systems. This article reviews aspects of the basic theory underlying the description of non-adiabatic transitions at conical intersections, with particular emphasis on the important case when the dynamics of the nuclei are treated classically. Within this classical nuclear motion framework, the main aspects of the surface hopping methodology in the conical intersection context are presented. The emerging picture from this treatment is that of electronic transitions around conical intersections dominated by the interplay of the nuclear velocity and the derivative non-adiabatic coupling vector field.
Collapse
Affiliation(s)
| | | | - James T Hynes
- Department of Chemistry and Biochemistry, University of Colorado, Boulder Boulder, CO, USA ; Département de Chimie, École Normale Supérieur, UMR ENS-CNRS-UPMC 8640 Paris, France
| |
Collapse
|
36
|
|
37
|
Blancafort L. Photochemistry and photophysics at extended seams of conical intersection. Chemphyschem 2014; 15:3166-81. [PMID: 25157686 DOI: 10.1002/cphc.201402359] [Citation(s) in RCA: 65] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2014] [Indexed: 11/07/2022]
Abstract
The role of extended seams of conical intersection in excited-state mechanisms is reviewed. Seams are crossings of the potential energy surface in many dimensions where the decay from the excited to the ground state can occur, and the extended seam is composed of different segments lying along a reaction coordinate. Every segment is associated with a different primary photoproduct, which gives rise to competing pathways. This idea is first illustrated for fulvene and ethylene, and then it is used to explain more complex cases such as the dependence of the isomerisation of retinal chromophore isomers on the protein environment, the dependence of the efficiency of the azobenzene photochemical switch on the wavelength of irradiation and the direction of the isomerisation, and the coexistence of different mechanisms in the photo-induced Wolff rearrangement of diazonaphthoquinone. The role of extended seams in the photophysics of the DNA nucleobases and the relationship between two-state seams and three-state crossings is also discussed. As an outlook, the design of optical control strategies based on the passage of the excited molecule through the seam is considered, and it is shown how the excited-state lifetime of fulvene can be modulated by shaping the energy of the seam.
Collapse
Affiliation(s)
- Lluís Blancafort
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Campus de Montilivi, 17071 Girona (Spain).
| |
Collapse
|
38
|
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]
|
39
|
Skouteris D, Barone V. A new Gaussian MCTDH program: implementation and validation on the levels of the water and glycine molecules. J Chem Phys 2014; 140:244104. [PMID: 24985615 DOI: 10.1063/1.4883677] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
We report the main features of a new general implementation of the Gaussian Multi-Configuration Time-Dependent Hartree model. The code allows effective computations of time-dependent phenomena, including calculation of vibronic spectra (in one or more electronic states), relative state populations, etc. Moreover, by expressing the Dirac-Frenkel variational principle in terms of an effective Hamiltonian, we are able to provide a new reliable estimate of the representation error. After validating the code on simple one-dimensional systems, we analyze the harmonic and anharmonic vibrational spectra of water and glycine showing that reliable and converged energy levels can be obtained with reasonable computing resources. The data obtained on water and glycine are compared with results of previous calculations using the vibrational second-order perturbation theory method. Additional features and perspectives are also shortly discussed.
Collapse
Affiliation(s)
- D Skouteris
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - V Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
40
|
Jenkins S, Blancafort L, Kirk SR, Bearpark MJ. The response of the electronic structure to electronic excitation and double bond torsion in fulvene: a combined QTAIM, stress tensor and MO perspective. Phys Chem Chem Phys 2014; 16:7115-26. [PMID: 24618666 DOI: 10.1039/c4cp00003j] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
New insights into the double bond isomerization of fulvene in the ground and excited electronic states are provided by newly developed QTAIM and stress tensor tools. The S0 and S1 states follow the 'biradical' torsion model, but the double bond is stiffer in the S0 state; by contrast, the S2 state follows the 'zwitterionic' torsion. Differences are explained in terms of the ellipticity and bond critical point (BCP) stiffness for both QTAIM and the stress tensor. Overall, the wave-function based analysis is found to be in agreement with the work of Bonačić-Koutecký and Michl that the bond-twisted species can have biradical or zwitterionic character, depending on the state. Using QTAIM and the stress tensor a new understanding of bond torsion is revealed; the electronic charge density around the twisted bond is found not to rotate in concert with the nuclei of the rotated -CH2 methylene group. The ability to visualize how the bond stiffness varies between individual electronic states and how this correlates with the QTAIM and stress tensor bond stiffness is highlighted. In addition, the most and least preferred morphologies of bond-path torsion are visualized. Briefly we discuss the prospects for using this new QTAIM and stress tensor analysis for excited state chemistry.
Collapse
Affiliation(s)
- Samantha Jenkins
- Key Laboratory of Chemical Biology and Traditional Chinese Medicine Research (Ministry of Education of China) and Key Laboratory of Resource Fine-Processing and Advanced Materials of Hunan Province, College of Chemistry and Chemical Engineering, Hunan Normal University, Changsha Hunan 410081, China.
| | | | | | | |
Collapse
|
41
|
Mendive-Tapia D, Vacher M, Bearpark MJ, Robb MA. Coupled electron-nuclear dynamics: Charge migration and charge transfer initiated near a conical intersection. J Chem Phys 2013; 139:044110. [DOI: 10.1063/1.4815914] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|