1
|
Song H, Freixas VM, Fernandez-Alberti S, White AJ, Zhang Y, Mukamel S, Govind N, Tretiak S. An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem. J Chem Theory Comput 2021; 17:3629-3643. [PMID: 34014085 DOI: 10.1021/acs.jctc.1c00131] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The recently developed ab initio multiple cloning (AIMC) approach based on the multiconfigurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to nonadiabatic dynamics with a particular strength in the proper description of decoherence effects because of the branching of vibrational wavepackets at a level crossing. Here, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-field theory to determine the equations of motion for classical trajectories. The multidimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based nonadiabatic coupling approach within the AIMC method, all observables are calculated on-the-fly in the nonadiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level nonadiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.
Collapse
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD, Bernal, Argentina
| | | | - Alexander J White
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Shaul Mukamel
- Departments of Chemistry, Physics, and Astronomy, University of California, Irvine, California 92697, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States.,Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
2
|
Mai S, González L. Molecular Photochemistry: Recent Developments in Theory. Angew Chem Int Ed Engl 2020; 59:16832-16846. [PMID: 32052547 PMCID: PMC7540682 DOI: 10.1002/anie.201916381] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Revised: 02/12/2020] [Indexed: 12/16/2022]
Abstract
Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.
Collapse
Affiliation(s)
- Sebastian Mai
- Photonics InstituteVienna University of TechnologyGusshausstrasse 27–291040ViennaAustria
| | - Leticia González
- Institute of Theoretical ChemistryFaculty of ChemistryUniversity of ViennaWähringer Strasse 171090ViennaAustria
| |
Collapse
|
3
|
Peters LDM, Kussmann J, Ochsenfeld C. A Fermi smearing variant of the Tamm–Dancoff approximation for nonadiabatic dynamics involving S 1–S 0 transitions: Validation and application to azobenzene. J Chem Phys 2020; 153:094104. [DOI: 10.1063/5.0016487] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Laurens D. M. Peters
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Jörg Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
4
|
Mai S, González L. Molekulare Photochemie: Moderne Entwicklungen in der theoretischen Chemie. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.201916381] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Sebastian Mai
- Institut für Photonik Technische Universität Wien Gußhausstraße 27–29 1040 Wien Österreich
| | - Leticia González
- Institut für theoretische Chemie Fakultät für Chemie Universität Wien Währinger Straße 17 1090 Wien Österreich
| |
Collapse
|
5
|
Peters LM, Kussmann J, Ochsenfeld C. Combining Graphics Processing Units, Simplified Time-Dependent Density Functional Theory, and Finite-Difference Couplings to Accelerate Nonadiabatic Molecular Dynamics. J Phys Chem Lett 2020; 11:3955-3961. [PMID: 32374606 PMCID: PMC7304892 DOI: 10.1021/acs.jpclett.0c00320] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Accepted: 04/23/2020] [Indexed: 06/11/2023]
Abstract
Starting from our recently published implementation of nonadiabatic molecular dynamics (NAMD) on graphics processing units (GPUs), we explore further approaches to accelerate ab initio NAMD calculations at the time-dependent density functional theory (TDDFT) level of theory. We employ (1) the simplified TDDFT schemes of Grimme et al. and (2) the Hammes-Schiffer-Tully approach to obtain nonadiabatic couplings from finite-difference calculations. The resulting scheme delivers an accurate physical picture while virtually eliminating the two computationally most demanding steps of the algorithm. Combined with our GPU-based integral routines for SCF, TDDFT, and TDDFT derivative calculations, NAMD simulations of systems of a few hundreds of atoms at a reasonable time scale become accessible on a single compute node. To demonstrate this and to present a first, illustrative example, we perform TDDFT/MM-NAMD simulations of the rhodopsin protein.
Collapse
Affiliation(s)
- Laurens
D. M. Peters
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Jörg Kussmann
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Max
Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
6
|
Lacombe L, Maitra NT. Developing new and understanding old approximations in TDDFT. Faraday Discuss 2020; 224:382-401. [PMID: 32926040 DOI: 10.1039/d0fd00049c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
When a system has evolved far from a ground-state, the adiabatic approximations commonly used in time-dependent density functional theory calculations completely fail in some applications, while giving qualitatively good predictions in others, and sometimes even quantitative predictions. It is not clearly understood why this is so, and developing practical approximations going beyond the adiabatic approximation remains a challenge. This paper explores three different lines of investigation. First, an expression for the exact time-dependent exchange-correlation potential suggests that the accuracy of an adiabatic approximation is intimately related to the deviation between the natural orbital occupation numbers of the physical system and those of the Kohn-Sham system, and we explore this on some exactly-solvable model systems. The exact expression further suggests a path to go beyond the adiabatic approximations, and in the second part we discuss a newly proposed class of memory-dependent approximations developed in this way. Finally, we derive a new expression for the exact exchange-correlation potential from a coupling-constant path integration.
Collapse
Affiliation(s)
- Lionel Lacombe
- Department of Physics, Rutgers University, Newark, New Jersey 07102, USA.
| | | |
Collapse
|