1
|
Ibele LM, Agostini F. Exploring Exact-Factorization-Based Trajectories for Low-Energy Dynamics near a Conical Intersection. J Phys Chem A 2024. [PMID: 38660710 DOI: 10.1021/acs.jpca.4c00555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
We study low-energy dynamics generated by a two-dimensional two-state Jahn-Teller Hamiltonian in the vicinity of a conical intersection using quantum wave packet and trajectory dynamics. Recently, these dynamics were studied by comparing the adiabatic representation and the exact factorization, with the purpose to highlight the different nature of topological-phase and geometric-phase effects arising in the two theoretical representations of the same problem. Here, we employ the exact factorization to understand how to accurately model low-energy dynamics in the vicinity of a conical intersection using an approximate description of the nuclear motion that uses trajectories. We find that since nonadiabatic effects are weak but non-negligible, the trajectory-based description that invokes the classical approximation struggles to capture the correct behavior.
Collapse
Affiliation(s)
- Lea M Ibele
- 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
|
2
|
Suchan J, Liang F, Durden AS, Levine BG. Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone. J Chem Phys 2024; 160:134310. [PMID: 38573851 DOI: 10.1063/5.0198333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 03/18/2024] [Indexed: 04/06/2024] Open
Abstract
Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum-classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?
Collapse
Affiliation(s)
- Jiří Suchan
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
Arribas EV, Ibele LM, Lauvergnat D, Maitra NT, Agostini F. Significance of Energy Conservation in Coupled-Trajectory Approaches to Nonadiabatic Dynamics. J Chem Theory Comput 2023; 19:7787-7800. [PMID: 37853509 DOI: 10.1021/acs.jctc.3c00845] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
Through approximating electron-nuclear correlation terms in the exact factorization approach, trajectory-based methods have been derived and successfully applied to the dynamics of a variety of light-induced molecular processes, capturing quantum (de)coherence effects rigorously. These terms account for the coupling among the trajectories, recovering the nonlocal nature of quantum nuclear dynamics that is completely overlooked in traditional independent-trajectory algorithms. Nevertheless, some of the approximations introduced in the derivation of some of these methods do not conserve the total energy. We analyze energy conservation in the coupled-trajectory mixed quantum-classical (CTMQC) algorithm and explore the performance of a modified algorithm, CTMQC-E, where some of the terms are redefined to restore energy conservation. A set of molecular models is used as a test, namely, 2-cis-penta-2,4-dienimium cation, bis(methylene) adamantyl radical cation, butatriene cation, uracil radical cation, and neutral pyrazine.
Collapse
Affiliation(s)
| | - Lea M Ibele
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
6
|
Villaseco Arribas E, Vindel-Zandbergen P, Roy S, Maitra NT. Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics. Phys Chem Chem Phys 2023; 25:26380-26395. [PMID: 37750820 DOI: 10.1039/d3cp03464j] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2023]
Abstract
The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.
Collapse
Affiliation(s)
| | - Patricia Vindel-Zandbergen
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
- Department of Chemistry, New York University, New York, New York 10003, USA
| | - Saswata Roy
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark 07102, New Jersey, USA.
| |
Collapse
|
7
|
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
|