1
|
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
|
2
|
Acosta-Matos JC, Meier C, Martínez-Mesa A, Uranga-Piña L. Effective Phase Space Representation of the Quantum Dynamics of Vibrational Predissociation of the ArBr 2(B,ν =16···25) Complex. J Phys Chem A 2022; 126:1805-1815. [PMID: 35285621 DOI: 10.1021/acs.jpca.1c08678] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We perform trajectory-based simulations of the vibrational predissociation of the ArBr2(B,ν=16···25) van der Waals triatomic complex, constrained to the T-shape geometry. To this aim, we employ a 2-fold mapping of the quantum dynamics into classical-like dynamics in an extended phase space. The effective phase space comprises two distinct sets of degrees of freedom, namely a collection of coupled harmonic oscillators and an ensemble of quantum trajectories. The time evolution of these variables represent bound and unbound motions of the quantum system, respectively. Quantum trajectories are propagated within the interacting trajectory representation. The comparison between the lifetimes of the predissociating complexes computed using the trajectory-based approach and the experimental results available for the target systems indicates that the present method is competitive with wavepacket propagation techniques. The competition between several simultaneous vibrational relaxation pathways was found to have a direct impact on the time scales of vibrational predissociation. Likewise, the analysis of the time evolution of the trajectories reveals the existence of regions in the effective phase space where transitions to vibrational states of higher energy are more likely to occur. The size and location of these regions influence the transient vibrational distributions and therefore the computed lifetimes. Furthermore, the mechanisms of energy redistribution along the dissociation coordinate are analyzed.
Collapse
Affiliation(s)
- Juan Carlos Acosta-Matos
- Department of Physics, Centre for Technological Applications and Nuclear Developments (CEADEN), Havana 11300, Cuba
| | - Christoph Meier
- Laboratoire Collisions Agrégats Réactivité (IRSAMC), UMR 5589, Université Toulouse III - Paul Sabatier, F-31062 Toulouse Cedex 09, France
| | - Aliezer Martínez-Mesa
- Laboratoire Collisions Agrégats Réactivité (IRSAMC), UMR 5589, Université Toulouse III - Paul Sabatier, F-31062 Toulouse Cedex 09, France.,DynAMoS (Dynamical Processes in Atomic and Molecular Systems), Facultad de Física, Universidad de la Habana, Havana 10400, Cuba
| | - Llinersy Uranga-Piña
- Laboratoire Collisions Agrégats Réactivité (IRSAMC), UMR 5589, Université Toulouse III - Paul Sabatier, F-31062 Toulouse Cedex 09, France.,DynAMoS (Dynamical Processes in Atomic and Molecular Systems), Facultad de Física, Universidad de la Habana, Havana 10400, Cuba
| |
Collapse
|
3
|
Lu CY, Lee TY, Chou CC. Moving boundary truncated grid method for electronic nonadiabatic dynamics. J Chem Phys 2022; 156:044107. [DOI: 10.1063/5.0078909] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chun-Yaung Lu
- Texas Advanced Computing Center, The University of Texas at Austin, Austin, Texas 78758, USA
| | - Tsung-Yen Lee
- Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| | - Chia-Chun Chou
- Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
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
|
7
|
Complex-valued derivative propagation method with approximate Bohmian trajectories: Application to electronic nonadiabatic dynamics. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.03.051] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
8
|
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
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
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: 80] [Impact Index Per Article: 10.0] [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
|
13
|
Cruz-Rodríguez L, Tremblay JC, Martínez-Mesa A, Uranga-Piña L. A Chebyshev expansion of hydrodynamical fields for ultrafast vibrational wave packet dynamics. COMPUT THEOR CHEM 2016. [DOI: 10.1016/j.comptc.2016.01.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
14
|
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: 89] [Impact Index Per Article: 9.9] [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
|
15
|
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: 78] [Impact Index Per Article: 8.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
|
16
|
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
|
17
|
|
18
|
|
19
|
Banerjee T, Makri N. Quantum-Classical Path Integral with Self-Consistent Solvent-Driven Reference Propagators. J Phys Chem B 2013; 117:13357-66. [DOI: 10.1021/jp4043123] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tuseeta Banerjee
- Department
of Chemistry, University of Illinois, 600 S. Goodwin Avenue, Urbana, Illinois 61801, United States
| | - Nancy Makri
- Department
of Chemistry, University of Illinois, 600 S. Goodwin Avenue, Urbana, Illinois 61801, United States
- Department
of Physics, University of Illinois, 1110 W. Green Street, Urbana, Illinois 61801, United States
| |
Collapse
|
20
|
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
|
21
|
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]
|
22
|
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
|
23
|
Lambert R, Makri N. Quantum-classical path integral. I. Classical memory and weak quantum nonlocality. J Chem Phys 2012; 137:22A552. [DOI: 10.1063/1.4767931] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
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
|
25
|
PARK KISAM, POIRIER BILL. QUANTUM TRAJECTORY CALCULATIONS FOR BIPOLAR WAVEPACKET DYNAMICS IN ONE DIMENSION: SYNTHETIC SINGLE-WAVEPACKET PROPAGATION. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2011. [DOI: 10.1142/s0219633610005967] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In a previous paper [Park K, Poirier B, Parlant G, J Chem Phys129:194112, 2008], a synthetic quantum trajectory method (QTM) was successfully implemented for wave-packet dynamics in a one-dimensional (1D) symmetric Eckart barrier system, utilizing a "double-wavepacket" version of the bipolar decomposition, ψ = ψ+ + ψ- = (ψ1+ + ψ2+) + (ψ1- + ψ2-), to avoid a technical difficulty involving negligible initial ψ- density. In this paper, we develop a new synthetic algorithm which overcomes this difficulty directly, utilizing the original "single-wavepacket" version of the bipolar decomposition, ψ =ψ+ + ψ-, and also show that the initial propagation of ψ- is mainly governed by probability transfer from ψ+, rather than by the given initial conditions for ψ-. The new algorithm makes it possible to apply the synthetic bipolar QTM to asymptotically asymmetric as well as symmetric potential systems. Successful application results for both symmetric and asymmetric Eckart barrier systems in 1D are presented.
Collapse
Affiliation(s)
- KISAM PARK
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409-1061, USA
| | - BILL POIRIER
- Department of Physics, Texas Tech University, Lubbock, Texas 79409-1061, USA
| |
Collapse
|
26
|
Lasorne B, Worth GA, Robb MA. Excited-state dynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.26] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
27
|
Curchod BFE, Tavernelli I, Rothlisberger U. Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly approach for molecular dynamics simulations. Phys Chem Chem Phys 2011; 13:3231-6. [DOI: 10.1039/c0cp02175j] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
28
|
Makri N. Forward–backward semiclassical and quantum trajectory methods for time correlation functions. Phys Chem Chem Phys 2011; 13:14442-52. [DOI: 10.1039/c0cp02374d] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
29
|
Garashchuk S, Rassolov V, Prezhdo O. Semiclassical Bohmian Dynamics. REVIEWS IN COMPUTATIONAL CHEMISTRY 2010. [DOI: 10.1002/9780470890905.ch6] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
30
|
Abstract
We present here a supersymmetric (SUSY) approach for determining excitation energies within the context of a quantum Monte Carlo scheme. By using the fact that SUSY quantum mechanics gives rises to a series of isospectral Hamiltonians, we show that Monte Carlo ground-state calculations in the SUSY partners can be used to reconstruct accurately both the spectrum and states of an arbitrary Schrodinger equation. Since the ground state of each partner potential is nodeless, we avoid any "node" problem typically associated with the Monte Carlo technique. Although we provide an example of using this approach to determine the tunneling states in a double-well potential, the method is applicable to any 1D potential problem. We conclude by discussing the extension to higher dimensions.
Collapse
Affiliation(s)
- Eric R Bittner
- Department of Chemistry, University of Houston, Houston, Texas 77204, USA.
| | | | | |
Collapse
|
31
|
Garashchuk S. Description of bound reactive dynamics within the approximate quantum trajectory framework. J Phys Chem A 2009; 113:4451-6. [PMID: 19290585 DOI: 10.1021/jp8110869] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The quantum trajectory framework incorporates quantum effects on dynamics through the quantum potential acting on a trajectory ensemble in addition to the classical potential. A global quadratic approximation to the quantum potential makes the method practical in many dimensions and captures dominant quantum effects in semiclassical systems. In this paper the approach is further developed to describe the "double well" dynamicsa prototype of the proton transfer reactionswhich exhibits the "hard" quantum effect of tunneling. Accurate description is achieved by combining the approximate quantum trajectory dynamics with the population amplitudes in the reactant and product wells. The quantum trajectory dynamics is defined by the asymptotic classical potentials. The population amplitudes represented in a small basis describe transfer between the wells. The method is exact if the reactant/product potentials are quadratic and the basis size is sufficiently large. In the semiclassical regime the trajectory dynamics is approximate and the basis size can be as small as two functions. The approach is fully compatible with the trajectory description of multidimensional systems capturing quantum tunneling along the reactive coordinate and zero-point energy flow among all degrees of freedom.
Collapse
Affiliation(s)
- Sophya Garashchuk
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA
| |
Collapse
|
32
|
Park K, Poirier B, Parlant G. Quantum trajectory calculations for bipolar wavepacket dynamics in one dimension. J Chem Phys 2008; 129:194112. [DOI: 10.1063/1.3013630] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Poirier B. Reconciling semiclassical and Bohmian mechanics. V. Wavepacket dynamics. J Chem Phys 2008; 128:164115. [DOI: 10.1063/1.2850207] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
|
34
|
Garashchuk S, Rassolov VA. Semiclassical nonadiabatic dynamics of NaFH with quantum trajectories. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2007.08.061] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
35
|
Poirier B, Parlant G. Reconciling semiclassical and Bohmian mechanics: IV. Multisurface dynamics. J Phys Chem A 2007; 111:10400-8. [PMID: 17887734 DOI: 10.1021/jp0731349] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
In previous articles (J. Chem. Phys. 2004, 121, 4501; 2006, 124, 034115; 2006, 124, 034116) a bipolar counter-propagating wave decomposition, Psi = Psi+ + Psi-, was presented for stationary states Psi of the one-dimensional Schrödinger equation, such that the components Psi+/- approach their semiclassical WKB analogs in the large action limit. The corresponding bipolar quantum trajectories are classical-like and well-behaved, even when Psi has many nodes or is wildly oscillatory. In this paper, the method is generalized for multisurface scattering applications and applied to several benchmark problems. A natural connection is established between intersurface transitions and (+ <--> -) transitions.
Collapse
Affiliation(s)
- Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, Box 41061, Lubbock, Texas 79409-1061, USA.
| | | |
Collapse
|
36
|
Gonzalez MF, Giménez X, Gonzalez-Aguilar J, Bofill JM. Quantum Trajectories from a Discrete−Variable Representation Method. J Phys Chem A 2007; 111:10226-33. [PMID: 17696411 DOI: 10.1021/jp072237o] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
A method for obtaining quantum trajectories from a discrete-variable representation computation of the quantum potential is presented. The method exploits the linearity of the Schrödinger equation, deals smoothly with the quantum potential singularities, and readily performs the time propagation up to fairly large total elapsed times. A one-dimensional test of the general n-dimensional formulation is included.
Collapse
Affiliation(s)
- María Fernanda Gonzalez
- Departament de Química Física i Centre especial de Recerca en Química Teorica, Universitat de Barcelona i Parc Científic de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain
| | | | | | | |
Collapse
|
37
|
Wu Y, Herman MF. On the properties of a primitive semiclassical surface hopping propagator for nonadiabatic quantum dynamics. J Chem Phys 2007; 127:044109. [PMID: 17672683 DOI: 10.1063/1.2756532] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
A previously developed nonadiabatic semiclassical surface hopping propagator [M. F. Herman J. Chem. Phys. 103, 8081 (1995)] is further studied. The propagator has been shown to satisfy the time-dependent Schrodinger equation (TDSE) through order h, and the O(h2) terms are treated as small errors, consistent with standard semiclassical analysis. Energy is conserved at each hopping point and the change in momentum accompanying each hop is parallel to the direction of the nonadiabatic coupling vector resulting in both transmission and reflection types of hops. Quantum mechanical analysis and numerical calculations presented in this paper show that the h2 terms involving the interstate coupling functions have significant effects on the quantum transition probabilities. Motivated by these data, the h2 terms are analyzed for the nonadiabatic semiclassical propagator. It is shown that the propagator can satisfy the TDSE for multidimensional systems by including another type of nonclassical trajectories that reflect on the same surfaces. This h2 analysis gives three conditions for these three types of trajectories so that their coefficients are uniquely determined. Besides the nonadiabatic semiclassical propagator, a numerically useful quantum propagator in the adiabatic representation is developed to describe nonadiabatic transitions.
Collapse
Affiliation(s)
- Yinghua Wu
- Department of Chemistry, Tulane University, New Orleans, Louisiana 70118, USA
| | | |
Collapse
|
38
|
|
39
|
Weeks DE, Niday TA, Yang SH. Inelastic scattering matrix elements for the nonadiabatic collision B(P1∕22)+H2(Σg+1,j)↔B(P3∕22)+H2(Σg+1,j′). J Chem Phys 2006; 125:164301. [PMID: 17092067 DOI: 10.1063/1.2222369] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Inelastic scattering matrix elements for the nonadiabatic collision B(2P1/2)+H2(1Sigmag+,j)<-->B(2P3/2)+H2(1Sigmag+,j') are calculated using the time dependent channel packet method (CPM). The calculation employs 1 2A', 2 2A', and 1 2A" adiabatic electronic potential energy surfaces determined by numerical computation at the multireference configuration-interaction level [M. H. Alexander, J. Chem. Phys. 99, 6041 (1993)]. The 1 2A' and 2 2A', adiabatic electronic potential energy surfaces are transformed to yield diabatic electronic potential energy surfaces that, when combined with the total B+H2 rotational kinetic energy, yield a set of effective potential energy surfaces [M. H. Alexander et al., J. Chem. Phys. 103, 7956 (1995)]. Within the framework of the CPM, the number of effective potential energy surfaces used for the scattering matrix calculation is then determined by the size of the angular momentum basis used as a representation. Twenty basis vectors are employed for these calculations, and the corresponding effective potential energy surfaces are identified in the asymptotic limit by the H2 rotor quantum numbers j=0, 2, 4, 6 and B electronic states 2Pja, ja=1/2, 3/2. Scattering matrix elements are obtained from the Fourier transform of the correlation function between channel packets evolving in time on these effective potential energy surfaces. For these calculations the H2 bond length is constrained to a constant value of req=1.402 a.u. and state to state scattering matrix elements corresponding to a total angular momentum of J=1/2 are discussed for j=0<-->j'=0,2,4 and 2P1/2<-->2P1/2, 2P3/2 over a range of total energy between 0.0 and 0.01 a.u.
Collapse
Affiliation(s)
- David E Weeks
- Department of Engineering Physics, Air Force Institute of Technology, Wright-Patterson AFB, Ohio 45433, USA.
| | | | | |
Collapse
|
40
|
Trahan C, Poirier B. Reconciling semiclassical and Bohmian mechanics. II. Scattering states for discontinuous potentials. J Chem Phys 2006; 124:034115. [PMID: 16438575 DOI: 10.1063/1.2145883] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In a previous paper [B. Poirier, J. Chem. Phys. 121, 4501 (2004)] a unique bipolar decomposition, psi = psi1 + psi2, was presented for stationary bound states Psi of the one-dimensional Schrodinger equation, such that the components psi1 and psi2 approach their semiclassical WKB analogs in the large action limit. Moreover, by applying the Madelung-Bohm ansatz to the components rather than to Psi itself, the resultant bipolar Bohmian mechanical formulation satisfies the correspondence principle. As a result, the bipolar quantum trajectories are classical-like and well behaved, even when psi has many nodes or is wildly oscillatory. In this paper, the previous decomposition scheme is modified in order to achieve the same desirable properties for stationary scattering states. Discontinuous potential systems are considered (hard wall, step potential, and square barrier/well), for which the bipolar quantum potential is found to be zero everywhere, except at the discontinuities. This approach leads to an exact numerical method for computing stationary scattering states of any desired boundary conditions, and reflection and transmission probabilities. The continuous potential case will be considered in a companion paper [C. Trahan and B. Poirier, J. Chem. Phys. 124, 034116 (2006), following paper].
Collapse
Affiliation(s)
- Corey Trahan
- Department of Chemistry and Biochemistry, Texas Tech University, P.O. Box 41061, Lubbock, Texas 79409-1061, USA
| | | |
Collapse
|
41
|
|
42
|
Gindensperger E, Meier C, Beswick JA, Parlant G. Combining fixed- and moving-grid methods to study direct dissociation processes involving nonadiabatic transitions. J Chem Phys 2005; 123:214107. [PMID: 16356039 DOI: 10.1063/1.2114807] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
We present a novel quantum-dynamics approach suitable for computing direct dissociation processes, including electronic transitions. This approach combines quantum trajectories in the Lagrangian reference frame with standard fixed-grid wave packets in order to overcome the limitations and difficulties of both techniques. As a model application, we consider the ultrafast photodissociation of H2 excited by a femtosecond extreme UV laser pulse.
Collapse
Affiliation(s)
- Etienne Gindensperger
- Laboratoire Collisions Agrégats Réactivité, Institut de Recherche sur les Systémes Atomiques et Moléculaires Complexes, Centre National de la Recherche Scientifique -Université Paul Sabatier, Toulouse, France.
| | | | | | | |
Collapse
|
43
|
Trahan CJ, Wyatt RE, Poirier B. Multidimensional quantum trajectories: Applications of the derivative propagation method. J Chem Phys 2005; 122:164104. [PMID: 15945669 DOI: 10.1063/1.1884606] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In a previous publication [J. Chem. Phys. 118, 9911 (2003)], the derivative propagation method (DPM) was introduced as a novel numerical scheme for solving the quantum hydrodynamic equations of motion (QHEM) and computing the time evolution of quantum mechanical wave packets. These equations are a set of coupled, nonlinear partial differential equations governing the time evolution of the real-valued functions C and S in the complex action, S=C(r,t) + iS(r,t)/Planck's over 2pi, where Psi(r,t)=exp(S). Past numerical solutions to the QHEM were obtained via ensemble trajectory propagation, where the required first- and second-order spatial derivatives were evaluated using fitting techniques such as moving least squares. In the DPM, however, equations of motion are developed for the derivatives themselves, and a truncated set of these are integrated along quantum trajectories concurrently with the original QHEM equations for C and S. Using the DPM quantum effects can be included at various orders of approximation; no spatial fitting is involved; there is no basis set expansion; and single, uncoupled quantum trajectories can be propagated (in parallel) rather than in correlated ensembles. In this study, the DPM is extended from previous one-dimensional (1D) results to calculate transmission probabilities for 2D and 3D wave packet evolution on coupled Eckart barrier/harmonic oscillator surfaces. In the 2D problem, the DPM results are compared to standard numerical integration of the time-dependent Schrodinger equation. Also in this study, the practicality of implementing the DPM for systems with many more degrees of freedom is discussed.
Collapse
Affiliation(s)
- Corey J Trahan
- Department of Chemistry and Biochemistry, Texas Tech University, Box 41061 Lubbock, Texas 79409-1061, USA.
| | | | | |
Collapse
|
44
|
Bonacić-Koutecký V, Mitrić R. Theoretical Exploration of Ultrafast Dynamics in Atomic Clusters: Analysis and Control. Chem Rev 2004; 105:11-66. [PMID: 15720151 DOI: 10.1021/cr0206925] [Citation(s) in RCA: 72] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Vlasta Bonacić-Koutecký
- Institut für Chemie, Humboldt-Universität zu Berlin, Brook-Taylor-Strasse 2, D-12489 Berlin, Germany.
| | | |
Collapse
|
45
|
Garashchuk S, Rassolov VA. Modified quantum trajectory dynamics using a mixed wave function representation. J Chem Phys 2004; 121:8711-5. [PMID: 15527334 DOI: 10.1063/1.1804177] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Dynamics of quantum trajectories provides an efficient framework for description of various quantum effects in large systems, but it is unstable near the wave function density nodes where the quantum potential becomes singular. A mixed coordinate space/polar representation of the wave function is used to circumvent this problem. The resulting modified trajectory dynamics associated with the polar representation is nonsingular and smooth. The interference structure and the nodes of the wave function density are described, in principle, exactly in the coordinate representation. The approximate version of this approach is consistent with the semiclassical linearized quantum force method [S. Garashchuk and V. A. Rassolov, J. Chem. Phys. 120, 1181 (2004)]. This approach is exact for general wave functions with the density nodes in a locally quadratic potential.
Collapse
Affiliation(s)
- Sophya Garashchuk
- Department of Chemistry and Biochemistry, University of South Carolina, South Carolina 29208, USA
| | | |
Collapse
|
46
|
Abstract
The semiclassical method is characterized by finite forces and smooth, well-behaved trajectories, but also by multivalued representational functions that are ill behaved at caustics. In contrast, quantum trajectory methods--based on Bohmian mechanics (quantum hydrodynamics)--are characterized by divergent forces and erratic trajectories near nodes, but also well-behaved, single-valued representational functions. In this paper, we unify these two approaches into a single method that captures the best features of both, and in addition, satisfies the correspondence principle. Stationary eigenstates in one degree of freedom are the primary focus, but more general applications are also anticipated.
Collapse
Affiliation(s)
- Bill Poirier
- Department of Chemistry and Biochemistry, and Department of Physics, Texas Tech University, Lubbock, Texas 79409-1061, USA.
| |
Collapse
|
47
|
Gonzalez J, Bofill JM, Gimenez X. A Bohmian total potential view to quantum effects. I. Methodology and simple model systems. J Chem Phys 2004; 120:10961-71. [PMID: 15268126 DOI: 10.1063/1.1747869] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The coherent-state wave packet dynamics of several model systems is analyzed in terms of Bohm's total potential. The quantum dynamics has been obtained by solving the time-dependent Schrodinger equation, and a method for obtaining the total potential from it, involving just matrix algebra, has been proposed. Contrary to what one may expect, it is shown that the time- and state-dependent features of the total potential admit a rationale, classical-like description of quantum effects, leading to a unified picture of them, which is not critically dependent, as for the key features, on the classical potential. An outstanding feature is found to be the relation of the state system's density amplitude and sharpness (in its dependence with position) with quantum effects. Sharp density profiles and low densities cause the total potential to strongly depart from the classical value, in both time regimes and position ranges, which provide a clearer, more deterministic view to quantum dynamics. Free motion as well as scattering processes by square and Eckart barriers have been analyzed by means of careful inspection of several time dependent snapshots. The result is an insightful picture of processes involving tunneling and antitunneling, including their dynamical variants, as well as resonances and quantization.
Collapse
Affiliation(s)
- Javier Gonzalez
- Departament de Quimica Fisica, Universitat de Barcelona, Marti i Franques 1Centre Especial de Recerca en Quimica Teorica, Parc Cientific de Barcelona, Baldiri Reixac 10-12, 08028 Barcelona, Spain
| | | | | |
Collapse
|
48
|
Rassolov VA, Garashchuk S. Bohmian dynamics on subspaces using linearized quantum force. J Chem Phys 2004; 120:6815-25. [PMID: 15267580 DOI: 10.1063/1.1669385] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In the de Broglie-Bohm formulation of quantum mechanics the time-dependent Schrodinger equation is solved in terms of quantum trajectories evolving under the influence of quantum and classical potentials. For a practical implementation that scales favorably with system size and is accurate for semiclassical systems, we use approximate quantum potentials. Recently, we have shown that optimization of the nonclassical component of the momentum operator in terms of fitting functions leads to the energy-conserving approximate quantum potential. In particular, linear fitting functions give the exact time evolution of a Gaussian wave packet in a locally quadratic potential and can describe the dominant quantum-mechanical effects in the semiclassical scattering problems of nuclear dynamics. In this paper we formulate the Bohmian dynamics on subspaces and define the energy-conserving approximate quantum potential in terms of optimized nonclassical momentum, extended to include the domain boundary functions. This generalization allows a better description of the non-Gaussian wave packets and general potentials in terms of simple fitting functions. The optimization is performed independently for each domain and each dimension. For linear fitting functions optimal parameters are expressed in terms of the first and second moments of the trajectory distribution. Examples are given for one-dimensional anharmonic systems and for the collinear hydrogen exchange reaction.
Collapse
Affiliation(s)
- Vitaly A Rassolov
- Department of Chemistry & Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA
| | | |
Collapse
|
49
|
Garashchuk S, Rassolov VA. Energy conserving approximations to the quantum potential: Dynamics with linearized quantum force. J Chem Phys 2004; 120:1181-90. [PMID: 15268241 DOI: 10.1063/1.1633263] [Citation(s) in RCA: 71] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Solution of the Schrodinger equation within the de Broglie-Bohm formulation is based on propagation of trajectories in the presence of a nonlocal quantum potential. We present a new strategy for defining approximate quantum potentials within a restricted trial function by performing the optimal fit to the log-derivatives of the wave function density. This procedure results in the energy-conserving dynamics for a closed system. For one particular form of the trial function leading to the linear quantum force, the optimization problem is solved analytically in terms of the first and second moments of the weighted trajectory distribution. This approach gives exact time-evolution of a correlated Gaussian wave function in a locally quadratic potential. The method is computationally cheap in many dimensions, conserves total energy and satisfies the criterion on the average quantum force. Expectation values are readily found by summing over trajectory weights. Efficient extraction of the phase-dependent quantities is discussed. We illustrate the efficiency and accuracy of the linear quantum force approximation by examining a one-dimensional scattering problem and by computing the wavepacket reaction probability for the hydrogen exchange reaction and the photodissociation spectrum of ICN in two dimensions.
Collapse
Affiliation(s)
- Sophya Garashchuk
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA.
| | | |
Collapse
|
50
|
Pauler DK, Kendrick BK. A new method for solving the quantum hydrodynamic equations of motion: Application to two-dimensional reactive scattering. J Chem Phys 2004; 120:603-11. [PMID: 15267894 DOI: 10.1063/1.1630302] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The de Broglie-Bohm hydrodynamic equations of motion are solved using a meshless method based on a moving least squares approach and an arbitrary Lagrangian-Eulerian frame of reference. A regridding algorithm adds and deletes computational points as needed in order to maintain a uniform interparticle spacing, and unitary time evolution is obtained by propagating the wave packet using averaged fields. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. The methodology is applied to a two-dimensional model collinear reaction with an activation barrier. Reaction probabilities are computed as a function of both time and energy, and are in excellent agreement with those based on the quantum trajectory method.
Collapse
Affiliation(s)
- Denise K Pauler
- Theoretical Division (T-12, MS-B268), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | |
Collapse
|