1
|
Kessing RK, Yang PY, Manmana SR, Cao J. Long-Range Nonequilibrium Coherent Tunneling Induced by Fractional Vibronic Resonances. J Phys Chem Lett 2022; 13:6831-6838. [PMID: 35857895 DOI: 10.1021/acs.jpclett.2c01455] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
We study the influence of a linear energy bias on a nonequilibrium excitation on a chain of molecules coupled to local vibrations (a tilted Holstein model) using both a random-walk rate kernel theory and a nonperturbative, massively parallelized adaptive-basis algorithm. We uncover structured and discrete vibronic resonance behavior fundamentally different from both linear response theory and homogeneous polaron dynamics. Remarkably, resonance between the phonon energy ℏω and the bias δϵ occurs not only at integer but also fractional ratios δϵ/(ℏω) = m/n, which effect long-range n-bond m-phonon tunneling. These observations are reproduced in a model calculation of a recently demonstrated Cy3 system, and the effect of dipole-dipole-type non-nearest-neighbor coupling and vibrationally relaxed initial states is also considered. Potential applications range from molecular electronics to optical lattices and artificial light harvesting via vibronic engineering of coherent quantum transport.
Collapse
Affiliation(s)
- R Kevin Kessing
- Institut für Theoretische Physik, Universität Ulm, Ulm, 89069, Germany
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Göttingen, 37077, Germany
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Pei-Yun Yang
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Physics, National Taiwan University, Taipei, 10617, Taiwan (R.O.C.)
- Beijing Computational Science Research Center, Beijing, 100193, China
| | - Salvatore R Manmana
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Göttingen, 37077, Germany
- Fachbereich Physik, Philipps-Universität Marburg, Marburg, 35032, Germany
| | - Jianshu Cao
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
2
|
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
|
3
|
Richings GW, Habershon S. Predicting Molecular Photochemistry Using Machine-Learning-Enhanced Quantum Dynamics Simulations. Acc Chem Res 2022; 55:209-220. [PMID: 34982533 DOI: 10.1021/acs.accounts.1c00665] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The processes which occur after molecules absorb light underpin an enormous range of fundamental technologies and applications, including photocatalysis to enable new chemical transformations, sunscreens to protect against the harmful effects of UV overexposure, efficient photovoltaics for energy generation from sunlight, and fluorescent probes to image the intricate details of complex biomolecular structures. Reflecting this broad range of applications, an enormously versatile set of experiments are now regularly used to interrogate light-driven chemical dynamics, ranging from the typical ultrafast transient absorption spectroscopy used in many university laboratories to the inspiring central facilities around the world, such as the next-generation of X-ray free-electron lasers.Computer simulations of light-driven molecular and material dynamics are an essential route to analyzing the enormous amount of transient electronic and structural data produced by these experimental sources. However, to date, the direct simulation of molecular photochemistry remains a frontier challenge in computational chemical science, simultaneously demanding the accurate treatment of molecular electronic structure, nuclear dynamics, and the impact of nonadiabatic couplings.To address these important challenges and to enable new computational methods which can be integrated with state-of-the-art experimental capabilities, the past few years have seen a burst of activity in the development of "direct" quantum dynamics methods, merging the machine learning of potential energy surfaces (PESs) and nonadiabatic couplings with accurate quantum propagation schemes such as the multiconfiguration time-dependent Hartree (MCTDH) method. The result of this approach is a new generation of direct quantum dynamics tools in which PESs are generated in tandem with wave function propagation, enabling accurate "on-the-fly" simulations of molecular photochemistry. These simulations offer an alternative route toward gaining quantum dynamics insights, circumventing the challenge of generating ab initio electronic structure data for PES fitting by instead only demanding expensive energy evaluations as and when they are needed.In this Account, we describe the chronological evolution of our own contributions to this field, focusing on describing the algorithmic developments that enable direct MCTDH simulations for complex molecular systems moving on multiple coupled electronic states. Specifically, we highlight active learning strategies for generating PESs during grid-based quantum chemical dynamics simulations, and we discuss the development and impact of novel diabatization schemes to enable direct grid-based simulations of photochemical dynamics; these developments are highlighted in a series of benchmark molecular simulations of systems containing multiple nuclear degrees of freedom moving on multiple coupled electronic states. We hope that the ongoing developments reported here represent a major step forward in tools for modeling excited-state chemistry such as photodissociation, proton and electron transfer, and ultrafast energy dissipation in complex molecular systems.
Collapse
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry, University of Warwick, Coventry, United Kingdom CV4 7AL
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry, United Kingdom CV4 7AL
| |
Collapse
|
4
|
Lee TY, Lu CY, Chou CC. Moving Boundary Truncated Grid Method: Application to the Time Evolution of Distribution Functions in Phase Space. J Phys Chem A 2021; 125:476-491. [PMID: 33372780 DOI: 10.1021/acs.jpca.0c09525] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The moving boundary truncated grid (TG) method, previously developed to integrate the time-dependent Schrödinger equation and the imaginary time Schrödinger equation, is extended to the time evolution of distribution functions in phase space. A variable number of phase space grid points in the Eulerian representation are used to integrate the equation of motion for the distribution function, and the boundaries of the TG are adaptively determined as the distribution function evolves in time. Appropriate grid points are activated and deactivated for propagation of the distribution function, and no advance information concerning the dynamics in phase space is required. The TG method is used to integrate the equations of motion for phase space distribution functions, including the Klein-Kramers, Wigner-Moyal, and modified Caldeira-Leggett equations. Even though the initial distribution function is nonnegative, the solutions to the Wigner-Moyal and modified Caldeira-Leggett equations may develop negative basins in phase space originating from interference effects. Trajectory-based methods for propagation of the distribution function do not permit the formation of negative regions. However, the TG method can correctly capture the negative basins. Comparisons between the computational results obtained from the full grid and TG calculations demonstrate that the TG method not only significantly reduces the computational effort but also permits accurate propagation of various distribution functions in phase space.
Collapse
Affiliation(s)
- Tsung-Yen Lee
- Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| | - Chun-Yaung Lu
- Texas Advanced Computing Center, The University of Texas at Austin, Austin, Texas 78758, United States
| | - Chia-Chun Chou
- Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| |
Collapse
|
5
|
Pandey A, Poirier B. An algorithm to find (and plug) “holes” in multi-dimensional surfaces. J Chem Phys 2020; 152:214102. [DOI: 10.1063/5.0005681] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Ankit Pandey
- Department of Chemistry and Biochemistry, Texas Tech University, P.O. Box 41061, Lubbock, Texas 79409-1061, USA
| | - Bill Poirier
- Department of Chemistry and Biochemistry, Texas Tech University, P.O. Box 41061, Lubbock, Texas 79409-1061, USA
| |
Collapse
|
6
|
Odunlami M, Le Bris V, Bégué D, Baraille I, Coulaud O. A-VCI: A flexible method to efficiently compute vibrational spectra. J Chem Phys 2018; 146:214108. [PMID: 28595393 DOI: 10.1063/1.4984266] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The adaptive vibrational configuration interaction algorithm has been introduced as a new method to efficiently reduce the dimension of the set of basis functions used in a vibrational configuration interaction process. It is based on the construction of nested bases for the discretization of the Hamiltonian operator according to a theoretical criterion that ensures the convergence of the method. In the present work, the Hamiltonian is written as a sum of products of operators. The purpose of this paper is to study the properties and outline the performance details of the main steps of the algorithm. New parameters have been incorporated to increase flexibility, and their influence has been thoroughly investigated. The robustness and reliability of the method are demonstrated for the computation of the vibrational spectrum up to 3000 cm-1 of a widely studied 6-atom molecule (acetonitrile). Our results are compared to the most accurate up to date computation; we also give a new reference calculation for future work on this system. The algorithm has also been applied to a more challenging 7-atom molecule (ethylene oxide). The computed spectrum up to 3200 cm-1 is the most accurate computation that exists today on such systems.
Collapse
Affiliation(s)
- Marc Odunlami
- Université de Pau et des Pays de l'Adour, CNRS, Institut des Sciences Analytiques et de Physico-Chimie pour l'Environnement et les Matériaux, UMR5254, 64000 Pau, France
| | - Vincent Le Bris
- Université de Pau et des Pays de l'Adour, CNRS, Institut des Sciences Analytiques et de Physico-Chimie pour l'Environnement et les Matériaux, UMR5254, 64000 Pau, France
| | - Didier Bégué
- Université de Pau et des Pays de l'Adour, CNRS, Institut des Sciences Analytiques et de Physico-Chimie pour l'Environnement et les Matériaux, UMR5254, 64000 Pau, France
| | - Isabelle Baraille
- Université de Pau et des Pays de l'Adour, CNRS, Institut des Sciences Analytiques et de Physico-Chimie pour l'Environnement et les Matériaux, UMR5254, 64000 Pau, France
| | - Olivier Coulaud
- HiePACS Project-Team, Inria Bordeaux Sud-Ouest, 200, Avenue de la Vieille Tour, 33405 Talence Cedex, France
| |
Collapse
|
7
|
Larsson HR, Riedel J, Wei J, Temps F, Hartke B. Resonance dynamics of DCO (X̃ A′2) simulated with the dynamically pruned discrete variable representation (DP-DVR). J Chem Phys 2018; 148:204309. [DOI: 10.1063/1.5026459] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Henrik R. Larsson
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Jens Riedel
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Jie Wei
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Friedrich Temps
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| | - Bernd Hartke
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, Olshausenstraße 40, 24098 Kiel, Germany
| |
Collapse
|
8
|
Lee TY, Chou CC. Moving Boundary Truncated Grid Method for Wave Packet Dynamics. J Phys Chem A 2018; 122:1451-1463. [DOI: 10.1021/acs.jpca.7b11932] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- 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
|
9
|
Avila G, Carrington T. Reducing the cost of using collocation to compute vibrational energy levels: Results for CH2NH. J Chem Phys 2017; 147:064103. [DOI: 10.1063/1.4994920] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gustavo Avila
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
10
|
Larsson HR, Tannor DJ. Dynamical pruning of the multiconfiguration time-dependent Hartree (DP-MCTDH) method: An efficient approach for multidimensional quantum dynamics. J Chem Phys 2017; 147:044103. [DOI: 10.1063/1.4993219] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- H. R. Larsson
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, 24098 Kiel, Germany
- Department of Chemical Physics, Weizmann Institute of Science, 76100 Rehovot, Israel
| | - D. J. Tannor
- Department of Chemical Physics, Weizmann Institute of Science, 76100 Rehovot, Israel
| |
Collapse
|
11
|
Carrington T. Perspective: Computing (ro-)vibrational spectra of molecules with more than four atoms. J Chem Phys 2017; 146:120902. [DOI: 10.1063/1.4979117] [Citation(s) in RCA: 67] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario
K7L 3N6, Canada
| |
Collapse
|
12
|
Larsson HR, Hartke B, Tannor DJ. Efficient molecular quantum dynamics in coordinate and phase space using pruned bases. J Chem Phys 2016; 145:204108. [DOI: 10.1063/1.4967432] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- H. R. Larsson
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, 24098 Kiel, Germany
| | - B. Hartke
- Institut für Physikalische Chemie, Christian-Albrechts-Universität zu Kiel, 24098 Kiel, Germany
| | - D. J. Tannor
- Department of Chemical Physics, Weizmann Institute of Science, 76100 Rehovot, Israel
| |
Collapse
|
13
|
Alborzpour JP, Tew DP, Habershon S. Efficient and accurate evaluation of potential energy matrix elements for quantum dynamics using Gaussian process regression. J Chem Phys 2016; 145:174112. [DOI: 10.1063/1.4964902] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jonathan P. Alborzpour
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - David P. Tew
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
14
|
Brown J, Carrington T. Using an expanding nondirect product harmonic basis with an iterative eigensolver to compute vibrational energy levels with as many as seven atoms. J Chem Phys 2016; 145:144104. [DOI: 10.1063/1.4963916] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- James Brown
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| | - Tucker Carrington
- Chemistry Department, Queen’s University, Kingston, Ontario K7L 3N6, Canada
| |
Collapse
|
15
|
Heaps CW, Mazziotti DA. Pseudospectral Gaussian quantum dynamics: Efficient sampling of potential energy surfaces. J Chem Phys 2016; 144:164108. [DOI: 10.1063/1.4946807] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
16
|
Phase Space Approach to Solving the Schrödinger Equation: Thinking Inside the Box. ADVANCES IN CHEMICAL PHYSICS 2014. [DOI: 10.1002/9781118949702.ch1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
17
|
Towards automated multi-dimensional quantum dynamical investigations of double-minimum potentials: Principles and example applications. Chem Phys 2011. [DOI: 10.1016/j.chemphys.2010.06.018] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
18
|
Thompson AL, Punwong C, Martínez TJ. Optimization of width parameters for quantum dynamics with frozen Gaussian basis sets. Chem Phys 2010. [DOI: 10.1016/j.chemphys.2010.03.020] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
19
|
Sielk J, von Horsten HF, Krüger F, Schneider R, Hartke B. Quantum-mechanical wavepacket propagation in a sparse, adaptive basis of interpolating Gaussians with collocation. Phys Chem Chem Phys 2009; 11:463-75. [DOI: 10.1039/b814315c] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
Pettey LR, Wyatt RE. Application of the Moving Boundary Truncation Method to Reactive Scattering: H + H2, O + H2, O + HD. J Phys Chem A 2008; 112:13335-42. [DOI: 10.1021/jp8067014] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Lucas R. Pettey
- Institute for Theoretical Chemistry and Department of Chemistry and Biochemistry, The University of Texas, Austin, Texas 78712
| | - Robert E. Wyatt
- Institute for Theoretical Chemistry and Department of Chemistry and Biochemistry, The University of Texas, Austin, Texas 78712
| |
Collapse
|
21
|
Ben-Nun M, Martínez TJ. A Continuous Spawning Method for Nonadiabatic Dynamics and Validation for the Zero-Temperature Spin-Boson Problem. Isr J Chem 2007. [DOI: 10.1560/ijc.47.1.75] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
22
|
von Horsten H, Hartke B. Links between potential energy structures and quantum cumulative reaction probabilities of double proton transfer reactions. Chem Phys 2007. [DOI: 10.1016/j.chemphys.2007.03.011] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
23
|
von Horsten HF, Rauhut G, Hartke B. Fingerprints of Delocalized Transition States in Quantum Dynamics. J Phys Chem A 2006; 110:13014-21. [PMID: 17134161 DOI: 10.1021/jp063051c] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Reactions with delocalized transition states (plateau reactions) can be characterized statically by their energy profile along the reaction path, where they exhibit a broad, flat region instead of one or several well-defined saddle points on the potential energy surface. Employing our new, highly flexible quantum dynamics code to perform two-dimensional and effective four-dimensional quantum wave packet propagations on ab initio based model potentials, we show that plateau reactions can also be discerned from the other standard reaction types by their dynamics.
Collapse
Affiliation(s)
- Hermann Frank von Horsten
- Institut für Physikalische Chemie, Christian-Albrechts-Universität, Olshausenstrasse 40, 24098 Kiel, Germany
| | | | | |
Collapse
|