1
|
Høyer NM, Christiansen O. Quasi-direct Quantum Molecular Dynamics: The Time-Dependent Adaptive Density-Guided Approach for Potential Energy Surface Construction. J Chem Theory Comput 2024; 20:558-579. [PMID: 38183272 DOI: 10.1021/acs.jctc.3c00962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2024]
Abstract
We present a new quasi-direct quantum molecular dynamics computational method which offers a compromise between quantum dynamics using a precomputed potential energy surface (PES) and fully direct quantum dynamics. This method is termed the time-dependent adaptive density-guided approach (TD-ADGA) and is a method for constructing a PES on the fly during a dynamics simulation. This is achieved by acquisition of new single-point (SP) calculations and refitting of the PES, depending on the need of the dynamics. The TD-ADGA is a further development of the adaptive density-guided approach (ADGA) for PES construction where the placement of SPs is guided by the density of the nuclear wave function. In TD-ADGA, the ADGA framework has been integrated into the time propagation of the time-dependent nuclear wave function and we use the reduced one-mode density of this wave function to guide when and where new SPs are placed. The PES is thus extended or updated if the wave function moves into new areas or if a certain area becomes more important. Here, we derive equations for the reduced one-mode density for the time-dependent Hartree (TDH) method and for multiconfiguration time-dependent Hartree (MCTDH) methods, but the TD-ADGA can be used with any time-dependent wave function method as long as a density is available. The TD-ADGA method has been investigated on molecular systems containing single- and double-minimum potentials and on single-mode and multi-mode systems. We explore different approaches to handle the fact that the TD-ADGA involves a PES that changes during the computation and show how results can be obtained that are in very good agreement with results obtained by using an accurate reference PES. Dynamics with TD-ADGA is essentially a black box procedure, where only the initialization of the system and how to compute SPs must be provided. The TD-ADGA thus makes it easier to carry out quantum molecular dynamics and the quasi-direct framework opens up the possibility to compute quantum dynamics accurately for larger molecular systems.
Collapse
Affiliation(s)
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
2
|
Miyazaki K, Ananth N. Nonadiabatic simulations of photoisomerization and dissociation in ethylene using ab initio classical trajectories. J Chem Phys 2023; 159:124110. [PMID: 38127384 DOI: 10.1063/5.0163371] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Accepted: 09/05/2023] [Indexed: 12/23/2023] Open
Abstract
We simulate the nonadiabatic dynamics of photo-induced isomerization and dissociation in ethylene using ab initio classical trajectories in an extended phase space of nuclear and electronic variables. This is achieved by employing the linearized semiclassical initial value representation method for nonadiabatic dynamics, where discrete electronic states are mapped to continuous classical variables using either the Meyer-Miller-Stock-Thoss representation or a more recently introduced spin mapping approach. Trajectory initial conditions are sampled by constraining electronic state variables to a single initial excited state and by drawing nuclear phase space configurations from a Wigner distribution at a finite temperature. An ensemble of classical ab initio trajectories is then generated to compute thermal population correlation functions and analyze the mechanisms of isomerization and dissociation. Our results serve as a demonstration that this parameter-free semiclassical approach is computationally efficient and accurate, identifying mechanistic pathways in agreement with previous theoretical studies and also uncovering dissociation pathways observed experimentally.
Collapse
Affiliation(s)
- K Miyazaki
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - N Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
3
|
Wang TY, Neville SP, Schuurman MS. Machine Learning Seams of Conical Intersection: A Characteristic Polynomial Approach. J Phys Chem Lett 2023; 14:7780-7786. [PMID: 37615964 PMCID: PMC10494228 DOI: 10.1021/acs.jpclett.3c01649] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Accepted: 08/15/2023] [Indexed: 08/25/2023]
Abstract
The machine learning of potential energy surfaces (PESs) has undergone rapid progress in recent years. The vast majority of this work, however, has been focused on the learning of ground state PESs. To reliably extend machine learning protocols to excited state PESs, the occurrence of seams of conical intersections between adiabatic electronic states must be correctly accounted for. This introduces a serious problem, for at such points, the adiabatic potentials are not differentiable to any order, complicating the application of standard machine learning methods. We show that this issue may be overcome by instead learning the coordinate-dependent coefficients of the characteristic polynomial of a simple decomposition of the potential matrix. We demonstrate that, through this approach, quantitatively accurate machine learning models of seams of conical intersection may be constructed.
Collapse
Affiliation(s)
- Tzu Yu Wang
- Department
of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Simon P. Neville
- National
Research Council Canada, 100 Sussex Dr., Ottawa, Ontario K1A 0R6, Canada
| | - Michael S. Schuurman
- Department
of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
- National
Research Council Canada, 100 Sussex Dr., Ottawa, Ontario K1A 0R6, Canada
| |
Collapse
|
4
|
Artiukhin DG, Godtliebsen IH, Schmitz G, Christiansen O. Gaussian process regression adaptive density-guided approach: Toward calculations of potential energy surfaces for larger molecules. J Chem Phys 2023; 159:024102. [PMID: 37428042 DOI: 10.1063/5.0152367] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Accepted: 06/15/2023] [Indexed: 07/11/2023] Open
Abstract
We present a new program implementation of the Gaussian process regression adaptive density-guided approach [Schmitz et al., J. Chem. Phys. 153, 064105 (2020)] for automatic and cost-efficient potential energy surface construction in the MidasCpp program. A number of technical and methodological improvements made allowed us to extend this approach toward calculations of larger molecular systems than those previously accessible and maintain the very high accuracy of constructed potential energy surfaces. On the methodological side, improvements were made by using a Δ-learning approach, predicting the difference against a fully harmonic potential, and employing a computationally more efficient hyperparameter optimization procedure. We demonstrate the performance of this method on a test set of molecules of growing size and show that up to 80% of single point calculations could be avoided, introducing a root mean square deviation in fundamental excitations of about 3 cm-1. A much higher accuracy with errors below 1 cm-1 could be achieved with tighter convergence thresholds still reducing the number of single point computations by up to 68%. We further support our findings with a detailed analysis of wall times measured while employing different electronic structure methods. Our results demonstrate that GPR-ADGA is an effective tool, which could be applied for cost-efficient calculations of potential energy surfaces suitable for highly accurate vibrational spectra simulations.
Collapse
Affiliation(s)
- Denis G Artiukhin
- Institut für Chemie und Biochemie, Freie Universität Berlin, Arnimallee 22, 14195 Berlin, Germany
| | - Ian H Godtliebsen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | - Gunnar Schmitz
- Lehrstuhl für Theoretische Chemie II, Ruhr-Universität Bochum, Universitätstraße 150, 44801 Bochum, Germany
| | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
5
|
Schneider M, Born D, Kästner J, Rauhut G. Positioning of grid points for spanning potential energy surfaces-How much effort is really needed? J Chem Phys 2023; 158:144118. [PMID: 37061506 DOI: 10.1063/5.0146020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
The positions of grid points for representing a multidimensional potential energy surface (PES) have a non-negligible impact on its accuracy and the associated computational effort for its generation. Six different positioning schemes were studied for PESs represented by n-mode expansions as needed for the accurate calculation of anharmonic vibrational frequencies by means of vibrational configuration interaction theory. A static approach, which has successfully been used in many applications, and five adaptive schemes based on Gaussian process regression have been investigated with respect to the number of necessary grid points and the accuracy of the fundamental modes for a small set of test molecules. A comparison with a related, more sophisticated, and consistent approach by Christiansen et al. is provided. The impact of the positions of the ab initio grid points is discussed for multilevel PESs, for which the computational effort of the individual electronic structure calculations decreases for increasing orders of the n-mode expansion. As a result of that, the ultimate goal is not the maximal reduction of grid points but rather the computational cost, which is not directly related.
Collapse
Affiliation(s)
- Moritz Schneider
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Daniel Born
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Johannes Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Guntram Rauhut
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
6
|
Miessen A, Ollitrault PJ, Tacchino F, Tavernelli I. Quantum algorithms for quantum dynamics. NATURE COMPUTATIONAL SCIENCE 2023; 3:25-37. [PMID: 38177956 DOI: 10.1038/s43588-022-00374-2] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/13/2021] [Accepted: 11/12/2022] [Indexed: 01/06/2024]
Abstract
Among the many computational challenges faced across different disciplines, quantum-mechanical systems pose some of the hardest ones and offer a natural playground for the growing field of quantum technologies. In this Perspective, we discuss quantum algorithmic solutions for quantum dynamics, reporting on the latest developments and offering a viewpoint on their potential and current limitations. We present some of the most promising areas of application and identify possible research directions for the coming years.
Collapse
Affiliation(s)
| | - Pauline J Ollitrault
- IBM Quantum, IBM Research - Zurich, Rüschlikon, Switzerland
- QC Ware, Palo Alto, CA, USA
| | | | | |
Collapse
|
7
|
Töpfer K, Upadhyay M, Meuwly M. Quantitative molecular simulations. Phys Chem Chem Phys 2022; 24:12767-12786. [PMID: 35593769 PMCID: PMC9158373 DOI: 10.1039/d2cp01211a] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 04/30/2022] [Indexed: 11/21/2022]
Abstract
All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.
Collapse
Affiliation(s)
- Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Meenu Upadhyay
- 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.
| |
Collapse
|
8
|
Asaad M, Joubert-Doriol L, Izmaylov AF. Controlling energy conservation in quantum dynamics with independently moving basis functions: Application to multi-configuration Ehrenfest. J Chem Phys 2022; 156:204121. [PMID: 35649883 DOI: 10.1063/5.0087797] [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
Application of the time-dependent variational principle to a linear combination of frozen-width Gaussians describing the nuclear wavefunction provides a formalism where the total energy is conserved. The computational downside of this formalism is that trajectories of individual Gaussians are solutions of a coupled system of differential equations, limiting implementation to serial propagation algorithms. To allow for parallelization and acceleration of the computation, independent trajectories based on simplified equations of motion were suggested. Unfortunately, within practical realizations involving finite Gaussian bases, this simplification leads to breaking the energy conservation. We offer a solution for this problem by using Lagrange multipliers to ensure the energy and norm conservation regardless of basis function trajectories or basis completeness. We illustrate our approach within the multi-configurational Ehrenfest method considering a linear vibronic coupling model.
Collapse
Affiliation(s)
- Mina Asaad
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| | - Loïc Joubert-Doriol
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| |
Collapse
|
9
|
Hino K, Kurashige Y. Matrix Product State Formulation of the MCTDH Theory in Local Mode Representations for Anharmonic Potentials. J Chem Theory Comput 2022; 18:3347-3356. [PMID: 35606892 DOI: 10.1021/acs.jctc.2c00243] [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
The matrix product state formulation of the multiconfiguration time-dependent Hartree theory, MPS-MCTDH, reported previously [Kurashige, J. Chem. Phys. 2018, 19, 194114] is extended to realistic anharmonic potentials with n-mode representations beyond the linear vibronic coupling model. For realistic vibrational potentials, the local mode representation should give a more compact representation of the potentials, i.e., lowering the dimensionality of the entanglements, than the normal coordinates, and the MPS-MCTDH formulation should work more efficiently and maintain the accuracy with a small bond dimension of the MPS ansatz. In fact, it was confirmed that the use of the local coordinates made the interaction matrices diagonal dominant and the number of terms in the n-body expansion of the potentials was significantly reduced. The method was applied to the IR spectrum of the CH2O molecule, the zero-point energies, and the vibrational energy redistribution dynamics of polyenes C2nH2n+2. The results showed that the efficiency of the MPS-MCTDH method is significantly accelerated by the use of local coordinates even if the long-range interactions are included in the potential.
Collapse
Affiliation(s)
- Kentaro Hino
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| | - Yuki Kurashige
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto 606-8502, Japan
| |
Collapse
|
10
|
Hammes-Schiffer S. Theoretical perspectives on non-Born-Oppenheimer effects in chemistry. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200377. [PMID: 35341306 DOI: 10.1098/rsta.2020.0377] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 07/12/2021] [Indexed: 06/14/2023]
Abstract
The Born-Oppenheimer approximation, which assumes that the electrons respond instantaneously to the motion of the nuclei, breaks down for a wide range of chemical and biological processes. The rate constants of such nonadiabatic processes can be calculated using analytical theories, and the real-time nonequilibrium dynamics can be described using numerical atomistic simulations. The selection of an approach depends on the desired balance between accuracy and efficiency. The computational expense of generating potential energy surfaces on-the-fly often favours the use of approximate, robust and efficient methods such as trajectory surface hopping for large, complex systems. The development of formally exact non-Born-Oppenheimer methods and the exploration of well-defined approximations to such methods are critical for providing benchmarks and preparing for the next generation of faster computers. Thus, the parallel development of rigorous but computationally expensive methods and more approximate but computationally efficient methods is optimal. This Perspective briefly summarizes the available theoretical and computational non-Born-Oppenheimer methods and presents examples illustrating how analytical theories and nonadiabatic dynamics simulations can elucidate the fundamental principles of chemical and biological processes. These examples also highlight how theoretical calculations are able to guide the interpretation of experimental data and provide experimentally testable predictions for nonadiabatic processes. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
|
11
|
Modelling Ultrafast Dynamics at a Conical Intersection with Regularized Diabatic States: An Approach Based on Multiplicative Neural Networks. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
12
|
Habershon S. Program Synthesis of Sparse Algorithms for Wave Function and Energy Prediction in Grid-Based Quantum Simulations. J Chem Theory Comput 2022; 18:2462-2478. [PMID: 35293216 PMCID: PMC9009083 DOI: 10.1021/acs.jctc.2c00035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We have recently shown how program synthesis (PS), or the concept of "self-writing code", can generate novel algorithms that solve the vibrational Schrödinger equation, providing approximations to the allowed wave functions for bound, one-dimensional (1-D) potential energy surfaces (PESs). The resulting algorithms use a grid-based representation of the underlying wave function ψ(x) and PES V(x), providing codes which represent approximations to standard discrete variable representation (DVR) methods. In this Article, we show how this inductive PS strategy can be improved and modified to enable prediction of both vibrational wave functions and energy eigenvalues of representative model PESs (both 1-D and multidimensional). We show that PS can generate algorithms that offer some improvements in energy eigenvalue accuracy over standard DVR schemes; however, we also demonstrate that PS can identify accurate numerical methods that exhibit desirable computational features, such as employing very sparse (tridiagonal) matrices. The resulting PS-generated algorithms are initially developed and tested for 1-D vibrational eigenproblems, before solution of multidimensional problems is demonstrated; we find that our new PS-generated algorithms can reduce calculation times for grid-based eigenvector computation by an order of magnitude or more. More generally, with further development and optimization, we anticipate that PS-generated algorithms based on effective Hamiltonian approximations, such as those proposed here, could be useful in direct simulations of quantum dynamics via wave function propagation and evaluation of molecular electronic structure.
Collapse
Affiliation(s)
- Scott Habershon
- Department of Chemistry, University of Warwick, Coventry, CV4 7AL, United Kingdom
| |
Collapse
|
13
|
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
|
14
|
Richings GW, Habershon S. Analyzing Grid-Based Direct Quantum Molecular Dynamics Using Non-Linear Dimensionality Reduction. Molecules 2021; 26:molecules26247418. [PMID: 34946499 PMCID: PMC8708769 DOI: 10.3390/molecules26247418] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Revised: 11/30/2021] [Accepted: 12/02/2021] [Indexed: 11/18/2022] Open
Abstract
Grid-based schemes for simulating quantum dynamics, such as the multi-configuration time-dependent Hartree (MCTDH) method, provide highly accurate predictions of the coupled nuclear and electronic dynamics in molecular systems. Such approaches provide a multi-dimensional, time-dependent view of the system wavefunction represented on a coordinate grid; in the case of non-adiabatic simulations, additional information about the state populations adds a further layer of complexity. As such, wavepacket motion on potential energy surfaces which couple many nuclear and electronic degrees-of-freedom can be extremely challenging to analyse in order to extract physical insight beyond the usual expectation-value picture. Here, we show that non-linear dimensionality reduction (NLDR) methods, notably diffusion maps, can be adapted to extract information from grid-based wavefunction dynamics simulations, providing insight into key nuclear motions which explain the observed dynamics. This approach is demonstrated for 2-D and 9-D models of proton transfer in salicylaldimine, as well as 8-D and full 12-D simulations of cis-trans isomerization in ethene; these simulations demonstrate how NLDR can provide alternative views of wavefunction dynamics, and also highlight future developments.
Collapse
|
15
|
Computational modeling of green hydrogen generation from photocatalytic H2S splitting: Overview and perspectives. JOURNAL OF PHOTOCHEMISTRY AND PHOTOBIOLOGY C: PHOTOCHEMISTRY REVIEWS 2021. [DOI: 10.1016/j.jphotochemrev.2021.100456] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
16
|
Wang C, Ulusoy IS, Aebersold LE, Wilson AK. Multi-configuration electron-nuclear dynamics: An open-shell approach. J Chem Phys 2021; 155:154103. [PMID: 34686063 DOI: 10.1063/5.0063478] [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
The multi-configuration electron-nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.
Collapse
Affiliation(s)
- Cong Wang
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Inga S Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Lucas E Aebersold
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| |
Collapse
|
17
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 176] [Impact Index Per Article: 58.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
18
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
19
|
Westermayr J, Gastegger M, Schütt KT, Maurer RJ. Perspective on integrating machine learning into computational chemistry and materials science. J Chem Phys 2021; 154:230903. [PMID: 34241249 DOI: 10.1063/5.0047760] [Citation(s) in RCA: 75] [Impact Index Per Article: 25.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Machine learning (ML) methods are being used in almost every conceivable area of electronic structure theory and molecular simulation. In particular, ML has become firmly established in the construction of high-dimensional interatomic potentials. Not a day goes by without another proof of principle being published on how ML methods can represent and predict quantum mechanical properties-be they observable, such as molecular polarizabilities, or not, such as atomic charges. As ML is becoming pervasive in electronic structure theory and molecular simulation, we provide an overview of how atomistic computational modeling is being transformed by the incorporation of ML approaches. From the perspective of the practitioner in the field, we assess how common workflows to predict structure, dynamics, and spectroscopy are affected by ML. Finally, we discuss how a tighter and lasting integration of ML methods with computational chemistry and materials science can be achieved and what it will mean for research practice, software development, and postgraduate training.
Collapse
Affiliation(s)
- Julia Westermayr
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Michael Gastegger
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
| | - Kristof T Schütt
- Machine Learning Group, Technische Universität Berlin, 10587 Berlin, Germany
| | - Reinhard J Maurer
- Department of Chemistry, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
20
|
Kanno M, Mignolet B, Remacle F, Kono H. Identification of an ultrafast internal conversion pathway of pyrazine by time-resolved vacuum ultraviolet photoelectron spectrum simulations. J Chem Phys 2021; 154:224304. [PMID: 34241214 DOI: 10.1063/5.0048900] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
The internal conversion from the optically bright S2 (1B2u, ππ*) state to the dark S1 (1B3u, nπ*) state in pyrazine is a standard benchmark for experimental and theoretical studies on ultrafast radiationless decay. Since 2008, a few theoretical groups have suggested significant contributions of other dark states S3 (1Au, nπ*) and S4 (1B2g, nπ*) to the decay of S2. We have previously reported the results of nuclear wave packet simulations [Kanno et al., Phys. Chem. Chem. Phys. 17, 2012 (2015)] and photoelectron spectrum calculations [Mignolet et al., Chem. Phys. 515, 704 (2018)] that support the conventional two-state picture. In this article, the two different approaches, i.e., wave packet simulation and photoelectron spectrum calculation, are combined: We computed the time-resolved vacuum ultraviolet photoelectron spectrum and photoelectron angular distribution for the ionization of the wave packet transferred from S2 to S1. The present results reproduce almost all the characteristic features of the corresponding experimental time-resolved spectrum [Horio et al., J. Chem. Phys. 145, 044306 (2016)], such as a rapid change from a three-band to two-band structure. This further supports the existence and character of the widely accepted pathway (S2 → S1) of ultrafast internal conversion in pyrazine.
Collapse
Affiliation(s)
- Manabu Kanno
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | - Benoît Mignolet
- Theoretical Physical Chemistry, UR MOLSYS, University of Liège, B4000 Liège, Belgium
| | - Françoise Remacle
- Theoretical Physical Chemistry, UR MOLSYS, University of Liège, B4000 Liège, Belgium
| | - Hirohiko Kono
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| |
Collapse
|
21
|
Abstract
Theoretical simulations of electronic excitations and associated processes in molecules are indispensable for fundamental research and technological innovations. However, such simulations are notoriously challenging to perform with quantum mechanical methods. Advances in machine learning open many new avenues for assisting molecular excited-state simulations. In this Review, we track such progress, assess the current state of the art and highlight the critical issues to solve in the future. We overview a broad range of machine learning applications in excited-state research, which include the prediction of molecular properties, improvements of quantum mechanical methods for the calculations of excited-state properties and the search for new materials. Machine learning approaches can help us understand hidden factors that influence photo-processes, leading to a better control of such processes and new rules for the design of materials for optoelectronic applications.
Collapse
|
22
|
Chowdhury SN, Huo P. Non-adiabatic Matsubara dynamics and non-adiabatic ring-polymer molecular dynamics. J Chem Phys 2021; 154:124124. [PMID: 33810665 DOI: 10.1063/5.0042136] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
We present the non-adiabatic Matsubara dynamics, a general framework for computing the time-correlation function (TCF) of electronically non-adiabatic systems. This new formalism is derived based on the generalized Kubo-transformed TCF using the Wigner representation for both the nuclear degrees of freedom and the electronic mapping variables. By dropping the non-Matsubara nuclear normal modes in the quantum Liouvillian and explicitly integrating these modes out from the expression of the TCF, we derived the non-adiabatic Matsubara dynamics approach. Further making the approximation to drop the imaginary part of the Matsubara Liouvillian and enforce the nuclear momentum integral to be real, we arrived at the non-adiabatic ring-polymer molecular dynamics (NRPMD) approach. We have further justified the capability of NRPMD for simulating the non-equilibrium TCF. This work provides the rigorous theoretical foundation for several recently proposed state-dependent RPMD approaches and offers a general framework for developing new non-adiabatic quantum dynamics methods in the future.
Collapse
Affiliation(s)
- Sutirtha N Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
23
|
Ibele LM, Lassmann Y, Martínez TJ, Curchod BFE. Comparing (stochastic-selection) ab initio multiple spawning with trajectory surface hopping for the photodynamics of cyclopropanone, fulvene, and dithiane. J Chem Phys 2021; 154:104110. [DOI: 10.1063/5.0045572] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Yorick Lassmann
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry, Stanford University, Stanford, California 94305, USA and PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
24
|
Panadés-Barrueta RL, Peláez D. Low-rank sum-of-products finite-basis-representation (SOP-FBR) of potential energy surfaces. J Chem Phys 2020; 153:234110. [PMID: 33353311 DOI: 10.1063/5.0027143] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
The sum-of-products finite-basis-representation (SOP-FBR) approach for the automated multidimensional fit of potential energy surfaces (PESs) is presented. In its current implementation, the method yields a PES in the so-called Tucker sum-of-products form, but it is not restricted to this specific ansatz. The novelty of our algorithm lies in the fact that the fit is performed in terms of a direct product of a Schmidt basis, also known as natural potentials. These encode in a non-trivial way all the physics of the problem and, hence, circumvent the usual extra ad hoc and a posteriori adjustments (e.g., damping functions) of the fitted PES. Moreover, we avoid the intermediate refitting stage common to other tensor-decomposition methods, typically used in the context of nuclear quantum dynamics. The resulting SOP-FBR PES is analytical and differentiable ad infinitum. Our ansatz is fully general and can be used in combination with most (molecular) dynamics codes. In particular, it has been interfaced and extensively tested with the Heidelberg implementation of the multiconfiguration time-dependent Hartree quantum dynamical software package.
Collapse
Affiliation(s)
- Ramón L Panadés-Barrueta
- Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM), Université Lille 1, Villeneuve d'Ascq Cedex, France
| | - Daniel Peláez
- Institut des Sciences Moléculaires d'Orsay (ISMO) - UMR 8214, Bât. 520, Université Paris-Saclay, 91405 Orsay Cedex, France
| |
Collapse
|
25
|
Dutta J, Mukherjee S, Naskar K, Ghosh S, Mukherjee B, Ravi S, Adhikari S. The role of electron-nuclear coupling on multi-state photoelectron spectra, scattering processes and phase transitions. Phys Chem Chem Phys 2020; 22:27496-27524. [PMID: 33283826 DOI: 10.1039/d0cp04052e] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present first principle based beyond Born-Oppenheimer (BBO) theory and its applications on various models as well as realistic spectroscopic and scattering processes, where the Jahn-Teller (JT) theory is brought in conjunction with the BBO approach on the phase transition of lanthanide complexes. Over one and half decades, our development of BBO theory is demonstrated with ab initio calculations on representative molecules of spectroscopic interest (NO2 radical, Na3 and K3 clusters, NO3 radical, C6H6+ and 1,3,5-C6H3F3+ radical cations) as well as triatomic reactive scattering processes (H+ + H2 and F + H2). Such an approach exhibits the effect of JT, Renner-Teller (RT) and pseudo Jahn-Teller (PJT) type of interactions. While implementing the BBO theory, we generate highly accurate diabatic potential energy surfaces (PESs) to carry out quantum dynamics calculation and find excellent agreement with experimental photoelectron spectra of spectroscopic systems and cross-sections/rate constants of scattering processes. On the other hand, such electron-nuclear couplings incorporated through JT theory play a crucial role in dictating higher energy satellite transitions in the dielectric function spectra of the LaMnO3 complex. Overall, this article thoroughly sketches the current perspective of the BBO approach and its connection with JT theory with various applications on physical and chemical processes.
Collapse
Affiliation(s)
- Joy Dutta
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India.
| | | | | | | | | | | | | |
Collapse
|
26
|
Dresselhaus T, Bungey CBA, Knowles PJ, Manby FR. Coupling electrons and vibrations in molecular quantum chemistry. J Chem Phys 2020; 153:214114. [PMID: 33291918 DOI: 10.1063/5.0032900] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
We derive an electron-vibration model Hamiltonian in a quantum chemical framework and explore the extent to which such a Hamiltonian can capture key effects of nonadiabatic dynamics. The model Hamiltonian is a simple two-body operator, and we make preliminary steps at applying standard quantum chemical methods to evaluate its properties, including mean-field theory, linear response, and a primitive correlated model. The Hamiltonian can be compared to standard vibronic Hamiltonians, but it is constructed without reference to potential energy surfaces through direct differentiation of the one- and two-electron integrals at a single reference geometry. The nature of the model Hamiltonian in the harmonic and linear-coupling regime is investigated for pyrazine, where a simple time-dependent calculation including electron-vibration correlation is demonstrated to exhibit the well-studied population transfer between the S2 and S1 excited states.
Collapse
Affiliation(s)
- Thomas Dresselhaus
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Callum B A Bungey
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
27
|
Madsen NK, Hansen MB, Christiansen O, Zoccante A. Time-dependent vibrational coupled cluster with variationally optimized time-dependent basis sets. J Chem Phys 2020; 153:174108. [DOI: 10.1063/5.0024428] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Niels Kristian Madsen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Mads Bøttger Hansen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Alberto Zoccante
- Department of Chemistry, University of Aarhus, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
28
|
Morzan UN, Videla PE, Soley MB, Nibbering ETJ, Batista VS. Vibronic Dynamics of Photodissociating ICN from Simulations of Ultrafast X‐Ray Absorption Spectroscopy. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202007192] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Condensed Matter Section The Abdus Salam International Center for Theoretical Physics Strada Costiera 11 34151 Trieste Italy
- Department of Chemistry Yale University P.O. Box 208107 New Haven CT 06520-8107 USA
| | - Pablo E. Videla
- Department of Chemistry Yale University P.O. Box 208107 New Haven CT 06520-8107 USA
- Energy Sciences Institute Yale University P.O. Box 27394 West Haven CT 06516-7394 USA
| | - Micheline B. Soley
- Department of Chemistry Yale University P.O. Box 208107 New Haven CT 06520-8107 USA
- Department of Chemistry and Chemical Biology Harvard University 12 Oxford Street Cambridge MA 02138 USA
- Yale Quantum Institute Yale University P.O. Box 208334 New Haven CT 06520-8263 USA
| | - Erik T. J. Nibbering
- Max Born Institute for Nonlinear Optics and Short Pulse Spectroscopy Max Born Strasse 2A 12489 Berlin Germany
| | - Victor S. Batista
- Department of Chemistry Yale University P.O. Box 208107 New Haven CT 06520-8107 USA
- Energy Sciences Institute Yale University P.O. Box 27394 West Haven CT 06516-7394 USA
| |
Collapse
|
29
|
Morzan UN, Videla PE, Soley MB, Nibbering ETJ, Batista VS. Vibronic Dynamics of Photodissociating ICN from Simulations of Ultrafast X-Ray Absorption Spectroscopy. Angew Chem Int Ed Engl 2020; 59:20044-20048. [PMID: 32691867 PMCID: PMC7693200 DOI: 10.1002/anie.202007192] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Revised: 06/23/2020] [Indexed: 11/07/2022]
Abstract
Ultrafast UV-pump/soft-X-ray-probe spectroscopy is a subject of great interest since it can provide detailed information about dynamical photochemical processes with ultrafast resolution and atomic specificity. Here, we focus on the photodissociation of ICN in the 1 Π1 excited state, with emphasis on the transient response in the soft-X-ray spectral region as described by the ab initio spectral lineshape averaged over the nuclear wavepacket probability density. We find that the carbon K-edge spectral region reveals a rich transient response that provides direct insights into the dynamics of frontier orbitals during the I-CN bond cleavage process. The simulated UV-pump/soft-X-ray-probe spectra exhibit detailed dynamical information, including a time-domain signature for coherent vibration associated with the photogenerated CN fragment.
Collapse
Affiliation(s)
- Uriel N. Morzan
- Condensed Matter SectionThe Abdus Salam International Center for Theoretical PhysicsStrada Costiera 1134151TriesteItaly
- Department of ChemistryYale UniversityP.O. Box 208107New HavenCT06520-8107USA
| | - Pablo E. Videla
- Department of ChemistryYale UniversityP.O. Box 208107New HavenCT06520-8107USA
- Energy Sciences InstituteYale UniversityP.O. Box 27394West HavenCT06516-7394USA
| | - Micheline B. Soley
- Department of ChemistryYale UniversityP.O. Box 208107New HavenCT06520-8107USA
- Department of Chemistry and Chemical BiologyHarvard University12 Oxford StreetCambridgeMA02138USA
- Yale Quantum InstituteYale UniversityP.O. Box 208334New HavenCT06520-8263USA
| | - Erik T. J. Nibbering
- Max Born Institute for Nonlinear Optics and Short Pulse SpectroscopyMax Born Strasse 2A12489BerlinGermany
| | - Victor S. Batista
- Department of ChemistryYale UniversityP.O. Box 208107New HavenCT06520-8107USA
- Energy Sciences InstituteYale UniversityP.O. Box 27394West HavenCT06516-7394USA
| |
Collapse
|
30
|
Richings GW, Habershon S. Direct Grid-Based Nonadiabatic Dynamics on Machine-Learned Potential Energy Surfaces: Application to Spin-Forbidden Processes. J Phys Chem A 2020; 124:9299-9313. [DOI: 10.1021/acs.jpca.0c06125] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
31
|
Westermayr J, Marquetand P. Deep learning for UV absorption spectra with SchNarc: First steps toward transferability in chemical compound space. J Chem Phys 2020; 153:154112. [DOI: 10.1063/5.0021915] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- J. Westermayr
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - P. Marquetand
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Faculty of Chemistry, Data Science @ Uni Vienna, University of Vienna, Währinger Str. 29, 1090 Vienna, Austria
| |
Collapse
|
32
|
Jung KA, Videla PE, Batista VS. Ring-polymer, centroid, and mean-field approximations to multi-time Matsubara dynamics. J Chem Phys 2020; 153:124112. [PMID: 33003707 DOI: 10.1063/5.0021843] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Based on a recently developed generalization of Matsubara dynamics to the multi-time realm, we present a formal derivation of multi-time generalizations of ring-polymer molecular dynamics, thermostatted ring-polymer molecular dynamics (TRPMD), centroid molecular dynamics (CMD), and mean-field Matsubara dynamics. Additionally, we analyze the short-time accuracy of each methodology. We find that for multi-time correlation functions of linear operators, (T)RPMD is accurate up to order t3, while CMD is only correct up to t, indicating a degradation in the accuracy of these methodologies with respect to the single-time counterparts. The present work provides a firm justification for the use of path-integral-based approximations for the calculation of multi-time correlation functions.
Collapse
Affiliation(s)
- Kenneth A Jung
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Pablo E Videla
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
33
|
Westermayr J, Marquetand P. Machine learning and excited-state molecular dynamics. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab9c3e] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
34
|
Schmitz G, Klinting EL, Christiansen O. A Gaussian process regression adaptive density guided approach for potential energy surface construction. J Chem Phys 2020; 153:064105. [DOI: 10.1063/5.0015344] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
35
|
Westermayr J, Gastegger M, Marquetand P. Combining SchNet and SHARC: The SchNarc Machine Learning Approach for Excited-State Dynamics. J Phys Chem Lett 2020; 11:3828-3834. [PMID: 32311258 PMCID: PMC7246974 DOI: 10.1021/acs.jpclett.0c00527] [Citation(s) in RCA: 97] [Impact Index Per Article: 24.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Accepted: 04/20/2020] [Indexed: 05/26/2023]
Abstract
In recent years, deep learning has become a part of our everyday life and is revolutionizing quantum chemistry as well. In this work, we show how deep learning can be used to advance the research field of photochemistry by learning all important properties-multiple energies, forces, and different couplings-for photodynamics simulations. We simplify such simulations substantially by (i) a phase-free training skipping costly preprocessing of raw quantum chemistry data; (ii) rotationally covariant nonadiabatic couplings, which can either be trained or (iii) alternatively be approximated from only ML potentials, their gradients, and Hessians; and (iv) incorporating spin-orbit couplings. As the deep-learning method, we employ SchNet with its automatically determined representation of molecular structures and extend it for multiple electronic states. In combination with the molecular dynamics program SHARC, our approach termed SchNarc is tested on two polyatomic molecules and paves the way toward efficient photodynamics simulations of complex systems.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Michael Gastegger
- Machine
Learning Group, Technical University of
Berlin, 10587 Berlin, Germany
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Str. 29, 1090 Vienna, Austria
| |
Collapse
|
36
|
Westermayr J, Faber FA, Christensen AS, von Lilienfeld OA, Marquetand P. Neural networks and kernel ridge regression for excited states dynamics of CH2NH$_2^+$: From single-state to multi-state representations and multi-property machine learning models. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab88d0] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
|
37
|
Richings GW, Habershon S. A new diabatization scheme for direct quantum dynamics: Procrustes diabatization. J Chem Phys 2020; 152:154108. [DOI: 10.1063/5.0003254] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Gareth W. Richings
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
38
|
Pavošević F, Culpitt T, Hammes-Schiffer S. Multicomponent Quantum Chemistry: Integrating Electronic and Nuclear Quantum Effects via the Nuclear–Electronic Orbital Method. Chem Rev 2020; 120:4222-4253. [DOI: 10.1021/acs.chemrev.9b00798] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Tanner Culpitt
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
39
|
Neville SP, Seidu I, Schuurman MS. Propagative block diagonalization diabatization of DFT/MRCI electronic states. J Chem Phys 2020; 152:114110. [DOI: 10.1063/1.5143126] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Simon P. Neville
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| | - Issaka Seidu
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
| | - Michael S. Schuurman
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, 10 Marie Curie, Ottawa, Ontario K1N 6N5, Canada
- National Research Council of Canada, 100 Sussex Drive, Ottawa, Ontario K1A 0R6, Canada
| |
Collapse
|
40
|
Naskar K, Mukherjee S, Mukherjee B, Ravi S, Mukherjee S, Sardar S, Adhikari S. ADT: A Generalized Algorithm and Program for Beyond Born–Oppenheimer Equations of “N” Dimensional Sub-Hilbert Space. J Chem Theory Comput 2020; 16:1666-1680. [DOI: 10.1021/acs.jctc.9b00948] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Soumya Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Bijit Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Satyam Ravi
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Saikat Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Subhankar Sardar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
- Department of Chemistry, Bhatter College, Dantan, Paschim Medinipur 721426, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
41
|
Dutra M, Wickramasinghe S, Garashchuk S. Quantum Dynamics with the Quantum Trajectory-Guided Adaptable Gaussian Bases. J Chem Theory Comput 2019; 16:18-34. [DOI: 10.1021/acs.jctc.9b00844] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Matthew Dutra
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, United States
| | - Sachith Wickramasinghe
- 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
|
42
|
Zhou W, Mandal A, Huo P. Quasi-Diabatic Scheme for Nonadiabatic On-the-Fly Simulations. J Phys Chem Lett 2019; 10:7062-7070. [PMID: 31665889 DOI: 10.1021/acs.jpclett.9b02747] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
We use the quasi-diabatic (QD) propagation scheme to perform on-the-fly nonadiabatic simulations of the photodynamics of ethylene. The QD scheme enables a seamless interface between accurate diabatic-based quantum dynamics approaches and adiabatic electronic structure calculations, explicitly avoiding any efforts to construct global diabatic states or reformulate the diabatic dynamics approach to the adiabatic representation. Using the partial linearized path-integral approach and the symmetrical quasi-classical approach as the diabatic dynamics methods, the QD propagation scheme enables direct nonadiabatic simulation with complete active space self-consistent field on-the-fly electronic structure calculations. The population dynamics obtained from both approaches are in a close agreement with the quantum wavepacket-based method and outperform the widely used trajectory surface-hopping approach. Further analysis of the ethylene photodeactivation pathways demonstrates the correct predictions of competing processes of nonradiative relaxation mechanism through various conical intersections. This work provides the foundation of using accurate diabatic dynamics approaches and on-the-fly adiabatic electronic structure information to perform ab initio nonadiabatic simulation.
Collapse
Affiliation(s)
- Wanghuai Zhou
- Advanced Functional Material and Photoelectric Technology Research Institution, School of Science , Hubei University of Automotive Technology , Shiyan , Hubei 442002 , People's Republic of China
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Arkajit Mandal
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| | - Pengfei Huo
- Department of Chemistry , University of Rochester , 120 Trustee Road , Rochester , New York 14627 , United States
| |
Collapse
|
43
|
Mukherjee B, Naskar K, Mukherjee S, Ghosh S, Sahoo T, Adhikari S. Beyond Born–Oppenheimer theory for spectroscopic and scattering processes. INT REV PHYS CHEM 2019. [DOI: 10.1080/0144235x.2019.1672987] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Bijit Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata, India
| | - Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata, India
| | - Soumya Mukherjee
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata, India
| | - Sandip Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata, India
| | - Tapas Sahoo
- Department of Chemistry, University of Waterloo, Waterloo, Ontario, Canada
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata, India
| |
Collapse
|
44
|
Kananenka AA, Yao K, Corcelli SA, Skinner JL. Machine Learning for Vibrational Spectroscopic Maps. J Chem Theory Comput 2019; 15:6850-6858. [DOI: 10.1021/acs.jctc.9b00698] [Citation(s) in RCA: 43] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Alexei A. Kananenka
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Kun Yao
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Steven A. Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - J. L. Skinner
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
45
|
Cho KH, Chung S, Rhee YM. Efficiently Transplanting Potential Energy Interpolation Database between Two Systems: Bacteriochlorophyll Case with FMO and LH2 Complexes. J Chem Inf Model 2019; 59:4228-4238. [PMID: 31487163 DOI: 10.1021/acs.jcim.9b00472] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Generating a reliable potential energy surface (PES) is an important issue for studying the dynamics of any system with computational simulations. Interpolation mechanics/molecular mechanics (IM/MM) based on a PES interpolation scheme is a useful tool in that regard as it provides an accuracy of a quantum chemistry (QC) level while maintaining its computational cost comparable to conventional MM force fields. Despite this benefit, constructing the database for interpolation itself is still challenging and time-consuming. Here, we present a method with which we can construct the IM database of one system based on a preexisting data set for another related system. We adopt the case of constructing bacteriochlorophyll PESs for the light-harvesting 2 (LH2) complex by utilizing already available IM database for the BChls from the Fenna-Matthews-Olson (FMO) complex. In this method, the IM database from FMO is first transplanted to LH2 by considering BChl displacement vectors that take into account the geometry differences induced by the protein scaffolds. From this transplanted primitive database entries, a relatively small number of effective ones are selected by a survival process based on a genetic algorithm such that the IM energies evaluated at geometries in a conveniently collected prediction set can closely match with the reference QC energies. The selection process is expedited by using two different levels of basis sets for the QC calculations. To demonstrate the utility of the PES thus constructed, we carry out 1 ns of IM/MM dynamics simulations with the finally optimized BChl database for LH2. Indeed, the energy profiles of the snapshots are found to be closely matching with the reference QC calculation data, with only ∼0.07 eV of errors in the ground- and excited-state energies and ∼0.008 eV of errors in the transition energies. We also show that properly selecting data points is actually quite important for generating an IM PES toward performing molecular dynamics simulations.
Collapse
Affiliation(s)
- Kwang Hyun Cho
- Department of Chemistry , Korea Advanced Institute of Science and Technology (KAIST) , Daejeon 34141 , Korea
| | - Seyoung Chung
- Department of Chemistry , Pohang University of Science and Technology (POSTECH) , Pohang 37673 , Korea
| | - Young Min Rhee
- Department of Chemistry , Korea Advanced Institute of Science and Technology (KAIST) , Daejeon 34141 , Korea
| |
Collapse
|
46
|
Panadés-Barrueta RL, Martínez-Núñez E, Peláez D. Specific Reaction Parameter Multigrid POTFIT (SRP-MGPF): Automatic Generation of Sum-of-Products Form Potential Energy Surfaces for Quantum Dynamical Calculations. Front Chem 2019; 7:576. [PMID: 31475138 PMCID: PMC6702682 DOI: 10.3389/fchem.2019.00576] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Accepted: 07/30/2019] [Indexed: 11/13/2022] Open
Abstract
We present Specific Reaction Parameter Multigrid POTFIT (SRP-MGPF), an automated methodology for the generation of global potential energy surfaces (PES), molecular properties surfaces, e.g., dipole, polarizabilities, etc. using a single random geometry as input. The SRP-MGPF workflow integrates: (i) a fully automated procedure for the global topographical characterization of a (intermolecular) PES based on the Transition State Search Using Chemical Dynamical Simulations (TSSCDS) family of methods;i (ii) the global optimization of the parameters of a semiempirical Hamiltonian in order to reproduce a given level of electronic structure theory; and (iii) a tensor decomposition algorithm which turns the resulting SRP-PES into sum of products (Tucker) form with the Multigrid POTFIT algorithm. The latter is necessary for quantum dynamical studies within the Multiconfiguration Time-Dependent Hartree (MCTDH) quantum dynamics method. To demonstrate our approach, we have applied our methodology to the cis-trans isomerization reaction in HONO in full dimensionality (6D). The resulting SRP-PES has been validated through the computation of classical on-the-fly dynamical calculations as well as calculations of the lowest vibrational eigenstates of HONO as well as high-energy wavepacket propagations.
Collapse
Affiliation(s)
- Ramón L. Panadés-Barrueta
- Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM), Université de Lille, Villeneuve-d'Ascq, France
| | - Emilio Martínez-Núñez
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela, Spain
| | - Daniel Peláez
- Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM), Université de Lille, Villeneuve-d'Ascq, France
| |
Collapse
|
47
|
Jung KA, Videla PE, Batista VS. Multi-time formulation of Matsubara dynamics. J Chem Phys 2019; 151:034108. [PMID: 31325942 DOI: 10.1063/1.5110427] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Matsubara dynamics has recently emerged as the most general form of a quantum-Boltzmann-conserving classical dynamics theory for the calculation of single-time correlation functions. Here, we present a generalization of Matsubara dynamics for the evaluation of multitime correlation functions. We first show that the Matsubara approximation can also be used to approximate the two-time symmetrized double Kubo transformed correlation function. By a straightforward extension of these ideas to the multitime realm, a multitime Matsubara dynamics approximation can be obtained for the multitime fully symmetrized Kubo transformed correlation function. Although not a practical method, due to the presence of a phase-term, this multitime formulation of Matsubara dynamics represents a benchmark theory for future development of Boltzmann preserving semiclassical approximations to general higher order multitime correlation functions. It also reveals a connection between imaginary time-ordering in the path integral and the classical dynamics of multitime correlation functions.
Collapse
Affiliation(s)
- Kenneth A Jung
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Pablo E Videla
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
48
|
Schmitz G, Godtliebsen IH, Christiansen O. Machine learning for potential energy surfaces: An extensive database and assessment of methods. J Chem Phys 2019; 150:244113. [DOI: 10.1063/1.5100141] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | | | - Ove Christiansen
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| |
Collapse
|
49
|
Chowdhury SN, Huo P. State dependent ring polymer molecular dynamics for investigating excited nonadiabatic dynamics. J Chem Phys 2019; 150:244102. [DOI: 10.1063/1.5096276] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sutirtha N. Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
50
|
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
|