1
|
Weike N, Fritsch F, Eisfeld W. Compensation States Approach in the Hybrid Diabatization Scheme: Extension to Multidimensional Data and Properties. J Phys Chem A 2024; 128:4353-4368. [PMID: 38748493 DOI: 10.1021/acs.jpca.4c01134] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2024]
Abstract
The diabatization of reactive systems for more than just a couple of states is a very demanding problem and generally requires advanced diabatization techniques. Especially for dissociative processes, the drastic changes in the adiabatic wave functions often would require large diabatic state bases, which quickly become impractical. Recently, we addressed this problem by the compensation states approach developed in the context of our hybrid diabatization scheme. This scheme utilizes wave function as well as energy data in combination with a diabatic potential model. In regions where the initial diabatic state basis becomes insufficient for an appropriate representation of the adiabatic states, new model states are generated. The new model states compensate for the state space not spanned by the initial diabatic basis. Such a compensation state is obtained by projecting the initial diabatic state space out of the adiabatic wave function. This yields a very efficient basis representation of the electronic Hamiltonian. The present work presents two new aspects. First, it is shown how other operators like the spin-orbit operator in the framework of the Effective Relativistic Coupling by Asymptotic Representation (ERCAR) can be evaluated in this compact model state space without losing the correct wave function information and accuracy. Second, the extension of the approach to multidimensional potential energy surface models is presented for methyl iodide including the C-I dissociation coordinate and the angular H3C-I bending coordinates.
Collapse
Affiliation(s)
- Nicole Weike
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Fabian Fritsch
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
2
|
Reddy SR, Coto PB, Thoss M. Intramolecular singlet fission: Quantum dynamical simulations including the effect of the laser field. J Chem Phys 2024; 160:194306. [PMID: 38767260 DOI: 10.1063/5.0209546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Accepted: 04/30/2024] [Indexed: 05/22/2024] Open
Abstract
In the previous work [Reddy et al., J. Chem. Phys. 151, 044307 (2019)], we have analyzed the dynamics of the intramolecular singlet fission process in a series of prototypical pentacene-based dimers, where the pentacene monomers are covalently bonded to a phenylene linker in ortho, meta, and para positions. The results obtained were qualitatively consistent with the experimental data available, showing an ultrafast population of the multiexcitonic state that mainly takes place via a mediated (superexchange-like) mechanism involving charge transfer and doubly excited states. Our results also highlighted the instrumental role of molecular vibrations in the process as a sizable population of the multiexcitonic state could only be obtained through vibronic coupling. Here, we extend these studies and investigate the effect of the laser field on the dynamics of intramolecular singlet fission by explicitly including the coupling to the laser field in our model. In this manner, and by selectively tuning the laser field to the different low-lying absorption bands of the systems investigated, we analyze the wavelength dependence of the intramolecular singlet fission process. In addition, we have also analyzed how the nature of the initially photoexcited electronic state (either localized or delocalized) affects its dynamics. Altogether, our results provide new insights into the design of intramolecular singlet fission-active molecules.
Collapse
Affiliation(s)
- S Rajagopala Reddy
- Department of Chemistry, School of Chemical Sciences and Pharmacy, Central University of Rajasthan, NH-8, Bandarsindri, Ajmer, Rajasthan 305817, India
| | - Pedro B Coto
- Materials Physics Center (CFM), Spanish National Research Council (CSIC) and Donostia International Physics Center (DIPC), Paseo Manuel de Lardizabal 5, 20018 Donostia-San Sebastián, Spain
| | - Michael Thoss
- Institute of Physics, Albert-Ludwigs University Freiburg, Hermann-Herder-Str. 3, 79104 Freiburg, Germany
| |
Collapse
|
3
|
Gu B. Nonadiabatic Conical Intersection Dynamics in the Local Diabatic Representation with Strang Splitting and Fourier Basis. J Chem Theory Comput 2024; 20:2711-2718. [PMID: 38536965 DOI: 10.1021/acs.jctc.3c01317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
We develop and implement an exact conical intersection nonadiabatic wave packet dynamics method that combines the local diabatic representation, Strang splitting for the total molecular propagator, and discrete variable representation with uniform grids. By employing the local diabatic representation, this method captures all nonadiabatic effects, including nonadiabatic transitions, electronic coherences, and geometric phase. Moreover, it is free of singularities in the first and second derivative couplings and does not require the electronic wave function to be continuous with respect to the nuclear coordinates. We further show that in contrast to the adiabatic representation, the split-operator method can be directly applied to the full molecular propagator with the locally diabatic ansatz. The Fourier series, employed as the primitive nuclear basis functions, is universal and can be applied to all types of reactive coordinates. The combination of local diabatic representation, Strang splitting, and Fourier basis allows numerically exact modeling of conical intersection quantum dynamics directly with adiabatic electronic states that can be obtained from standard electronic structure computations.
Collapse
Affiliation(s)
- Bing Gu
- Department of Chemistry and Department of Physics, Westlake University, Hangzhou, Zhejiang 310030, China
- Institute of Natural Sciences, Westlake Institute for Advanced Study, Hangzhou, Zhejiang 310024, China
| |
Collapse
|
4
|
Brady RP, Drury C, Yurchenko SN, Tennyson J. Numerical Equivalence of Diabatic and Adiabatic Representations in Diatomic Molecules. J Chem Theory Comput 2024; 20:2127-2139. [PMID: 38171539 PMCID: PMC10938500 DOI: 10.1021/acs.jctc.3c01150] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Revised: 12/19/2023] [Accepted: 12/20/2023] [Indexed: 01/05/2024]
Abstract
The (time-independent) Schrödinger equation for atomistic systems is solved by using the adiabatic potential energy curves (PECs) and the associated adiabatic approximation. In cases where interactions between electronic states become important, the associated nonadiabatic effects are taken into account via derivative couplings (DDRs), also known as nonadiabatic couplings (NACs). For diatomic molecules, the corresponding PECs in the adiabatic representation are characterized by avoided crossings. The alternative to the adiabatic approach is the diabatic representation obtained via a unitary transformation of the adiabatic states by minimizing the DDRs. For diatomics, the diabatic representation has zero DDR and nondiagonal diabatic couplings ensue. The two representations are fully equivalent and so should be the rovibronic energies and wave functions, which result from the solution of the corresponding Schrödinger equations. We demonstrate (for the first time) the numerical equivalence between the adiabatic and diabatic rovibronic calculations of diatomic molecules using the ab initio curves of yttrium oxide (YO) and carbon monohydride (CH) as examples of two-state systems, where YO is characterized by a strong NAC, while CH has a strong diabatic coupling. Rovibronic energies and wave functions are computed using a new diabatic module implemented in the variational rovibronic code Duo. We show that it is important to include both the diagonal Born-Oppenheimer correction and nondiagonal DDRs. We also show that the convergence of the vibronic energy calculations can strongly depend on the representation of nuclear motion used and that no one representation is best in all cases.
Collapse
Affiliation(s)
- Ryan P. Brady
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, U.K.
| | - Charlie Drury
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, U.K.
| | - Sergei N. Yurchenko
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, U.K.
| | - Jonathan Tennyson
- Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, U.K.
| |
Collapse
|
5
|
Han S, Xie C, Hu X, Yarkony DR, Guo H, Xie D. Quantum Dynamics of Photodissociation: Recent Advances and Challenges. J Phys Chem Lett 2023; 14:10517-10530. [PMID: 37970789 DOI: 10.1021/acs.jpclett.3c02735] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
Recent advances in constructing accurate potential energy surfaces and nonadiabatic couplings from high-level ab initio data have revealed detailed potential landscapes in not only the ground electronic state but also excited ones. They enabled quantitatively accurate characterization of photoexcited reactive systems using quantum mechanical methods. In this Perspective, we survey the recent progress in quantum mechanical studies of adiabatic and nonadiabatic photodissociation dynamics, focusing on initial state control and product energy disposal. These new insights helped to understand quantum effects in small prototypical systems, and the results serve as benchmarks for developing more approximate theoretical methods.
Collapse
Affiliation(s)
- Shanyu Han
- International Center for Isotope Effects Research, School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China
| | - Changjian Xie
- Institute of Modern Physics, Shaanxi Key Laboratory for Theoretical Physics Frontiers, Northwest University, Xi'an 710127, China
| | - Xixi Hu
- Kuang Yaming Honors School, Institute for Brain Sciences, Nanjing University, Nanjing 210023, China
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, Center for Computational Chemistry, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|
6
|
Mutsuji A, Saita K, Maeda S. An energy decomposition and extrapolation scheme for evaluating electron transfer rate constants: a case study on electron self-exchange reactions of transition metal complexes. RSC Adv 2023; 13:32097-32103. [PMID: 37920761 PMCID: PMC10619204 DOI: 10.1039/d3ra05784d] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2023] [Accepted: 10/23/2023] [Indexed: 11/04/2023] Open
Abstract
A simple approach to the analysis of electron transfer (ET) reactions based on energy decomposition and extrapolation schemes is proposed. The present energy decomposition and extrapolation-based electron localization (EDEEL) method represents the diabatic energies for the initial and final states using the adiabatic energies of the donor and acceptor species and their complex. A scheme for the efficient estimation of ET rate constants is also proposed. EDEEL is semi-quantitative by directly evaluating the seam-of-crossing region of two diabatic potentials. In a numerical test, EDEEL successfully provided ET rate constants for electron self-exchange reactions of thirteen transition metal complexes with reasonable accuracy. In addition, its energy decomposition and extrapolation schemes provide all the energy values required for activation-strain model (ASM) analysis. The ASM analysis using EDEEL provided rational interpretations of the variation of the ET rate constants as a function of the transition metal complexes. These results suggest that EDEEL is useful for efficiently evaluating ET rate constants and obtaining a rational understanding of their magnitudes.
Collapse
Affiliation(s)
- Akihiro Mutsuji
- Graduate School of Chemical Sciences and Engineering, Hokkaido University Sapporo Hokkaido 060-8628 Japan
| | - Kenichiro Saita
- Department of Chemistry, Graduate School of Science, Hokkaido University Sapporo Hokkaido 060-0810 Japan
| | - Satoshi Maeda
- Department of Chemistry, Graduate School of Science, Hokkaido University Sapporo Hokkaido 060-0810 Japan
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University Sapporo Hokkaido 001-0021 Japan
- ERATO Maeda Artificial Intelligence for Chemical Reaction Design and Discovery Project, Hokkaido University Sapporo Hokkaido 060-0810 Japan
- Research and Services Division of Materials Data and Integrated System (MaDIS), National Institute for Materials Science (NIMS) Tsukuba Ibaraki 305-0044 Japan
| |
Collapse
|
7
|
Hennefarth MR, Hermes MR, Truhlar DG, Gagliardi L. Linearized Pair-Density Functional Theory. J Chem Theory Comput 2023. [PMID: 37207365 DOI: 10.1021/acs.jctc.3c00207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Multiconfiguration pair-density functional theory (MC-PDFT) is a post-SCF multireference method that has been successful at computing ground- and excited-state energies. However, MC-PDFT is a single-state method in which the final MC-PDFT energies do not come from diagonalization of a model-space Hamiltonian matrix, and this can lead to inaccurate topologies of potential energy surfaces near locally avoided crossings and conical intersections. Therefore, in order to perform physically correct ab initio molecular dynamics with electronically excited states or to treat Jahn-Teller instabilities, it is necessary to develop a PDFT method that recovers the correct topology throughout the entire nuclear configuration space. Here we construct an effective Hamiltonian operator, called the linearized PDFT (L-PDFT) Hamiltonian, by expanding the MC-PDFT energy expression to first order in a Taylor series of the wave function density. Diagonalization of the L-PDFT Hamiltonian gives the correct potential energy surface topology near conical intersections and locally avoided crossings for a variety of challenging cases including phenol, methylamine, and the spiro cation. Furthermore, L-PDFT outperforms MC-PDFT and previous multistate PDFT methods for predicting vertical excitations from a variety of representative organic chromophores.
Collapse
Affiliation(s)
- Matthew R Hennefarth
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute, and Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R Hermes
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute, and Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Argonne National Laboratory, 9700 South Cass Avenue, Lemont, Illinois 60439, United States
| |
Collapse
|
8
|
Merritt ICD, Jacquemin D, Vacher M. Nonadiabatic Coupling in Trajectory Surface Hopping: How Approximations Impact Excited-State Reaction Dynamics. J Chem Theory Comput 2023; 19:1827-1842. [PMID: 36897995 DOI: 10.1021/acs.jctc.2c00968] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023]
Abstract
Photochemical reactions are widely modeled using the popular trajectory surface hopping (TSH) method, an affordable mixed quantum-classical approximation to the full quantum dynamics of the system. TSH is able to account for nonadiabatic effects using an ensemble of trajectories, which are propagated on a single potential energy surface at a time and which can hop from one electronic state to another. The occurrences and locations of these hops are typically determined using the nonadiabatic coupling between electronic states, which can be assessed in a number of ways. In this work, we benchmark the impact of some approximations to the coupling term on the TSH dynamics for several typical isomerization and ring-opening reactions. We have identified that two of the schemes tested, the popular local diabatization scheme and a scheme based on biorthonormal wave function overlap implemented in the OpenMOLCAS code as part of this work, reproduce at a much reduced cost the dynamics obtained using the explicitly calculated nonadiabatic coupling vectors. The other two schemes tested can give different results, and in some cases, even entirely incorrect dynamics. Of these two, the scheme based on configuration interaction vectors gives unpredictable failures, while the other scheme based on the Baeck-An approximation systematically overestimates hopping to the ground state as compared to the reference approaches.
Collapse
Affiliation(s)
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
9
|
Rani VJ, Kanakati AK, Mahapatra S. Multi-state and Multi-mode Vibronic Coupling Effects in the Photoionization Spectroscopy of Acetaldehyde. J Phys Chem A 2022; 126:6581-6593. [PMID: 36126257 DOI: 10.1021/acs.jpca.2c03073] [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
Multi-state and multi-mode vibronic dynamics in the seven energetically low-lying (X~2A', A~2A″, B~2A', C~2A', D~2A″, E~2A', and F~2A') electronic states of the acetaldehyde radical cation is theoretically studied in this article. Adiabatic energies of these electronic states are calculated by ab initio quantum chemistry methods. A vibronic coupling model of seven electronic states is constructed in a diabatic electronic basis to carry out the first-principles nuclear dynamics study. The vibronic spectrum is calculated and compared with the experimental findings reported in the literature. The progressions of vibrational modes found in the spectrum are assigned. The findings reveal that the X~2A' and F~2A' electronic states are energetically well-separated from the other electronic states and the remaining states (A~2A″ to E~2A') are energetically very close or even quasi-degenerate at the equilibrium geometry of the reference electronic ground state of acetaldehyde. The energetic proximity of A~2A″ to E~2A' electronic states results in multiple multi-state conical intersections. The impact of electronic nonadiabatic interactions due to conical intersections on the vibronic structure of the photoionization band and nonradiative internal conversion dynamics is discussed.
Collapse
Affiliation(s)
| | | | - S Mahapatra
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| |
Collapse
|
10
|
Ren M, Liu X, Zhang L, Lin X, Wu W, Chen Z. Compact and accurate ab initio valence bond wave functions for electron transfer: the classic but challenging covalent-ionic interaction in LiF. J Chem Phys 2022; 157:084106. [DOI: 10.1063/5.0097614] [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 paper combines the valence bond block diabatization approach (VBBDA) with the idea of orbital breathing. With highly compact wave functions, the breathing orbital valence bond (BOVB) method is applied to investigate several atomic and molecular properties including the electron affinity of F, the adiabatic and diabatic potential energy curves and the dipole moment curves of the two lowest-lying 1Σ+ states, the electronic coupling curve and the crossing distance of the two diabatic states, and the spectroscopic constants of the ground states for LiF. The configuration selection scheme proposed in this work is quite general, requiring only the selection of several de-excitation and excitation orbitals in a sense like the restricted active space self-consistent field method. Practically, this is also the first time that BOVB results are extrapolated to complete basis set limit. Armed with the chemical intuition provided by VB theory, the classic but challenging covalent-ionic interaction in the title molecule is not only conceptually interpreted but is also accurately computed.
Collapse
Affiliation(s)
- Mingxing Ren
- College of Chemistry and Chemical Engineering, Xiamen University, China
| | | | - Lina Zhang
- College of Chemistry and Chemical Engineering, Xiamen University, China
| | | | - Wei Wu
- Chemistry, Xiamen University, China
| | | |
Collapse
|
11
|
Shu Y, Varga Z, Kanchanakungwankul S, Zhang L, Truhlar DG. Diabatic States of Molecules. J Phys Chem A 2022; 126:992-1018. [PMID: 35138102 DOI: 10.1021/acs.jpca.1c10583] [Citation(s) in RCA: 31] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zoltan Varga
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
12
|
Han S, de Oliveira-Filho AGS, Shu Y, Truhlar DG, Guo H. Semiclassical Trajectory Studies of Reactive and Nonreactive Scattering of OH(A 2S+) by H2 Based on an Improved Full-Dimensional Ab Initio Diabatic Potential Energy Matrix. Chemphyschem 2022; 23:e202200039. [PMID: 35179813 DOI: 10.1002/cphc.202200039] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2022] [Revised: 02/17/2022] [Indexed: 11/06/2022]
Abstract
We present a new full-dimensional diabatic potential energy matrix (DPEM) for electronically nonadiabatic collisions of OH( A 2 Σ + ) with H 2 , and we calculate the probabilities of electronically adiabatic inelastic collisions, nonreactive quenching, and reactive quenching to form H 2 O + H. The DPEM was fitted using a many-body expansion with permutationally invariant polynomials in bond-order functions to represent the many-body part. The dynamics calculations were carried out with the fewest-switches with time uncertainty and stochastic decoherence (FSTU/SD) semiclassical trajectory method. We present results both for head-on collisions (impact parameter b equal to zero) and for a full range of impact parameters. The results are compared to experiment and to earlier FSTU/SD and quantum dynamics calculations with a previously published DPEM. The various theoretical results all agree that nonreactive quenching dominates reactive quenching, but there are quantitative differences between the two DPEMs and between the b = 0 results and the all- b results, especially for the probability of reactive quenching.
Collapse
Affiliation(s)
- Shanyu Han
- University of New Mexico - Albuquerque: The University of New Mexico, Chemistry and Chemical Biology, UNITED STATES
| | | | - Yinan Shu
- University of Minnesota Twin Cities: University of Minnesota Twin Cities, Chemistry, UNITED STATES
| | - Donald G Truhlar
- University of Minnesota Twin Cities Campus: University of Minnesota Twin Cities, Chemistry, UNITED STATES
| | - Hua Guo
- University of New Mexico, Department of Chemistry, Clark Hall, 87131, Albuquerque, UNITED STATES
| |
Collapse
|
13
|
Ananth N. Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics. Annu Rev Phys Chem 2022; 73:299-322. [PMID: 35081325 DOI: 10.1146/annurev-physchem-082620-021809] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This review focuses on a recent class of path-integral-based methods that simulate nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. The remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, USA;
| |
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
|
Guan Y, Xie C, Yarkony DR, Guo H. High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics. Phys Chem Chem Phys 2021; 23:24962-24983. [PMID: 34473156 DOI: 10.1039/d1cp03008f] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Nonadiabatic dynamics, which goes beyond the Born-Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.
Collapse
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Changjian Xie
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China.
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico, 87131, USA.
| |
Collapse
|
16
|
Ren M, Zhang L, Jiao Y, Chen Z, Wu W. Extended Mulliken-Hush Method with Applications to the Theoretical Study of Electron Transfer. J Chem Theory Comput 2021; 17:6861-6875. [PMID: 34605634 DOI: 10.1021/acs.jctc.1c00603] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A novel adiabatic-to-diabatic (ATD) transformation strategy, namely, the extended Mulliken-Hush (XMH) method, is proposed to evaluate diabatic properties including electronic couplings, potential energy surfaces, and their crossings. The XMH method is developed by adopting our recently proposed ATD transformation formula of a general vectorial physical observable, in which a useful ATD transformation is further determined by using an auxiliary dipole between localized frontier orbitals as a simple approximation of the diabatic transition dipole. The XMH method is simple and practical that provides a flexible way to construct diabatic states. To some extent, it can be regarded as an extension of the generalized Mulliken-Hush (GMH) method since the latter takes a stronger approximation, in which the diabatic transition dipole is assumed to be vanishing. Test calculations on the HeH2+ system show that the electronic couplings predicted by the XMH method are closer to the ones calculated by the valence bond block-diagonalization approach than the GMH ones since the XMH method takes into account both the magnitude and direction of the diabatic transition dipole, which is consistent with the properties of this molecule. In the study of electron transfer in the two kinds of donor-bridge-acceptor systems, the XMH method maintains the simplicity of the GMH method and gives reasonable results even when the latter fails, wherein the diabatic transition dipole is nearly perpendicular to the difference of the initial and final adiabatic dipoles. More importantly, the XMH method can be easily combined with high-level electronic structure methods, in which the properties of the ground and excited states may be more accurately calculated, and hence, one may expect that further development of the XMH method would result in a general computational model for studying electron transfer reactions.
Collapse
Affiliation(s)
- Mingxing Ren
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Lina Zhang
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Yang Jiao
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Zhenhua Chen
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
17
|
Bai X, Guo X, Wang L. Machine Learning Approach to Calculate Electronic Couplings between Quasi-diabatic Molecular Orbitals: The Case of DNA. J Phys Chem Lett 2021; 12:10457-10464. [PMID: 34672582 DOI: 10.1021/acs.jpclett.1c03053] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Diabatization of one-electron states in flexible molecular aggregates is a great challenge due to the presence of surface crossings between molecular orbital (MO) levels and the complex interaction between MOs of neighboring molecules. In this work, we present an efficient machine learning approach to calculate electronic couplings between quasi-diabatic MOs without the need of nonadiabatic coupling calculations. Using MOs of rigid molecules as references, the MOs that can be directly regarded to be quasi-diabatic in molecular dynamics are selected out, state tracked, and phase corrected. On the basis of this information, artificial neural networks are trained to characterize the structure-dependent onsite energies of quasi-diabatic MOs and the intermolecular electronic couplings. A representative sequence of DNA is systematically studied as an illustration. Smooth time evolution of electronic couplings in all base pairs is obtained with quasi-diabatic MOs. In particular, our method can calculate electronic couplings between different quasi-diabatic MOs independently, and thus, this possesses unique advantages in many applications.
Collapse
Affiliation(s)
- Xin Bai
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Xin Guo
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
18
|
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: 162] [Impact Index Per Article: 54.0] [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
|
19
|
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
|
20
|
Fumanal M, Corminboeuf C. Pushing the Limits of the Donor-Acceptor Copolymer Strategy for Intramolecular Singlet Fission. J Phys Chem Lett 2021; 12:7270-7277. [PMID: 34318679 DOI: 10.1021/acs.jpclett.1c01986] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Donor-acceptor (D-A) copolymers have shown great potential for intramolecular singlet fission (iSF). Nonetheless, very few design principles exist for optimizing these systems for iSF, with very little knowledge about how to engineer them for this purpose. In recent work, a fundamental trade-off between the main electronic ingredients required for iSF capable D-A coplanar copolymers was revealed. Still, further investigations are needed to understand these limitations and learn how to bypass them. In this work, we propose to induce torsion as an effective way to circumvent the limits of the coplanar approach. We disclose the potential of noncoplanar copolymers with inherently low triplet energies that encompass all the characteristics required for iSF beyond the limiting values associated with fully coplanar systems. Our findings shed some light on the electronic structure aspects of D-A copolymers for iSF and offer a new avenue for the rational design of novel promising candidates.
Collapse
Affiliation(s)
- Maria Fumanal
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Clémence Corminboeuf
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
21
|
Yin Z, Braams BJ, Fu B, Zhang DH. Neural Network Representation of Three-State Quasidiabatic Hamiltonians Based on the Transformation Properties from a Valence Bond Model: Three Singlet States of H3+. J Chem Theory Comput 2021; 17:1678-1690. [DOI: 10.1021/acs.jctc.0c01336] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zhengxi Yin
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Bastiaan J. Braams
- Centrum Wiskunde & Informatica (CWI), the Dutch National Center for Mathematics and Computer Science, 1098 XG Amsterdam, Netherlands
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, P. R. China
| |
Collapse
|
22
|
Paz ASP, Glover WJ. Diabatic Many-Body Expansion: Development and Application to Charge-Transfer Reactions. J Chem Theory Comput 2021; 17:1497-1511. [DOI: 10.1021/acs.jctc.0c01231] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Amiel S. P. Paz
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshang Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - William J. Glover
- NYU Shanghai, 1555 Century Avenue, Shanghai 200122, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshang Road North, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| |
Collapse
|
23
|
Wang YC, Feng S, Liang W, Zhao Y. Electronic Couplings for Photoinduced Charge Transfer and Excitation Energy Transfer Based on Fragment Particle-Hole Densities. J Phys Chem Lett 2021; 12:1032-1039. [PMID: 33470827 DOI: 10.1021/acs.jpclett.0c03514] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
A new scheme is proposed to calculate the electronic couplings for photoinduced charge transfer and excitation energy transfer for both singlet and triplet states. In this scheme, the locally excited and charge-transfer states are constructed from the adiabatic ones by maximally localizing the particle (i.e., electron) and hole densities in terms of predefined molecular fragments. The construction process, after which the electronic couplings are directly obtained, is highly efficient and can be combined with various kinds of preliminary electronic structure calculations as long as the adiabatic excitation energies and transition densities are available. The method also applies to the systems with multiple charge or excitation centers. Its validity is demonstrated by the applications to the 6,13-dichloropentacene dimer and tetramer and the C60-Zn porphyrin dyad. The results reveal that the environment has a strong impact on the electronic couplings and can even enlarge those for long-range charge transfer by several orders of magnitude.
Collapse
Affiliation(s)
- Yu-Chen Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Shishi Feng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| | - Yi Zhao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People's Republic of China
| |
Collapse
|
24
|
Mao Y, Montoya-Castillo A, Markland TE. Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer. J Chem Phys 2020; 153:244111. [PMID: 33380087 DOI: 10.1063/5.0035593] [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/14/2022] Open
Abstract
Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time-dependent density functional theory, and exorbitant computational cost when combined with high-level wavefunction-based methods. Here, we introduce a density functional theory (DFT)-based diabatization scheme that directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs), which we denote as Δ-ALMO(MSDFT2). We demonstrate that our method, which combines ALMO calculations with the ΔSCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.
Collapse
Affiliation(s)
- Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | | | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
25
|
Fumanal M, Corminboeuf C. Direct, Mediated, and Delayed Intramolecular Singlet Fission Mechanism in Donor-Acceptor Copolymers. J Phys Chem Lett 2020; 11:9788-9794. [PMID: 33147966 DOI: 10.1021/acs.jpclett.0c03076] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Donor-acceptor (D-A) extended copolymers have shown great potential to be exploited for intramolecular singlet fission (iSF) because of their modular tunability and intrinsic ability to incorporate low-lying charge-transfer (CT) and a triplet-pair (1TT) states. While the SF mechanism has been widely debated in homo- and heterodimers, little is known about the singlet splitting process in A-D-A copolymer trimers. Unlike traditional two-site SF, the process of iSF in D-A copolymers involves three molecular units consisting of two A's and one D following an A-D-A polymeric chain. This scenario is, therefore, different from that of the homodimer analogues in terms of which states (if any) may drive the SF process. In this work, we identify how singlet splitting occurs in prototypical iSF D-A copolymer poly(benzodithiophene-alt-thiophene-1,1-dioxide) by means of wave packet propagations on the basis of the linear vibronic coupling model Hamiltonian. Our results reveal that three different mechanisms drive the S1 → 1TT population transfer via antisymmetric and symmetric vibrational motion, including two favorable mechanisms of direct and mediated interactions, as well as a parasitic decay pathway that potentially delays the process. Remarkably, we uncover the interplay between an upper state of marked multiexcitonic character and a low-lying CT state in balancing the splitting efficiency, which anticipates their major role in defining future guidelines for the molecular design of D-A copolymers for iSF.
Collapse
Affiliation(s)
- Maria Fumanal
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Clémence Corminboeuf
- Laboratory for Computational Molecular Design, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
26
|
Li X, Parrish RM, Martínez TJ. An ab initio exciton model for singlet fission. J Chem Phys 2020; 153:184116. [PMID: 33187442 DOI: 10.1063/5.0028605] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
We present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet-multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.
Collapse
Affiliation(s)
- Xin Li
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Robert M Parrish
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
27
|
Sasmal S, Vendrell O. Non-adiabatic quantum dynamics without potential energy surfaces based on second-quantized electrons: Application within the framework of the MCTDH method. J Chem Phys 2020; 153:154110. [PMID: 33092359 DOI: 10.1063/5.0028116] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
A first principles quantum formalism to describe the non-adiabatic dynamics of electrons and nuclei based on a second quantization representation (SQR) of the electronic motion combined with the usual representation of the nuclear coordinates is introduced. This procedure circumvents the introduction of potential energy surfaces and non-adiabatic couplings, providing an alternative to the Born-Oppenheimer approximation. An important feature of the molecular Hamiltonian in the mixed first quantized representation for the nuclei and the SQR representation for the electrons is that all degrees of freedom, nuclear positions and electronic occupations, are distinguishable. This makes the approach compatible with various tensor decomposition Ansätze for the propagation of the nuclear-electronic wavefunction. Here, we describe the application of this formalism within the multi-configuration time-dependent Hartree framework and its multilayer generalization, corresponding to Tucker and hierarchical Tucker tensor decompositions of the wavefunction, respectively. The approach is applied to the calculation of the photodissociation cross section of the HeH+ molecule under extreme ultraviolet irradiation, which features non-adiabatic effects and quantum interferences between the two possible fragmentation channels, He + H+ and He+ + H. These calculations are compared with the usual description based on ab initio potential energy surfaces and non-adiabatic coupling matrix elements, which fully agree. The proof-of-principle calculations serve to illustrate the advantages and drawbacks of this formalism, which are discussed in detail, as well as possible ways to overcome them. We close with an outlook of possible application domains where the formalism might outperform the usual approach, for example, in situations that combine a strong static correlation of the electrons with non-adiabatic electronic-nuclear effects.
Collapse
Affiliation(s)
- Sudip Sasmal
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuneheimer Feld 229, 69120 Heidelberg, Germany
| | - Oriol Vendrell
- Theoretische Chemie, Physikalisch-Chemisches Institut, Universität Heidelberg, Im Neuneheimer Feld 229, 69120 Heidelberg, Germany
| |
Collapse
|
28
|
Han S, Wang Y, Guan Y, Yarkony DR, Guo H. Impact of Diabolical Singular Points on Nonadiabatic Dynamics and a Remedy: Photodissociation of Ammonia in the First Band. J Chem Theory Comput 2020; 16:6776-6784. [DOI: 10.1021/acs.jctc.0c00811] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Shanyu Han
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yucheng Wang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
29
|
Abstract
Understanding nonadiabatic dynamics is important for chemical and physical processes involving multiple electronic states. Direct nonadiabatic dynamics simulations are often employed to observe such processes on a femtosecond time scale. One often needs to do the simulation on a longer time scale, but direct simulation based on electronic structure calculations of the surfaces and couplings is expensive due to the large number of electronic structure calculations needed for ensemble averaging or simulation of longer-time processes. An alternative approach is to construct an analytical representation of potential energy surfaces (PESs) and couplings, which allows for faster dynamics calculations. Diabatic representations are preferred for such purposes because of the smoothness of the surfaces and couplings and the scalar nature of the couplings. However, many diabatization procedures are complicated by the need to consider orbitals or vector coupling elements, and these can make the process very labor-intensive. To circumvent these difficulties, we here propose diabatization by a deep neural network (DDNN) based on a new architecture for a deep neural network that requires neither orbital input nor vector input. The DDNN method allows convenient and semiautomatic diabatization, and it is demonstrated here for a model problem and for producing diabatic potential energy matrices for thiophenol.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
30
|
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]
|
31
|
Hong Y, Yin Z, Guan Y, Zhang Z, Fu B, Zhang DH. Exclusive Neural Network Representation of the Quasi-Diabatic Hamiltonians Including Conical Intersections. J Phys Chem Lett 2020; 11:7552-7558. [PMID: 32835486 DOI: 10.1021/acs.jpclett.0c02173] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
We propose a numerically simple and straightforward, yet accurate and efficient neural networks-based fitting strategy to construct coupled potential energy surfaces (PESs) in a quasi-diabatic representation. The fundamental invariants are incorporated to account for the complete nuclear permutation inversion symmetry. Instead of derivative couplings or interstate couplings, a so-called modified derivative coupling term is fitted by neural networks, resulting in accurate description of near degeneracy points, such as the conical intersections. The adiabatic energies, energy gradients, and derivative couplings are well reproduced, and the vanishing of derivative couplings as well as the isotropic topography of adiabatic and diabatic energies in asymptotic regions are automatically satisfied. All of these features of the coupled global PESs are requisite for accurate dynamics simulations. Our approach is expected to be very useful in developing highly accurate coupled PESs in a quasi-diabatic representation in an efficient machine learning-based way.
Collapse
Affiliation(s)
- Yingyue Hong
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
- University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
| | - Zhengxi Yin
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
- University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
| | - Yafu Guan
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Zhaojun Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
| | - Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
| | - Dong H Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian, P.R. China 116023
| |
Collapse
|
32
|
Conti I, Cerullo G, Nenov A, Garavelli M. Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going. J Am Chem Soc 2020; 142:16117-16139. [PMID: 32841559 PMCID: PMC7901644 DOI: 10.1021/jacs.0c04952] [Citation(s) in RCA: 32] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Computational spectroscopy is becoming a mandatory tool for the interpretation of the
complex, and often congested, spectral maps delivered by modern non-linear multi-pulse
techniques. The fields of Electronic Structure Methods,
Non-Adiabatic Molecular Dynamics, and Theoretical
Spectroscopy represent the three pillars of the virtual ultrafast
optical spectrometer, able to deliver transient spectra in
silico from first principles. A successful simulation strategy requires a
synergistic approach that balances between the three fields, each one having its very
own challenges and bottlenecks. The aim of this Perspective is to demonstrate that,
despite these challenges, an impressive agreement between theory and experiment is
achievable now regarding the modeling of ultrafast photoinduced processes in complex
molecular architectures. Beyond that, some key recent developments in the three fields
are presented that we believe will have major impacts on spectroscopic simulations in
the very near future. Potential directions of development, pending challenges, and
rising opportunities are illustrated.
Collapse
Affiliation(s)
- Irene Conti
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Giulio Cerullo
- Dipartimento di Fisica, Politecnico di Milano, IFN-CNR, Piazza Leonardo da Vinci 32, I-20133 Milano, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| |
Collapse
|
33
|
Williams DMG, Eisfeld W. Complete Nuclear Permutation Inversion Invariant Artificial Neural Network (CNPI-ANN) Diabatization for the Accurate Treatment of Vibronic Coupling Problems. J Phys Chem A 2020; 124:7608-7621. [DOI: 10.1021/acs.jpca.0c05991] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David M. G. Williams
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
34
|
Shen Y, Yarkony DR. Construction of Quasi-diabatic Hamiltonians That Accurately Represent ab Initio Determined Adiabatic Electronic States Coupled by Conical Intersections for Systems on the Order of 15 Atoms. Application to Cyclopentoxide Photoelectron Detachment in the Full 39 Degrees of Freedom. J Phys Chem A 2020; 124:4539-4548. [PMID: 32374600 DOI: 10.1021/acs.jpca.0c02763] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present, for systems of moderate dimension, a fitting framework to construct quasi-diabatic Hamiltonians that accurately represent ab initio adiabatic electronic structure data including the effects of conical intersections. The framework introduced here minimizes the difference between the fit prediction and the ab initio data obtained in the adiabatic representation, which is singular at a conical intersection seam. We define a general and flexible merit function to allow arbitrary representations and propose a representation to measure the fit-ab initio difference at geometries near electronic degeneracies. A fit Hamiltonian may behave poorly in insufficiently sampled regions, in which case a machine learning theory analysis of the fit representation suggests a regularization to address the deficiency. Our fitting framework including the regularization is used to construct the full 39-dimensional coupled diabatic potential energy surfaces for cyclopentoxy relevant to cyclopentoxide photoelectron detachment.
Collapse
Affiliation(s)
- Yifan Shen
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
35
|
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
|
36
|
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
|
37
|
Jin Z, Dou W, Subotnik JE. Configuration interaction approaches for solving quantum impurity models. J Chem Phys 2020; 152:064105. [PMID: 32061216 DOI: 10.1063/1.5131624] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We develop several configuration interaction approaches for characterizing the electronic structure of an adsorbate on a metal surface (at least in model form). When one can separate the adsorbate from the substrate, these methods can achieve a reasonable description of adsorbate on-site electron-electron correlation in the presence of a continuum of states. While the present paper is restricted to the Anderson impurity model, there is hope that these methods can be extended to ab initio Hamiltonians and provide insight into the structure and dynamics of molecule-metal surface interactions.
Collapse
Affiliation(s)
- Zuxin Jin
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Wenjie Dou
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
38
|
Guan Y, Guo H, Yarkony DR. Extending the Representation of Multistate Coupled Potential Energy Surfaces To Include Properties Operators Using Neural Networks: Application to the 1,21A States of Ammonia. J Chem Theory Comput 2019; 16:302-313. [DOI: 10.1021/acs.jctc.9b00898] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
39
|
Sapunar M, Domcke W, Došlić N. UV absorption spectra of DNA bases in the 350-190 nm range: assignment and state specific analysis of solvation effects. Phys Chem Chem Phys 2019; 21:22782-22793. [PMID: 31595896 DOI: 10.1039/c9cp04662c] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.
Collapse
Affiliation(s)
- Marin Sapunar
- Department of Physical Chemistry, Ruder Bošković Institute, 10000 Zagreb, Croatia.
| | | | | |
Collapse
|
40
|
Williams DMG, Viel A, Eisfeld W. Diabatic neural network potentials for accurate vibronic quantum dynamics—The test case of planar NO3. J Chem Phys 2019; 151:164118. [DOI: 10.1063/1.5125851] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David M. G. Williams
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| | - Alexandra Viel
- Univ Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, F-35000 Rennes, France
| | - Wolfgang Eisfeld
- Theoretische Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany
| |
Collapse
|
41
|
Zhang L, Truhlar DG, Sun S. Full-dimensional three-state potential energy surfaces and state couplings for photodissociation of thiophenol. J Chem Phys 2019; 151:154306. [DOI: 10.1063/1.5124870] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People’s Republic of China
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People’s Republic of China
| |
Collapse
|
42
|
Shu Y, Kryven J, Sampaio de Oliveira-Filho AG, Zhang L, Song GL, Li SL, Meana-Pañeda R, Fu B, Bowman JM, Truhlar DG. Direct diabatization and analytic representation of coupled potential energy surfaces and couplings for the reactive quenching of the excited 2Σ+ state of OH by molecular hydrogen. J Chem Phys 2019; 151:104311. [DOI: 10.1063/1.5111547] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Joanna Kryven
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Antonio Gustavo Sampaio de Oliveira-Filho
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
- Departamento de Química, Laboratório Computacional de Espectroscopia e Cinética, Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto, Universidade de São Paulo, 14040-901 Ribeirão Preto-SP, Brazil
| | - Linyao Zhang
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Guo-Liang Song
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Shaohong L. Li
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Rubén Meana-Pañeda
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Bina Fu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Joel M. Bowman
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| |
Collapse
|
43
|
Wang L, Qiu J, Bai X, Xu J. Surface hopping methods for nonadiabatic dynamics in extended systems. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1435] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Linjun Wang
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Jing Qiu
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Xin Bai
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Jiabo Xu
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| |
Collapse
|
44
|
Reddy SR, Coto PB, Thoss M. Quantum dynamical simulation of intramolecular singlet fission in covalently coupled pentacene dimers. J Chem Phys 2019; 151:044307. [DOI: 10.1063/1.5109897] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- S. Rajagopala Reddy
- Institute of Physics, Albert-Ludwigs University Freiburg, Hermann-Herder-Str. 3, 79104 Freiburg, Germany
| | - Pedro B. Coto
- Department of Physical and Analytical Chemistry, University of Oviedo, Avda. Julián Clavería 8, 33006, Oviedo, Spain
| | - Michael Thoss
- Institute of Physics, Albert-Ludwigs University Freiburg, Hermann-Herder-Str. 3, 79104 Freiburg, Germany
| |
Collapse
|
45
|
Guan Y, Guo H, Yarkony DR. Neural network based quasi-diabatic Hamiltonians with symmetry adaptation and a correct description of conical intersections. J Chem Phys 2019; 150:214101. [DOI: 10.1063/1.5099106] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yafu Guan
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - David R. Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
46
|
Roulet J, Choi S, Vaníček J. Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation. J Chem Phys 2019; 150:204113. [PMID: 31153180 DOI: 10.1063/1.5094046] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Exact nonadiabatic quantum evolution preserves many geometric properties of the molecular Hilbert space. In the first paper of this series ["Paper I," S. Choi and J. Vaníček, J. Chem. Phys. 150, 204112 (2019)], we presented numerical integrators of arbitrary-order of accuracy that preserve these geometric properties exactly even in the adiabatic representation, in which the molecular Hamiltonian is not separable into kinetic and potential terms. Here, we focus on the separable Hamiltonian in diabatic representation, where the split-operator algorithm provides a popular alternative because it is explicit and easy to implement, while preserving most geometric invariants. Whereas the standard version has only second-order accuracy, we implemented, in an automated fashion, its recursive symmetric compositions, using the same schemes as in Paper I, and obtained integrators of arbitrary even order that still preserve the geometric properties exactly. Because the automatically generated splitting coefficients are redundant, we reduce the computational cost by pruning these coefficients and lower memory requirements by identifying unique coefficients. The order of convergence and preservation of geometric properties are justified analytically and confirmed numerically on a one-dimensional two-surface model of NaI and a three-dimensional three-surface model of pyrazine. As for efficiency, we find that to reach a convergence error of 10-10, a 600-fold speedup in the case of NaI and a 900-fold speedup in the case of pyrazine are obtained with the higher-order compositions instead of the second-order split-operator algorithm. The pyrazine results suggest that the efficiency gain survives in higher dimensions.
Collapse
Affiliation(s)
- Julien Roulet
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| | - Seonghoon Choi
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| | - Jiří Vaníček
- Laboratory of Theoretical Physical Chemistry, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| |
Collapse
|
47
|
Sapunar M, Piteša T, Davidović D, Došlić N. Highly Efficient Algorithms for CIS Type Excited State Wave Function Overlaps. J Chem Theory Comput 2019; 15:3461-3469. [DOI: 10.1021/acs.jctc.9b00235] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Marin Sapunar
- Centre for Informatics and Computing, Rudjer Bošković Institute, Zagreb 10000, Croatia
| | - Tomislav Piteša
- Centre for Informatics and Computing, Rudjer Bošković Institute, Zagreb 10000, Croatia
| | - Davor Davidović
- Centre for Informatics and Computing, Rudjer Bošković Institute, Zagreb 10000, Croatia
| | - Nadja Došlić
- Centre for Informatics and Computing, Rudjer Bošković Institute, Zagreb 10000, Croatia
| |
Collapse
|
48
|
Yarkony DR, Xie C, Zhu X, Wang Y, Malbon CL, Guo H. Diabatic and adiabatic representations: Electronic structure caveats. COMPUT THEOR CHEM 2019. [DOI: 10.1016/j.comptc.2019.01.020] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
49
|
Seidu I, Goel P, Wang XG, Chen B, Wang XB, Zeng T. Vibronic interaction in CO3− photo-detachment: Jahn–Teller effects beyond structural distortion and general formalisms for vibronic Hamiltonians in trigonal symmetries. Phys Chem Chem Phys 2019; 21:8679-8690. [DOI: 10.1039/c9cp01352k] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Expansion formalisms for trigonal Jahn–Teller and pseudo-Jahn–Teller vibronic Hamiltonians are developed and used to study and correctly interpret the photoelectron spectrum of CO3−.
Collapse
Affiliation(s)
- Issaka Seidu
- Department of Chemistry
- Carleton University
- Ottawa
- Canada
- Department of Chemistry and Biomolecular Sciences
| | - Prateek Goel
- Department of Chemistry
- University of Florida
- Gainesville
- USA
| | | | - Bo Chen
- Department of Chemistry
- Pennsylvania State University
- State College
- USA
| | - Xue-Bin Wang
- Pacific Northwest National Laboratory
- Richland
- USA
| | - Tao Zeng
- Department of Chemistry
- Carleton University
- Ottawa
- Canada
| |
Collapse
|
50
|
Robertson C, González-Vázquez J, Corral I, Díaz-Tendero S, Díaz C. Nonadiabatic scattering of NO off Au3
clusters: A simple and robust diabatic state manifold generation method for multiconfigurational wavefunctions. J Comput Chem 2018; 40:794-810. [DOI: 10.1002/jcc.25764] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2018] [Revised: 10/28/2018] [Accepted: 11/01/2018] [Indexed: 01/19/2023]
Affiliation(s)
- Christopher Robertson
- Department of Chemistry and Centre for Scientific Computing; University Of Warwick, CV4 7AL; Coventry United Kingdom
| | - Jesús González-Vázquez
- Departamento de Química Módulo 13; Universidad Autónoma de Madrid, 28049; Madrid Spain
- Institute for Advanced Research in Chemistry (IAdChem); Universidad Autónoma de Madrid, 28049; Madrid Spain
| | - Ines Corral
- Institute for Advanced Research in Chemistry (IAdChem); Universidad Autónoma de Madrid, 28049; Madrid Spain
- Departamento de Química Módulo 13; Universidad Autónoma de Madrid, 28049; Madrid Spain
| | - Sergio Díaz-Tendero
- Condensed Matter Physics Center (IFIMAC); Universidad Autónoma de Madrid, 28049; Madrid Spain
- Departamento de Química Módulo 13; Universidad Autónoma de Madrid, 28049; Madrid Spain
- Institute for Advanced Research in Chemistry (IAdChem); Universidad Autónoma de Madrid, 28049; Madrid Spain
| | - Cristina Díaz
- Departamento de Química Módulo 13; Universidad Autónoma de Madrid, 28049; Madrid Spain
- Institute for Advanced Research in Chemistry (IAdChem); Universidad Autónoma de Madrid, 28049; Madrid Spain
- Condensed Matter Physics Center (IFIMAC); Universidad Autónoma de Madrid, 28049; Madrid Spain
| |
Collapse
|