1
|
Holzer C, Franzke YJ. Beyond Electrons: Correlation and Self-Energy in Multicomponent Density Functional Theory. Chemphyschem 2024; 25:e202400120. [PMID: 38456204 DOI: 10.1002/cphc.202400120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Revised: 03/07/2024] [Accepted: 03/08/2024] [Indexed: 03/09/2024]
Abstract
Post-Kohn-Sham methods are used to evaluate the ground-state correlation energy and the orbital self-energy of systems consisting of multiple flavors of different fermions. Starting from multicomponent density functional theory, suitable ways to arrive at the corresponding multicomponent random-phase approximation and the multicomponent Green's functionG W ${GW}$ approximation, including relativistic effects, are outlined. Given the importance of both of this methods in the development of modern Kohn-Sham density functional approximations, this work will provide a foundation to design advanced multicomponent density functional approximations. Additionally, theG W ${GW}$ quasiparticle energies are needed to study light-matter interactions with the Bethe-Salpeter equation.
Collapse
Affiliation(s)
- Christof Holzer
- Karlsruhe Institute of Technology (KIT), Institute of Theoretical Solid State Physics, Kaiserstraße 12, 76131, Karlsruhe, Germany
| | - Yannick J Franzke
- Friedrich Schiller University Jena, Otto Schott Institute of Materials Research, Löbdergraben 32, 07743, Jena, Germany
| |
Collapse
|
2
|
Garner SM, Upadhyay S, Li X, Hammes-Schiffer S. Nuclear-Electronic Orbital Time-Dependent Configuration Interaction Method. J Phys Chem Lett 2024; 15:6017-6023. [PMID: 38815051 DOI: 10.1021/acs.jpclett.4c00805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Combining real-time electronic structure with the nuclear-electronic orbital (NEO) method has enabled the simulation of complex nonadiabatic chemical processes. However, accurate descriptions of hydrogen tunneling and double excitations require multiconfigurational treatments. Herein, we develop and implement the real-time NEO time-dependent configuration interaction (NEO-TDCI) approach. Comparison to NEO-full CI calculations of absorption spectra for a molecular system shows that the NEO-TDCI approach can accurately capture the tunneling splitting associated with the electronic ground state as well as vibronic progressions corresponding to double electron-proton excitations associated with excited electronic states. Both of these features are absent from spectra obtained with single reference real-time NEO methods. Our simulations of hydrogen tunneling dynamics illustrate the oscillation of the proton density from one side to the other via a delocalized, bilobal proton wave function. These results indicate that the NEO-TDCI approach is highly suitable for studying hydrogen tunneling and other inherently multiconfigurational systems.
Collapse
Affiliation(s)
- Scott M Garner
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Shiv Upadhyay
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
3
|
Hasecke L, Mata RA. Optimization of Quantum Nuclei Positions with the Adaptive Nuclear-Electronic Orbital Approach. J Phys Chem A 2024; 128:3205-3211. [PMID: 38619054 PMCID: PMC11056972 DOI: 10.1021/acs.jpca.4c00096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/08/2024] [Accepted: 03/08/2024] [Indexed: 04/16/2024]
Abstract
The use of multicomponent methods has become increasingly popular over the last years. Under this framework, nuclei (commonly protons) are treated quantum mechanically on the same footing as the electronic structure problem. Under the use of atomic-centered orbitals, this can lead to some complications as the ideal location of the nuclear basis centers must be optimized. In this contribution, we propose a straightforward approach to determine the position of such centers within the self-consistent cycle of a multicomponent calculation, making use of individual proton charge centroids. We test the method on model systems including the water dimer, a protonated water tetramer, and a porphine system. Comparing to numerical gradient calculations, the adaptive nuclear-electronic orbital (NEO) procedure is able to converge the basis centers to within a few cents of an Ångström and with less than 0.1 kcal/mol differences in absolute energies. This is achieved in one single calculation and with a small added computational effort of up to 80% compared to a regular NEO- self-consistent field run. An example application for the human transketolase proton wire is also provided.
Collapse
Affiliation(s)
- Lukas Hasecke
- Institute of Physical Chemistry, University
of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Ricardo A. Mata
- Institute of Physical Chemistry, University
of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| |
Collapse
|
4
|
Lambros E, Fetherolf JH, Hammes-Schiffer S, Li X. A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters. J Phys Chem Lett 2024; 15:4070-4075. [PMID: 38587257 DOI: 10.1021/acs.jpclett.4c00439] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Abstract
Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. By performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Jonathan H Fetherolf
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
5
|
Xu J, Carney TE, Zhou R, Shepard C, Kanai Y. Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics. J Am Chem Soc 2024; 146:5011-5029. [PMID: 38362887 DOI: 10.1021/jacs.3c08226] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2024]
Abstract
The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. In this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.
Collapse
Affiliation(s)
- Jianhang Xu
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Thomas E Carney
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Christopher Shepard
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
- Department of Physics and Astronomy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
6
|
Ramos P, Friedman H, Li BY, Garcia C, Sletten E, Caram JR, Jang SJ. Nonadiabatic Derivative Couplings through Multiple Franck-Condon Modes Dictate the Energy Gap Law for Near and Short-Wave Infrared Dye Molecules. J Phys Chem Lett 2024; 15:1802-1810. [PMID: 38329913 DOI: 10.1021/acs.jpclett.3c02629] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2024]
Abstract
Near infrared (NIR, 700-1000 nm) and short-wave infrared (SWIR, 1000-2000 nm) dye molecules exhibit significant nonradiative decay rates from the first singlet excited state to the ground state. While these trends can be empirically explained by a simple energy gap law, detailed mechanisms of nearly universal behavior have remained unsettled for many cases. Theoretical and experimental results for two representative NIR/SWIR dye molecules reported here clarify the key mechanism for the observed energy gap law behavior. It is shown that the first derivative nonadiabatic coupling terms serve as major coupling pathways for nonadiabatic decay processes from the first excited singlet state to the ground state for these NIR and SWIR dye molecules and that vibrational modes other than the highest frequency modes also make significant contributions to the rate. This assessment is corroborated by further theoretical comparison with possible alternative mechanisms of intersystem crossing to triplet states and also by comparison with experimental data for deuterated molecules.
Collapse
Affiliation(s)
- Pablo Ramos
- Department of Chemistry and Biochemistry, Queens College, City University of New York, 65-30 Kissena Boulevard, New York, New York 11367, United States
| | - Hannah Friedman
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Barry Y Li
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Cesar Garcia
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Ellen Sletten
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Justin R Caram
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, United States
| | - Seogjoo J Jang
- Department of Chemistry and Biochemistry, Queens College, City University of New York, 65-30 Kissena Boulevard, New York, New York 11367, United States
- Chemistry and Physics PhD programs, Graduate Center, City University of New York, New York, New York 10016, United States
| |
Collapse
|
7
|
Li TE, Paenurk E, Hammes-Schiffer S. Squeezed Protons and Infrared Plasmonic Resonance Energy Transfer. J Phys Chem Lett 2024; 15:751-757. [PMID: 38226772 DOI: 10.1021/acs.jpclett.3c03112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2024]
Abstract
Unusual nuclear quantum effects may emerge near noble metal nanostructures such as squeezed vibrational states in molecular junctions and plasmonic resonance energy transfer in the infrared domain. Herein, nuclear quantum effects near heavy metals are studied by nuclear-electronic orbital density functional theory (NEO-DFT) with an effective core potential. For a quantum proton sandwiched between a pair of gold tips modeled by two Au6 clusters, NEO-DFT calculations suggest that the quantum proton density can be squeezed as the tip distance decreases. For an HF molecule placed near a one-dimensional Au nanowire composed of up to 34 Au atoms, real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) shows that the infrared plasmonic motion within the Au nanowire may resonantly transfer electronic energy to the HF proton vibrational stretch mode. Overall, these calculations illustrate the advantages of the NEO approach for probing nuclear quantum effects, such as squeezed proton vibrational states and infrared plasmonic resonance energy transfer.
Collapse
Affiliation(s)
- Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Eno Paenurk
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
8
|
Xu J, Zhou R, Blum V, Li TE, Hammes-Schiffer S, Kanai Y. First-Principles Approach for Coupled Quantum Dynamics of Electrons and Protons in Heterogeneous Systems. PHYSICAL REVIEW LETTERS 2023; 131:238002. [PMID: 38134781 DOI: 10.1103/physrevlett.131.238002] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Accepted: 11/01/2023] [Indexed: 12/24/2023]
Abstract
The coupled quantum dynamics of electrons and protons is ubiquitous in many dynamical processes involving light-matter interaction, such as solar energy conversion in chemical systems and photosynthesis. A first-principles description of such nuclear-electronic quantum dynamics requires not only the time-dependent treatment of nonequilibrium electron dynamics but also that of quantum protons. Quantum mechanical correlation between electrons and protons adds further complexity to such coupled dynamics. Here we extend real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) to periodic systems and perform first-principles simulations of coupled quantum dynamics of electrons and protons in complex heterogeneous systems. The process studied is an electronically excited-state intramolecular proton transfer of o-hydroxybenzaldehyde in water and at a silicon (111) semiconductor-molecule interface. These simulations illustrate how environments such as hydrogen-bonding water molecules and an extended material surface impact the dynamical process on the atomistic level. Depending on how the molecule is chemisorbed on the surface, excited-state electron transfer from the molecule to the semiconductor surface can inhibit ultrafast proton transfer within the molecule. This Letter elucidates how heterogeneous environments influence the balance between the quantum mechanical proton transfer and excited electron dynamics. The periodic RT-NEO-TDDFT approach is applicable to a wide range of other photoinduced heterogeneous processes.
Collapse
Affiliation(s)
- Jianhang Xu
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Volker Blum
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina, USA and Department of Chemistry, Duke University, Durham, North Carolina, USA
| | - Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut, USA
| | | | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA and Department of Physics and Astronomy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| |
Collapse
|
9
|
Hasecke L, Mata RA. Nuclear Quantum Effects Made Accessible: Local Density Fitting in Multicomponent Methods. J Chem Theory Comput 2023; 19:8223-8233. [PMID: 37920900 PMCID: PMC10687858 DOI: 10.1021/acs.jctc.3c01055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2023] [Revised: 10/18/2023] [Accepted: 10/19/2023] [Indexed: 11/04/2023]
Abstract
The simulation of nuclear quantum effects (NQEs) is crucial for an accurate description of systems and processes involving light nuclei, such as hydrogen atoms. Within the last years, the importance of those effects has been highlighted for a vast range of systems with tremendous implications in chemistry, biology, physics, and materials sciences. However, while electronic structure theory methods have become routine tools for quantum chemical investigations, there is still a lack of approaches to address NQEs that are computationally accessible and straightforward to use. To address this, we present the first combination of the nuclear-electronic orbital Hartree-Fock approach with both local and density fitting approximations (LDF-NEO-HF). This results in a low-order scaling approach that enables the inclusion of NQEs for large systems within a fraction of a day and for small to medium size systems in minutes. Moreover, we demonstrate the qualitative accuracy and robustness of our approach to retrieve NQEs for three real-use cases motivated by chemical, biological, and materials science applications.
Collapse
Affiliation(s)
- Lukas Hasecke
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Ricardo A. Mata
- Institute of Physical Chemistry, University of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| |
Collapse
|
10
|
Lambros E, Link B, Chow M, Lipparini F, Hammes-Schiffer S, Li X. Assessing Implicit and Explicit Polarizable Solvation Models for Nuclear-Electronic Orbital Systems: Quantum Proton Polarization and Solvation Energetics. J Phys Chem A 2023; 127:9322-9333. [PMID: 37889479 DOI: 10.1021/acs.jpca.3c03153] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2023]
Abstract
Accurate simulations of many chemical processes require the inclusion of both nuclear quantum effects and a solvent environment. The nuclear-electronic orbital (NEO) approach, which treats electrons and select nuclei quantum mechanically on the same level, combined with a polarizable continuum model (PCM) for the solvent environment, addresses this challenge in a computationally practical manner. In this work, the NEO-PCM approach is extended beyond the IEF-PCM (integral equation formalism PCM) and C-PCM (conductor PCM) approaches to the SS(V)PE (surface and simulation of volume polarization for electrostatics) and ddCOSMO (domain decomposed conductor-like screening model) approaches. IEF-PCM, SS(V)PE, C-PCM, and ddCOSMO all exhibit similar solvation energies as well as comparable nuclear polarization within the NEO framework. The calculations show that the nuclear density does not leak out of the molecular cavity because it is much more localized than the electronic density. Finally, the polarization of quantized protons is analyzed in both continuum solvent and explicit solvent environments described by the polarizable MB-pol model, illustrating the impact of specific hydrogen-bonding interactions captured only by explicit solvation. These calculations highlight the relationship among solvation formalism, nuclear polarization, and energetics.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Benjamin Link
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Mathew Chow
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | | | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
11
|
Chow M, Li TE, Hammes-Schiffer S. Nuclear-Electronic Orbital Quantum Mechanical/Molecular Mechanical Real-Time Dynamics. J Phys Chem Lett 2023; 14:9556-9562. [PMID: 37857272 DOI: 10.1021/acs.jpclett.3c02275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2023]
Abstract
Simulating the nuclear-electronic quantum dynamics of large-scale molecular systems in the condensed phase is key for studying biologically and chemically important processes such as proton transfer and proton-coupled electron transfer reactions. Herein, the real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) approach is combined with a hybrid quantum mechanical/molecular mechanical (QM/MM) strategy to enable the accurate description of coupled nuclear-electronic quantum dynamics in the presence of heterogeneous environments such as solvent or proteins. The densities of the electrons and quantum protons are propagated in real time, while the other nuclei are propagated classically on the instantaneous electron-proton vibronic surface. This approach is applied to phenol bound to lysozyme, intramolecular proton transfer in malonaldehyde, and nonequilibrium excited-state intramolecular proton transfer in o-hydroxybenzaldehyde. These examples illustrate that the RT-NEO-TDDFT framework, coupled with an atomistic representation of the environment, allows the simulation of condensed-phase systems that exhibit significant nuclear quantum effects.
Collapse
Affiliation(s)
- Mathew Chow
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
12
|
Di Felice R, Mayes ML, Richard RM, Williams-Young DB, Chan GKL, de Jong WA, Govind N, Head-Gordon M, Hermes MR, Kowalski K, Li X, Lischka H, Mueller KT, Mutlu E, Niklasson AMN, Pederson MR, Peng B, Shepard R, Valeev EF, van Schilfgaarde M, Vlaisavljevich B, Windus TL, Xantheas SS, Zhang X, Zimmerman PM. A Perspective on Sustainable Computational Chemistry Software Development and Integration. J Chem Theory Comput 2023; 19:7056-7076. [PMID: 37769271 PMCID: PMC10601486 DOI: 10.1021/acs.jctc.3c00419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Indexed: 09/30/2023]
Abstract
The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.
Collapse
Affiliation(s)
- Rosa Di Felice
- Departments
of Physics and Astronomy and Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, United States
- CNR-NANO
Modena, Modena 41125, Italy
| | - Maricris L. Mayes
- Department
of Chemistry and Biochemistry, University
of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747, United States
| | | | | | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Wibe A. de Jong
- Lawrence
Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Niranjan Govind
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Martin Head-Gordon
- Pitzer Center
for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Karol Kowalski
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409, United States
| | - Karl T. Mueller
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced
Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Anders M. N. Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mark R. Pederson
- Department
of Physics, The University of Texas at El
Paso, El Paso, Texas 79968, United States
| | - Bo Peng
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| | - Edward F. Valeev
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - Bess Vlaisavljevich
- Department
of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Theresa L. Windus
- Department
of Chemistry, Iowa State University and
Ames Laboratory, Ames, Iowa 50011, United States
| | - Sotiris S. Xantheas
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced
Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Xing Zhang
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Paul M. Zimmerman
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
13
|
Su R, Huang Z. A brand-new type of excited-state proton transfer (ESIPT) molecule based on sulfoxide/sulfenic acid tautomerism. Phys Chem Chem Phys 2023; 25:27566-27573. [PMID: 37807837 DOI: 10.1039/d3cp02624h] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/10/2023]
Abstract
The excited-state proton transfer (ESIPT) behavior of organic fluorophores has attracted much attention due to their unique photophysical properties. So far, ESIPT studies have mainly focused on the transfer of hydrogen atoms between N-N, N-O, or O-O. In this work, a brand-new type of ESIPT molecule based on sulfoxide/sulfenic acid tautomerism has been thoroughly investigated. The sulfoxide/sulfenic acid tautomerization process requires one step and two steps in the ground and first excited singlet states, respectively. A range of density functional theory and time-dependent density functional theory methods have been employed to investigate these structures, and the changes in aromaticity may be responsible for obtaining the ESIPT process. This work presents a novel ESIPT process, showcasing molecules that exhibit distinctive properties compared to conventional ESIPT compounds. These findings are expected to expand the horizons of experimental research in ESIPT.
Collapse
Affiliation(s)
- Rongchuan Su
- Department of Pharmacology, North Sichuan Medical College, Nanchong, 637100, China.
| | - Zhenmei Huang
- College of Chemistry, Sichuan University, Chengdu, Sichuan, 610064, China
| |
Collapse
|
14
|
Li TE, Hammes-Schiffer S. Nuclear-Electronic Orbital Quantum Dynamics of Plasmon-Driven H 2 Photodissociation. J Am Chem Soc 2023; 145:18210-18214. [PMID: 37555733 DOI: 10.1021/jacs.3c04927] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2023]
Abstract
Leveraging localized surface plasmon resonances of metal nanoparticles to trigger chemical reactions is a promising approach for heterogeneous catalysis. First-principles modeling of such processes is challenging due to the large number of electrons and electronic excited states as well as the significance of nuclear quantum effects when hydrogen is involved. Herein, the nonadiabatic nuclear-electronic quantum dynamics of plasmon-induced H2 photodissociation near an Al13- cluster is simulated with real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT). This approach propagates the nonequilibrium quantum dynamics of both electrons and protons. The plasmonic oscillations are shown to inject hot electrons into the antibonding orbital of H2, thereby inducing H2 dissociation. The quantum mechanical treatment of the hydrogen nuclei leads to faster H2 photodissociation and slightly larger isotope effects. Analysis of the nonequilibrium electronic density suggests that these findings stem from enhanced excited-state electronic coupling between the plasmonic mode and the H2 antibonding orbital due to proton delocalization or zero-point energy effects. Given the low computational overhead for including nuclear quantum effects with the RT-NEO-TDDFT approach, this work paves the way for simulating nonadiabatic nuclear-electronic quantum dynamics in other plasmonic systems.
Collapse
Affiliation(s)
- Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
15
|
Kovyrshin A, Skogh M, Tornberg L, Broo A, Mensa S, Sahin E, Symons BCB, Crain J, Tavernelli I. Nonadiabatic Nuclear-Electron Dynamics: A Quantum Computing Approach. J Phys Chem Lett 2023; 14:7065-7072. [PMID: 37527463 DOI: 10.1021/acs.jpclett.3c01589] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/03/2023]
Abstract
Coupled quantum electron-nuclear dynamics is often associated with the Born-Huang expansion of the molecular wave function and the appearance of nonadiabatic effects as a perturbation. On the other hand, native multicomponent representations of electrons and nuclei also exist, which do not rely on any a priori approximation. However, their implementation is hampered by prohibitive scaling. Consequently, quantum computers offer a unique opportunity for extending their use to larger systems. Here, we propose a quantum algorithm for simulating the time-evolution of molecular systems and apply it to proton transfer dynamics in malonaldehyde, described as a rigid scaffold. The proposed quantum algorithm can be easily generalized to include the explicit dynamics of the classically described molecular scaffold. We show how entanglement between electronic and nuclear degrees of freedom can persist over long times if electrons do not follow the nuclear displacement adiabatically. The proposed quantum algorithm may become a valid candidate for the study of such phenomena when sufficiently powerful quantum computers become available.
Collapse
Affiliation(s)
- Arseny Kovyrshin
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca Gothenburg, Pepparedsleden 1, Molndal SE-431 83, Sweden
| | - Mårten Skogh
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca Gothenburg, Pepparedsleden 1, Molndal SE-431 83, Sweden
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology, 412 96 Gothenburg, Sweden
| | - Lars Tornberg
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca Gothenburg, Pepparedsleden 1, Molndal SE-431 83, Sweden
| | - Anders Broo
- Data Science and Modelling, Pharmaceutical Sciences, R&D, AstraZeneca Gothenburg, Pepparedsleden 1, Molndal SE-431 83, Sweden
| | - Stefano Mensa
- The Hartree Centre, STFC, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Emre Sahin
- The Hartree Centre, STFC, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Benjamin C B Symons
- The Hartree Centre, STFC, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Jason Crain
- IBM Research Europe, Hartree Centre STFC Laboratory, Sci-Tech Daresbury, Warrington WA4 4AD, United Kingdom
- Department of Biochemistry, University of Oxford, Oxford OX1 3QU, U.K
| | - Ivano Tavernelli
- IBM Quantum, IBM Research Europe-Zurich, Säumerstrasse 4, 8803 Rüschlikon, Switzerland
| |
Collapse
|
16
|
Xu X. Constrained Nuclear-Electronic Orbital Density Functional Theory with a Dielectric Continuum Solvent Model. J Phys Chem A 2023. [PMID: 37470267 DOI: 10.1021/acs.jpca.3c02507] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/21/2023]
Abstract
Solvent effects are crucial for simulating chemical and biological processes in solutions. The continuum solvation model is widely used for incorporating solvent effects with different levels of theoretical descriptions of solutes. For solutes and solutions containing hydrogen atoms, nuclear quantum effects can also be nonnegligible for reliable simulations. In this work, we couple our recently developed constrained nuclear-electronic orbital density functional theory with a dielectric continuum solvation model to cover nuclear quantum effects and solvent effects simultaneously. This approach is applied to the formate ion, where an anomalous solvatochromic shift in C-H stretch frequency was reported in experiments. By using this new approach to account for nuclear quantum effects and solvent effects, we show that the vibrational frequency of the C-H stretch and the solvatochromic shift are accurately described.
Collapse
Affiliation(s)
- Xi Xu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| |
Collapse
|
17
|
Dickinson JA, Yu Q, Hammes-Schiffer S. Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory for Multiple Proton Transfer Processes. J Phys Chem Lett 2023:6170-6178. [PMID: 37379485 DOI: 10.1021/acs.jpclett.3c01422] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Proton transfer and hydrogen tunneling play pivotal roles in many chemical and biological processes. The nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) approach was developed to describe hydrogen tunneling systems within the multicomponent NEO framework, where the transferring proton is quantized and treated with molecular orbital techniques on the same level as the electrons. Herein, the NEO-MSDFT framework is generalized to an arbitrary number of quantum protons to allow applications to systems involving the transfer and tunneling of multiple protons. The generalized NEO-MSDFT approach is shown to produce delocalized, bilobal proton densities and accurate tunneling splittings for fixed geometries of the formic acid dimer and asymmetric substituted variants, as well as the porphycene molecule. Investigation of a protonated water chain highlights the applicability of this approach to proton relay systems. This work provides the foundation for nuclear-electronic quantum dynamics simulations of a wide range of multiple proton transfer processes.
Collapse
Affiliation(s)
- Joseph A Dickinson
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
18
|
Chow M, Lambros E, Li X, Hammes-Schiffer S. Nuclear-Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics. J Chem Theory Comput 2023. [PMID: 37329317 DOI: 10.1021/acs.jctc.3c00361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear-electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. These developments and initial examples provide the foundation for future studies of nuclear-electronic quantum dynamics in complex chemical and biological environments.
Collapse
Affiliation(s)
- Mathew Chow
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Eleftherios Lambros
- 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
|
19
|
Hammes-Schiffer S. Exploring Proton-Coupled Electron Transfer at Multiple Scales. NATURE COMPUTATIONAL SCIENCE 2023; 3:291-300. [PMID: 37577057 PMCID: PMC10416817 DOI: 10.1038/s43588-023-00422-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Accepted: 02/23/2023] [Indexed: 08/15/2023]
Abstract
The coupling of electron and proton transfer is critical for chemical and biological processes spanning a wide range of length and time scales and often occurring in complex environments. Thus, diverse modeling strategies, including analytical theories, quantum chemistry, molecular dynamics, and kinetic modeling, are essential for a comprehensive understanding of such proton-coupled electron transfer reactions. Each of these computational methods provides one piece of the puzzle, and all these pieces must be viewed together to produce the full picture.
Collapse
|
20
|
Lambros E, Link B, Chow M, Hammes-Schiffer S, Li X. Solvent Induced Proton Polarization within the Nuclear-Electronic Orbital Framework. J Phys Chem Lett 2023; 14:2990-2995. [PMID: 36942912 DOI: 10.1021/acs.jpclett.3c00471] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
To explicitly account for nuclear quantum effects and solvent environments in simulations of chemical processes, the nuclear-electronic orbital approach is coupled with a polarizable continuum model (PCM). This NEO-PCM approach is used to explore the influence of solvation on nuclear polarization through applications to a water dimer and a set of protonated water tetramers. Nuclear polarization in these species is analyzed in terms of changes in proton density and oxygen-hydrogen bond length. Solvation is shown to enhance nuclear polarization with increasing dielectric constant. For the water dimer, the internal, hydrogen-bonded proton is shown to polarize more than the external, free proton. Moreover, proton quantization leads to greater solvent polarization through their mutual polarization. These calculations highlight the complex interplay among electronic, nuclear, and solvent polarization in chemical systems.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Benjamin Link
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Mathew Chow
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
21
|
Li TE, Hammes-Schiffer S. Electronic Born-Oppenheimer approximation in nuclear-electronic orbital dynamics. J Chem Phys 2023; 158:114118. [PMID: 36948810 DOI: 10.1063/5.0142007] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/24/2023] Open
Abstract
Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born-Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.
Collapse
Affiliation(s)
- Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | | |
Collapse
|
22
|
Culpitt T, Peters LDM, Tellgren EI, Helgaker T. Time-dependent nuclear-electronic orbital Hartree-Fock theory in a strong uniform magnetic field. J Chem Phys 2023; 158:114115. [PMID: 36948801 DOI: 10.1063/5.0139675] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2023] Open
Abstract
In an ultrastrong magnetic field, with field strength B ≈ B0 = 2.35 × 105 T, molecular structure and dynamics differ strongly from that observed on the Earth. Within the Born-Oppenheimer (BO) approximation, for example, frequent (near) crossings of electronic energy surfaces are induced by the field, suggesting that nonadiabatic phenomena and processes may play a more important role in this mixed-field regime than in the weak-field regime on Earth. To understand the chemistry in the mixed regime, it therefore becomes important to explore non-BO methods. In this work, the nuclear-electronic orbital (NEO) method is employed to study protonic vibrational excitation energies in the presence of a strong magnetic field. The NEO generalized Hartree-Fock theory and time-dependent Hartree-Fock (TDHF) theory are derived and implemented, accounting for all terms that result as a consequence of the nonperturbative treatment of molecular systems in a magnetic field. The NEO results for HCN and FHF- with clamped heavy nuclei are compared against the quadratic eigenvalue problem. Each molecule has three semi-classical modes owing to the hydrogen-two precession modes that are degenerate in the absence of a field and one stretching mode. The NEO-TDHF model is found to perform well; in particular, it automatically captures the screening effects of the electrons on the nuclei, which are quantified through the difference in energy of the precession modes.
Collapse
Affiliation(s)
- Tanner Culpitt
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Laurens D M Peters
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Erik I Tellgren
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
23
|
Samsonova I, Tucker GB, Alaal N, Brorsen KR. Hydrogen-Atom Electronic Basis Sets for Multicomponent Quantum Chemistry. ACS OMEGA 2023; 8:5033-5041. [PMID: 36777583 PMCID: PMC9910068 DOI: 10.1021/acsomega.2c07782] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 01/13/2023] [Indexed: 06/18/2023]
Abstract
Multicomponent methods are a conceptually simple way to include nuclear quantum effects into quantum chemistry calculations. In multicomponent methods, the electronic molecular orbitals are described using the linear combination of atomic orbitals approximation. This requires the selection of a one-particle electronic basis set which, in practice, is commonly a correlation-consistent basis set. In multicomponent method studies, it has been demonstrated that large electronic basis sets are required for quantum hydrogen nuclei to accurately describe electron-nuclear correlation. However, as we show in this study, much of the need for large electronic basis sets is due to the correlation-consistent electronic basis sets not being optimized to describe nuclear properties and electron-nuclear correlation. Herein, we introduce a series of correlation-consistent electronic basis sets for hydrogen atoms called cc-pVnZ-mc with additional basis functions optimized to reproduce multicomponent density functional theory protonic densities. These new electronic basis sets are shown to yield better protonic densities with fewer electronic basis functions than the standard correlation-consistent basis sets and reproduce other protonic properties such as proton affinities and protonic excitation energies, even though they were not optimized for these purposes. The cc-pVnZ-mc basis sets should enable multicomponent many-body calculations on larger systems due to the improved computational efficiency they provide for a given level of accuracy.
Collapse
|
24
|
Chen Z, Yang Y. Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface. J Phys Chem Lett 2023; 14:279-286. [PMID: 36595586 DOI: 10.1021/acs.jpclett.2c02905] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The accurate incorporation of nuclear quantum effects in large-scale molecular dynamics (MD) simulations remains a significant challenge. Recently, we combined constrained nuclear-electronic orbital (CNEO) theory with classical MD and obtained a new approach (CNEO-MD) that can accurately and efficiently incorporate nuclear quantum effects into classical simulations. In this Letter, we provide the theoretical foundation for CNEO-MD by developing an alternative formulation of the equations of motion for MD. In this new formulation, the expectation values of quantum nuclear positions evolve classically on an effective energy surface that is obtained from a constrained energy minimization procedure when solving for the quantum nuclear wave function, thus enabling the incorporation of nuclear quantum effects in classical MD simulations. For comparison with other existing approaches, we examined a series of model systems and found that this new MD approach is significantly more accurate than the conventional way of performing classical MD and generally outperforms centroid MD and ring-polymer MD in describing vibrations in these model systems.
Collapse
Affiliation(s)
- Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| |
Collapse
|
25
|
Yu Q, Roy S, Hammes-Schiffer S. Nonadiabatic Dynamics of Hydrogen Tunneling with Nuclear-Electronic Orbital Multistate Density Functional Theory. J Chem Theory Comput 2022; 18:7132-7141. [PMID: 36378867 DOI: 10.1021/acs.jctc.2c00938] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Proton transfer reactions play a critical role in many chemical and biological processes. The development of computationally efficient approaches to describe the quantum dynamics of proton transfer, which often involves hydrogen tunneling, is challenging. Herein, the nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) method is combined with both Ehrenfest and surface hopping nonadiabatic dynamics methods to describe hydrogen tunneling. The NEO-MSDFT method treats the transferring hydrogen nucleus quantum mechanically on the same level as the electrons and incorporates both static and dynamical correlation by mixing localized NEO-DFT solutions with a nonorthogonal configuration interaction scheme. The other nuclei are propagated on the NEO-MSDFT vibronic surfaces during the Ehrenfest or surface hopping dynamics. These methods are applied to proton transfer in malonaldehyde as a prototypical hydrogen tunneling system. The inclusion of vibronically nonadiabatic effects is found to significantly impact the proton transfer time and tunneling dynamics. This approach is applicable to a wide range of other proton transfer reactions.
Collapse
Affiliation(s)
- Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Saswata Roy
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
26
|
Liu A, Chow M, Wildman A, Frisch MJ, Hammes-Schiffer S, Li X. Simultaneous Optimization of Nuclear-Electronic Orbitals. J Phys Chem A 2022; 126:7033-7039. [PMID: 36154137 DOI: 10.1021/acs.jpca.2c05172] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Accurate modeling of important nuclear quantum effects, such as nuclear delocalization, zero-point energy, and tunneling, as well as non-Born-Oppenheimer effects, requires treatment of both nuclei and electrons quantum mechanically. The nuclear-electronic orbital (NEO) method provides an elegant framework to treat specified nuclei, typically protons, on the same level as the electrons. In conventional electronic structure theory, finding a converged ground state can be a computationally demanding task; converging NEO wavefunctions, due to their coupled electronic and nuclear nature, is even more demanding. Herein, we present an efficient simultaneous optimization method that uses the direct inversion in the iterative subspace method to simultaneously converge wavefunctions for both the electrons and quantum nuclei. Benchmark studies show that the simultaneous optimization method can significantly reduce the computational cost compared to the conventional stepwise method for optimizing NEO wavefunctions for multicomponent systems.
Collapse
Affiliation(s)
- Aodong Liu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Mathew Chow
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Andrew Wildman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J Frisch
- Gaussian Incorporated, 340 Quinnipiac Street, Bldg 40, Wallingford, Connecticut 06492, United States
| | | | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States.,Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
27
|
Hammes-Schiffer S. Theoretical perspectives on non-Born-Oppenheimer effects in chemistry. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200377. [PMID: 35341306 DOI: 10.1098/rsta.2020.0377] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 07/12/2021] [Indexed: 06/14/2023]
Abstract
The Born-Oppenheimer approximation, which assumes that the electrons respond instantaneously to the motion of the nuclei, breaks down for a wide range of chemical and biological processes. The rate constants of such nonadiabatic processes can be calculated using analytical theories, and the real-time nonequilibrium dynamics can be described using numerical atomistic simulations. The selection of an approach depends on the desired balance between accuracy and efficiency. The computational expense of generating potential energy surfaces on-the-fly often favours the use of approximate, robust and efficient methods such as trajectory surface hopping for large, complex systems. The development of formally exact non-Born-Oppenheimer methods and the exploration of well-defined approximations to such methods are critical for providing benchmarks and preparing for the next generation of faster computers. Thus, the parallel development of rigorous but computationally expensive methods and more approximate but computationally efficient methods is optimal. This Perspective briefly summarizes the available theoretical and computational non-Born-Oppenheimer methods and presents examples illustrating how analytical theories and nonadiabatic dynamics simulations can elucidate the fundamental principles of chemical and biological processes. These examples also highlight how theoretical calculations are able to guide the interpretation of experimental data and provide experimentally testable predictions for nonadiabatic processes. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
|
28
|
Xu J, Zhou R, Tao Z, Malbon C, Blum V, Hammes-Schiffer S, Kanai Y. Nuclear-Electronic Orbital Approach to Quantization of Protons in Periodic Electronic Structure Calculations. J Chem Phys 2022; 156:224111. [DOI: 10.1063/5.0088427] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The nuclear-electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born-Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems such as two-dimensional materials and liquid-semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide-water interface. The zero-point energy effects of the protons, as well as electron-proton correlation, are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.
Collapse
Affiliation(s)
- Jianhang Xu
- Chemistry, The University of North Carolina at Chapel Hill, United States of America
| | | | - Zhen Tao
- Yale University, United States of America
| | | | - Volker Blum
- Duke University Department of Mechanical Engineering and Materials Science, United States of America
| | | | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, United States of America
| |
Collapse
|
29
|
Li TE, Tao Z, Hammes-Schiffer S. Semiclassical Real-Time Nuclear-Electronic Orbital Dynamics for Molecular Polaritons: Unified Theory of Electronic and Vibrational Strong Couplings. J Chem Theory Comput 2022; 18:2774-2784. [PMID: 35420037 DOI: 10.1021/acs.jctc.2c00096] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Molecular polaritons have become an emerging platform for remotely controlling molecular properties through strong light-matter interactions. Herein, a semiclassical approach is developed for describing molecular polaritons by self-consistently propagating the real-time dynamics of classical cavity modes and a quantum molecular subsystem described by the nuclear-electronic orbital (NEO) method, where electrons and specified nuclei are treated quantum mechanically on the same level. This semiclassical real-time NEO approach provides a unified description of electronic and vibrational strong couplings and describes the impact of the cavity on coupled nuclear-electronic dynamics while including nuclear quantum effects. For a single o-hydroxybenzaldehyde molecule under electronic strong coupling, this approach shows that the cavity suppression of excited state intramolecular proton transfer is influenced not only by the polaritonic potential energy surface but also by the time scale of the chemical reaction. This work provides the foundation for exploring collective strong coupling in nuclear-electronic quantum dynamical systems within optical cavities.
Collapse
Affiliation(s)
- Tao E Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Zhen Tao
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
30
|
Wildman A, Tao Z, Zhao L, Hammes-Schiffer S, Li X. Solvated Nuclear-Electronic Orbital Structure and Dynamics. J Chem Theory Comput 2022; 18:1340-1346. [PMID: 35179376 DOI: 10.1021/acs.jctc.1c01285] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Nonadiabatic dynamical processes such as proton-coupled electron transfer and excited state intramolecular proton transfer have been the subject of much research. One of the promising theoretical methods to describe these processes is the nuclear-electronic orbital (NEO) approach. This approach inherently accounts for nuclear quantum effects within quantum chemistry calculations, and it has recently been extended to directly simulate nonadiabatic processes with the development of real-time NEO methods. These processes can also be significantly dependent on the surrounding chemical environment, however, and capturing the effects of the environment is often necessary for analyzing experimentally relevant systems. This work couples the NEO density functional theory and real-time time-dependent density functional theory approaches with solvation through the polarizable continuum model. The effects of this coupling are investigated for ground state properties, solvent-dependent vibrational frequencies, and direct excited state intramolecular proton transfer dynamics.
Collapse
Affiliation(s)
- Andrew Wildman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Zhen Tao
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Luning Zhao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States.,Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
31
|
Feldmann R, Muolo A, Baiardi A, Reiher M. Quantum Proton Effects from Density Matrix Renormalization Group Calculations. J Chem Theory Comput 2022; 18:234-250. [PMID: 34978441 DOI: 10.1021/acs.jctc.1c00913] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
We recently introduced [J. Chem. Phys. 2020, 152, 204103] the nuclear-electronic all-particle density matrix renormalization group (NEAP-DMRG) method to solve the molecular Schrödinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG method with the nuclear-electronic Hartree-Fock (NEHF-DMRG) approach, treating nuclei and electrons on the same footing. Inter- and intraspecies correlations are described within the DMRG method without truncating the excitation degree of the full configuration interaction wave function. We extend the concept of orbital entanglement and mutual information to nuclear-electronic wave functions and demonstrate that they are reliable metrics to detect strong correlation effects. We apply the NEHF-DMRG method to the HeHHe+ molecular ion, to obtain accurate proton densities, ground-state total energies, and vibrational transition frequencies by comparison with state-of-the-art data obtained with grid-based approaches and modern configuration interaction methods. For HCN, we improve on the accuracy of the latter approaches with respect to both the ground-state absolute energy and proton density, which is a major challenge for multireference nuclear-electronic state-of-the-art methods.
Collapse
Affiliation(s)
- Robin Feldmann
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Andrea Muolo
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, 91904 Jerusalem, Israel
| | - Alberto Baiardi
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
32
|
Tao Z, Yu Q, Roy S, Hammes-Schiffer S. Direct Dynamics with Nuclear-Electronic Orbital Density Functional Theory. Acc Chem Res 2021; 54:4131-4141. [PMID: 34726895 DOI: 10.1021/acs.accounts.1c00516] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Direct dynamics simulations of chemical reactions typically require the selection of a method for generating the potential energy surfaces and a method for the dynamical propagation of the nuclei on these surfaces. The nuclear-electronic orbital (NEO) framework avoids this Born-Oppenheimer separation by treating specified nuclei on the same level as the electrons with wave function methods or density functional theory (DFT). The NEO approach is particularly applicable to proton, hydride, and proton-coupled electron transfer reactions, where the transferring proton(s) and all electrons are treated quantum mechanically. In this manner, the zero-point energy, density delocalization, and anharmonicity of the transferring protons are inherently and efficiently included in the energies, optimized geometries, and dynamics.This Account describes how various NEO methods can be used for direct dynamics simulations on electron-proton vibronic surfaces. The strengths and limitations of these approaches are discussed, and illustrative examples are presented. The NEO-DFT method can be used to simulate chemical reactions on the ground state vibronic surface, as illustrated by the application to hydride transfer in C4H9+. The NEO multistate DFT (NEO-MSDFT) method is useful for simulating ground state reactions in which the proton density becomes bilobal during the dynamics, a characteristic of hydrogen tunneling, as illustrated by proton transfer in malonaldehyde. The NEO time-dependent DFT (NEO-TDDFT) method produces excited electronic, vibrational, and vibronic surfaces. The application of linear-response NEO-TDDFT to H2 and H3+, as well as the partially and fully deuterated counterparts, shows that this approach produces accurate fundamental vibrational excitation energies when all nuclei and all electrons are treated quantum mechanically. Moreover, when only specified nuclei are treated quantum mechanically, this approach can be used to optimize geometries on excited state vibronic surfaces, as illustrated by photoinduced single and double proton transfer systems, and to conduct adiabatic dynamics on these surfaces. The real-time NEO-TDDFT method provides an alternative approach for simulating nonequilibrium nuclear-electronic dynamics of such systems. These various NEO methods can be combined with nonadiabatic dynamics methods such as Ehrenfest and surface hopping dynamics to include the nonadiabatic effects between the quantum and classical subsystems. The real-time NEO-TDDFT Ehrenfest dynamics simulation of excited state intramolecular proton transfer in o-hydroxybenzaldehyde illustrates the power of this type of combined approach. The field of multicomponent quantum chemistry is in the early stages, and the methods discussed herein provide the foundation for a wide range of promising future directions to be explored. An appealing future direction is the expansion of the real-time NEO-TDDFT method to describe the dynamics of all nuclei and electrons on the same level. Direct dynamics simulations using NEO wave function methods such as equation-of-motion coupled cluster or multiconfigurational approaches are also attractive but computationally expensive options. The further development of NEO direct dynamics methods will enable the simulation of the nuclear-electronic dynamics for a vast array of chemical and biological processes that extend beyond the Born-Oppenheimer approximation.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Qi Yu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Saswata Roy
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | | |
Collapse
|
33
|
Wang C, Ulusoy IS, Aebersold LE, Wilson AK. Multi-configuration electron-nuclear dynamics: An open-shell approach. J Chem Phys 2021; 155:154103. [PMID: 34686063 DOI: 10.1063/5.0063478] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The multi-configuration electron-nuclear dynamics for open shell systems with a spin-unrestricted formalism is described. The mean fields are evaluated using second-order reduced density matrices for electronic and nuclear degrees of freedom. Applications to light-element diatomics including equilibrium geometries, electronic energies, dipole moments, and absorption spectra are presented. The von Neumann entropies for different spin states of a LiH molecule are compared.
Collapse
Affiliation(s)
- Cong Wang
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Inga S Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Lucas E Aebersold
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, USA
| |
Collapse
|
34
|
Shepard C, Zhou R, Yost DC, Yao Y, Kanai Y. Simulating electronic excitation and dynamics with real-time propagation approach to TDDFT within plane-wave pseudopotential formulation. J Chem Phys 2021; 155:100901. [PMID: 34525811 DOI: 10.1063/5.0057587] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We give a perspective on simulating electronic excitation and dynamics using the real-time propagation approach to time-dependent density functional theory (RT-TDDFT) in the plane-wave pseudopotential formulation. RT-TDDFT is implemented in various numerical formalisms in recent years, and its practical application often dictates the most appropriate implementation of the theory. We discuss recent developments and challenges, emphasizing numerical aspects of studying real systems. Several applications of RT-TDDFT simulation are discussed to highlight how the approach is used to study interesting electronic excitation and dynamics phenomena in recent years.
Collapse
Affiliation(s)
- Christopher Shepard
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Dillon C Yost
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02142, USA
| | - Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| |
Collapse
|
35
|
Epifanovsky E, Gilbert ATB, Feng X, Lee J, Mao Y, Mardirossian N, Pokhilko P, White AF, Coons MP, Dempwolff AL, Gan Z, Hait D, Horn PR, Jacobson LD, Kaliman I, Kussmann J, Lange AW, Lao KU, Levine DS, Liu J, McKenzie SC, Morrison AF, Nanda KD, Plasser F, Rehn DR, Vidal ML, You ZQ, Zhu Y, Alam B, Albrecht BJ, Aldossary A, Alguire E, Andersen JH, Athavale V, Barton D, Begam K, Behn A, Bellonzi N, Bernard YA, Berquist EJ, Burton HGA, Carreras A, Carter-Fenk K, Chakraborty R, Chien AD, Closser KD, Cofer-Shabica V, Dasgupta S, de Wergifosse M, Deng J, Diedenhofen M, Do H, Ehlert S, Fang PT, Fatehi S, Feng Q, Friedhoff T, Gayvert J, Ge Q, Gidofalvi G, Goldey M, Gomes J, González-Espinoza CE, Gulania S, Gunina AO, Hanson-Heine MWD, Harbach PHP, Hauser A, Herbst MF, Hernández Vera M, Hodecker M, Holden ZC, Houck S, Huang X, Hui K, Huynh BC, Ivanov M, Jász Á, Ji H, Jiang H, Kaduk B, Kähler S, Khistyaev K, Kim J, Kis G, Klunzinger P, Koczor-Benda Z, Koh JH, Kosenkov D, Koulias L, Kowalczyk T, Krauter CM, Kue K, Kunitsa A, Kus T, Ladjánszki I, Landau A, Lawler KV, Lefrancois D, Lehtola S, Li RR, Li YP, Liang J, Liebenthal M, Lin HH, Lin YS, Liu F, Liu KY, Loipersberger M, Luenser A, Manjanath A, Manohar P, Mansoor E, Manzer SF, Mao SP, Marenich AV, Markovich T, Mason S, Maurer SA, McLaughlin PF, Menger MFSJ, Mewes JM, Mewes SA, Morgante P, Mullinax JW, Oosterbaan KJ, Paran G, Paul AC, Paul SK, Pavošević F, Pei Z, Prager S, Proynov EI, Rák Á, Ramos-Cordoba E, Rana B, Rask AE, Rettig A, Richard RM, Rob F, Rossomme E, Scheele T, Scheurer M, Schneider M, Sergueev N, Sharada SM, Skomorowski W, Small DW, Stein CJ, Su YC, Sundstrom EJ, Tao Z, Thirman J, Tornai GJ, Tsuchimochi T, Tubman NM, Veccham SP, Vydrov O, Wenzel J, Witte J, Yamada A, Yao K, Yeganeh S, Yost SR, Zech A, Zhang IY, Zhang X, Zhang Y, Zuev D, Aspuru-Guzik A, Bell AT, Besley NA, Bravaya KB, Brooks BR, Casanova D, Chai JD, Coriani S, Cramer CJ, Cserey G, DePrince AE, DiStasio RA, Dreuw A, Dunietz BD, Furlani TR, Goddard WA, Hammes-Schiffer S, Head-Gordon T, Hehre WJ, Hsu CP, Jagau TC, Jung Y, Klamt A, Kong J, Lambrecht DS, Liang W, Mayhall NJ, McCurdy CW, Neaton JB, Ochsenfeld C, Parkhill JA, Peverati R, Rassolov VA, Shao Y, Slipchenko LV, Stauch T, Steele RP, Subotnik JE, Thom AJW, Tkatchenko A, Truhlar DG, Van Voorhis T, Wesolowski TA, Whaley KB, Woodcock HL, Zimmerman PM, Faraji S, Gill PMW, Head-Gordon M, Herbert JM, Krylov AI. Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package. J Chem Phys 2021; 155:084801. [PMID: 34470363 PMCID: PMC9984241 DOI: 10.1063/5.0055522] [Citation(s) in RCA: 480] [Impact Index Per Article: 160.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.
Collapse
Affiliation(s)
- Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alec F. White
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marc P. Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Adrian L. Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zhengting Gan
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Diptarka Hait
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Paul R. Horn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Leif D. Jacobson
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Jörg Kussmann
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Adrian W. Lange
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Daniel S. Levine
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Simon C. McKenzie
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Kaushik D. Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Dirk R. Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Marta L. Vidal
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | - Ying Zhu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Bushra Alam
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Benjamin J. Albrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | - Ethan Alguire
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Josefine H. Andersen
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | - Vishikh Athavale
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Dennis Barton
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Khadiza Begam
- Department of Physics, Kent State University, Kent, Ohio 44242, USA
| | - Andrew Behn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Yves A. Bernard
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Hugh G. A. Burton
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Abel Carreras
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Alan D. Chien
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | - Vale Cofer-Shabica
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Marc de Wergifosse
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jia Deng
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Hainam Do
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Beringstr. 4, 53115 Bonn, Germany
| | - Po-Tung Fang
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Qingguo Feng
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Triet Friedhoff
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - James Gayvert
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Qinghui Ge
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, USA
| | - Matthew Goldey
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joe Gomes
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Sahil Gulania
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Anastasia O. Gunina
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Phillip H. P. Harbach
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Hauser
- Institute of Experimental Physics, Graz University of Technology, Graz, Austria
| | | | - Mario Hernández Vera
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Manuel Hodecker
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zachary C. Holden
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Shannon Houck
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Xunkun Huang
- Department of Chemistry, Xiamen University, Xiamen 361005, China
| | - Kerwin Hui
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Bang C. Huynh
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Maxim Ivanov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Ádám Jász
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Hyunjun Ji
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Hanjie Jiang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Benjamin Kaduk
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Sven Kähler
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Kirill Khistyaev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jaehoon Kim
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Gergely Kis
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | | | - Zsuzsanna Koczor-Benda
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Joong Hoon Koh
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Dimitri Kosenkov
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Laura Koulias
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | | | - Caroline M. Krauter
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Karl Kue
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Alexander Kunitsa
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Thomas Kus
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Arie Landau
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Keith V. Lawler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | | | - Run R. Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Yi-Pei Li
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jiashu Liang
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marcus Liebenthal
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Hung-Hsuan Lin
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - You-Sheng Lin
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Fenglai Liu
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Arne Luenser
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Aaditya Manjanath
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Prashant Manohar
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Erum Mansoor
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Sam F. Manzer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Shan-Ping Mao
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Thomas Markovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Stephen Mason
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Simon A. Maurer
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Peter F. McLaughlin
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | - Jan-Michael Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Stefanie A. Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Pierpaolo Morgante
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - J. Wayne Mullinax
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | | | | | - Alexander C. Paul
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Suranjan K. Paul
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Zheng Pei
- School of Electrical and Computer Engineering, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Stefan Prager
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Emil I. Proynov
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Ádám Rák
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Eloy Ramos-Cordoba
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Bhaskar Rana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Alan E. Rask
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Ryan M. Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fazle Rob
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Elliot Rossomme
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Tarek Scheele
- Institute for Physical and Theoretical Chemistry, University of Bremen, Bremen, Germany
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Matthias Schneider
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Nickolai Sergueev
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Shaama M. Sharada
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Wojciech Skomorowski
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Christopher J. Stein
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yu-Chuan Su
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Eric J. Sundstrom
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Zhen Tao
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Jonathan Thirman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gábor J. Tornai
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Takashi Tsuchimochi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Norm M. Tubman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Oleg Vydrov
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Jan Wenzel
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Jon Witte
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Atsushi Yamada
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Kun Yao
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Sina Yeganeh
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Shane R. Yost
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Alexander Zech
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - Igor Ying Zhang
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Xing Zhang
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Yu Zhang
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Dmitry Zuev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Alexis T. Bell
- Department of Chemical Engineering, University of California, Berkeley, California 94720, USA
| | - Nicholas A. Besley
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Ksenia B. Bravaya
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Bernard R. Brooks
- Laboratory of Computational Biophysics, National Institute of Health, Bethesda, Maryland 20892, USA
| | - David Casanova
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | | | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | | | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Thomas R. Furlani
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - William A. Goddard
- Materials and Process Simulation Center, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | | | | | - Yousung Jung
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Andreas Klamt
- COSMOlogic GmbH & Co. KG, Imbacher Weg 46, D-51379 Leverkusen, Germany
| | - Jing Kong
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Daniel S. Lambrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | | - C. William McCurdy
- Department of Chemistry, University of California, Davis, California 95616, USA
| | - Jeffrey B. Neaton
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Christian Ochsenfeld
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - John A. Parkhill
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Roberto Peverati
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - Vitaly A. Rassolov
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA
| | | | | | | | - Ryan P. Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Donald G. Truhlar
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - K. Birgitta Whaley
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - H. Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, USA
| | - Paul M. Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Shirin Faraji
- Zernike Institute for Advanced Materials, University of Groningen, 9774AG Groningen, The Netherlands
| | | | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA,Author to whom correspondence should be addressed:
| |
Collapse
|
36
|
Hammes-Schiffer S. Nuclear-electronic orbital methods: Foundations and prospects. J Chem Phys 2021; 155:030901. [PMID: 34293877 DOI: 10.1063/5.0053576] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
The incorporation of nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations and molecular dynamics simulations is a longstanding challenge. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons with wave function and density functional theory methods. This approach inherently includes nuclear delocalization and zero-point energy in molecular energy calculations, geometry optimizations, reaction paths, and dynamics. It can also provide accurate descriptions of excited electronic, vibrational, and vibronic states as well as nuclear tunneling and nonadiabatic dynamics. Nonequilibrium nuclear-electronic dynamics simulations beyond the Born-Oppenheimer approximation can be used to investigate a wide range of excited state processes. This Perspective provides an overview of the foundational NEO methods and enumerates the prospects for using these methods as building blocks for future developments. The conceptual simplicity and computational efficiency of the NEO approach will enhance its accessibility and applicability to diverse chemical and biological systems.
Collapse
|
37
|
Udagawa T, Murphy RB, Darwish TA, Tachikawa M, Mori S. H/D Isotope Effects in Keto-Enol Tautomerism of β-Dicarbonyl Compounds —Importance of Nuclear Quantum Effects of Hydrogen Nuclei—. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2021. [DOI: 10.1246/bcsj.20210083] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Taro Udagawa
- Department of Chemistry and Biomolecular Science, Faculty of Engineering, Gifu University, 1-1 Yanagido, Gifu 501-1193, Japan
| | - Rhys B. Murphy
- Australian Nuclear Science and Technology Organisation, New Illawarra Rd, Lucas Heights NSW 2234, Australia
| | - Tamim A. Darwish
- Australian Nuclear Science and Technology Organisation, New Illawarra Rd, Lucas Heights NSW 2234, Australia
| | - Masanori Tachikawa
- Graduate School of NanobioScience, Yokohama City University, 22-2 Seto, Kanazawa-ku, Yokohama, Kanagawa 236-0027, Japan
| | - Seiji Mori
- Institute of Quantum Beam Science, Graduate School of Science and Engineering, Ibaraki University, 2-1-1 Bunkyo, Mito, Ibaraki 310-8512, Japan
- Frontier Research Center for Applied Atomic Sciences, Ibaraki University, 162-1 Shirakata, Tokai, Ibaraki 319-1106, Japan
| |
Collapse
|
38
|
Tao Z, Roy S, Schneider PE, Pavošević F, Hammes-Schiffer S. Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited-State Geometry Optimizations and Adiabatic Excitation Energies. J Chem Theory Comput 2021; 17:5110-5122. [PMID: 34260237 DOI: 10.1021/acs.jctc.1c00454] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The computational investigation of photochemical processes often entails the calculation of excited-state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients as well as the NEO-DFT analytical Hessian are provided. The NEO approach includes the anharmonic zero-point energy (ZPE) and density delocalization associated with the quantum protons as well as vibronic mixing in geometry optimizations and energy calculations of ground and excited states. The harmonic ZPE associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited-state intramolecular proton-transfer systems, [2,2'-bipyridyl]-3-ol and [2,2'-bipyridyl]-3,3'-diol, are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited-state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Saswata Roy
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Patrick E Schneider
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
39
|
Barone V, Alessandrini S, Biczysko M, Cheeseman JR, Clary DC, McCoy AB, DiRisio RJ, Neese F, Melosso M, Puzzarini C. Computational molecular spectroscopy. ACTA ACUST UNITED AC 2021. [DOI: 10.1038/s43586-021-00034-1] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
40
|
Pavošević F, Hammes-Schiffer S. Multicomponent Unitary Coupled Cluster and Equation-of-Motion for Quantum Computation. J Chem Theory Comput 2021; 17:3252-3258. [PMID: 33945684 DOI: 10.1021/acs.jctc.1c00220] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The variational quantum eigensolver (VQE) algorithm combined with the unitary coupled cluster (UCC) ansatz has been developed for the quantum computation of molecular energies and wave functions within the Born-Oppenheimer approximation. Herein, this approach is extended to multicomponent systems to enable the quantum mechanical treatment of more than one type of particle, such as electrons and positrons or electrons and nuclei, without invoking the Born-Oppenheimer approximation. Specifically, we introduce the multicomponent unitary coupled cluster (mcUCC) method combined with the VQE algorithm for the calculation of ground-state energies and wave functions as well as the multicomponent equation-of-motion (mcEOM) method for the calculation of excitation energies. These methods are developed within the nuclear-electronic orbital (NEO) framework and are formulated in the qubit basis to enable implementations on quantum computers. Moreover, these methods are used to calculate the ground-state energy and excitation energies of positronium hydride, where both electrons and the positron are treated quantum mechanically, as well as the H2 molecule, where both electrons and one proton are treated quantum mechanically. These applications validate the implementation and provide benchmark data for future calculations. The errors due to Trotterization of the mcUCC ansatz are also analyzed. This formalism, as well as the accompanying computer code, will serve as the basis for applications to more complex multicomponent systems, such as simulations of photoinduced nonadiabatic molecular processes, on both classical and quantum computers.
Collapse
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
41
|
Omar KA, Hasnaoui K, de la Lande A. First-Principles Simulations of Biological Molecules Subjected to Ionizing Radiation. Annu Rev Phys Chem 2021; 72:445-465. [PMID: 33878897 DOI: 10.1146/annurev-physchem-101419-013639] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Ionizing rays cause damage to genomes, proteins, and signaling pathways that normally regulate cell activity, with harmful consequences such as accelerated aging, tumors, and cancers but also with beneficial effects in the context of radiotherapies. While the great pace of research in the twentieth century led to the identification of the molecular mechanisms for chemical lesions on the building blocks of biomacromolecules, the last two decades have brought renewed questions, for example, regarding the formation of clustered damage or the rich chemistry involving the secondary electrons produced by radiolysis. Radiation chemistry is now meeting attosecond science, providing extraordinary opportunities to unravel the very first stages of biological matter radiolysis. This review provides an overview of the recent progress made in this direction, focusing mainly on the atto- to femto- to picosecond timescales. We review promising applications of time-dependent density functional theory in this context.
Collapse
Affiliation(s)
- Karwan Ali Omar
- Institut de Chimie Physique, CNRS UMR 8000, Université Paris-Saclay, 91405 Orsay, France; .,Department of Chemistry, College of Education, University of Sulaimani, 41005 Kurdistan, Iraq
| | - Karim Hasnaoui
- High Performance Computing User Support Team, Institut du Développement et des Ressources en Informatique Scientifique (IDRIS), 91403 Orsay, France.,Maison de la Simulation, CNRS, Commissariat à l'Energie Atomique et aux Énergies Alternatives (CEA), Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Aurélien de la Lande
- Institut de Chimie Physique, CNRS UMR 8000, Université Paris-Saclay, 91405 Orsay, France;
| |
Collapse
|
42
|
Zhao L, Wildman A, Pavošević F, Tully JC, Hammes-Schiffer S, Li X. Excited State Intramolecular Proton Transfer with Nuclear-Electronic Orbital Ehrenfest Dynamics. J Phys Chem Lett 2021; 12:3497-3502. [PMID: 33792317 DOI: 10.1021/acs.jpclett.1c00564] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The recent development of the Ehrenfest dynamics approach in the nuclear-electronic orbital (NEO) framework provides a promising way to simulate coupled nuclear-electronic dynamics. Our previous study showed that the NEO-Ehrenfest approach with a semiclassical traveling proton basis method yields accurate predictions of molecular vibrational frequencies. In this work, we provide a more thorough analysis of the semiclassical traveling proton basis method to elucidate its validity and convergence behavior. We also conduct NEO-Ehrenfest dynamics simulations to study an excited state intramolecular proton transfer process. These simulations reveal that nuclear quantum effects influence the predictions of proton transfer reaction rates and kinetic isotope effects due to the intrinsic delocalized nature of the quantum nuclear wave function. This work illustrates the importance of nuclear quantum effects in coupled nuclear-electronic dynamical processes and shows that the NEO-Ehrenfest approach can be a powerful tool for providing insights and predictions for these processes.
Collapse
Affiliation(s)
- Luning Zhao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew Wildman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - John C Tully
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
43
|
Zhao L, Wildman A, Tao Z, Schneider P, Hammes-Schiffer S, Li X. Nuclear–electronic orbital Ehrenfest dynamics. J Chem Phys 2020; 153:224111. [DOI: 10.1063/5.0031019] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Luning Zhao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Andrew Wildman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Zhen Tao
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Patrick Schneider
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
44
|
Paz JL, Ruiz-Hinojosa E, Alvarado Y, Lascano L, Fernández L, Espinoza-Montero P, Costa-Vera C. Ecuaciones de Bloch Ópticas en Sistemas Complejos con Acoplamiento Intramolecular. REVISTA POLITÉCNICA 2020. [DOI: 10.33333/rp.vol46n2.03] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
Proponemos modificaciones de las ecuaciones de Bloch ópticas convencionales para un sistema molecular, cuando consideramos los efectos de acoplamiento intramolecular. Modelamos la molécula aislada como curvas de energía de Born-Oppenheimer que consisten en dos estados electrónicos cruzados descritos como potenciales armónicos, con los mínimos desplazados en coordenadas nucleares y energía. Consideramos dos estados vibracionales y una perturbación, que puede surgir de una correlación residual electrón-electrón y/o términos de acoplamiento spin-órbita en el Hamiltoniano del sistema, causando la separación de las dos curvas según la regla del cruce evitado. Las ecuaciones extendidas de Bloch ópticas permiten establecer la dinámica de un sistema molecular de estados adiabáticos sujeto a la interacción de un reservorio térmico e interactuando con campos electromagnéticos. En las ecuaciones se observa que haciendo nulo el factor de acoplamiento intramolecular, se recuperan las ecuaciones de Bloch ópticas convencionales para un sistema molecular con potenciales armónicos cuyos mínimos se encuentran exactamente en la misma coordenada nuclear sujeto a la aproximación Born-Oppenheimer.
Collapse
Affiliation(s)
| | | | - Ysaias Alvarado
- Instituto Venezolano de Investigaciones Científicas (IVIC) Región Zulia
| | | | | | | | | |
Collapse
|
45
|
DiRisio RJ, Jones CM, Ma H, Rousseau BJG. Viewpoints on the 2020 Virtual Conference on Theoretical Chemistry. J Phys Chem A 2020; 124:8875-8883. [PMID: 33054223 DOI: 10.1021/acs.jpca.0c08955] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Ryan J DiRisio
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Chey M Jones
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - He Ma
- Institute for Molecular engineering, University of Chicago, 5640 S. Ellis Avenue, Chicago, Illinois 60637, United States
| | - Benjamin J G Rousseau
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
46
|
Sibaev M, Polyak I, Manby FR, Knowles PJ. Molecular second-quantized Hamiltonian: Electron correlation and non-adiabatic coupling treated on an equal footing. J Chem Phys 2020; 153:124102. [DOI: 10.1063/5.0018930] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Marat Sibaev
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Iakov Polyak
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| | - Frederick R. Manby
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| | - Peter J. Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
47
|
Goings JJ, Hammes-Schiffer S. Nonequilibrium Dynamics of Proton-Coupled Electron Transfer in Proton Wires: Concerted but Asynchronous Mechanisms. ACS CENTRAL SCIENCE 2020; 6:1594-1601. [PMID: 32999935 PMCID: PMC7517869 DOI: 10.1021/acscentsci.0c00756] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/09/2020] [Indexed: 05/29/2023]
Abstract
The coupling between electrons and protons and the long-range transport of protons play important roles throughout biology. Biomimetic systems derived from benzimidazole-phenol (BIP) constructs have been designed to undergo proton-coupled electron transfer (PCET) upon electrochemical or photochemical oxidation. Moreover, these systems can transport protons along hydrogen-bonded networks or proton wires through multiproton PCET. Herein, the nonequilibrium dynamics of both single and double proton transfer in BIP molecules initiated by oxidation are investigated with first-principles molecular dynamics simulations. Although these processes are concerted in that no thermodynamically stable intermediate is observed, the simulations predict that they are predominantly asynchronous on the ultrafast time scale. For both systems, the first proton transfer typically occurs ∼100 fs after electron transfer. For the double proton transfer system, typically the second proton transfer occurs hundreds of femtoseconds after the initial proton transfer. A machine learning algorithm was used to identify the key molecular vibrational modes essential for proton transfer: a slow, in-plane bending mode that dominates the overall inner-sphere reorganization, the proton donor-acceptor motion that leads to vibrational coherence, and the faster donor-hydrogen stretching mode. The asynchronous double proton transfer mechanism can be understood in terms of a significant mode corresponding to the two anticorrelated proton donor-acceptor motions, typically decreasing only one donor-acceptor distance at a time. Although these PCET processes appear concerted on the time scale of typical electrochemical experiments, attaching these BIP constructs to photosensitizers may enable the detection of the asynchronicity of the electron and multiple proton transfers with ultrafast two-dimensional spectroscopy. Understanding the fundamental PCET mechanisms at this level will guide the design of PCET systems for catalysis and energy conversion processes.
Collapse
Affiliation(s)
- Joshua J. Goings
- Department of Chemistry, Yale University, 225
Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225
Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
48
|
McLaurin PM, Morales JA. Electron nuclear dynamics of H + + CO 2 (000) → H + + CO 2 ( v1v2v3) at ELab = 20.5-30 eV with coherent-states quantum reconstruction procedure. Phys Chem Chem Phys 2020; 22:19549-19559. [PMID: 32869775 PMCID: PMC7501211 DOI: 10.1039/d0cp03722b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2023]
Abstract
The simplest-level electron nuclear dynamics (SLEND) method with the coherent-states (CSs) quantum reconstruction procedure (CSQRP) is applied to the scattering system H+ + CO2 (000) → H+ + CO2 (v1v2v3) at ELab = 20.5-30 eV. Relevant for astrophysics, atmospheric chemistry and proton cancer therapy, this system undergoes collision-induced vibrational excitations in CO2. SLEND is a time-dependent, variational, direct, and non-adiabatic method that adopts a classical-mechanics description for nuclei and a single-determinantal wavefunction for electrons. The CSQRP employs the canonical CS to reconstruct quantum state-to-state vibrational properties from the SLEND classical nuclear dynamics. Overall, the calculated collision-induced vibrational properties agree well with experimental data. SLEND total differential cross sections (DCSs) agree remarkably well with their experimental counterparts and accurately display rainbow scattering angles structures. SLEND averaged target excitation energies for vibrational + rotational and rotational motions exhibit reasonable and good agreements with experimental data, respectively. These properties show that rotational excitation is low and that the asymmetric stretch normal mode of CO2 is much more excited than the others. SLEND/CSQRP state-to-state vibrational DCSs agree reasonably well with the sparse experimental data for final states v1v2v3 = 000-002, but less satisfactorily for 003. These DCSs also accurately display rainbow scattering angles structures. Finally, SLEND/CSQRP vibrational proton energy loss spectra agree remarkably well with their experimental counterparts for various final vibrational states of CO2, collisions energies and scattering angles. Present results demonstrate the accuracy of SLEND/CSQRP to predict state-to-state vibrational properties in scattering systems with multiple normal modes.
Collapse
Affiliation(s)
- Patrick M McLaurin
- Department of Chemistry and Biochemistry, Texas Tech University, PO Box 41061, Lubbock, TX 79409-1061, USA.
| | - Jorge A Morales
- Department of Chemistry and Biochemistry, Texas Tech University, PO Box 41061, Lubbock, TX 79409-1061, USA.
| |
Collapse
|
49
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
50
|
Pavošević F, Tao Z, Culpitt T, Zhao L, Li X, Hammes-Schiffer S. Frequency and Time Domain Nuclear-Electronic Orbital Equation-of-Motion Coupled Cluster Methods: Combination Bands and Electronic-Protonic Double Excitations. J Phys Chem Lett 2020; 11:6435-6442. [PMID: 32658486 DOI: 10.1021/acs.jpclett.0c01891] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The accurate description of excited vibronic states is important for modeling a wide range of photoinduced processes. The nuclear-electronic orbital (NEO) approach, which treats specified protons on the same level as the electrons, can describe excited electronic-protonic states. Herein the multicomponent equation-of-motion coupled cluster with singles and doubles (NEO-EOM-CCSD) method and its time-domain counterpart, TD-NEO-EOM-CCSD, are developed and implemented. The application of these methods to the HCN molecule highlights their capabilities. These methods predict qualitatively reasonable energies and intensities for a combination band corresponding to simultaneous excitation of two vibrational modes, as well as an overtone. These methods also describe states with double excitation character, such as excited electronic-protonic states corresponding to the simultaneous excitation of an electron and a proton. The ability of the NEO-EOM-CCSD method and its time-dependent counterpart to describe combination bands, overtones, and double excitations will enable a wide range of photochemical applications.
Collapse
Affiliation(s)
- Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Zhen Tao
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Tanner Culpitt
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Luning Zhao
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|