1
|
Chen Z, Yang Y. Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface. J Phys Chem Lett 2023; 14:279-286. [PMID: 36595586 DOI: 10.1021/acs.jpclett.2c02905] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The accurate incorporation of nuclear quantum effects in large-scale molecular dynamics (MD) simulations remains a significant challenge. Recently, we combined constrained nuclear-electronic orbital (CNEO) theory with classical MD and obtained a new approach (CNEO-MD) that can accurately and efficiently incorporate nuclear quantum effects into classical simulations. In this Letter, we provide the theoretical foundation for CNEO-MD by developing an alternative formulation of the equations of motion for MD. In this new formulation, the expectation values of quantum nuclear positions evolve classically on an effective energy surface that is obtained from a constrained energy minimization procedure when solving for the quantum nuclear wave function, thus enabling the incorporation of nuclear quantum effects in classical MD simulations. For comparison with other existing approaches, we examined a series of model systems and found that this new MD approach is significantly more accurate than the conventional way of performing classical MD and generally outperforms centroid MD and ring-polymer MD in describing vibrations in these model systems.
Collapse
Affiliation(s)
- Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| |
Collapse
|
2
|
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
|
3
|
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
|
4
|
Bouakline F, Saalfrank P. Seemingly asymmetric atom-localized electronic densities following laser-dissociation of homonuclear diatomics. J Chem Phys 2021; 154:234305. [PMID: 34241262 DOI: 10.1063/5.0049710] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Recent experiments on laser-dissociation of aligned homonuclear diatomic molecules show an asymmetric forward-backward (spatial) electron-localization along the laser polarization axis. Most theoretical models attribute this asymmetry to interference effects between gerade and ungerade vibronic states. Presumably due to alignment, these models neglect molecular rotations and hence infer an asymmetric (post-dissociation) charge distribution over the two identical nuclei. In this paper, we question the equivalence that is made between spatial electron-localization, observed in experiments, and atomic electron-localization, alluded by these theoretical models. We show that (seeming) agreement between these models and experiments is due to an unfortunate omission of nuclear permutation symmetry, i.e., quantum statistics. Enforcement of the latter requires mandatory inclusion of the molecular rotational degree of freedom, even for perfectly aligned molecules. Unlike previous interpretations, we ascribe spatial electron-localization to the laser creation of a rovibronic wavepacket that involves field-free molecular eigenstates with opposite space-inversion symmetry i.e., even and odd parity. Space-inversion symmetry breaking would then lead to an asymmetric distribution of the (space-fixed) electronic density over the forward and backward hemisphere. However, owing to the simultaneous coexistence of two indistinguishable molecular orientational isomers, our analytical and computational results show that the post-dissociation electronic density along a specified space-fixed axis is equally shared between the two identical nuclei-a result that is in perfect accordance with the principle of the indistinguishability of identical particles.
Collapse
Affiliation(s)
- Foudhil Bouakline
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| | - Peter Saalfrank
- Institut für Chemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, D-14476 Potsdam-Golm, Germany
| |
Collapse
|
5
|
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
|
6
|
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: 10] [Impact Index Per Article: 2.5] [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
|
7
|
Zhou Z, Chen HT, Nitzan A, Subotnik JE. Nonadiabatic Dynamics in a Laser Field: Using Floquet Fewest Switches Surface Hopping To Calculate Electronic Populations for Slow Nuclear Velocities. J Chem Theory Comput 2020; 16:821-834. [PMID: 31951404 DOI: 10.1021/acs.jctc.9b00950] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We investigate two well-known approaches for extending the fewest switches surface hopping (FSSH) algorithm to periodic time-dependent couplings. The first formalism acts as if the instantaneous adiabatic electronic states were standard adiabatic states, which just happen to evolve in time. The second formalism replaces the role of the usual adiabatic states by the time-independent adiabatic Floquet states. For a set of modified Tully model problems, the Floquet FSSH (F-FSSH) formalism gives a better estimate for both transmission and reflection probabilities than the instantaneous adiabatic FSSH (IA-FSSH) formalism, especially for slow nuclear velocities. More importantly, only F-FSSH predicts the correct final scattering momentum. Finally, in order to use Floquet theory accurately, we find that it is crucial to account for the interference between wavepackets on different Floquet states. Our results should be of interest to all those interested in laser-induced molecular dynamics.
Collapse
Affiliation(s)
- Zeyu Zhou
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| | - Hsing-Ta Chen
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| | - Abraham Nitzan
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| | - Joseph Eli Subotnik
- Department of Chemistry , University of Pennsylvania , Philadelphia , Pennsylvania 19104 , United States
| |
Collapse
|
8
|
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
|
9
|
Santana-Bonilla A, Medrano Sandonas L, Gutierrez R, Cuniberti G. Exploring the write-in process in molecular quantum cellular automata: a combined modelingand first-principle approach. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:405502. [PMID: 31195387 DOI: 10.1088/1361-648x/ab29c1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The molecular quantum cellular automata paradigm (m-QCA) offers a promising alternative framework to current CMOS implementations. A crucial aspect for implementing this technology concerns the construction of a device which effectively controls intramolecular charge-transfer processes. Tentative experimental implementations have been developed in which a voltage drop is created generating the forces that drive a molecule into a logic state. However, important factors such as the electric field profile, its possible time-dependency and the influence of temperature in the overall success of charge-transfer are relevant issues to be considered in the design of a reliable device. In this work, we theoretically study the role played by these processes in the overall intramolecular charge-transfer process. We have used a Landau-Zener (LZ) model, where different time-dependent electric field profiles have been simulated. The results have been further corroborated employing density functional tight-binding method. The role played by the nuclear motions in the electron-transfer process has been investigated beyond the Born-Oppenheimer approximation by computing the effect of the external electric field in the behavior of the potential energy surface. Hence, we demonstrate that the intramolecular charge-transfer process is a direct consequence of the coherent LZ nonadiabatic tunneling and the hybridization of the diabatic vibronic states which effectively reduces the trapping of the itinerant electron at the donor group.
Collapse
Affiliation(s)
- Alejandro Santana-Bonilla
- Institute for Materials Science and Max Bergmann Center of Biomaterials, Dresden University of Technology, 01062 Dresden, Germany. Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany
| | | | | | | |
Collapse
|
10
|
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
|
11
|
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]
|
12
|
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
|
13
|
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
|
14
|
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
|
15
|
Khosravi E, Abedi A, Rubio A, Maitra NT. Electronic non-adiabatic dynamics in enhanced ionization of isotopologues of hydrogen molecular ions from the exact factorization perspective. Phys Chem Chem Phys 2017; 19:8269-8281. [DOI: 10.1039/c6cp08539c] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
An exact-factorization perspective of enhanced ionization in isotopologues of H2+ demonstrates the concept of the exact potential driving the electrons in non-adiabatic motion of molecules in strong fields, and sets a new platform for introducing various approximations.
Collapse
Affiliation(s)
- Elham Khosravi
- Nano-Bio Spectroscopy Group and ETSF
- Universidad del País Vasco
- CFM CSIC-UPV/EHU
- 20018 San Sebastián
- Spain
| | - Ali Abedi
- Nano-Bio Spectroscopy Group and ETSF
- Universidad del País Vasco
- CFM CSIC-UPV/EHU
- 20018 San Sebastián
- Spain
| | - Angel Rubio
- Max Planck Institute for the Structure and Dynamics of Matter
- 22761 Hamburg
- Germany
- Nano-Bio Spectroscopy Group and ETSF
- Universidad del País Vasco
| | - Neepa T. Maitra
- Department of Physics and Astronomy
- Hunter College and the Graduate Center of the City University of New York
- New York
- USA
| |
Collapse
|
16
|
Eich FG, Agostini F. The adiabatic limit of the exact factorization of the electron-nuclear wave function. J Chem Phys 2016; 145:054110. [DOI: 10.1063/1.4959962] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- F. G. Eich
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, D-22761 Hamburg, Germany
| | - Federica Agostini
- Max Planck Institute of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| |
Collapse
|
17
|
Curchod BFE, Agostini F, Gross EKU. An exact factorization perspective on quantum interferences in nonadiabatic dynamics. J Chem Phys 2016; 145:034103. [DOI: 10.1063/1.4958637] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
Agostini F, Min SK, Abedi A, Gross EKU. Quantum-Classical Nonadiabatic Dynamics: Coupled- vs Independent-Trajectory Methods. J Chem Theory Comput 2016; 12:2127-43. [PMID: 27030209 DOI: 10.1021/acs.jctc.5b01180] [Citation(s) in RCA: 84] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Trajectory-based mixed quantum-classical approaches to coupled electron-nuclear dynamics suffer from well-studied problems such as the lack of (or incorrect account for) decoherence in the trajectory surface hopping method and the inability of reproducing the spatial splitting of a nuclear wave packet in Ehrenfest-like dynamics. In the context of electronic nonadiabatic processes, these problems can result in wrong predictions for quantum populations and in unphysical outcomes for the nuclear dynamics. In this paper, we propose a solution to these issues by approximating the coupled electronic and nuclear equations within the framework of the exact factorization of the electron-nuclear wave function. We present a simple quantum-classical scheme based on coupled classical trajectories and test it against the full quantum mechanical solution from wave packet dynamics for some model situations which represent particularly challenging problems for the above-mentioned traditional methods.
Collapse
Affiliation(s)
- Federica Agostini
- Max-Planck Institut für Mikrostrukturphysik , Weinberg 2, D-06120 Halle, Germany
| | - Seung Kyu Min
- Department of Chemistry, School of Natural Science, Ulsan National Institute of Science and Technology (UNIST) , 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Ali Abedi
- Nano-Bio Spectroscopy group and European Theoretical Spectroscopy Facility (ETSF), Dpto. Física de Materiales, Universidad del País Vasco, Centro de Física de Materiales CSIC-UPV/EHU-MPC and DIPC , Av. Tolosa 72, E-20018 San Sebastián, Spain
| | - E K U Gross
- Max-Planck Institut für Mikrostrukturphysik , Weinberg 2, D-06120 Halle, Germany
| |
Collapse
|
19
|
Akimov AV. Libra: An open-Source “methodology discovery” library for quantum and classical dynamics simulations. J Comput Chem 2016; 37:1626-49. [DOI: 10.1002/jcc.24367] [Citation(s) in RCA: 51] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2015] [Revised: 02/03/2016] [Accepted: 03/02/2016] [Indexed: 01/10/2023]
Affiliation(s)
- Alexey V. Akimov
- Department of Chemistry; University at Buffalo, the State University of New York; Buffalo, New York 14260-3000
| |
Collapse
|
20
|
Schild A, Agostini F, Gross EKU. Electronic Flux Density beyond the Born–Oppenheimer Approximation. J Phys Chem A 2016; 120:3316-25. [DOI: 10.1021/acs.jpca.5b12657] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Axel Schild
- Max-Planck Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - Federica Agostini
- Max-Planck Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - E. K. U. Gross
- Max-Planck Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| |
Collapse
|
21
|
Khosravi E, Abedi A, Maitra NT. Exact Potential Driving the Electron Dynamics in Enhanced Ionization of H(2)(+). PHYSICAL REVIEW LETTERS 2015; 115:263002. [PMID: 26764989 DOI: 10.1103/physrevlett.115.263002] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/06/2015] [Indexed: 06/05/2023]
Abstract
It was recently shown that the exact factorization of the electron-nuclear wave function allows the construction of a Schrödinger equation for the electronic system, in which the potential contains exactly the effect of coupling to the nuclear degrees of freedom and any external fields. Here we study the exact potential acting on the electron in charge-resonance enhanced ionization in a model one-dimensional H(2)(+) molecule. We show there can be significant differences between the exact potential and that used in the traditional quasistatic analyses, arising from nonadiabatic coupling to the nuclear system, and that these are crucial to include for accurate simulations of time-resolved ionization dynamics and predictions of the ionization yield.
Collapse
Affiliation(s)
- Elham Khosravi
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
- Nano-Bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Universidad del País Vasco CFM CSIC-UPV/EHU-MPC and DIPC, Av. Tolosa 72, 20018 San Sebastián, Spain
| | - Ali Abedi
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
- Nano-Bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Universidad del País Vasco CFM CSIC-UPV/EHU-MPC and DIPC, Av. Tolosa 72, 20018 San Sebastián, Spain
| | - Neepa T Maitra
- Department of Physics and Astronomy, Hunter College and the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| |
Collapse
|