1
|
Giarrusso S, Gori-Giorgi P, Agostini F. Electronic Vector Potential from the Exact Factorization of a Complex Wavefunction. Chemphyschem 2024; 25: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
|
Sangiogo Gil E, Lauvergnat D, Agostini F. Exact factorization of the photon-electron-nuclear wavefunction: Formulation and coupled-trajectory dynamics. J Chem Phys 2024; 161:084112. [PMID: 39189656 DOI: 10.1063/5.0224779] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2024] [Accepted: 07/29/2024] [Indexed: 08/28/2024] Open
Abstract
We employ the exact-factorization formalism to study the coupled dynamics of photons, electrons, and nuclei at the quantum mechanical level, proposing illustrative examples of model situations of nonadiabatic dynamics and spontaneous emission of electron-nuclear systems in the regime of strong light-matter coupling. We make a particular choice of factorization for such a multi-component system, where the full wavefunction is factored as a conditional electronic amplitude and a marginal photon-nuclear amplitude. Then, we apply the coupled-trajectory mixed quantum-classical (CTMQC) algorithm to perform trajectory-based simulations, by treating photonic and nuclear degrees of freedom on equal footing in terms of classical-like trajectories. The analysis of the time-dependent potentials of the theory along with the assessment of the performance of CTMQC allows us to point out some limitations of the current approximations used in CTMQC. Meanwhile, comparing CTMQC with other trajectory-based algorithms, namely multi-trajectory Ehrenfest and Tully surface hopping, demonstrates the better quality of CTMQC predictions.
Collapse
Affiliation(s)
- Eduarda Sangiogo Gil
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
- Institute of Theoretical Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - David Lauvergnat
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| | - Federica Agostini
- CNRS, Institut de Chimie Physique UMR8000, Université Paris-Saclay, 91405 Orsay, France
| |
Collapse
|
3
|
Karak P, Manna P, Banerjee A, Ruud K, Chakrabarti S. Reverse Intersystem Crossing Dynamics in Vibronically Modulated Inverted Singlet-Triplet Gap System: A Wigner Phase Space Study. J Phys Chem Lett 2024; 15:7603-7609. [PMID: 39028946 DOI: 10.1021/acs.jpclett.4c01437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/21/2024]
Abstract
We inspect the origin of the inverted singlet-triplet gap (INVEST) and slow change in the reverse intersystem crossing (rISC) rate with temperature, as recently observed. A Wigner phase space study reveals that, though INVEST is found at equilibrium geometry, variation in the exchange interaction and the doubles-excitation for other geometries in the harmonic region leads to non-INVEST behavior. This highlights the importance of nuclear degrees of freedom for the INVEST phenomenon, and in this case, geometric puckering of the studied molecule determines INVEST and the associated rISC dynamics.
Collapse
Affiliation(s)
- Pijush Karak
- Department of Chemistry, University of Calcutta, 92 A. P. C. Road, Kolkata 700009, West Bengal, India
| | - Pradipta Manna
- Department of Chemistry, University of Calcutta, 92 A. P. C. Road, Kolkata 700009, West Bengal, India
| | - Ambar Banerjee
- Research Institute for Sustainable Energy (RISE), TCG Centres for Research and Education in Science and Technology (TCG-CREST), Kolkata 700091, India
| | - Kenneth Ruud
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Tromsø - The Arctic University of Norway, 9037 Tromsø, Norway
- Norwegian Defence Research Establishment, P.O. Box 25, 2027 Kjeller, Norway
| | - Swapan Chakrabarti
- Department of Chemistry, University of Calcutta, 92 A. P. C. Road, Kolkata 700009, West Bengal, India
| |
Collapse
|
4
|
Toldo JM, do Casal MT, Ventura E, do Monte SA, Barbatti M. Surface hopping modeling of charge and energy transfer in active environments. Phys Chem Chem Phys 2023; 25:8293-8316. [PMID: 36916738 PMCID: PMC10034598 DOI: 10.1039/d3cp00247k] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 03/01/2023] [Indexed: 03/05/2023]
Abstract
An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.
Collapse
Affiliation(s)
| | | | - Elizete Ventura
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Silmar A do Monte
- Departamento de Química, CCEN, Universidade Federal da Paraíba, 58059-900, João Pessoa, Brazil.
| | - Mario Barbatti
- Aix-Marseille University, CNRS, ICR, Marseille, France.
- Institut Universitaire de France, 75231, Paris, 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
|
Liu W, Wang Z, Chen Z, Luo J, Li S, Wang L. Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1577] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Wen‐Hao Liu
- State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors Chinese Academy of Sciences Beijing China
- Center of Materials Science and Optoelectronics Engineering University of Chinese Academy of Sciences Beijing China
| | - Zhi Wang
- State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors Chinese Academy of Sciences Beijing China
| | - Zhang‐Hui Chen
- Materials Science Division Lawrence Berkeley National Laboratory Berkeley California USA
| | - Jun‐Wei Luo
- State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors Chinese Academy of Sciences Beijing China
- Center of Materials Science and Optoelectronics Engineering University of Chinese Academy of Sciences Beijing China
- Beijing Academy of Quantum Information Sciences Beijing China
| | - Shu‐Shen Li
- State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors Chinese Academy of Sciences Beijing China
- Center of Materials Science and Optoelectronics Engineering University of Chinese Academy of Sciences Beijing China
- Beijing Academy of Quantum Information Sciences Beijing China
| | - Lin‐Wang Wang
- Materials Science Division Lawrence Berkeley National Laboratory Berkeley California USA
| |
Collapse
|
7
|
Mukherjee S, Varganov SA. Intersystem crossing and internal conversion dynamics with GAIMS-TeraChem: Excited state relaxation in 2-cyclopentenone. J Chem Phys 2021; 155:174107. [PMID: 34742200 DOI: 10.1063/5.0068040] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Excited states relaxation in complex molecules often involves two types of nonradiative transitions, internal conversion (IC) and intersystem crossing (ISC). In the situations when the timescales of IC and ISC are comparable, an interplay between these two types of transitions can lead to complex nonadiabatic dynamics on multiple electronic states of different characters and spin multiplicities. We demonstrate that the generalized ab initio multiple spawning (GAIMS) method interfaced with the fast graphics processing unit-based TeraChem electronic structure code can be used to model such nonadiabatic dynamics involving both the IC and ISC transitions in molecules of moderate size. We carried out 1500 fs GAIMS simulations leading to the creation of up to 2500 trajectory basis functions to study the excited states relaxation in 2-cyclopentenone. After a vertical excitation from the ground state to the bright S2 state, the molecule quickly relaxes to the S1 state via conical intersection. The following relaxation proceeds along two competing pathways: one involves IC to the ground state, and the other is dominated by ISC to the low-lying triplet states. The time constants describing the population transfer between the six lowest singlet and triplet states predicted by the GAIMS dynamics are in good agreement with the characteristic times of IC and ISC obtained from the analysis of the time-resolved photoelectron spectrum.
Collapse
Affiliation(s)
- Saikat Mukherjee
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA
| | - Sergey A Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, USA
| |
Collapse
|
8
|
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
|
9
|
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
|
10
|
Yang JJ, Liu XY, Li ZW, Frauenheim T, Yam C, Fang WH, Cui G. The spin-orbit interaction controls photoinduced interfacial electron transfer in fullerene-perovskite heterojunctions: C 60versus C 70. Phys Chem Chem Phys 2021; 23:6536-6543. [PMID: 33690742 DOI: 10.1039/d0cp06579j] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Here, we used collinear and noncollinear density functional theory (DFT) methods to explore the interfacial properties of two heterojunctions between a fullerene (C60 and C70) and the MAPbI3(110) surface. Methodologically, consideration of the spin-orbit interaction has been proven to be required to obtain accurate energy-level alignment and interfacial carrier dynamics between fullerenes and perovskites in hybrid perovskite solar cells including heavy atoms (such as Pb atoms). Both heterojunctions are predicted to be the same type-II heterojunction, but the interfacial electron transfer process from MAPbI3 to C60 is completely distinct from that to C70. In the former, the interfacial electron transfer is slow because of the associated large energy gap, and the excited electrons are thus trapped in MAPbI3 for a while. In contrast, in the latter, the smaller energy gap induces ultrafast electron transfer from MAPbI3 to C70. These points are further supported by DFT-based nonadiabatic dynamics simulations including the spin-orbit coupling (SOC) effects. These gained insights could help rationally design superior fullerene-perovskite interfaces to achieve high power conversion efficiencies of fullerene-perovskite solar cells.
Collapse
Affiliation(s)
- Jia-Jia Yang
- College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | | | | | | | | | | | | |
Collapse
|
11
|
Liu XY, Yang JJ, Chen WK, Akimov AV, Fang WH, Cui G. Spin-Orbit Coupling Accelerates the Photoinduced Interfacial Electron Transfer in a Fullerene-Based Perovskite Heterojunction. J Phys Chem Lett 2021; 12:1131-1137. [PMID: 33475363 DOI: 10.1021/acs.jpclett.0c03782] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In this work, we explore the interfacial properties of the C60-Py@MAPbI3 heterojunction of the PbI-terminated MAPbI3(001) surface and pyridine-functionalized C60-Py fullerene derivative through both collinear and noncollinear density functional theory calculations with and without spin-orbit coupling (SOC) effects. C60-Py is bound to the MAPbI3 surface through interfacial Pb-O and Pb-N bonds. Although C60-Py@MAPbI3 is predicted to be the same type II heterojunction at all of the computational levels considered, the SOC effects largely decrease the energy gap of the first conduction bands of C60-Py and MAPbI3, thereby accelerating the interfacial electron transfer. Further dynamics simulations show that the inclusion of the SOC effects induces the transfer of approximately 80% of electrons from MAPbI3 to C60-Py within 1 ps. The work demonstrates that the SOC effects are indispensable for the interfacial properties of C60-Py@MAPbI3 and could also play a non-negligible role in tuning the optoelectronic properties of fullerene-based or similar perovskite devices.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Jia-Jia Yang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
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
|
Marsili E, Olivucci M, Lauvergnat D, Agostini F. Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model. J Chem Theory Comput 2020; 16:6032-6048. [PMID: 32931266 DOI: 10.1021/acs.jctc.0c00679] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.
Collapse
Affiliation(s)
- Emanuele Marsili
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France.,Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Via A. Moro 2, I-53100 Siena, Italy.,Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| |
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
|