51
|
Hanasaki K, Kanno M, Niehaus TA, Kono H. An efficient approximate algorithm for nonadiabatic molecular dynamics. J Chem Phys 2019; 149:244117. [PMID: 30599729 DOI: 10.1063/1.5046757] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
We propose a modification to the nonadiabatic surface hopping calculation method formulated in a paper by Yu et al. [Phys. Chem. Chem. Phys. 16, 25883 (2014)], which is a multidimensional extension of the Zhu-Nakamura theory with a practical diabatic gradient estimation algorithm. In our modification, their diabatic gradient estimation algorithm, which is based on a simple interpolation of the adiabatic potential energy surfaces, is replaced by an algorithm using the numerical derivatives of the adiabatic gradients. We then apply the algorithm to several models of nonadiabatic dynamics, both analytic and ab initio models, to numerically demonstrate that our method indeed widens the applicability and robustness of their method. We also discuss the validity and limitations of our new nonadiabatic surface hopping method while considering in mind potential applications to excited-state dynamics of biomolecules or unconventional nonadiabatic dynamics such as radiation decay processes in ultraintense X-ray fields.
Collapse
Affiliation(s)
- Kota Hanasaki
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | - Manabu Kanno
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| | - Thomas A Niehaus
- University of Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeubanne, France
| | - Hirohiko Kono
- Department of Chemistry, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan
| |
Collapse
|
52
|
Parker SM, Roy S, Furche F. Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation. Phys Chem Chem Phys 2019; 21:18999-19010. [DOI: 10.1039/c9cp03127h] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions.
Collapse
Affiliation(s)
- Shane M. Parker
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| | - Saswata Roy
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| | - Filipp Furche
- Department of Chemistry
- University of California, Irvine
- Irvine
- USA
| |
Collapse
|
53
|
Zhang X, Lu G. First-order nonadiabatic couplings in extended systems by time-dependent density functional theory. J Chem Phys 2018; 149:244103. [DOI: 10.1063/1.5065504] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xu Zhang
- Department of Physics and Astronomy, California State University Northridge, Northridge, California 91330, USA
| | - Gang Lu
- Department of Physics and Astronomy, California State University Northridge, Northridge, California 91330, USA
| |
Collapse
|
54
|
Liu M, Chen X, Grofe A, Gao J. Diabatic States at Construction (DAC) through Generalized Singular Value Decomposition. J Phys Chem Lett 2018; 9:6038-6046. [PMID: 30277783 DOI: 10.1021/acs.jpclett.8b02472] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
A procedure, called generalized diabatic-at-construction (GDAC), is presented to transform adiabatic potential energy surfaces into a diabatic representation by generalized singular value decomposition. First, we use a set of localized, valence bond-like configuration state functions, called DAC, as the basis states. Then, the adiabatic ground and relevant excited states are determined using multistate density functional theory (MSDFT). GDAC differs in the opposite direction from traditional approaches based on adiabatic-to-diabatic transformation with certain property restraints. The method is illustrated with applications to a model first-order bond dissociation reaction of CH3OCH2Cl polarized by a solvent molecule, the ground- and first-excited-state potential energy surfaces near the minimum conical intersection for the ammonia dimer photodissociation, and the multiple avoided curve crossings in the dissociation of lithium hydride. The GDAC diabatization method may be useful for defining charge-localized states in studies of electron transfer and proton-coupled electron transfer reactions in proteins.
Collapse
Affiliation(s)
- Meiyi Liu
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry , Jilin University , Changchun , Jilin Province 130023 , China
| | - Xin Chen
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry , Jilin University , Changchun , Jilin Province 130023 , China
| | - Adam Grofe
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry , Jilin University , Changchun , Jilin Province 130023 , China
| | - Jiali Gao
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry , Jilin University , Changchun , Jilin Province 130023 , China
- Department of Chemistry and Supercomputing Institute , University of Minnesota , Minneapolis , Minnesota 55455 , United States
| |
Collapse
|
55
|
Véril M, Romaniello P, Berger JA, Loos PF. Unphysical Discontinuities in GW Methods. J Chem Theory Comput 2018; 14:5220-5228. [DOI: 10.1021/acs.jctc.8b00745] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
56
|
COBRAMM 2.0 — A software interface for tailoring molecular electronic structure calculations and running nanoscale (QM/MM) simulations. J Mol Model 2018; 24:271. [DOI: 10.1007/s00894-018-3769-6] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2018] [Accepted: 07/27/2018] [Indexed: 01/04/2023]
|
57
|
Mandal A, Shakib FA, Huo P. Investigating photoinduced proton coupled electron transfer reaction using quasi diabatic dynamics propagation. J Chem Phys 2018; 148:244102. [DOI: 10.1063/1.5030634] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Arkajit Mandal
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Farnaz A. Shakib
- 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
|
58
|
Zhang X. Large-scaleab initiocalculations of Raman scattering spectra within time-dependent density functional perturbation theory. J Chem Phys 2018; 148:244103. [DOI: 10.1063/1.5038112] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xu Zhang
- Department of Physics and Astronomy, California State University Northridge, Northridge, California 91330, USA
| |
Collapse
|
59
|
Gonon B, Perveaux A, Gatti F, Lauvergnat D, Lasorne B. On the applicability of a wavefunction-free, energy-based procedure for generating first-order non-adiabatic couplings around conical intersections. J Chem Phys 2018; 147:114114. [PMID: 28938825 DOI: 10.1063/1.4991635] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The primal definition of first-order non-adiabatic couplings among electronic states relies on the knowledge of how electronic wavefunctions vary with nuclear coordinates. However, the non-adiabatic coupling between two electronic states can be obtained in the vicinity of a conical intersection from energies only, as this vector spans the branching plane along which degeneracy is lifted to first order. The gradient difference and derivative coupling are responsible of the two-dimensional cusp of a conical intersection between both potential-energy surfaces and can be identified to the non-trivial eigenvectors of the second derivative of the square energy difference, as first pointed out in Köppel and Schubert [Mol. Phys. 104(5-7), 1069 (2006)]. Such quantities can always be computed in principle for the cost of two numerical Hessians in the worst-case scenario. Analytic-derivative techniques may help in terms of accuracy and efficiency but also raise potential traps due to singularities and ill-defined derivatives at degeneracies. We compare here two approaches, one fully numerical, the other semianalytic, where analytic gradients are available but Hessians are not, and investigate their respective conditions of applicability. Benzene and 3-hydroxychromone are used as illustrative application cases. It is shown that non-adiabatic couplings can thus be estimated with decent accuracy in regions of significant size around conical intersections. This procedure is robust and could be useful in the context of on-the-fly non-adiabatic dynamics or be used for producing model representations of intersecting potential energy surfaces with complete obviation of the electronic wavefunctions.
Collapse
Affiliation(s)
- Benjamin Gonon
- Institut Charles Gerhardt Montpellier (UMR 5253), CNRS, Université de Montpellier, F-34095 Montpellier, France
| | - Aurelie Perveaux
- Institut Charles Gerhardt Montpellier (UMR 5253), CNRS, Université de Montpellier, F-34095 Montpellier, France
| | - Fabien Gatti
- Institut Charles Gerhardt Montpellier (UMR 5253), CNRS, Université de Montpellier, F-34095 Montpellier, France
| | - David Lauvergnat
- Laboratoire de Chimie Physique (UMR 8000), CNRS, Université Paris-Sud/Paris-Saclay, F-91405 Orsay, France
| | - Benjamin Lasorne
- Institut Charles Gerhardt Montpellier (UMR 5253), CNRS, Université de Montpellier, F-34095 Montpellier, France
| |
Collapse
|
60
|
Crespo-Otero R, Barbatti M. Recent Advances and Perspectives on Nonadiabatic Mixed Quantum–Classical Dynamics. Chem Rev 2018; 118:7026-7068. [DOI: 10.1021/acs.chemrev.7b00577] [Citation(s) in RCA: 301] [Impact Index Per Article: 50.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rachel Crespo-Otero
- School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | | |
Collapse
|
61
|
Yang CH, Yam C, Wang H. Approximate DFT-based methods for generating diabatic states and calculating electronic couplings: models of two and more states. Phys Chem Chem Phys 2018; 20:2571-2584. [PMID: 29318238 DOI: 10.1039/c7cp06660k] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Four types of density functional theory (DFT)-based approaches are assessed in this work for the approximate construction of diabatic states and the evaluation of electronic couplings between these states. These approaches include the constrained DFT (CDFT) method, the constrained noninteracting electron (CNE) model to post-process Kohn-Sham operators, the approximate block-diagonalization (BD) of the Kohn-Sham operators, and the generalized Mulliken-Hush method. It is shown that the first three approaches provide a good description for long-distance intermolecular electron transfer (ET) reactions. On the other hand, inconsistent results were found when applying these approaches to intramolecular ET in strongly coupled, mixed-valence systems. Model analysis shows that this discrepancy is caused by the inappropriate use of the two-state model rather than the defects of the approaches themselves. The situation is much improved when more states are included in the model electronic Hamiltonian. The CNE and BD approaches can thus serve as efficient and robust alternatives for building ET models based on DFT calculations.
Collapse
Affiliation(s)
- Chou-Hsun Yang
- Beijing Computational Science Research Center, Haidian District, Beijing 100193, China
| | | | | |
Collapse
|
62
|
Abstract
The foundations, formalisms, technicalities, and practicalities of relativistic time-dependent density functional theories (R-TD-DFT) for spinor excited states of molecular systems containing heavy elements are critically reviewed.
Collapse
Affiliation(s)
- Wenjian Liu
- Beijing National Center for Molecular Sciences
- Institute of Theoretical and Computational Chemistry
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
| | - Yunlong Xiao
- Beijing National Center for Molecular Sciences
- Institute of Theoretical and Computational Chemistry
- College of Chemistry and Molecular Engineering
- Peking University
- Beijing 100871
| |
Collapse
|
63
|
Maitra NT. Charge transfer in time-dependent density functional theory. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:423001. [PMID: 28766507 DOI: 10.1088/1361-648x/aa836e] [Citation(s) in RCA: 58] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Charge transfer plays a crucial role in many processes of interest in physics, chemistry, and bio-chemistry. In many applications the size of the systems involved calls for time-dependent density functional theory (TDDFT) to be used in their computational modeling, due to its unprecedented balance between accuracy and efficiency. However, although exact in principle, in practise approximations must be made for the exchange-correlation functional in this theory, and the standard functional approximations perform poorly for excitations which have a long-range charge-transfer component. Intense progress has been made in developing more sophisticated functionals for this problem, which we review. We point out an essential difference between the properties of the exchange-correlation kernel needed for an accurate description of charge-transfer between open-shell fragments and between closed-shell fragments. We then turn to charge-transfer dynamics, which, in contrast to the excitation problem, is a highly non-equilibrium, non-perturbative, process involving a transfer of one full electron in space. This turns out to be a much more challenging problem for TDDFT functionals. We describe dynamical step and peak features in the exact functional evolving over time, that are missing in the functionals currently used. The latter underestimate the amount of charge transferred and manifest a spurious shift in the charge transfer resonance position. We discuss some explicit examples.
Collapse
Affiliation(s)
- Neepa T Maitra
- Department of Physics and Astronomy, Hunter College and the Physics Program at the Graduate Center of the City University of New York, 695 Park Avenue, NY 10065, United States of America
| |
Collapse
|
64
|
Morrison AF, Herbert JM. Analytic derivative couplings and first-principles exciton/phonon coupling constants for an ab initio Frenkel-Davydov exciton model: Theory, implementation, and application to compute triplet exciton mobility parameters for crystalline tetracene. J Chem Phys 2017; 146:224110. [DOI: 10.1063/1.4985607] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adrian F. Morrison
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
65
|
Petrone A, Williams-Young DB, Lingerfelt DB, Li X. Ab Initio Excited-State Transient Raman Analysis. J Phys Chem A 2017; 121:3958-3965. [DOI: 10.1021/acs.jpca.7b02905] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Alessio Petrone
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | | | - David B. Lingerfelt
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
66
|
Grofe A, Qu Z, Truhlar DG, Li H, Gao J. Diabatic-At-Construction Method for Diabatic and Adiabatic Ground and Excited States Based on Multistate Density Functional Theory. J Chem Theory Comput 2017; 13:1176-1187. [PMID: 28135420 PMCID: PMC5793876 DOI: 10.1021/acs.jctc.6b01176] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We describe a diabatic-at-construction (DAC) strategy for defining diabatic states to determine the adiabatic ground and excited electronic states and their potential energy surfaces using the multistate density functional theory (MSDFT). The DAC approach differs in two fundamental ways from the adiabatic-to-diabatic (ATD) procedures that transform a set of preselected adiabatic electronic states to a new representation. (1) The DAC states are defined in the first computation step to form an active space, whose configuration interaction produces the adiabatic ground and excited states in the second step of MSDFT. Thus, they do not result from a similarity transformation of the adiabatic states as in the ATD procedure; they are the basis for producing the adiabatic states. The appropriateness and completeness of the DAC active space can be validated by comparison with experimental observables of the ground and excited states. (2) The DAC diabatic states are defined using the valence bond characters of the asymptotic dissociation limits of the adiabatic states of interest, and they are strictly maintained at all molecular geometries. Consequently, DAC diabatic states have specific and well-defined physical and chemical meanings that can be used for understanding the nature of the adiabatic states and their energetic components. Here we present results for the four lowest singlet states of LiH and compare them to a well-tested ATD diabatization method, namely the 3-fold way; the comparison reveals both similarities and differences between the ATD diabatic states and the orthogonalized DAC diabatic states. Furthermore, MSDFT can provide a quantitative description of the ground and excited states for LiH with multiple strongly and weakly avoided curve crossings spanning over 10 Å of interatomic separation.
Collapse
Affiliation(s)
- Adam Grofe
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Zexing Qu
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Hui Li
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
| | - Jiali Gao
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
67
|
Gao X, Bai S, Fazzi D, Niehaus T, Barbatti M, Thiel W. Evaluation of Spin-Orbit Couplings with Linear-Response Time-Dependent Density Functional Methods. J Chem Theory Comput 2017; 13:515-524. [DOI: 10.1021/acs.jctc.6b00915] [Citation(s) in RCA: 171] [Impact Index Per Article: 24.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Xing Gao
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz, D-45470, Mülheim an der Ruhr, Germany
| | - Shuming Bai
- Aix Marseille Univ, CNRS, ICR, Marseille, France
| | - Daniele Fazzi
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz, D-45470, Mülheim an der Ruhr, Germany
| | - Thomas Niehaus
- Univ Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, F-69622, Villeurbanne, France
| | | | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz, D-45470, Mülheim an der Ruhr, Germany
| |
Collapse
|
68
|
Cisneros C, Thompson T, Baluyot N, Smith AC, Tapavicza E. The role of tachysterol in vitamin D photosynthesis – a non-adiabatic molecular dynamics study. Phys Chem Chem Phys 2017; 19:5763-5777. [DOI: 10.1039/c6cp08064b] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
To investigate the role of tachysterol in the regulation of vitamin D photosynthesis, we studied its absorption properties and photodynamics by ab initio methods and non-adiabatic molecular dynamics.
Collapse
Affiliation(s)
- Cecilia Cisneros
- Department of Chemistry and Biochemistry
- California State University
- Long Beach
- Long Beach
- USA
| | - Travis Thompson
- Department of Chemistry and Biochemistry
- California State University
- Long Beach
- Long Beach
- USA
| | - Noel Baluyot
- Department of Chemistry and Biochemistry
- California State University
- Long Beach
- Long Beach
- USA
| | - Adam C. Smith
- Department of Chemistry and Biochemistry
- California State University
- Long Beach
- Long Beach
- USA
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry
- California State University
- Long Beach
- Long Beach
- USA
| |
Collapse
|
69
|
Petrone A, Lingerfelt DB, Williams-Young DB, Li X. Ab Initio Transient Vibrational Spectral Analysis. J Phys Chem Lett 2016; 7:4501-4508. [PMID: 27788583 DOI: 10.1021/acs.jpclett.6b02292] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Pump probe spectroscopy techniques have enabled the direct observation of a variety of transient molecular species in both ground and excited electronic states. Time-resolved vibrational spectroscopy is becoming an indispensable tool for investigating photoinduced nuclear dynamics of chemical systems of all kinds. On the other hand, a complete picture of the chemical dynamics encoded in these spectra cannot be achieved without a full temporal description of the structural relaxation, including the explicit time-dependence of vibrational coordinates that are substantially displaced from equilibrium by electronic excitation. Here we present a transient vibrational analysis protocol combining ab initio direct molecular dynamics and time-integrated normal modes introduced in this work, relying on the recent development of analytic time-dependent density functional theory (TDDFT) second derivatives for excited states. Prototypical molecules will be used as test cases, showing the evolution of the vibrational signatures that follow electronic excitation. This protocol provides a direct route to assigning the vibrations implicated in the (photo)dynamics of several (photoactive) systems.
Collapse
Affiliation(s)
- Alessio Petrone
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - David B Lingerfelt
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - David B Williams-Young
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington , Seattle, Washington 98195, United States
| |
Collapse
|
70
|
Vincent JC, Muuronen M, Pearce KC, Mohanam LN, Tapavicza E, Furche F. That Little Extra Kick: Nonadiabatic Effects in Acetaldehyde Photodissociation. J Phys Chem Lett 2016; 7:4185-4190. [PMID: 27704839 DOI: 10.1021/acs.jpclett.6b02037] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The effect of nonadiabatic transitions on branching ratios, kinetic and internal energy distribution of fragments, and reaction mechanisms observed in acetaldehyde photodissociation is investigated by nonadiabatic molecular dynamics (NAMD) simulations using time-dependent hybrid density functional theory and Tully surface hopping. Homolytic bond breaking is approximately captured by allowing spin symmetry to break. The NAMD simulations reveal that nonadiabatic transitions selectively enhance the kinetic energy of certain internal degrees of freedom within approximately 50 fs. Branching ratios from NAMD and conventional "hot" Born-Oppenheimer molecular dynamics (BOMD) are similar and qualitatively agree with experiment. However, as opposed to the BOMD simulations, NAMD captures the high-energy tail of the experimental kinetic energy distribution. The extra "kick" of the nuclei in the direction of the nonadiabatic coupling vector results from nonadiabatic transitions close to conical intersections. From a mechanistic perspective, the nonadiabatic effects favor asynchronous over synchronous fragmentation and tend to suppress roaming.
Collapse
Affiliation(s)
- Jordan C Vincent
- Department of Chemistry, University of California, Irvine , 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Mikko Muuronen
- Department of Chemistry, University of California, Irvine , 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Kirk C Pearce
- Department of Chemistry, University of California, Irvine , 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Luke N Mohanam
- Department of Chemistry, University of California, Irvine , 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach , 1250 Bellflower Boulevard, Long Beach, California 90840-9507, United States
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine , 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
71
|
Parker SM, Roy S, Furche F. Unphysical divergences in response theory. J Chem Phys 2016; 145:134105. [DOI: 10.1063/1.4963749] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Shane M. Parker
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Saswata Roy
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
72
|
Arbelo-González W, Crespo-Otero R, Barbatti M. Steady and Time-Resolved Photoelectron Spectra Based on Nuclear Ensembles. J Chem Theory Comput 2016; 12:5037-5049. [DOI: 10.1021/acs.jctc.6b00704] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Wilmer Arbelo-González
- Max-Planck-Institut für Kohlenforschung, Kaizer-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Rachel Crespo-Otero
- School
of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | - Mario Barbatti
- Max-Planck-Institut für Kohlenforschung, Kaizer-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
- Aix Marseille Univ, CNRS, ICR, Marseille, France
| |
Collapse
|
73
|
Maitra NT. Perspective: Fundamental aspects of time-dependent density functional theory. J Chem Phys 2016; 144:220901. [DOI: 10.1063/1.4953039] [Citation(s) in RCA: 228] [Impact Index Per Article: 28.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Neepa T. Maitra
- Department of Physics and Astronomy, Hunter College and the Physics Program at the Graduate Center of the City University of New York, 695 Park Avenue, New York, New York 10065, USA
| |
Collapse
|
74
|
Cao Y, Hughes T, Giesen D, Halls MD, Goldberg A, Vadicherla TR, Sastry M, Patel B, Sherman W, Weisman AL, Friesner RA. Highly efficient implementation of pseudospectral time-dependent density-functional theory for the calculation of excitation energies of large molecules. J Comput Chem 2016; 37:1425-41. [DOI: 10.1002/jcc.24350] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2015] [Revised: 02/02/2016] [Accepted: 02/03/2016] [Indexed: 12/25/2022]
Affiliation(s)
- Yixiang Cao
- Schrödinger Inc; 120 West 45th Street, Tower 45, 17th Floor New York New York 10036
| | - Thomas Hughes
- Schrödinger Inc; 120 West 45th Street, Tower 45, 17th Floor New York New York 10036
| | - Dave Giesen
- Schrödinger Inc; 120 West 45th Street, Tower 45, 17th Floor New York New York 10036
| | - Mathew D. Halls
- Schrödinger Inc; 120 West 45th Street, Tower 45, 17th Floor New York New York 10036
| | - Alexander Goldberg
- Schrödinger Inc; 8910 University Lane, Suite 270 San Diego California 92122
| | - Tati Reddy Vadicherla
- Schrödinger; Plot No. 573, B & C, Road No. 1 Jubilee Hills, Hyderabad Telangana 5000096 India
| | - Madhavi Sastry
- Schrödinger; Sanali Infopark, 8-2-120/113 Banjara Hills, Hyderabad Andhra Pradesh 500034 India
| | - Bhargav Patel
- Schrödinger; Sanali Infopark, 8-2-120/113 Banjara Hills, Hyderabad Andhra Pradesh 500034 India
| | - Woody Sherman
- Schrödinger Inc; 245 First St. Riverview II, 18th Floor Cambridge Massachusetts 02142
| | - Andrew L. Weisman
- Department of Applied Physics and Applied Mathematics; Columbia University; New York New York 10027
| | | |
Collapse
|
75
|
Hu Z, Autschbach J, Jensen L. Simulating Third-Order Nonlinear Optical Properties Using Damped Cubic Response Theory within Time-Dependent Density Functional Theory. J Chem Theory Comput 2016; 12:1294-304. [DOI: 10.1021/acs.jctc.5b01060] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Zhongwei Hu
- Department
of Chemistry, The Pennsylvania State University, 104 Chemistry Building, University Park, Pennsylvania 16802-4615, United States
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Lasse Jensen
- Department
of Chemistry, The Pennsylvania State University, 104 Chemistry Building, University Park, Pennsylvania 16802-4615, United States
| |
Collapse
|
76
|
Ryabinkin IG, Nagesh J, Izmaylov AF. Fast Numerical Evaluation of Time-Derivative Nonadiabatic Couplings for Mixed Quantum-Classical Methods. J Phys Chem Lett 2015; 6:4200-4203. [PMID: 26538034 DOI: 10.1021/acs.jpclett.5b02062] [Citation(s) in RCA: 69] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
We have developed a numerical differentiation scheme that eliminates evaluation of overlap determinants in calculating the time-derivative nonadiabatic couplings (TDNACs). Evaluation of these determinants was the bottleneck in previous implementations of mixed quantum-classical methods using numerical differentiation of electronic wave functions in the Slater determinant representation. The central idea of our approach is, first, to reduce the analytic time derivatives of Slater determinants to time derivatives of molecular orbitals and then to apply a finite-difference formula. Benchmark calculations prove the efficiency of the proposed scheme showing impressive several-order-of-magnitude speedups of the TDNAC calculation step for midsize molecules.
Collapse
Affiliation(s)
- Ilya G Ryabinkin
- Department of Physical and Environmental Sciences, University of Toronto Scarborough , Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto , Toronto, Ontario M5S 3H6, Canada
| | - Jayashree Nagesh
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto , Toronto, Ontario M5S 3H6, Canada
| | - Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough , Toronto, Ontario M1C 1A4, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto , Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
77
|
Fischer SA, Cramer CJ, Govind N. Excited State Absorption from Real-Time Time-Dependent Density Functional Theory. J Chem Theory Comput 2015; 11:4294-303. [DOI: 10.1021/acs.jctc.5b00473] [Citation(s) in RCA: 56] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Sean A. Fischer
- Environmental
Molecular Sciences Laboratory, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| | - Christopher J. Cramer
- Department
of Chemistry, Supercomputing Institute and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Niranjan Govind
- Environmental
Molecular Sciences Laboratory, Pacific Northwest National Laboratory, P.O. Box 999, Richland, Washington 99352, United States
| |
Collapse
|
78
|
Hammes-Schiffer S. Proton-Coupled Electron Transfer: Moving Together and Charging Forward. J Am Chem Soc 2015; 137:8860-71. [PMID: 26110700 PMCID: PMC4601483 DOI: 10.1021/jacs.5b04087] [Citation(s) in RCA: 291] [Impact Index Per Article: 32.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2015] [Indexed: 12/24/2022]
Abstract
Proton-coupled electron transfer (PCET) is ubiquitous throughout chemistry and biology. This Perspective discusses recent advances and current challenges in the field of PCET, with an emphasis on the role of theory and computation. The fundamental theoretical concepts are summarized, and expressions for rate constants and kinetic isotope effects are provided. Computational methods for calculating reduction potentials and pKa's for molecular electrocatalysts, as well as insights into linear correlations and non-innocent ligands, are also described. In addition, computational methods for simulating the nonadiabatic dynamics of photoexcited PCET are discussed. Representative applications to PCET in solution, proteins, electrochemistry, and photoinduced processes are presented, highlighting the interplay between theoretical and experimental studies. The current challenges and suggested future directions are outlined for each type of application, concluding with an overall view to the future.
Collapse
Affiliation(s)
- Sharon Hammes-Schiffer
- Department of Chemistry, University of
Illinois at Urbana−Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United
States
| |
Collapse
|
79
|
Subotnik JE, Alguire EC, Ou Q, Landry BR, Fatehi S. The requisite electronic structure theory to describe photoexcited nonadiabatic dynamics: nonadiabatic derivative couplings and diabatic electronic couplings. Acc Chem Res 2015; 48:1340-50. [PMID: 25932499 DOI: 10.1021/acs.accounts.5b00026] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Electronically photoexcited dynamics are complicated because there are so many different relaxation pathways: fluorescence, phosphorescence, radiationless decay, electon transfer, etc. In practice, to model photoexcited systems is a very difficult enterprise, requiring accurate and very efficient tools in both electronic structure theory and nonadiabatic chemical dynamics. Moreover, these theoretical tools are not traditional tools. On the one hand, the electronic structure tools involve couplings between electonic states (rather than typical single state energies and gradients). On the other hand, the dynamics tools involve propagating nuclei on multiple potential energy surfaces (rather than the usual ground state dynamics). In this Account, we review recent developments in electronic structure theory as directly applicable for modeling photoexcited systems. In particular, we focus on how one may evaluate the couplings between two different electronic states. These couplings come in two flavors. If we order states energetically, the resulting adiabatic states are coupled via derivative couplings. Derivative couplings capture how electronic wave functions change as a function of nuclear geometry and can usually be calculated with straightforward tools from analytic gradient theory. One nuance arises, however, in the context of time-dependent density functional theory (TD-DFT): how do we evaluate derivative couplings between TD-DFT excited states (which are tricky, because no wave function is available)? This conundrum was recently solved, and we review the solution below. We also discuss the solution to a second, pesky problem of origin dependence, whereby the derivative couplings do not (strictly) satisfy translation variance, which can lead to a lack of momentum conservation. Apart from adiabatic states, if we order states according to their electronic character, the resulting diabatic states are coupled via electronic or diabatic couplings. The couplings between diabatic states |ΞA⟩ and |ΞB⟩ are just the simple matrix elements, ⟨ΞA|H|ΞB⟩. A difficulty arises, however, because constructing exactly diabatic states is formally impossible and constructing quasi-diabatic states is not unique. To that end, we review recent advances in localized diabatization, which is one approach for generating adiabatic-to-diabatic (ATD) transformations. We also highlight outstanding questions in the arena of diabatization, especially how to generate multiple globally stable diabatic surfaces.
Collapse
Affiliation(s)
- Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, 231
South 34th Street, Philadelphia, Pennsylvania 19104, United States
| | - Ethan C. Alguire
- Department of Chemistry, University of Pennsylvania, 231
South 34th Street, Philadelphia, Pennsylvania 19104, United States
| | - Qi Ou
- Department of Chemistry, University of Pennsylvania, 231
South 34th Street, Philadelphia, Pennsylvania 19104, United States
| | - Brian R. Landry
- Department of Chemistry, University of Pennsylvania, 231
South 34th Street, Philadelphia, Pennsylvania 19104, United States
| | - Shervin Fatehi
- Department of Chemistry, University of Utah, 315 South 1400
East, Room 2020, Salt Lake
City, Utah 84112, United States
| |
Collapse
|