1
|
Giarrusso S, Gori-Giorgi P, Agostini F. Electronic Vector Potential from the Exact Factorization of a Complex Wavefunction. Chemphyschem 2024:e202400127. [PMID: 38837609 DOI: 10.1002/cphc.202400127] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2024] [Revised: 06/01/2024] [Accepted: 06/03/2024] [Indexed: 06/07/2024]
Abstract
We generalize the definitions of local scalar potentials namedυ kin ${\upsilon _{{\rm{kin}}} }$ andυ N - 1 ${\upsilon _{N - 1} }$ , which are relevant to properly describe phenomena such as molecular dissociation with density-functional theory, to the case in which the electronic wavefunction corresponds to a complex current-carrying state. In such a case, an extra term in the form of a vector potential appears which cannot be gauged away. Both scalar and vector potentials are introduced via the exact factorization formalism which allows us to express the given Schrödinger equation as two coupled equations, one for the marginal and one for the conditional amplitude. The electronic vector potential is directly related to the paramagnetic current density carried by the total wavefunction and to the diamagnetic current density in the equation for the marginal amplitude. An explicit example of this vector potential in a triplet state of two non-interacting electrons is showcased together with its associated circulation, giving rise to a non-vanishing geometric phase. Some connections with the exact factorization for the full molecular wavefunction beyond the Born-Oppenheimer approximation are also discussed.
Collapse
Affiliation(s)
- Sara Giarrusso
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| | - Paola Gori-Giorgi
- Department of Chemistry & Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV, Amsterdam, The Netherlands
- Microsoft Research AI4Science, Evert van de Beekstraat 354, 1118CZ, Schiphol, The Netherlands
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| |
Collapse
|
2
|
Villaseco Arribas E, Maitra NT, Agostini F. Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states. J Chem Phys 2024; 160:054102. [PMID: 38310471 DOI: 10.1063/5.0186984] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Accepted: 01/11/2024] [Indexed: 02/05/2024] Open
Abstract
Advances in coherent light sources and development of pump-probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron-nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. This work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.
Collapse
Affiliation(s)
- Evaristo Villaseco Arribas
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
3
|
Pieroni C, Sangiogo Gil E, Ibele LM, Persico M, Granucci G, Agostini F. Investigating the Photodynamics of trans-Azobenzene with Coupled Trajectories. J Chem Theory Comput 2024; 20:580-596. [PMID: 38177105 DOI: 10.1021/acs.jctc.3c00978] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2024]
Abstract
In this work, we present the first implementation of coupled-trajectory Tully surface hopping (CT-TSH) suitable for applications to molecular systems. We combine CT-TSH with the semiempirical floating occupation molecular orbital-configuration interaction electronic structure method to investigate the photoisomerization dynamics of trans-azobenzene. Our study shows that CT-TSH can capture correctly decoherence effects in this system, yielding consistent electronic and nuclear dynamics in agreement with (standard) decoherence-corrected TSH. Specifically, CT-TSH is derived from the exact factorization and the electronic coefficients' evolution is directly influenced by the coupling of trajectories, resulting in the improvement of internal consistency if compared to standard TSH.
Collapse
Affiliation(s)
- Carlotta Pieroni
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Eduarda Sangiogo Gil
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| | - Lea M Ibele
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| | - Maurizio Persico
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Federica Agostini
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
4
|
Dupuy L, Talotta F, Agostini F, Lauvergnat D, Poirier B, Scribano Y. Adiabatic and Nonadiabatic Dynamics with Interacting Quantum Trajectories. J Chem Theory Comput 2022; 18:6447-6462. [DOI: 10.1021/acs.jctc.2c00744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Lucien Dupuy
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095Montpellier, France
| | - Francesco Talotta
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR-CNRS 8000, 91405Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR-CNRS 8000, 91405Orsay, France
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR-CNRS 8000, 91405Orsay, France
| | - Bill Poirier
- Department of Chemistry and Biochemistry, and Department of Physics, Texas Tech University, Box 41061, 79409-1061Lubbock, Texas, United States
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095Montpellier, France
| |
Collapse
|
5
|
Villaseco Arribas E, Agostini F, Maitra NT. Exact Factorization Adventures: A Promising Approach for Non-Bound States. Molecules 2022; 27:molecules27134002. [PMID: 35807246 PMCID: PMC9267945 DOI: 10.3390/molecules27134002] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Revised: 06/07/2022] [Accepted: 06/09/2022] [Indexed: 11/29/2022] Open
Abstract
Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.
Collapse
Affiliation(s)
| | - Federica Agostini
- Institut de Chimie Physique UMR8000, Université Paris-Saclay, CNRS, 91405 Orsay, France;
| | - Neepa T. Maitra
- Department of Physics, Rutgers University, Newark, NJ 07102, USA;
- Correspondence:
| |
Collapse
|
6
|
Talotta F, Lauvergnat D, Agostini F. Describing the photo-isomerization of a retinal chromophore model with coupled and quantum trajectories. J Chem Phys 2022; 156:184104. [DOI: 10.1063/5.0089415] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The exact factorization of the electron-nuclear wavefunction is applied to the study of the photo- isomerization of a retinal chromophore model. We describe such an ultrafast nonadiabatic process by analyzing the time-dependent potentials of the theory and by mimicking nuclear dynamics with quantum and coupled trajectories. The time-dependent vector and scalar potentials are the signature of the exact factorization, as they guide nuclear dynamics by encoding the complete electronic dynamics and including excited-state effects. Analysis of the potentials is, thus, essential - when possible - to predict the time-dependent behavior of the system of interest. In this work, we employ the exact time-dependent potentials, available for the numerically-exactly solvable model used here, to propagate quantum nuclear trajectories representing the isomerization reaction of the retinal chromophore. The quantum trajectories are the best possible trajectory-based description of the reaction when using the exact-factorization formalism, and thus allow us to assess the performance of the coupled-trajectory, fully approximate, schemes derived from the exact-factorization equations.
Collapse
Affiliation(s)
| | - David Lauvergnat
- Institut de Chimie Physique, UMR 8000, CNRS Délégation Ile-de-France Sud, France
| | | |
Collapse
|
7
|
Ibele LM, Curchod BFE, Agostini F. A Photochemical Reaction in Different Theoretical Representations. J Phys Chem A 2022; 126:1263-1281. [PMID: 35157450 PMCID: PMC8883471 DOI: 10.1021/acs.jpca.1c09604] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
![]()
The Born–Oppenheimer
picture has forged our representation
and interpretation of photochemical processes, from photoexcitation
down to the passage through a conical intersection, a funnel connecting
different electronic states. In this work, we analyze a full in silico
photochemical experiment, from the explicit electronic excitation
by a laser pulse to the formation of photoproducts following a nonradiative
decay through a conical intersection, by contrasting the picture offered
by Born–Oppenheimer and that proposed by the exact factorization.
The exact factorization offers an alternative understanding of photochemistry
that does not rely on concepts such as electronic states, nonadiabatic
couplings, and conical intersections. On the basis of nonadiabatic
quantum dynamics performed for a two-state 2D model system, this work
allows us to compare Born–Oppenheimer and exact factorization
for (i) an explicit photoexcitation with and without the Condon approximation,
(ii) the passage of a nuclear wavepacket through a conical intersection,
(iii) the formation of excited stationary states in the Franck–Condon
region, and (iv) the use of classical and quantum trajectories in
the exact factorization picture to capture nonadiabatic processes
triggered by a laser pulse.
Collapse
Affiliation(s)
- Lea M Ibele
- Department of Chemistry, Durham University, Durham DH1 3LE, United Kingdom
| | - Basile F E Curchod
- Department of Chemistry, Durham University, Durham DH1 3LE, United Kingdom
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
8
|
Vindel-Zandbergen P, Matsika S, Maitra NT. Exact-Factorization-Based Surface Hopping for Multistate Dynamics. J Phys Chem Lett 2022; 13:1785-1790. [PMID: 35170972 DOI: 10.1021/acs.jpclett.1c04132] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
A surface-hopping algorithm recently derived from the exact factorization approach, SHXF [Ha et al. J. Phys. Chem. Lett. 2018, 9, 1097], introduces an additional term in the electronic equation of surface hopping that couples electronic states through the quantum momentum. This term not only provides a first-principles description of decoherence, but here we show it is crucial to accurately capture nonadiabatic dynamics when more than two states are occupied at any given time. Using a vibronic coupling model of the uracil cation, we show that the lack of this term in traditional surface-hopping methods, including those with decoherence corrections, leads to failure to predict the dynamics through a three-state intersection, while SHXF performs similarly to the multiconfiguration time-dependent Hartree quantum dynamics benchmark.
Collapse
Affiliation(s)
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
9
|
Zhang J, Peng J, Hu D, Lan Z. Investigation of nonadiabatic dynamics in the photolysis of methyl nitrate (CH 3ONO 2) by on-the-fly surface hopping simulation. Phys Chem Chem Phys 2021; 23:25597-25611. [PMID: 34546246 DOI: 10.1039/d1cp03226g] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The photolysis mechanism of methyl nitrate (CH3ONO2) was studied using the on-the-fly surface hopping dynamics at the XMS-CASPT2 level. Several critical geometries, including electronic state minima and conical intersections, were obtained, which play essential roles in the nonadiabatic dynamics of CH3ONO2. The ultrafast nonadiabatic decay dynamics to the ground state were simulated, which gives a proper explanation on the broad and structureless absorption spectra of CH3ONO2. The photodissociation channels, including CH3O + NO2, CH3O + NO + O, and others, as well as their branching ratios, were identified. When the dynamics starts from the lowest two electronic states (S1 and S2), the CH3O + NO2 channel is the dominant photolysis pathway, although we observed the minor contributions of other channels. In contrast, when the trajectories start from the third excited state S3, both CH3O + NO2 and CH3O + NO + O channels become important. Here the CH3O-NO2 bond dissociation takes place first, and then for some trajectories, the N-O bond of the NO2 part breaks successively. The quasi-degeneracy of electronic states may exist in the dissociation limits of both CH3O + NO2 and CH3O + NO + O channels. The current work provides valuable information in the understanding of experimental findings of the wavelength-dependent photolysis mechanism of CH3ONO2.
Collapse
Affiliation(s)
- Juanjuan Zhang
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China. .,School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
10
|
Abstract
In this paper, we discuss coupled-trajectory schemes for molecular-dynamics simulations of excited-state processes. New coupled-trajectory strategies to capture decoherence effects, revival of coherence and nonadiabatic interferences in long-time dynamics are proposed, and compared to independent-trajectory schemes. The working framework is provided by the exact factorization of the electron-nuclear wave function, and it exploits ideas emanating from various surface-hopping schemes. The new coupled-trajectory algorithms are tested on a one-dimensional two-state system using different model parameters which allow one to induce different dynamics. The benchmark is provided by the numerically exact solution of the time-dependent Schrödinger equation.
Collapse
Affiliation(s)
- Carlotta Pieroni
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France.,Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Federica Agostini
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
11
|
Schirò M, Eich FG, Agostini F. Quantum-classical nonadiabatic dynamics of Floquet driven systems. J Chem Phys 2021; 154:114101. [PMID: 33752379 DOI: 10.1063/5.0043790] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop a trajectory-based approach for excited-state molecular dynamics simulations of systems subject to an external periodic drive. We combine the exact-factorization formalism, allowing us to treat electron-nuclear systems in nonadiabatic regimes, with the Floquet formalism for time-periodic processes. The theory is developed starting with the molecular time-dependent Schrödinger equation with the inclusion of an external periodic drive that couples to the system dipole moment. With the support of the Floquet formalism, quantum dynamics is approximated by combining classical-like, trajectory-based, nuclear evolution with electronic dynamics represented in the Floquet basis. The resulting algorithm, which is an extension of the coupled-trajectory mixed quantum-classical scheme for periodically driven systems, is applied to a model study, exactly solvable, with different field intensities.
Collapse
Affiliation(s)
- Marco Schirò
- JEIP, USR 3573 CNRS, Collège de France, PSL Research University, 11 Place Marcelin Berthelot, 75321 Paris Cedex 05, France
| | - Florian G Eich
- HQS Quantum Simulations GmbH, Haid-und-Neu-Straße 7, D-76131 Karlsruhe, Germany
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
12
|
Pieroni C, Marsili E, Lauvergnat D, Agostini F. Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies. J Chem Phys 2021; 154:034104. [PMID: 33499611 DOI: 10.1063/5.0036726] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.
Collapse
Affiliation(s)
- Carlotta Pieroni
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Emanuele Marsili
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
13
|
Martinez P, Rosenzweig B, Hoffmann NM, Lacombe L, Maitra NT. Case studies of the time-dependent potential energy surface for dynamics in cavities. J Chem Phys 2021; 154:014102. [PMID: 33412864 PMCID: PMC7968936 DOI: 10.1063/5.0033386] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2020] [Accepted: 12/10/2020] [Indexed: 11/14/2022] Open
Abstract
The exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed.
Collapse
Affiliation(s)
- Phillip Martinez
- Department of Physics and Astronomy, Hunter College of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| | | | - Norah M. Hoffmann
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA
| | - Lionel Lacombe
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA
| | - Neepa T. Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA
| |
Collapse
|
14
|
Talotta F, Agostini F, Ciccotti G. Quantum Trajectories for the Dynamics in the Exact Factorization Framework: A Proof-of-Principle Test. J Phys Chem A 2020; 124:6764-6777. [PMID: 32786992 DOI: 10.1021/acs.jpca.0c03969] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation. Those characteristics are referred to as quantum trajectories, since they are generated via ordinary differential equations similar to Hamilton's equations, but including the so-called quantum potential, and they can be used to reconstruct exactly the quantum-mechanical nuclear wave function, provided infinite initial conditions are propagated in time.
Collapse
Affiliation(s)
- Francesco Talotta
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France.,Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405, Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France.,Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405, Orsay, France
| | - Giovanni Ciccotti
- CNR, Institute for Applied Computing "Mauro Picone" (IAC), Via dei Taurini 19, 00185 Rome, Italy.,School of Physics, University College of Dublin UCD - Belfield, Dublin 4, Ireland.,Dipartimento di Fisica, Università di Roma La Sapienza, P. le A. Moro 5, 00185 Roma, Italy
| |
Collapse
|
15
|
Talotta F, Morisset S, Rougeau N, Lauvergnat D, Agostini F. Internal Conversion and Intersystem Crossing with the Exact Factorization. J Chem Theory Comput 2020; 16:4833-4848. [PMID: 32633509 DOI: 10.1021/acs.jctc.0c00493] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present a detailed derivation of the generalized coupled-trajectory mixed quantum-classical (G-CT-MQC) algorithm based on the exact-factorization equations. The ultimate goal is to propose an algorithm that can be employed for molecular dynamics simulations of nonradiative phenomena, as the spin-allowed internal conversions and the spin-forbidden intersystem crossings. Internal conversions are nonadiabatic processes driven by the kinetic coupling between electronic states, whereas intersystem crossings are mediated by the spin-orbit coupling. In this paper, we discuss computational issues related to the suitable representation for electronic dynamics and the different natures of kinetic and spin-orbit coupling. Numerical studies on model systems allow us to test the performance of the G-CT-MQC algorithm in different situations.
Collapse
Affiliation(s)
- Francesco Talotta
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France.,Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Sabine Morisset
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Nathalie Rougeau
- Institut des Sciences Moléculaires d'Orsay, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - David Lauvergnat
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Federica Agostini
- Institut de Chimie Physique, UMR8000, CNRS, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
16
|
Talotta F, Morisset S, Rougeau N, Lauvergnat D, Agostini F. Spin-Orbit Interactions in Ultrafast Molecular Processes. PHYSICAL REVIEW LETTERS 2020; 124:033001. [PMID: 32031839 DOI: 10.1103/physrevlett.124.033001] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Indexed: 06/10/2023]
Abstract
We investigate spin-orbit interactions in ultrafast molecular processes employing the exact factorization of the electron-nuclear wave function. We revisit the original derivation by including spin-orbit coupling, and show how the dynamics driven by the time-dependent potential energy surface alleviates inconsistencies arising from different electronic representations. We propose a novel trajectory-based scheme to simulate spin-forbidden non-radiative processes, and we show its performance in the treatment of excited-state dynamics where spin-orbit effects couple different spin multiplets.
Collapse
Affiliation(s)
- Francesco Talotta
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Sabine Morisset
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Nathalie Rougeau
- Institut de Sciences Moléculaires d'Orsay, UMR 8214 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - David Lauvergnat
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| | - Federica Agostini
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
17
|
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
|
18
|
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
|
19
|
Agostini F, Curchod BFE. Different flavors of nonadiabatic molecular dynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1417] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Federica Agostini
- Laboratoire de Chimie Physique UMR 8000 CNRS/University Paris‐Sud Orsay France
| | | |
Collapse
|
20
|
Gossel GH, Lacombe L, Maitra NT. On the numerical solution of the exact factorization equations. J Chem Phys 2019; 150:154112. [PMID: 31005081 DOI: 10.1063/1.5090802] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The exact factorization (EF) approach to coupled electron-ion dynamics recasts the time-dependent molecular Schrödinger equation as two coupled equations, one for the nuclear wavefunction and one for the conditional electronic wavefunction. The potentials appearing in these equations have provided insight into non-adiabatic processes, and new practical non-adiabatic dynamics methods have been formulated starting from these equations. Here, we provide a first demonstration of a self-consistent solution of the exact equations, with a preliminary analysis of their stability and convergence properties. The equations have an unprecedented mathematical form, involving a Hamiltonian outside the class of Hermitian Hamiltonians usually encountered in time-propagation, and so the usual numerical methods for time-dependent Schrödinger fail when applied in a straightforward way to the EF equations. We find an approach that enables stable propagation long enough to witness non-adiabatic behavior in a model system before non-trivial instabilities take over. Implications for the development and analysis of EF-based methods are discussed.
Collapse
Affiliation(s)
- Graeme H Gossel
- Department of Physics and Astronomy, Hunter College and the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| | - Lionel Lacombe
- Department of Physics and Astronomy, Hunter College and the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| | - Neepa T Maitra
- Department of Physics and Astronomy, Hunter College and the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| |
Collapse
|
21
|
Agostini F, Gross E, Curchod BF. Electron-nuclear entanglement in the time-dependent molecular wavefunction. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.01.021] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
22
|
Gossel GH, Agostini F, Maitra NT. Coupled-Trajectory Mixed Quantum-Classical Algorithm: A Deconstruction. J Chem Theory Comput 2018; 14:4513-4529. [DOI: 10.1021/acs.jctc.8b00449] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Graeme H. Gossel
- Department of Physics and Astronomy, Hunter College and the City University of New York, 695 Park Avenue, New York, New York 10065, United States
| | - Federica Agostini
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, 91405 Orsay, France
| | - Neepa T. Maitra
- Department of Physics and Astronomy, Hunter College and the City University of New York, 695 Park Avenue, New York, New York 10065, United States
- The Physics Program and the Chemistry Program of the Graduate Center, City University of New York, 365 Fifth Avenue, New York, United States
| |
Collapse
|
23
|
Crespo-Otero R, Barbatti M. Recent Advances and Perspectives on Nonadiabatic Mixed Quantum–Classical Dynamics. Chem Rev 2018; 118:7026-7068. [DOI: 10.1021/acs.chemrev.7b00577] [Citation(s) in RCA: 301] [Impact Index Per Article: 50.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | | |
Collapse
|
24
|
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
|
25
|
Min SK, Agostini F, Tavernelli I, Gross EKU. Ab Initio Nonadiabatic Dynamics with Coupled Trajectories: A Rigorous Approach to Quantum (De)Coherence. J Phys Chem Lett 2017; 8:3048-3055. [PMID: 28618782 DOI: 10.1021/acs.jpclett.7b01249] [Citation(s) in RCA: 93] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We report the first nonadiabatic molecular dynamics study based on the exact factorization of the electron-nuclear wave function. Our approach (a coupled-trajectory mixed quantum-classical, CT-MQC, scheme) is based on the quantum-classical limit derived from systematic and controlled approximations to the full quantum-mechanical problem formulated in the exact-factorization framework. Its strength is the ability to correctly capture quantum (de)coherence effects in a trajectory-based approach to excited-state dynamics. We show this by benchmarking CT-MQC dynamics against a revised version of the popular fewest-switches surface-hopping scheme that is able to fix its well-documented overcoherence issue. The CT-MQC approach is successfully applied to investigation of the photochemistry (ring-opening) of oxirane in the gas phase, analyzing in detail the role of decoherence. This work represents a significant step forward in the establishment of the exact factorization as a powerful tool to study excited-state dynamics, not only for interpretation purposes but mainly for nonadiabatic ab initio molecular dynamics simulations.
Collapse
Affiliation(s)
- Seung Kyu Min
- Department of Chemistry, School of Natural Science, Ulsan National Institute of Science and Technology (UNIST) , Ulsan 44919, Korea
| | - Federica Agostini
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay , 91405 Orsay, France
| | - Ivano Tavernelli
- IBM Research GmbH, Zürich Research Laboratory , 8803 Rüschlikon, Switzerland
| | - E K U Gross
- Max-Planck Institut für Mikrostrukturphysik , Weinberg 2, D-06120 Halle, Germany
| |
Collapse
|
26
|
Abstract
Conical intersections represent critical topological features of potential energy surfaces and open ultrafast nonradiative deactivation channels for photoexcited molecules. In the following, we investigate how this funneling picture is transposed in the eyes of the exact factorization formalism for a 2D model system. The exact factorization of the total molecular wave function leads to the fundamental concept of time-dependent potential energy surface and time-dependent vector potential, whose behavior during a dynamics through a conical intersection has up to now remained unexplored. Despite the fact that these quantities might be viewed as time-dependent generalizations of the adiabatic potential energy surfaces and the nonadiabatic coupling vectors, characteristic quantities appearing in the Born-Oppenheimer framework, we observe that they do not exhibit particular topological features in the region of conical intersection but still reflect the complex dynamics of the nuclear wavepacket.
Collapse
Affiliation(s)
- Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, United Kingdom
| | - Federica Agostini
- Laboratoire de Chimie Physique, UMR 8000 CNRS/University Paris-Sud, University Paris-Saclay , 91405 Orsay, France
| |
Collapse
|