1
|
Gay-Balmaz F, Tronci C. Complex Fluid Models of Mixed Quantum-Classical Dynamics. JOURNAL OF NONLINEAR SCIENCE 2024; 34:81. [PMID: 38966687 PMCID: PMC11219418 DOI: 10.1007/s00332-024-10044-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Accepted: 05/08/2024] [Indexed: 07/06/2024]
Abstract
Several methods in nonadiabatic molecular dynamics are based on Madelung's hydrodynamic description of nuclear motion, while the electronic component is treated as a finite-dimensional quantum system. In this context, the quantum potential leads to severe computational challenges and one often seeks to neglect its contribution, thereby approximating nuclear motion as classical. The resulting model couples classical hydrodynamics for the nuclei to the quantum motion of the electronic component, leading to the structure of a complex fluid system. This type of mixed quantum-classical fluid models has also appeared in solvation dynamics to describe the coupling between liquid solvents and the quantum solute molecule. While these approaches represent a promising direction, their mathematical structure requires a certain care. In some cases, challenging higher-order gradients make these equations hardly tractable. In other cases, these models are based on phase-space formulations that suffer from well-known consistency issues. Here, we present a new complex fluid system that resolves these difficulties. Unlike common approaches, the current system is obtained by applying the fluid closure at the level of the action principle of the original phase-space model. As a result, the system inherits a Hamiltonian structure and retains energy/momentum balance. After discussing some of its structural properties and dynamical invariants, we illustrate the model in the case of pure-dephasing dynamics. We conclude by presenting some invariant planar models.
Collapse
Affiliation(s)
- François Gay-Balmaz
- Division of Mathematical Sciences, Nanyang Technological University, Singapore, Singapore
| | - Cesare Tronci
- Department of Mathematics, University of Surrey, Guildford, UK
| |
Collapse
|
2
|
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
|
3
|
Santiago J, de Faria JC, San-Miguel M, Bernal MA. A TD-DFT-Based Study on the Attack of the OH· Radical on a Guanine Nucleotide. Int J Mol Sci 2022; 23:ijms231710007. [PMID: 36077404 PMCID: PMC9456168 DOI: 10.3390/ijms231710007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2022] [Revised: 08/24/2022] [Accepted: 08/26/2022] [Indexed: 11/16/2022] Open
Abstract
Heavy charged particles induce severe damage in DNA, which is a radiobiological advantage when treating radioresistant tumors. However, these particles can also induce cancer in humans exposed to them, such as astronauts in space missions. This damage can be directly induced by the radiation or indirectly by the attack of free radicals mainly produced by water radiolysis. We previously studied the impact of a proton on a DNA base pair, using the Time Dependent-Density Functional Theory (TD-DFT). In this work, we go a step further and study the attack of the OH· radical on the Guanine nucleotide to unveil how this molecule subsequently dissociates. The OH· attack on the H1', H2', H3', and H5' atoms in the guanine was investigated using the Ehrenfest dynamics within the TD-DFT framework. In all cases, the hydrogen abstraction succeeded, and the subsequent base pair dissociation was observed. The DNA dissociates in three major fragments: the phosphate group, the deoxyribose sugar, and the nitrogenous base, with slight differences, no matter which hydrogen atom was attacked. Hydrogen abstraction occurs at about 6 fs, and the nucleotide dissociation at about 100 fs, which agrees with our previous result for the direct proton impact on the DNA. These calculations may be a reference for adjusting reactive force fields so that more complex DNA structures can be studied using classical molecular dynamics, including both direct and indirect DNA damage.
Collapse
Affiliation(s)
- João Santiago
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, Campinas 13083-859, SP, Brazil
| | - Jhaison C. de Faria
- Instituto de Química, Universidade Estadual de Campinas, Campinas 13083-970, SP, Brazil
| | - Miguel San-Miguel
- Instituto de Química, Universidade Estadual de Campinas, Campinas 13083-970, SP, Brazil
| | - Mario A. Bernal
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, Campinas 13083-859, SP, Brazil
- Correspondence:
| |
Collapse
|
4
|
Sun J, Sasmal S, Vendrell O. A bosonic perspective on the classical mapping of fermionic quantum dynamics. J Chem Phys 2021; 155:134110. [PMID: 34624967 DOI: 10.1063/5.0066740] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We consider the application of the original Meyer-Miller (MM) Hamiltonian to mapping fermionic quantum dynamics to classical equations of motion. Non-interacting fermionic and bosonic systems share the same one-body density dynamics when evolving from the same initial many-body state. The MM classical mapping is exact for non-interacting bosons, and therefore, it yields the exact time-dependent one-body density for non-interacting fermions as well. Starting from this observation, the MM mapping is compared to different mappings specific for fermionic systems, namely, the spin mapping with and without including a Jordan-Wigner transformation and the Li-Miller mapping (LMM). For non-interacting systems, the inclusion of fermionic anti-symmetry through the Jordan-Wigner transform does not lead to any improvement in the performance of the mappings, and instead, it worsens the classical description. For an interacting impurity model and for models of excitonic energy transfer, the MM and LMM mappings perform similarly, and in some cases, the former outperforms the latter when compared to a full quantum description. The classical mappings are able to capture interference effects, both constructive and destructive, that originate from equivalent energy transfer pathways in the models.
Collapse
Affiliation(s)
- Jing Sun
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| | - Sudip Sasmal
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| | - Oriol Vendrell
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Heidelberg, Germany
| |
Collapse
|
5
|
Rassolov V, Garashchuk S. Local Measure of Quantum Effects in Quantum Dynamics. J Phys Chem A 2021; 125:4653-4667. [PMID: 34014096 DOI: 10.1021/acs.jpca.1c02533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The Madelung-de Broglie-Bohm formulation of the Schrödinger equation casts the time-evolution of a wave function as dynamics of an ensemble of quantum, or Bohmian, trajectories, interacting via the nonlocal quantum potential. This trajectory perspective gives insight into the quantumness (or classicality) of a given system due to clear partitioning of the energy into classical and quantum components. Here, we propose a system-independent measure of the quantumness of dynamics, based on the energy time-change, referred to as "quantum power". This measure is local in the coordinate space. Based on applications to model chemical systems, we argue that during the transition from the quantum to classical regime, defined as compression of quantization, the quantum features in dynamics do not "disappear" but are pushed forward in time. This feature may be used to gauge the validity of the semiclassical and other approximate dynamics approaches in applications to anharmonic systems.
Collapse
Affiliation(s)
- Vitaly Rassolov
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Sophya Garashchuk
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| |
Collapse
|
6
|
Farfan CA, Turner DB. A systematic model study quantifying how conical intersection topography modulates photochemical reactions. Phys Chem Chem Phys 2020; 22:20265-20283. [PMID: 32966428 DOI: 10.1039/d0cp03464a] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.
Collapse
Affiliation(s)
- Camille A Farfan
- Department of Chemistry, New York University, New York, NY 10003, USA
| | - Daniel B Turner
- Department of Chemistry, New York University, New York, NY 10003, USA
| |
Collapse
|
7
|
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
|
8
|
You P, Chen D, Lian C, Zhang C, Meng S. First‐principles dynamics of photoexcited molecules and materials towards a quantum description. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1492] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Peiwei You
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics Chinese Academy of Sciences Beijing China
- School of Physical Sciences University of Chinese Academy of Sciences Beijing China
| | - Daqiang Chen
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics Chinese Academy of Sciences Beijing China
- School of Physical Sciences University of Chinese Academy of Sciences Beijing China
| | - Chao Lian
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics Chinese Academy of Sciences Beijing China
| | - Cui Zhang
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics Chinese Academy of Sciences Beijing China
- Songshan Lake Materials Laboratory Dongguan China
| | - Sheng Meng
- Beijing National Laboratory for Condensed Matter Physics and Institute of Physics Chinese Academy of Sciences Beijing China
- School of Physical Sciences University of Chinese Academy of Sciences Beijing China
- Songshan Lake Materials Laboratory Dongguan China
| |
Collapse
|
9
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
10
|
Prlj A, Begušić T, Zhang ZT, Fish GC, Wehrle M, Zimmermann T, Choi S, Roulet J, Moser JE, Vaníček J. Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene. J Chem Theory Comput 2020; 16:2617-2626. [PMID: 32119547 DOI: 10.1021/acs.jctc.0c00079] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Azulene is a prototypical molecule with an anomalous fluorescence from the second excited electronic state, thus violating Kasha's rule, and with an emission spectrum that cannot be understood within the Condon approximation. To better understand the photophysics and spectroscopy of azulene and other nonconventional molecules, we developed a systematic, general, and efficient computational approach combining the semiclassical dynamics of nuclei with ab initio electronic structure. First, to analyze the nonadiabatic effects, we complement the standard population dynamics by a rigorous measure of adiabaticity, estimated with the multiple-surface dephasing representation. Second, we propose a new semiclassical method for simulating non-Condon spectra, which combines the extended thawed Gaussian approximation with the efficient single-Hessian approach. S1 ← S0 and S2 ← S0 absorption and S2 → S0 emission spectra of azulene, recorded in a new set of experiments, agree very well with our calculations. We find that accuracy of the evaluated spectra requires the treatment of anharmonicity, Herzberg-Teller, and mode-mixing effects.
Collapse
Affiliation(s)
- Antonio Prlj
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Tomislav Begušić
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Zhan Tong Zhang
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - George Cameron Fish
- Photochemical Dynamics Group, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Marius Wehrle
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Tomáš Zimmermann
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Julien Roulet
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jacques-Edouard Moser
- Photochemical Dynamics Group, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
11
|
Ghosh S, Giannini S, Lively K, Blumberger J. Nonadiabatic dynamics with quantum nuclei: simulating charge transfer with ring polymer surface hopping. Faraday Discuss 2020; 221:501-525. [DOI: 10.1039/c9fd00046a] [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/31/2022]
Abstract
Exploring effects of quantizing nuclei in non-adiabatic dynamics for simulating charge transfer in a dimer of “ethylene-like-molecules” at different temperatures.
Collapse
Affiliation(s)
- Soumya Ghosh
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
| | - Samuele Giannini
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
| | - Kevin Lively
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
| | - Jochen Blumberger
- Department of Physics and Astronomy
- University College London
- London WC1E 6BT
- UK
- Institute for Advanced Study
| |
Collapse
|
12
|
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
|
13
|
Carof A, Giannini S, Blumberger J. How to calculate charge mobility in molecular materials from surface hopping non-adiabatic molecular dynamics - beyond the hopping/band paradigm. Phys Chem Chem Phys 2019; 21:26368-26386. [PMID: 31793569 DOI: 10.1039/c9cp04770k] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Charge transport in high mobility organic semiconductors is in an intermediate regime between small polaron hopping and band transport limits. We have recently shown that surface hopping non-adiabatic molecular dynamics is a powerful method for prediction of charge transport mechanisms in organic materials and for near-quantitative prediction of charge mobilities at room temperature where the effects of nuclear zero-point motion and tunneling are still relatively small [S. Giannini et al., Nat. Commun., 2019, 10, 3843]. Here we assess and critically discuss the extensions to Tully's original method that have led to this success: (i) correction for missing electronic decoherence, (ii) detection of trivial crossings and (iii) removal of decoherence correction-induced spurious charge transfer. If any one of these corrections is not included, the charge mobility diverges with system size, each for different physical reasons. Yet if they are included, convergence with system size, detailed balance and good internal consistency are achieved.
Collapse
Affiliation(s)
- Antoine Carof
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK.
| | - Samuele Giannini
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK.
| | - Jochen Blumberger
- Department of Physics and Astronomy, University College London, London WC1E 6BT, UK. and Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 2 a, D-85748 Garching, Germany
| |
Collapse
|
14
|
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
|
15
|
Li X, Hu D, Xie Y, Lan Z. Analysis of trajectory similarity and configuration similarity in on-the-fly surface-hopping simulation on multi-channel nonadiabatic photoisomerization dynamics. J Chem Phys 2018; 149:244104. [DOI: 10.1063/1.5048049] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Xusong Li
- CAS Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
- Sino-Danish Center for Education and Research/Sino-Danish College, University of Chinese Academy of Sciences, Beijing 100049, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Deping Hu
- CAS Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yu Xie
- The Environmental Research Institute, MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- CAS Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
- Sino-Danish Center for Education and Research/Sino-Danish College, University of Chinese Academy of Sciences, Beijing 100049, China
- The Environmental Research Institute, MOE Key Laboratory of Theoretical Chemistry of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
16
|
Mignolet B, Curchod BFE. A walk through the approximations of ab initio multiple spawning. J Chem Phys 2018; 148:134110. [PMID: 29626896 DOI: 10.1063/1.5022877] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Full multiple spawning offers an in principle exact framework for excited-state dynamics, where nuclear wavefunctions in different electronic states are represented by a set of coupled trajectory basis functions that follow classical trajectories. The couplings between trajectory basis functions can be approximated to treat molecular systems, leading to the ab initio multiple spawning method which has been successfully employed to study the photochemistry and photophysics of several molecules. However, a detailed investigation of its approximations and their consequences is currently missing in the literature. In this work, we simulate the explicit photoexcitation and subsequent excited-state dynamics of a simple system, LiH, and we analyze (i) the effect of the ab initio multiple spawning approximations on different observables and (ii) the convergence of the ab initio multiple spawning results towards numerically exact quantum dynamics upon a progressive relaxation of these approximations. We show that, despite the crude character of the approximations underlying ab initio multiple spawning for this low-dimensional system, the qualitative excited-state dynamics is adequately captured, and affordable corrections can further be applied to ameliorate the coupling between trajectory basis functions.
Collapse
Affiliation(s)
- Benoit Mignolet
- Theoretical Physical Chemistry, UR MolSYS, B6c, University of Liège, B4000 Liège, Belgium
| | - Basile F E Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
17
|
Tzemos AC, Efthymiopoulos C, Contopoulos G. Origin of chaos near three-dimensional quantum vortices: A general Bohmian theory. Phys Rev E 2018; 97:042201. [PMID: 29758605 DOI: 10.1103/physreve.97.042201] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2017] [Indexed: 11/07/2022]
Abstract
We provide a general theory for the structure of the quantum flow near three-dimensional (3D) nodal lines, i.e., one-dimensional loci where the 3D wave function becomes equal to zero. In suitably defined coordinates (comoving with the nodal line) the generic structure of the flow implies the formation of 3D quantum vortices. We show that such vortices are accompanied by nearby invariant lines of the comoving quantum flow, called X lines, which are normally hyperbolic. Furthermore, the stable and unstable manifolds of the X lines produce chaotic scatterings of nearby quantum (Bohmian) trajectories, thus inducing an intricate form of the quantum current in the neighborhood of each 3D quantum vortex. Generic formulas describing the structure around 3D quantum vortices are provided, applicable to an arbitrary choice of 3D wave function. We also give specific numerical examples as well as a discussion of the physical consequences of chaos near 3D quantum vortices.
Collapse
Affiliation(s)
- Athanasios C Tzemos
- Research Center for Astronomy and Applied Mathematics, Academy of Athens, 4 Soranou Efesiou Street, Athens 11527, Greece
| | - Christos Efthymiopoulos
- Research Center for Astronomy and Applied Mathematics, Academy of Athens, 4 Soranou Efesiou Street, Athens 11527, Greece
| | - George Contopoulos
- Research Center for Astronomy and Applied Mathematics, Academy of Athens, 4 Soranou Efesiou Street, Athens 11527, Greece
| |
Collapse
|
18
|
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
|
19
|
Ha JK, Lee IS, Min SK. Surface Hopping Dynamics beyond Nonadiabatic Couplings for Quantum Coherence. J Phys Chem Lett 2018; 9:1097-1104. [PMID: 29439572 DOI: 10.1021/acs.jpclett.8b00060] [Citation(s) in RCA: 64] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
Description of correct electron-nuclear couplings is crucial in modeling of nonadiabatic dynamics. Within traditional semiclassical or mixed quantum-classical dynamics, the coupling between quantum electronic states and classical nuclear trajectories is governed by nonadiabatic coupling vectors coupled to classical nuclear momenta. This enables us to develop a very powerful nonadiabatic dynamics algorithm, namely, surface hopping dynamics, which can describe the splitting of nuclear wave packets and detailed balance. Despite its efficiency and practicality, it suffers from the lack of quantum decoherence due to incorrect accounts for the electron-nuclear coupling. Here we present a new surface hopping algorithm based on the exact electron-nuclear correlation from the exact factorization of molecular wave functions. This algorithm demands comparable computational costs to existing surface hopping methods. Numerical simulations with two-state models and a multidimensional multistate realistic molecule show that the electron-nuclear coupling beyond the nonadiabatic coupling terms can describe the quantum coherence properly.
Collapse
Affiliation(s)
- Jong-Kwon Ha
- Department of Chemistry, School of Natural Science, Ulsan National Institute of Science and Technology (UNIST) , 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - In Seong Lee
- Department of Chemistry, School of Natural Science, Ulsan National Institute of Science and Technology (UNIST) , 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - 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
| |
Collapse
|
20
|
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
|
21
|
Antipov SV, Bhattacharyya S, El Hage K, Xu ZH, Meuwly M, Rothlisberger U, Vaníček J. Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2017; 4:061509. [PMID: 29376107 PMCID: PMC5758379 DOI: 10.1063/1.4996559] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2017] [Accepted: 10/11/2017] [Indexed: 06/07/2023]
Abstract
Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.
Collapse
Affiliation(s)
- Sergey V Antipov
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Swarnendu Bhattacharyya
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Krystel El Hage
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Zhen-Hao Xu
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
22
|
Li X, Xie Y, Hu D, Lan Z. Analysis of the Geometrical Evolution in On-the-Fly Surface-Hopping Nonadiabatic Dynamics with Machine Learning Dimensionality Reduction Approaches: Classical Multidimensional Scaling and Isometric Feature Mapping. J Chem Theory Comput 2017; 13:4611-4623. [DOI: 10.1021/acs.jctc.7b00394] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Xusong Li
- CAS
Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
- Sino-Danish
Center for Education and Research/Sino-Danish College, University of Chinese Academy of Sciences, Beijing 100049, China
- University of Chinese Academy of Sciences, Beijing 100049, China
- Department
of Chemistry, Technical University of Denmark, Kemitorvet 207, 2800, Kgs. Lyngby, Denmark
| | - Yu Xie
- CAS
Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
| | - Deping Hu
- CAS
Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zhenggang Lan
- CAS
Key Laboratory of Biobased Materials, Qingdao Institute of Bioenergy
and Bioprocess Technology, Chinese Academy of Sciences, Qingdao 266101, China
- Sino-Danish
Center for Education and Research/Sino-Danish College, University of Chinese Academy of Sciences, Beijing 100049, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
23
|
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
|
24
|
Gu B, Franco I. Partial hydrodynamic representation of quantum molecular dynamics. J Chem Phys 2017; 146:194104. [PMID: 28527433 PMCID: PMC5648576 DOI: 10.1063/1.4983495] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2017] [Accepted: 05/03/2017] [Indexed: 11/14/2022] Open
Abstract
A hybrid method is proposed to propagate system-bath quantum dynamics that use both basis functions and coupled quantum trajectories. In it, the bath is represented with an ensemble of Bohmian trajectories while the system degrees of freedom are accounted through reduced density matrices. By retaining the Hilbert space structure for the system, the method is able to capture interference processes that are challenging to describe in Bohmian dynamics due to singularities that these processes introduce in the quantum potential. By adopting quantum trajectories to represent the bath, the method beats the exponential scaling of the computational cost with the bath size. This combination makes the method suitable for large-scale ground and excited state fully quantum molecular dynamics simulations. Equations of motion for the quantum trajectories and reduced density matrices are derived from the Schrödinger equation and a computational algorithm to solve these equations is proposed. Through computations in two-dimensional model systems, the method is shown to offer an accurate description of subsystem observables and of quantum decoherence, which is difficult to obtain when the quantum nature of the bath is ignored. The scaling of the method is demonstrated using a model with 21 degrees of freedom. The limit of independent trajectories is recovered when the mass of bath degrees of freedom is much larger than the one of the system, in agreement with mixed quantum-classical descriptions.
Collapse
Affiliation(s)
- Bing Gu
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Ignacio Franco
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
25
|
Heaps CW, Mazziotti DA. Accurate non-adiabatic quantum dynamics from pseudospectral sampling of time-dependent Gaussian basis sets. J Chem Phys 2016. [DOI: 10.1063/1.4959872] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Charles W. Heaps
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - David A. Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
26
|
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
|
27
|
Min SK, Agostini F, Gross EKU. Coupled-Trajectory Quantum-Classical Approach to Electronic Decoherence in Nonadiabatic Processes. PHYSICAL REVIEW LETTERS 2015; 115:073001. [PMID: 26317716 DOI: 10.1103/physrevlett.115.073001] [Citation(s) in RCA: 92] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/01/2015] [Indexed: 06/04/2023]
Abstract
We present a novel quantum-classical approach to nonadiabatic dynamics, deduced from the coupled electronic and nuclear equations in the framework of the exact factorization of the electron-nuclear wave function. The method is based on the quasiclassical interpretation of the nuclear wave function, whose phase is related to the classical momentum and whose density is represented in terms of classical trajectories. In this approximation, electronic decoherence is naturally induced as an effect of the coupling to the nuclei and correctly reproduces the expected quantum behavior. Moreover, the splitting of the nuclear wave packet is captured as a consequence of the correct approximation of the time-dependent potential of the theory. This new approach offers a clear improvement over Ehrenfest-like dynamics. The theoretical derivation presented in this Letter is supported by numerical results that are compared to quantum mechanical calculations.
Collapse
Affiliation(s)
- Seung Kyu Min
- Max-Planck Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
- Department of Chemistry, School of Natural Science, Ulsan National Institute of Science and Technology (UNIST), Ulsan 689-798, Korea
| | - 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
|
28
|
Agostini F, Abedi A, Suzuki Y, Min SK, Maitra NT, Gross EKU. The exact forces on classical nuclei in non-adiabatic charge transfer. J Chem Phys 2015; 142:084303. [PMID: 25725727 DOI: 10.1063/1.4908133] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The decomposition of electronic and nuclear motion presented in Abedi et al. [Phys. Rev. Lett. 105, 123002 (2010)] yields a time-dependent potential that drives the nuclear motion and fully accounts for the coupling to the electronic subsystem. Here, we show that propagation of an ensemble of independent classical nuclear trajectories on this exact potential yields dynamics that are essentially indistinguishable from the exact quantum dynamics for a model non-adiabatic charge transfer problem. We point out the importance of step and bump features in the exact potential that are critical in obtaining the correct splitting of the quasiclassical nuclear wave packet in space after it passes through an avoided crossing between two Born-Oppenheimer surfaces and analyze their structure. Finally, an analysis of the exact potentials in the context of trajectory surface hopping is presented, including preliminary investigations of velocity-adjustment and the force-induced decoherence effect.
Collapse
Affiliation(s)
- Federica Agostini
- Max-Planck-Institut of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| | - Ali Abedi
- Max-Planck-Institut of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| | - Yasumitsu Suzuki
- Max-Planck-Institut of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| | - Seung Kyu Min
- Max-Planck-Institut of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| | - 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
| | - E K U Gross
- Max-Planck-Institut of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| |
Collapse
|
29
|
Albareda G, Bofill JM, Tavernelli I, Huarte-Larrañaga F, Illas F, Rubio A. Conditional Born-Oppenheimer Dynamics: Quantum Dynamics Simulations for the Model Porphine. J Phys Chem Lett 2015; 6:1529-1535. [PMID: 26263307 DOI: 10.1021/acs.jpclett.5b00422] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
We report a new theoretical approach to solve adiabatic quantum molecular dynamics halfway between wave function and trajectory-based methods. The evolution of a N-body nuclear wave function moving on a 3N-dimensional Born-Oppenheimer potential-energy hyper-surface is rewritten in terms of single-nuclei wave functions evolving nonunitarily on a 3-dimensional potential-energy surface that depends parametrically on the configuration of an ensemble of generally defined trajectories. The scheme is exact and, together with the use of trajectory-based statistical techniques, can be exploited to circumvent the calculation and storage of many-body quantities (e.g., wave function and potential-energy surface) whose size scales exponentially with the number of nuclear degrees of freedom. As a proof of concept, we present numerical simulations of a 2-dimensional model porphine where switching from concerted to sequential double proton transfer (and back) is induced quantum mechanically.
Collapse
Affiliation(s)
- Guillermo Albareda
- †Institut de Química Teòrica i Computacional, Universitat de Barcelona, Barcelona 08028, Spain
- ‡Departament de Química Física, Universitat de Barcelona, Barcelona 08028, Spain
| | - Josep Maria Bofill
- †Institut de Química Teòrica i Computacional, Universitat de Barcelona, Barcelona 08028, Spain
- §Departament de Química Orgànica, Universitat de Barcelona, Barcelona 08028, Spain
| | - Ivano Tavernelli
- ∥IBM Research Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| | - Fermin Huarte-Larrañaga
- †Institut de Química Teòrica i Computacional, Universitat de Barcelona, Barcelona 08028, Spain
- ‡Departament de Química Física, Universitat de Barcelona, Barcelona 08028, Spain
| | - Francesc Illas
- †Institut de Química Teòrica i Computacional, Universitat de Barcelona, Barcelona 08028, Spain
- ‡Departament de Química Física, Universitat de Barcelona, Barcelona 08028, Spain
| | - Angel Rubio
- ⊥Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
30
|
Brunk E, Rothlisberger U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem Rev 2015; 115:6217-63. [PMID: 25880693 DOI: 10.1021/cr500628b] [Citation(s) in RCA: 301] [Impact Index Per Article: 33.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Elizabeth Brunk
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,‡Joint BioEnergy Institute, Lawrence Berkeley National Laboratory, Emeryville, California 94618, United States
| | - Ursula Rothlisberger
- †Laboratory of Computational Chemistry and Biochemistry, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland.,§National Competence Center of Research (NCCR) MARVEL-Materials' Revolution: Computational Design and Discovery of Novel Materials, 1015 Lausanne, Switzerland
| |
Collapse
|
31
|
Agostini F, Abedi A, Gross EKU. Classical nuclear motion coupled to electronic non-adiabatic transitions. J Chem Phys 2014; 141:214101. [DOI: 10.1063/1.4902225] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Federica Agostini
- Max-Planck Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - Ali Abedi
- 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
|
32
|
An overview of nonadiabatic dynamics simulations methods, with focus on the direct approach versus the fitting of potential energy surfaces. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1526-1] [Citation(s) in RCA: 102] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
33
|
Franco de Carvalho F, Curchod BFE, Penfold TJ, Tavernelli I. Derivation of spin-orbit couplings in collinear linear-response TDDFT: A rigorous formulation. J Chem Phys 2014; 140:144103. [DOI: 10.1063/1.4870010] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
|
34
|
González L, Marquetand P, Richter M, González-Vázquez J, Sola I. Ultrafast Laser-Induced Processes Described by Ab Initio Molecular Dynamics. ACTA ACUST UNITED AC 2014. [DOI: 10.1007/978-3-319-02051-8_7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/20/2023]
|
35
|
|
36
|
Curchod BFE, Tavernelli I. On trajectory-based nonadiabatic dynamics: Bohmian dynamics versus trajectory surface hopping. J Chem Phys 2013; 138:184112. [DOI: 10.1063/1.4803835] [Citation(s) in RCA: 57] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
37
|
Curchod BFE, Rothlisberger U, Tavernelli I. Trajectory-Based Nonadiabatic Dynamics with Time-Dependent Density Functional Theory. Chemphyschem 2013; 14:1314-40. [DOI: 10.1002/cphc.201200941] [Citation(s) in RCA: 153] [Impact Index Per Article: 13.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2012] [Indexed: 11/11/2022]
|
38
|
Zimmermann T, Vaníček J. Evaluation of the importance of spin-orbit couplings in the nonadiabatic quantum dynamics with quantum fidelity and with its efficient "on-the-fly" ab initio semiclassical approximation. J Chem Phys 2013; 137:22A516. [PMID: 23249053 DOI: 10.1063/1.4738878] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose to measure the importance of spin-orbit couplings (SOCs) in the nonadiabatic molecular quantum dynamics rigorously with quantum fidelity. To make the criterion practical, quantum fidelity is estimated efficiently with the multiple-surface dephasing representation (MSDR). The MSDR is a semiclassical method that includes nuclear quantum effects through interference of mixed quantum-classical trajectories without the need for the Hessian of potential energy surfaces. Two variants of the MSDR are studied, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics. The fidelity criterion and MSDR are first tested on one-dimensional model systems amenable to numerically exact quantum dynamics. Then, the MSDR is combined with "on-the-fly" computed electronic structure to measure the importance of SOCs and nonadiabatic couplings in the photoisomerization dynamics of CH(2)NH(2)(+) considering 20 electronic states and in the collision of F + H(2) considering six electronic states.
Collapse
Affiliation(s)
- Tomáš Zimmermann
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | | |
Collapse
|
39
|
|
40
|
Zimmermann T, Vaníček J. Measuring nonadiabaticity of molecular quantum dynamics with quantum fidelity and with its efficient semiclassical approximation. J Chem Phys 2012; 136:094106. [PMID: 22401428 DOI: 10.1063/1.3690458] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose to measure nonadiabaticity of molecular quantum dynamics rigorously with the quantum fidelity between the Born-Oppenheimer and fully nonadiabatic dynamics. It is shown that this measure of nonadiabaticity applies in situations where other criteria, such as the energy gap criterion or the extent of population transfer, fail. We further propose to estimate this quantum fidelity efficiently with a generalization of the dephasing representation to multiple surfaces. Two variants of the multiple-surface dephasing representation (MSDR) are introduced, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics (LMFD). The LMFD can be interpreted as the Ehrenfest dynamics of an ensemble of nuclear trajectories, and has been used previously in the nonadiabatic semiclassical initial value representation. In addition to propagating an ensemble of classical trajectories, the MSDR requires evaluating nonadiabatic couplings and solving the Schrödinger (or more generally, the quantum Liouville-von Neumann) equation for a single discrete degree of freedom. The MSDR can be also used in the diabatic basis to measure the importance of the diabatic couplings. The method is tested on three model problems introduced by Tully and on a two-surface model of dissociation of NaI.
Collapse
Affiliation(s)
- Tomáš Zimmermann
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| | | |
Collapse
|
41
|
Tavernelli I, Curchod BF, Rothlisberger U. Nonadiabatic molecular dynamics with solvent effects: A LR-TDDFT QM/MM study of ruthenium (II) tris (bipyridine) in water. Chem Phys 2011. [DOI: 10.1016/j.chemphys.2011.03.021] [Citation(s) in RCA: 89] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|