1
|
Singh K, Lee KH, Peláez D, Bande A. Accelerating wavepacket propagation with machine learning. J Comput Chem 2024; 45:2360-2373. [PMID: 39031712 DOI: 10.1002/jcc.27443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 05/13/2024] [Accepted: 05/16/2024] [Indexed: 07/22/2024]
Abstract
In this work, we discuss the use of a recently introduced machine learning (ML) technique known as Fourier neural operators (FNO) as an efficient alternative to the traditional solution of the time-dependent Schrödinger equation (TDSE). FNOs are ML models which are employed in the approximated solution of partial differential equations. For a wavepacket propagating in an anharmonic potential and for a tunneling system, we show that the FNO approach can accurately and faithfully model wavepacket propagation via the density. Additionally, we demonstrate that FNOs can be a suitable replacement for traditional TDSE solvers in cases where the results of the quantum dynamical simulation are required repeatedly such as in the case of parameter optimization problems (e.g., control). The speed-up from the FNO method allows for its combination with the Markov-chain Monte Carlo approach in applications that involve solving inverse problems such as optimal and coherent laser control of the outcome of dynamical processes.
Collapse
Affiliation(s)
- Kanishka Singh
- Theory of Electron Dynamics and Spectroscopy, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Berlin, Germany
- Institute of Chemistry and Biochemistry, Freie Universität Berlin, Berlin, Germany
| | - Ka Hei Lee
- Theory of Electron Dynamics and Spectroscopy, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Berlin, Germany
- Fachbereich Physik, Freie Universität Berlin, Berlin, Germany
| | - Daniel Peláez
- CNRS, Institut des Sciences Moléculaires d'Orsay, Université Paris-Saclay, Orsay, France
| | - Annika Bande
- Theory of Electron Dynamics and Spectroscopy, Helmholtz-Zentrum Berlin für Materialien und Energie GmbH, Berlin, Germany
- Institute of Inorganic Chemistry, Leibniz University Hannover, Hannover, Germany
- Cluster of Excellence PhoenixD, Leibniz University Hannover, Hannover, Germany
| |
Collapse
|
2
|
Shi L, Schröder M, Meyer HD, Peláez D, Wodtke AM, Golibrzuch K, Schönemann AM, Kandratsenka A, Gatti F. Quantum and classical molecular dynamics for H atom scattering from graphene. J Chem Phys 2023; 159:194102. [PMID: 37965999 DOI: 10.1063/5.0176655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2023] [Accepted: 10/25/2023] [Indexed: 11/16/2023] Open
Abstract
This work presents systematic comparisons between classical molecular dynamics (cMD) and quantum dynamics (QD) simulations of 15-dimensional and 75-dimensional models in their description of H atom scattering from graphene. We use an experimentally validated full-dimensional neural network potential energy surface of a hydrogen atom interacting with a large cell of graphene containing 24 carbon atoms. For quantum dynamics simulations, we apply Monte Carlo canonical polyadic decomposition to transform the original potential energy surface (PES) into a sum of products form and use the multi-layer multi-configuration time-dependent Hartree method to simulate the quantum scattering of a hydrogen or deuterium atom with an initial kinetic energy of 1.96 or 0.96 eV and an incident angle of 0°, i.e., perpendicular to the graphene surface. The cMD and QD initial conditions have been carefully chosen in order to be as close as possible. Our results show little differences between cMD and QD simulations when the incident energy of the H atom is equal to 1.96 eV. However, a large difference in sticking probability is observed when the incident energy of the H atom is equal to 0.96 eV, indicating the predominance of quantum effects. To the best of our knowledge, our work provides the first benchmark of quantum against classical simulations for a system of this size with a realistic PES. Additionally, new projectors are implemented in the Heidelberg multi-configuration time-dependent Hartree package for the calculation of the atom scattering energy transfer distribution as a function of outgoing angles.
Collapse
Affiliation(s)
- Lei Shi
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay UMR 8214, 91405 Orsay, France
| | - Markus Schröder
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Hans-Dieter Meyer
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuenheimer Feld 229, 69120 Heidelberg, Germany
| | - Daniel Peláez
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay UMR 8214, 91405 Orsay, France
| | - Alec M Wodtke
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, am Faßberg 11, 37077 Göttingen, Germany
| | - Kai Golibrzuch
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, am Faßberg 11, 37077 Göttingen, Germany
| | - Anna-Maria Schönemann
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, am Faßberg 11, 37077 Göttingen, Germany
| | - Alexander Kandratsenka
- Department of Dynamics at Surfaces, Max Planck Institute for Multidisciplinary Sciences, am Faßberg 11, 37077 Göttingen, Germany
| | - Fabien Gatti
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay UMR 8214, 91405 Orsay, France
| |
Collapse
|
3
|
Nadoveza N, Panadés-Barrueta RL, Shi L, Gatti F, Peláez D. Analytical high-dimensional operators in canonical polyadic finite basis representation (CP-FBR). J Chem Phys 2023; 158:114109. [PMID: 36948837 DOI: 10.1063/5.0139224] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023] Open
Abstract
In the present work, we introduce a simple means of obtaining an analytical (i.e., grid-free) canonical polyadic (CP) representation of a multidimensional function that is expressed in terms of a set of discrete data. For this, we make use of an initial CP guess, even not fully converged, and a set of auxiliary basis functions [finite basis representation (FBR)]. The resulting CP-FBR expression constitutes the CP counterpart of our previous Tucker sum-of-products-FBR approach. However, as is well-known, CP expressions are much more compact. This has obvious advantages in high-dimensional quantum dynamics. The power of CP-FBR lies in the fact that it requires a grid much coarser than the one needed for the dynamics. In a subsequent step, the basis functions can be interpolated to any desired density of grid points. This is useful, for instance, when different initial conditions (e.g., energy content) of a system are to be considered. We show the application of the method to bound systems of increased dimensionality: H2 (3D), HONO (6D), and CH4 (9D).
Collapse
Affiliation(s)
- Nataša Nadoveza
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405 Orsay, France
| | | | - Lei Shi
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405 Orsay, France
| | - Fabien Gatti
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405 Orsay, France
| | - Daniel Peláez
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405 Orsay, France
| |
Collapse
|
4
|
Song Q, Zhang X, Peláez D, Meng Q. Direct Canonical-Polyadic-Decomposition of the Potential Energy Surface from Discrete Data by Decoupled Gaussian Process Regression. J Phys Chem Lett 2022; 13:11128-11135. [PMID: 36442084 DOI: 10.1021/acs.jpclett.2c03080] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
A Gaussian process regression (GPR) approach for directly constructing the canonical polyadic decomposition (CPD) of a multidimensional potential energy surface (PES) by discrete training energies is proposed and denoted by CPD-GPR. The present CPD-GPR method requires the kernel function in a product of a series of one-dimensional functions. To test CPD-GPR, the reactive probabilities of H + H2 as a function of kinetics energy are performed. Comparing the dynamics results computed by the CPD-GPR PES with those by the original PES, a good agreement between these results can be clearly found. Discussions on the previous algorithms for building the decomposed form are also given. We further show that the CPD-GPR method might be the general algorithm for building the decomposed form. However, further development is needed to reduce the CPD rank. Therefore, the present CPD-GPR method might be helpful to inspire ideas for developing new tools in building decomposed potential functions.
Collapse
Affiliation(s)
- Qingfei Song
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072Xi'an, China
- Institut des Sciences Moléculaires d'Orsay, CNRS-UMR 8214, Université Paris-Saclay, Bâtiment 520, F-91405Orsay, France
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072Xi'an, China
| | - Daniel Peláez
- Institut des Sciences Moléculaires d'Orsay, CNRS-UMR 8214, Université Paris-Saclay, Bâtiment 520, F-91405Orsay, France
| | - Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072Xi'an, China
| |
Collapse
|
5
|
Papp D, Tajti V, Avila G, Mátyus E, Czakó G. CH 4·F − revisited: full-dimensional ab initio potential energy surface and variational vibrational states. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2113565] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Dóra Papp
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary
| | - Viktor Tajti
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary
| | - Gustavo Avila
- ELTE, Eötvös Loránd University, Institute of Chemistry, Budapest, Hungary
| | - Edit Mátyus
- ELTE, Eötvös Loránd University, Institute of Chemistry, Budapest, Hungary
| | - Gábor Czakó
- MTA-SZTE Lendület Computational Reaction Dynamics Research Group, Interdisciplinary Excellence Centre and Department of Physical Chemistry and Materials Science, Institute of Chemistry, University of Szeged, Szeged, Hungary
| |
Collapse
|
6
|
Aerts A, Schaefer MR, Brown A. Adaptive Fitting of Potential Energy Surfaces of Small to Medium-Sized Molecules in Sum-of-Product Form: Application to Vibrational Spectroscopy. J Chem Phys 2022; 156:164106. [DOI: 10.1063/5.0089570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A semi-automatic sampling and fitting procedure for generating sum-of-product (Born-Oppenheimer) potential energy surfaces based on a high-dimensional model representation is presented. The adaptive sampling procedure and subsequent fitting relies on energies only and can be used for re-fitting existing analytic potential energy surfaces in sum-of-product form or for direct fits from ab initio computa- tions. The method is tested by fitting ground electronic state potential energy surfaces for small to medium sized semi-rigid molecules, i.e., HFCO, HONO, and HCOOH, based upon ab initio computations at the CCSD(T)-F12/cc-pVTZ-F12 or MP2/aug-cc-pVTZ levels of theory. Vibrational eigenstates are computed using block improved relaxation in the Heidelberg MCTDH package and compared to available experimental and theoretical data. The new potential energy surfaces are compared to the best ones currently available for these molecules, in terms of accuracy, including of resulting vibrational states, required numbers of sampling points, and number of fitting parameters. The present procedure leads to compact expansions and scales well with the number of dimensions for simple potentials such as single or double wells.
Collapse
Affiliation(s)
| | | | - Alex Brown
- Department of Chemistry, University of Alberta, Canada
| |
Collapse
|
7
|
Guo H, Worth G, Domcke W. Quantum dynamics with ab initio potentials. J Chem Phys 2021; 155:080401. [PMID: 34470339 DOI: 10.1063/5.0066234] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Graham Worth
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| |
Collapse
|