1
|
Suchan J, Liang F, Durden AS, Levine BG. Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone. J Chem Phys 2024; 160:134310. [PMID: 38573851 DOI: 10.1063/5.0198333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 03/18/2024] [Indexed: 04/06/2024] Open
Abstract
Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum-classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?
Collapse
Affiliation(s)
- Jiří Suchan
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
2
|
Rana B, Hohenstein EG, Martínez TJ. Simulating the Excited-State Dynamics of Polaritons with Ab Initio Multiple Spawning. J Phys Chem A 2024; 128:139-151. [PMID: 38110364 DOI: 10.1021/acs.jpca.3c06607] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2023]
Abstract
Over the past decade, there has been a growth of interest in polaritonic chemistry, where the formation of hybrid light-matter states (polaritons) can alter the course of photochemical reactions. These hybrid states are created by strong coupling between molecules and photons in resonant optical cavities and can even occur in the absence of light when the molecule is strongly coupled with the electromagnetic fluctuations of the vacuum field. We present a first-principles model to simulate nonadiabatic dynamics of such polaritonic states inside optical cavities by leveraging graphical processing units (GPUs). Our first implementation of this model is specialized for a single molecule coupled to a single-photon mode confined inside the optical cavity but with any number of excited states computed using complete active space configuration interaction (CASCI) and a Jaynes-Cummings-type Hamiltonian. Using this model, we have simulated the excited-state dynamics of a single salicylideneaniline (SA) molecule strongly coupled to a cavity photon with the ab initio multiple spawning (AIMS) method. We demonstrate how the branching ratios of the photodeactivation pathways for this molecule can be manipulated by coupling to the cavity. We also show how one can stop the photoreaction from happening inside of an optical cavity. Finally, we also investigate cavity-based control of the ordering of two excited states (one optically bright and the other optically dark) inside a cavity for a set of molecules, where the dark and bright states are close in energy.
Collapse
Affiliation(s)
- Bhaskar Rana
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
3
|
Janoš J, Slavíček P. What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm. J Chem Theory Comput 2023; 19:8273-8284. [PMID: 37939301 PMCID: PMC10688183 DOI: 10.1021/acs.jctc.3c00908] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Revised: 10/18/2023] [Accepted: 10/19/2023] [Indexed: 11/10/2023]
Abstract
The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.
Collapse
Affiliation(s)
- Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague 6, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague 6, Czech Republic
| |
Collapse
|
4
|
Hohenstein EG, Oumarou O, Al-Saadon R, Anselmetti GLR, Scheurer M, Gogolin C, Parrish RM. Efficient quantum analytic nuclear gradients with double factorization. J Chem Phys 2023; 158:114119. [PMID: 36948843 DOI: 10.1063/5.0137167] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2023] Open
Abstract
Efficient representations of the Hamiltonian, such as double factorization, drastically reduce the circuit depth or the number of repetitions in error corrected and noisy intermediate-scale quantum (NISQ) algorithms for chemistry. We report a Lagrangian-based approach for evaluating relaxed one- and two-particle reduced density matrices from double factorized Hamiltonians, unlocking efficiency improvements in computing the nuclear gradient and related derivative properties. We demonstrate the accuracy and feasibility of our Lagrangian-based approach to recover all off-diagonal density matrix elements in classically simulated examples with up to 327 quantum and 18 470 total atoms in QM/MM simulations with modest-sized quantum active spaces. We show this in the context of the variational quantum eigensolver in case studies, such as transition state optimization, ab initio molecular dynamics simulation, and energy minimization of large molecular systems.
Collapse
|
5
|
Bannwarth C, Martínez TJ. SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods. J Chem Phys 2023; 158:074109. [PMID: 36813714 DOI: 10.1063/5.0132776] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. We present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn-Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.
Collapse
Affiliation(s)
- Christoph Bannwarth
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, 2575 Sand Hill Road, 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, 2575 Sand Hill Road, Menlo Park, California 94025, USA
| |
Collapse
|
6
|
Cruzeiro VWD, Wang Y, Pieri E, Hohenstein EG, Martínez TJ. TeraChem protocol buffers (TCPB): Accelerating QM and QM/MM simulations with a client-server model. J Chem Phys 2023; 158:044801. [PMID: 36725506 DOI: 10.1063/5.0130886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
The routine use of electronic structures in many chemical simulation applications calls for efficient and easy ways to access electronic structure programs. We describe how the graphics processing unit (GPU) accelerated electronic structure program TeraChem can be set up as an electronic structure server, to be easily accessed by third-party client programs. We exploit Google's protocol buffer framework for data serialization and communication. The client interface, called TeraChem protocol buffers (TCPB), has been designed for ease of use and compatibility with multiple programming languages, such as C++, Fortran, and Python. To demonstrate the ease of coupling third-party programs with electronic structures using TCPB, we have incorporated the TCPB client into Amber for quantum mechanics/molecular mechanics (QM/MM) simulations. The TCPB interface saves time with GPU initialization and I/O operations, achieving a speedup of more than 2× compared to a prior file-based implementation for a QM region with ∼250 basis functions. We demonstrate the practical application of TCPB by computing the free energy profile of p-hydroxybenzylidene-2,3-dimethylimidazolinone (p-HBDI-)-a model chromophore in green fluorescent proteins-on the first excited singlet state using Hamiltonian replica exchange for enhanced sampling. All calculations in this work have been performed with the non-commercial freely-available version of TeraChem, which is sufficient for many QM region sizes in common use.
Collapse
Affiliation(s)
| | - Yuanheng Wang
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Elisa Pieri
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
7
|
Nottoli T, Burger S, Stopkowicz S, Gauss J, Lipparini F. Computation of NMR shieldings at the CASSCF level using gauge-including atomic orbitals and Cholesky decomposition. J Chem Phys 2022; 157:084122. [DOI: 10.1063/5.0101838] [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
We present an implementation of coupled-perturbed complete active space self-consistent field (CP-CASSCF) theory for the computation of nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals and Cholesky decomposed two-electron integrals. The CP-CASSCF equations are solved using a direct algorithm where the magnetic Hessian matrix-vector product is expressed in terms of one-index transformed quantities. Numerical tests on systems with up to about 1300 basis functions provide information regarding both the computational efficiency and limitations of our implementation.
Collapse
Affiliation(s)
- Tommaso Nottoli
- Università degli Studi di Pisa Dipartimento di Chimica e Chimica Industriale, Italy
| | - Sophia Burger
- Chemistry, Johannes Gutenberg Universität Mainz, Germany
| | - Stella Stopkowicz
- Department of Chemistry, Johannes Gutenberg Universität Mainz, Germany
| | - Jürgen Gauss
- Johannes Gutenberg Universität Mainz, Department Chemie, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università degli Studi di Pisa Dipartimento di Chimica e Chimica Industriale, Italy
| |
Collapse
|
8
|
Kratz T, Steinbach P, Breitenlechner S, Storch G, Bannwarth C, Bach T. Photochemical Deracemization of Chiral Alkenes via Triplet Energy Transfer. J Am Chem Soc 2022; 144:10133-10138. [PMID: 35658423 DOI: 10.1021/jacs.2c02511] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
A visible-light-mediated, enantioselective approach to axially chiral alkenes is described. Starting from a racemic mixture, a major alkene enantiomer is formed due to selective triplet energy transfer from a catalytically active chiral sensitizer. A catalyst loading of 2 mol % was sufficient to guarantee consistently high enantioselectivities and yields (16 examples, 51%-quant., 81-96% ee). NMR studies and DFT computations revealed that triplet energy transfer is more rapid within the substrate-catalyst complex of the minor alkene enantiomer. Since this enantiomer is continuously racemized, the major enantiomer is enriched in the photostationary state.
Collapse
Affiliation(s)
- Thilo Kratz
- School of Natural Sciences, Department Chemie, and Catalysis Research Center (CRC), Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Pit Steinbach
- Institut für Physikalische Chemie, RWTH Aachen University, 52074 Aachen, Germany
| | - Stefan Breitenlechner
- School of Natural Sciences, Department Chemie, and Catalysis Research Center (CRC), Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Golo Storch
- School of Natural Sciences, Department Chemie, and Catalysis Research Center (CRC), Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Christoph Bannwarth
- Institut für Physikalische Chemie, RWTH Aachen University, 52074 Aachen, Germany
| | - Thorsten Bach
- School of Natural Sciences, Department Chemie, and Catalysis Research Center (CRC), Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| |
Collapse
|
9
|
Wang Y, Seritan S, Lahana D, Ford JE, Valentini A, Hohenstein EG, Martínez TJ. InteraChem: Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics. J Chem Theory Comput 2022; 18:3308-3317. [PMID: 35649124 DOI: 10.1021/acs.jctc.2c00005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree-Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.
Collapse
Affiliation(s)
- Yuanheng Wang
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Stefan Seritan
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Dean Lahana
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Jason E Ford
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Alessio Valentini
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
10
|
Poštulka J, Slavíček P, Pysanenko A, Poterya V, Fárník M. Bimolecular reactions on sticky and slippery clusters: Electron-induced reactions of hydrogen peroxide. J Chem Phys 2022; 156:054306. [DOI: 10.1063/5.0079283] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Affiliation(s)
- Jan Poštulka
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Andriy Pysanenko
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| | - Viktoriya Poterya
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| | - Michal Fárník
- J. Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 2155/3, Prague 8, Czech Republic
| |
Collapse
|
11
|
Matsika S. Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections. Chem Rev 2021; 121:9407-9449. [PMID: 34156838 DOI: 10.1021/acs.chemrev.1c00074] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.
Collapse
Affiliation(s)
- Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
12
|
Plaza M, Großkopf J, Breitenlechner S, Bannwarth C, Bach T. Photochemical Deracemization of Primary Allene Amides by Triplet Energy Transfer: A Combined Synthetic and Theoretical Study. J Am Chem Soc 2021; 143:11209-11217. [PMID: 34279085 DOI: 10.1021/jacs.1c05286] [Citation(s) in RCA: 38] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
The photochemical deracemization of 2,4-disubstituted 2,3-butadienamides (allene amides) was investigated both experimentally and theoretically. The reaction was catalyzed by a thioxanthone which is covalently linked to a chiral 1,5,7-trimethyl-3-azabicyclo[3.3.1]nonan-2-one skeleton providing a U-shaped arrangement of the sensitizing unit relative to a potential hydrogen-bonding site. Upon irradiation at λ = 420 nm in the presence of the sensitizer (2.5 mol %), the amides reached at -10 °C a photostationary state in which one enantiomer prevailed. The enantioenriched allene amides (70-93% ee) were isolated in 74% to quantitative yield (19 examples). Based on luminescence data and DFT calculations, energy transfer from the thioxanthone to the allene amides is thermodynamically feasible, and the achiral triplet allene intermediate was structurally characterized. Hydrogen bonding of the amide enantiomers to the sensitizer was monitored by NMR titration. The experimental association constants (Ka) were similar (59.8 vs 25.7 L·mol-1). DFT calculations, however, revealed a significant difference in the binding properties of the two enantiomers. The major product enantiomer exhibits a noncovalent dispersion interaction of its arylmethyl group to the external benzene ring of the thioxanthone, thus moving away the allene from the carbonyl chromophore. The minor enantiomer displays a CH-π interaction of the hydrogen atom at the terminal allene carbon atom to the same benzene ring, thus forcing the allene into close proximity to the chromophore. The binding behavior explains the observed enantioselectivity which, as corroborated by additional calculations, is due to a rapid triplet energy transfer within the substrate-catalyst complex of the minor enantiomer.
Collapse
Affiliation(s)
- Manuel Plaza
- Department Chemie and Catalysis Research Center (CRC), Technische Universität München, D-85747 Garching, Germany
| | - Johannes Großkopf
- Department Chemie and Catalysis Research Center (CRC), Technische Universität München, D-85747 Garching, Germany
| | - Stefan Breitenlechner
- Department Chemie and Catalysis Research Center (CRC), Technische Universität München, D-85747 Garching, Germany
| | - Christoph Bannwarth
- Institute of Physical Chemistry, RWTH Aachen University, D-52056 Aachen, Germany
| | - Thorsten Bach
- Department Chemie and Catalysis Research Center (CRC), Technische Universität München, D-85747 Garching, Germany
| |
Collapse
|
13
|
Levine BG, Durden AS, Esch MP, Liang F, Shu Y. CAS without SCF-Why to use CASCI and where to get the orbitals. J Chem Phys 2021; 154:090902. [PMID: 33685182 DOI: 10.1063/5.0042147] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree-Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Finally, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.
Collapse
Affiliation(s)
- Benjamin G Levine
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Michael P Esch
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute for Advanced Computational Science and Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
14
|
Song C, Neaton JB, Martínez TJ. Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method. J Chem Phys 2021; 154:014103. [PMID: 33412861 DOI: 10.1063/5.0035233] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Physics, University of California Berkeley, Berkeley, California 94720, USA
| | - Jeffrey B Neaton
- Department of Physics, University of California Berkeley, Berkeley, California 94720, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
15
|
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
|
16
|
Yu JK, Bannwarth C, Hohenstein EG, Martínez TJ. Ab Initio Nonadiabatic Molecular Dynamics with Hole–Hole Tamm–Dancoff Approximated Density Functional Theory. J Chem Theory Comput 2020; 16:5499-5511. [DOI: 10.1021/acs.jctc.0c00644] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Jimmy K. Yu
- Biophysics Program, Stanford University, Stanford, California 94305, United States
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Christoph Bannwarth
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Edward G. Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Todd J. Martínez
- Biophysics Program, Stanford University, Stanford, California 94305, United States
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
17
|
Seritan S, Bannwarth C, Fales BS, Hohenstein EG, Isborn CM, Kokkila‐Schumacher SIL, Li X, Liu F, Luehr N, Snyder JW, Song C, Titov AV, Ufimtsev IS, Wang L, Martínez TJ. TeraChem
: A graphical processing unit
‐accelerated
electronic structure package for
large‐scale
ab initio molecular dynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1494] [Citation(s) in RCA: 53] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Stefan Seritan
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Christoph Bannwarth
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Bryan S. Fales
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Edward G. Hohenstein
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| | - Christine M. Isborn
- Department of Chemistry University of California Merced Merced California USA
| | | | - Xin Li
- Division of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Fang Liu
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts USA
| | | | | | - Chenchen Song
- Department of Physics University of California Berkeley Berkeley California USA
- Molecular Foundry Lawrence Berkeley National Laboratory Berkeley California USA
| | | | - Ivan S. Ufimtsev
- Department of Structural Biology Stanford University School of Medicine Stanford California USA
| | - Lee‐Ping Wang
- Department of Chemistry University of California Davis Davis California USA
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute Stanford University Stanford California USA
- SLAC National Accelerator Laboratory Menlo Park California USA
| |
Collapse
|
18
|
Suchan J, Janoš J, Slavíček P. Pragmatic Approach to Photodynamics: Mixed Landau–Zener Surface Hopping with Intersystem Crossing. J Chem Theory Comput 2020; 16:5809-5820. [DOI: 10.1021/acs.jctc.0c00512] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Jiří Suchan
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| |
Collapse
|
19
|
Song C, Martínez TJ. Reduced scaling extended multi-state CASPT2 (XMS-CASPT2) using supporting subspaces and tensor hyper-contraction. J Chem Phys 2020; 152:234113. [PMID: 32571032 DOI: 10.1063/5.0007417] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
We present a reduced scaling formulation of the extended multi-state CASPT2 (XMS-CASPT2) method, which is based on our recently developed state-specific CASPT2 (SS-CASPT2) formulation using supporting subspaces and tensor hyper-contraction. By using these two techniques, the off-diagonal elements of the effective Hamiltonian can be computed with only O(N3) operations and O(N2) memory, where N is the number of basis functions. This limits the overall computational scaling to O(N4) operations and O(N2) memory. Thus, excited states can now be obtained at the same reduced (relative to previous algorithms) scaling we achieved for SS-CASPT2. In addition, we also investigate how the energy denominators can be factorized with the Laplace quadrature when some of the denominators are negative, which is critical for excited state calculations. An efficient implementation of the method has been developed using graphical processing units while also exploiting spatial sparsity in tensor operations. We benchmark the accuracy of the new method by comparison to non-THC formulated XMS-CASPT2 for the excited states of various molecules. In our tests, the THC approximation introduces negligible errors (≈0.01 eV) compared to the non-THC reference method. Scaling behavior and computational timings are presented to demonstrate performance. The new method is also interfaced with quantum mechanics/molecular mechanics (QM/MM). In an example study of green fluorescent protein, we show how the XMS-CASPT2 potential energy surfaces and excitation energies are affected by increasing the size of the QM region up to 278 QM atoms with more than 2300 basis functions.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
20
|
Fales BS, Martínez TJ. Efficient Treatment of Large Active Spaces through Multi-GPU Parallel Implementation of Direct Configuration Interaction. J Chem Theory Comput 2020; 16:1586-1596. [DOI: 10.1021/acs.jctc.9b01165] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- B. Scott Fales
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Todd J. Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
21
|
Mullinax JW, Maradzike E, Koulias LN, Mostafanejad M, Epifanovsky E, Gidofalvi G, DePrince AE. Heterogeneous CPU + GPU Algorithm for Variational Two-Electron Reduced-Density Matrix-Driven Complete Active-Space Self-Consistent Field Theory. J Chem Theory Comput 2019; 15:6164-6178. [DOI: 10.1021/acs.jctc.9b00768] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- J. Wayne Mullinax
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Elvis Maradzike
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Lauren N. Koulias
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Mohammad Mostafanejad
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| |
Collapse
|
22
|
Valleau S, Martínez TJ. Reaction Dynamics of Cyanohydrins with Hydrosulfide in Water. J Phys Chem A 2019; 123:7210-7217. [DOI: 10.1021/acs.jpca.9b05735] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Stéphanie Valleau
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Todd J. Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
23
|
Thongyod W, Buranachai C, Pengpan T, Punwong C. Fluorescence quenching by photoinduced electron transfer between 7-methoxycoumarin and guanine base facilitated by hydrogen bonds: an in silico study. Phys Chem Chem Phys 2019; 21:16258-16269. [PMID: 31304496 DOI: 10.1039/c9cp02037c] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this study, the effects of hydrogen bond (H-bond) formation on fluorescence quenching of 7-methoxycoumarin (7MC) via photo-induced electron transfer from a guanine base (Gua) are investigated using a combined quantum mechanics/molecular mechanics simulation. The electronic structure is calculated by the floating occupation molecular orbital complete active space configuration interaction modification on a semiempirical method. Then the full multiple spawning method is employed for the dynamics simulations on multiple electronic states. The methods employed here are validated by simulating direct dynamics of 7MC (without Gua) and compared with available experimental results. Our computational results are in good agreement with the previously reported experimental results in terms of spectroscopic properties of 7MC. In the case of a H-bonded 7MC-Gua complex, the results from constrained dynamics simulations and single-point calculations suggest that the electron transfer occurs on the second excited state and it depends not only on the H-bond length but also on the intermolecular planarity between 7MC and Gua. Moreover, a proton coupled electron transfer can occur at ≈1 Å of H-bond length, where a proton from Gua is also transferred together with the electron to 7MC. The obtained simulations are expected to be greatly beneficial for designing effective fluorescently labeled nucleotide probes as well as providing information for precise fluorescence signal interpretation.
Collapse
Affiliation(s)
- Wutthinan Thongyod
- Department of Physics, Faculty of Science, Prince of Songkla University, Songkhla 90112, Thailand. and Center of Excellence for Trace Analysis and Biosensor, Prince of Songkla University, Songkhla 90112, Thailand
| | - Chittanon Buranachai
- Department of Physics, Faculty of Science, Prince of Songkla University, Songkhla 90112, Thailand. and Center of Excellence for Trace Analysis and Biosensor, Prince of Songkla University, Songkhla 90112, Thailand
| | - Teparksorn Pengpan
- Department of Physics, Faculty of Science, Prince of Songkla University, Songkhla 90112, Thailand.
| | - Chutintorn Punwong
- Department of Physics, Faculty of Science, Prince of Songkla University, Songkhla 90112, Thailand.
| |
Collapse
|
24
|
Park JW. Single-State Single-Reference and Multistate Multireference Zeroth-Order Hamiltonians in MS-CASPT2 and Conical Intersections. J Chem Theory Comput 2019; 15:3960-3973. [DOI: 10.1021/acs.jctc.9b00067] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University (CBNU), Cheongju 28644, Korea
| |
Collapse
|
25
|
Esch MP, Shu Y, Levine BG. A Conical Intersection Perspective on the Low Nonradiative Recombination Rate in Lead Halide Perovskites. J Phys Chem A 2019; 123:2661-2673. [DOI: 10.1021/acs.jpca.9b00952] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Michael P. Esch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
26
|
Ranković M, Chalabala J, Zawadzki M, Kočišek J, Slavíček P, Fedor J. Dissociative ionization dynamics of dielectric gas C3F7CN. Phys Chem Chem Phys 2019; 21:16451-16458. [PMID: 31312828 DOI: 10.1039/c9cp02188d] [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
Fluoronitrile C3F7CN is a promising candidate for the replacement of SF6 dielectric gas in high-voltage insulation.
Collapse
Affiliation(s)
- M. Ranković
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
| | - J. Chalabala
- Department of Physical Chemistry
- University of Chemistry and Technology
- 16628 Prague
- Czech Republic
| | - M. Zawadzki
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
- Department of Atomic, Molecular, and Optical Physics
| | - J. Kočišek
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
| | - P. Slavíček
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
- Department of Physical Chemistry
| | - J. Fedor
- J. Heyrovský Institute of Physical Chemistry
- Czech Academy of Sciences
- 18223 Prague
- Czech Republic
| |
Collapse
|
27
|
|
28
|
Yang J, Zhu X, Wolf TJA, Li Z, Nunes JPF, Coffee R, Cryan JP, Gühr M, Hegazy K, Heinz TF, Jobe K, Li R, Shen X, Veccione T, Weathersby S, Wilkin KJ, Yoneda C, Zheng Q, Martinez TJ, Centurion M, Wang X. Imaging CF 3I conical intersection and photodissociation dynamics with ultrafast electron diffraction. Science 2018; 361:64-67. [PMID: 29976821 DOI: 10.1126/science.aat0049] [Citation(s) in RCA: 117] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2018] [Accepted: 05/15/2018] [Indexed: 11/02/2022]
Abstract
Conical intersections play a critical role in excited-state dynamics of polyatomic molecules because they govern the reaction pathways of many nonadiabatic processes. However, ultrafast probes have lacked sufficient spatial resolution to image wave-packet trajectories through these intersections directly. Here, we present the simultaneous experimental characterization of one-photon and two-photon excitation channels in isolated CF3I molecules using ultrafast gas-phase electron diffraction. In the two-photon channel, we have mapped out the real-space trajectories of a coherent nuclear wave packet, which bifurcates onto two potential energy surfaces when passing through a conical intersection. In the one-photon channel, we have resolved excitation of both the umbrella and the breathing vibrational modes in the CF3 fragment in multiple nuclear dimensions. These findings benchmark and validate ab initio nonadiabatic dynamics calculations.
Collapse
Affiliation(s)
- Jie Yang
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA. .,Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Xiaolei Zhu
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Department of Chemistry, Stanford University, Stanford, CA, USA
| | - Thomas J A Wolf
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Zheng Li
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron (DESY), Hamburg, Germany.,Max Planck Institute for the Structure and Dynamics of Matter, Hamburg, Germany
| | - J Pedro F Nunes
- Department of Chemistry, University of York, Heslington, York, UK
| | - Ryan Coffee
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Department of Physics, Stanford University, Stanford, CA, USA
| | - James P Cryan
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Markus Gühr
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Institut für Physik und Astronomie, Universität Potsdam, Potsdam, Germany
| | - Kareem Hegazy
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Department of Physics, Stanford University, Stanford, CA, USA
| | - Tony F Heinz
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA.,Department of Applied Physics, Stanford University, Stanford, CA, USA
| | - Keith Jobe
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Renkai Li
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Xiaozhe Shen
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | | | | | - Kyle J Wilkin
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Lincoln, NE, USA
| | - Charles Yoneda
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Qiang Zheng
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA
| | - Todd J Martinez
- Stanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, CA, USA. .,Department of Chemistry, Stanford University, Stanford, CA, USA
| | - Martin Centurion
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Lincoln, NE, USA.
| | - Xijie Wang
- SLAC National Accelerator Laboratory, Menlo Park, CA, USA.
| |
Collapse
|
29
|
Pijeau S, Foster D, Hohenstein EG. Effect of Nonplanarity on Excited-State Proton Transfer and Internal Conversion in Salicylideneaniline. J Phys Chem A 2018; 122:5555-5562. [DOI: 10.1021/acs.jpca.8b02426] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Shiela Pijeau
- Department of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
| | - Donneille Foster
- Department of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
| | - Edward G. Hohenstein
- Department of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| |
Collapse
|
30
|
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
|
31
|
Affiliation(s)
- Basile F. E. Curchod
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
32
|
Peng WT, Fales BS, Shu Y, Levine BG. Dynamics of recombination via conical intersection in a semiconductor nanocrystal. Chem Sci 2018; 9:681-687. [PMID: 29629136 PMCID: PMC5869574 DOI: 10.1039/c7sc04221c] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2017] [Accepted: 11/13/2017] [Indexed: 12/17/2022] Open
Abstract
Conical intersections are well known to introduce nonradiative decay pathways in molecules, but have only recently been implicated in nonradiative recombination processes in materials. Here we apply excited state ab initio molecular dynamics simulations based on a multireference description of the electronic structure to defective silicon nanocrystals up to 1.7 nm in diameter to search for accessible nonradiative recombination pathways. Dangling bond defects are found to induce conical intersections between the ground and first excited electronic states of five systems of various sizes. These defect-induced conical intersections are accessible at energies that are in the visible range (2.4-2.7 eV) and very weakly dependent on particle size. The dynamic simulations suggest that these intersections are accessed 40-60 fs after creation of a defect-localized excitation. This ultrafast recombination is attributed to the fact that Jahn-Teller distortion on the first excited state drives the defect directly towards a conical intersection with the ground electronic state.
Collapse
Affiliation(s)
- Wei-Tao Peng
- Department of Chemistry , Michigan State University , East Lansing , MI 48824 , USA .
| | - B Scott Fales
- Department of Chemistry , The PULSE Institute , Stanford University , Stanford , CA 94305 , USA.,SLAC National Accelerator Laboratory , Menlo Park , CA 94025 , USA
| | - Yinan Shu
- Department of Chemistry , University of Minnesota , Minneapolis , MN 55455 , USA
| | - Benjamin G Levine
- Department of Chemistry , Michigan State University , East Lansing , MI 48824 , USA .
| |
Collapse
|
33
|
Hollas D, Šištík L, Hohenstein EG, Martínez TJ, Slavíček P. Nonadiabatic Ab Initio Molecular Dynamics with the Floating Occupation Molecular Orbital-Complete Active Space Configuration Interaction Method. J Chem Theory Comput 2017; 14:339-350. [DOI: 10.1021/acs.jctc.7b00958] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daniel Hollas
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Lukáš Šištík
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
| | - Edward G. Hohenstein
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
- PhD
Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Todd J. Martínez
- Department
of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Petr Slavíček
- Department
of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague 6, Czech Republic
- J. Heyrovský Institute of Physical Chemistry, Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
34
|
Fales BS, Shu Y, Levine BG, Hohenstein EG. Complete active space configuration interaction from state-averaged configuration interaction singles natural orbitals: Analytic first derivatives and derivative coupling vectors. J Chem Phys 2017; 147:094104. [DOI: 10.1063/1.5000476] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Affiliation(s)
- B. Scott Fales
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Yinan Shu
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Benjamin G. Levine
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Edward G. Hohenstein
- Department of Chemistry and Biochemistry, City College of New York, New York, New York 10031, USA
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, USA
| |
Collapse
|
35
|
Pijeau S, Foster D, Hohenstein EG. Excited-State Dynamics of a Benzotriazole Photostabilizer: 2-(2′-Hydroxy-5′-methylphenyl)benzotriazole. J Phys Chem A 2017; 121:6377-6387. [DOI: 10.1021/acs.jpca.7b04504] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Shiela Pijeau
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
| | - Donneille Foster
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
| | - Edward G. Hohenstein
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
- Ph.D.
Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| |
Collapse
|
36
|
Pijeau S, Foster D, Hohenstein EG. Excited-State Dynamics of 2-(2′-Hydroxyphenyl)benzothiazole: Ultrafast Proton Transfer and Internal Conversion. J Phys Chem A 2017; 121:4595-4605. [DOI: 10.1021/acs.jpca.7b01215] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Shiela Pijeau
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
| | - Donneille Foster
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
| | - Edward G. Hohenstein
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
- Ph.D.
Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| |
Collapse
|
37
|
Nishimoto Y. Analytic hyperpolarizability and polarizability derivative with fractional occupation numbers for large extended systems. J Chem Phys 2017; 146:084101. [DOI: 10.1063/1.4976551] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
38
|
Pijeau S, Hohenstein EG. Improved Complete Active Space Configuration Interaction Energies with a Simple Correction from Density Functional Theory. J Chem Theory Comput 2017; 13:1130-1146. [DOI: 10.1021/acs.jctc.6b00893] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Shiela Pijeau
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
| | - Edward G. Hohenstein
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
- Ph.D.
Program in Chemistry, The City University of New York, New York, New York 10016, United States
| |
Collapse
|
39
|
Hohenstein EG. Analytic formulation of derivative coupling vectors for complete active space configuration interaction wavefunctions with floating occupation molecular orbitals. J Chem Phys 2016; 145:174110. [DOI: 10.1063/1.4966235] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Edward G. Hohenstein
- Department of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, USA and Ph.D. Program in Chemistry, The City University of New York, New York, New York 10016, USA
| |
Collapse
|
40
|
Curchod BFE, Rauer C, Marquetand P, González L, Martínez TJ. Communication: GAIMS—Generalized Ab Initio Multiple Spawning for both internal conversion and intersystem crossing processes. J Chem Phys 2016; 144:101102. [DOI: 10.1063/1.4943571] [Citation(s) in RCA: 77] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Basile F. E. Curchod
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Clemens Rauer
- Institute of Theoretical Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Leticia González
- Institute of Theoretical Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
41
|
Wang LP, McGibbon RT, Pande VS, Martinez TJ. Automated Discovery and Refinement of Reactive Molecular Dynamics Pathways. J Chem Theory Comput 2016; 12:638-49. [PMID: 26683346 DOI: 10.1021/acs.jctc.5b00830] [Citation(s) in RCA: 81] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We describe a flexible and broadly applicable energy refinement method, "nebterpolation," for identifying and characterizing the reaction events in a molecular dynamics (MD) simulation. The new method is applicable to ab initio simulations with hundreds of atoms containing complex and multimolecular reaction events. A key aspect of nebterpolation is smoothing of the reactive MD trajectory in internal coordinates to initiate the search for the reaction path on the potential energy surface. We apply nebterpolation to analyze the reaction events in an ab initio nanoreactor simulation that discovers new molecules and mechanisms, including a C-C coupling pathway for glycolaldehyde synthesis. We find that the new method, which incorporates information from the MD trajectory that connects reactants with products, produces a dramatically distinct set of minimum energy paths compared to existing approaches that start from information for the reaction end points alone. The energy refinement method described here represents a key component of an emerging simulation paradigm where molecular dynamics simulations are applied to discover the possible reaction mechanisms.
Collapse
Affiliation(s)
| | | | | | - Todd J Martinez
- SLAC Linear Accelerator Laboratory , Menlo Park, California 94025, United States
| |
Collapse
|