101
|
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
|
102
|
Talbot JJ, Head-Gordon M, Cotton SJ. The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2153761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Justin J. Talbot
- Department of Chemistry, University of California, Berkeley, CA, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| | - Stephen J. Cotton
- Department of Chemistry, University of California, Berkeley, CA, USA
| |
Collapse
|
103
|
Zhang W, Huang W, Tan J, Guo Q, Wu B. Heterogeneous catalysis mediated by light, electricity and enzyme via machine learning: Paradigms, applications and prospects. CHEMOSPHERE 2022; 308:136447. [PMID: 36116627 DOI: 10.1016/j.chemosphere.2022.136447] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/09/2022] [Revised: 09/08/2022] [Accepted: 09/11/2022] [Indexed: 06/15/2023]
Abstract
Energy crisis and environmental pollution have become the bottleneck of human sustainable development. Therefore, there is an urgent need to develop new catalysts for energy production and environmental remediation. Due to the high cost caused by blind screening and limited valuable computing resources, the traditional experimental methods and theoretical calculations are difficult to meet with the requirements. In the past decades, computer science has made great progress, especially in the field of machine learning (ML). As a new research paradigm, ML greatly accelerates the theoretical calculation methods represented by first principal calculation and molecular dynamics, and establish the physical picture of heterogeneous catalytic processes for energy and environment. This review firstly summarized the general research paradigms of ML in the discovery of catalysts. Then, the latest progresses of ML in light-, electricity- and enzyme-mediated heterogeneous catalysis were reviewed from the perspective of catalytic performance, operating conditions and reaction mechanism. The general guidelines of ML for heterogeneous catalysis were proposed. Finally, the existing problems and future development trend of ML in heterogeneous catalysis mediated by light, electricity and enzyme were summarized. We highly expect that this review will facilitate the interaction between ML and heterogeneous catalysis, and illuminate the development prospect of heterogeneous catalysis.
Collapse
Affiliation(s)
- Wentao Zhang
- Tsinghua Shenzhen International Graduate School, Tsinghua University, Shenzhen, 518055, People's Republic of China
| | - Wenguang Huang
- South China Institute of Environmental Sciences, Ministry of Ecology and Environment of PRC, Guangzhou, 510655, People's Republic of China.
| | - Jie Tan
- South China Institute of Environmental Sciences, Ministry of Ecology and Environment of PRC, Guangzhou, 510655, People's Republic of China
| | - Qingwei Guo
- South China Institute of Environmental Sciences, Ministry of Ecology and Environment of PRC, Guangzhou, 510655, People's Republic of China
| | - Bingdang Wu
- School of Environmental Science and Engineering, Suzhou University of Science and Technology, Suzhou, 215009, People's Republic of China; Key Laboratory of Suzhou Sponge City Technology, Suzhou, 215002, People's Republic of China.
| |
Collapse
|
104
|
Rajpurohit S, Simoni J, Tan LZ. Photo-induced phase-transitions in complex solids. NANOSCALE ADVANCES 2022; 4:4997-5008. [PMID: 36504738 PMCID: PMC9680828 DOI: 10.1039/d2na00481j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Accepted: 11/01/2022] [Indexed: 06/17/2023]
Abstract
Photo-induced phase-transitions (PIPTs) driven by highly cooperative interactions are of fundamental interest as they offer a way to tune and control material properties on ultrafast timescales. Due to strong correlations and interactions, complex quantum materials host several fascinating PIPTs such as light-induced charge density waves and ferroelectricity and have become a desirable setting for studying these PIPTs. A central issue in this field is the proper understanding of the underlying mechanisms driving the PIPTs. As these PIPTs are highly nonlinear processes and often involve multiple time and length scales, different theoretical approaches are often needed to understand the underlying mechanisms. In this review, we present a brief overview of PIPTs realized in complex materials, followed by a discussion of the available theoretical methods with selected examples of recent progress in understanding of the nonequilibrium pathways of PIPTs.
Collapse
Affiliation(s)
| | - Jacopo Simoni
- Molecular Foundry, Lawrence Berkeley National Laboratory USA
| | - Liang Z Tan
- Molecular Foundry, Lawrence Berkeley National Laboratory USA
| |
Collapse
|
105
|
Schnappinger T, Jadoun D, Gudem M, Kowalewski M. Time-resolved X-ray and XUV based spectroscopic methods for nonadiabatic processes in photochemistry. Chem Commun (Camb) 2022; 58:12763-12781. [PMID: 36317595 PMCID: PMC9671098 DOI: 10.1039/d2cc04875b] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Accepted: 10/21/2022] [Indexed: 11/03/2023]
Abstract
The photochemistry of numerous molecular systems is influenced by conical intersections (CIs). These omnipresent nonadiabatic phenomena provide ultra-fast radiationless relaxation channels by creating degeneracies between electronic states and decide over the final photoproducts. In their presence, the Born-Oppenheimer approximation breaks down, and the timescales of the electron and nuclear dynamics become comparable. Due to the ultra-fast dynamics and the complex interplay between nuclear and electronic degrees of freedom, the direct experimental observation of nonadiabatic processes close to CIs remains challenging. In this article, we give a theoretical perspective on novel spectroscopic techniques capable of observing clear signatures of CIs. We discuss methods that are based on ultra-short laser pulses in the extreme ultraviolet and X-ray regime, as their spectral and temporal resolution allow for resolving the ultra-fast dynamics near CIs.
Collapse
Affiliation(s)
- Thomas Schnappinger
- Department of Physics, Stockholm University, Albanova University Centre, SE-106 91 Stockholm, Sweden.
| | - Deependra Jadoun
- Department of Physics, Stockholm University, Albanova University Centre, SE-106 91 Stockholm, Sweden.
| | - Mahesh Gudem
- Department of Physics, Stockholm University, Albanova University Centre, SE-106 91 Stockholm, Sweden.
| | - Markus Kowalewski
- Department of Physics, Stockholm University, Albanova University Centre, SE-106 91 Stockholm, Sweden.
| |
Collapse
|
106
|
Tang D, Jia L, Shen L, Fang WH. Fewest-Switches Surface Hopping with Long Short-Term Memory Networks. J Phys Chem Lett 2022; 13:10377-10387. [PMID: 36317657 DOI: 10.1021/acs.jpclett.2c02299] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The mixed quantum-classical dynamical simulation is essential for studying nonadiabatic phenomena in photophysics and photochemistry. In recent years, many machine learning models have been developed to accelerate the time evolution of the nuclear subsystem. Herein, we implement long short-term memory (LSTM) networks as a propagator to accelerate the time evolution of the electronic subsystem during the fewest-switches surface hopping (FSSH) simulations. A small number of reference trajectories are generated using the original FSSH method, and then the LSTM networks can be built, accompanied by careful examination of typical LSTM-FSSH trajectories that employ the same initial condition and random numbers as the corresponding reference. The constructed network is applied to FSSH to further produce a trajectory ensemble to reveal the mechanism of nonadiabatic processes. Taking Tully's three models as test systems, we qualitatively reproduced the collective results. This work demonstrates that LSTM can be applied to the most popular surface hopping simulations.
Collapse
Affiliation(s)
- Diandong Tang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Luyang Jia
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, China
| |
Collapse
|
107
|
Nelson TR, Fernandez-Alberti S, Tretiak S. Modeling excited-state molecular dynamics beyond the Born-Oppenheimer regime. NATURE COMPUTATIONAL SCIENCE 2022; 2:689-692. [PMID: 38177357 DOI: 10.1038/s43588-022-00357-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2024]
Affiliation(s)
- Tammie R Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA
| | | | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, USA.
- Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, NM, USA.
| |
Collapse
|
108
|
Anderson MC, Schile AJ, Limmer DT. Nonadiabatic transition paths from quantum jump trajectories. J Chem Phys 2022; 157:164105. [DOI: 10.1063/5.0102891] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. We detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.
Collapse
Affiliation(s)
- Michelle C. Anderson
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Addison J. Schile
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - David T. Limmer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Kavli Energy NanoSciences Institute, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
109
|
Green JA, Gómez S, Worth G, Santoro F, Improta R. Solvent Effects on Ultrafast Charge Transfer Population: Insights from the Quantum Dynamics of Guanine-Cytosine in Chloroform. Chemistry 2022; 28:e202201731. [PMID: 35950519 PMCID: PMC9828530 DOI: 10.1002/chem.202201731] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Indexed: 01/12/2023]
Abstract
We study the ultrafast photoactivated dynamics of the hydrogen bonded dimer Guanine-Cytosine in chloroform solution, focusing on the population of the Guanine→Cytosine charge transfer state (GC-CT), an important elementary process for the photophysics and photochemistry of nucleic acids. We integrate a quantum dynamics propagation scheme, based on a linear vibronic model parameterized through time dependent density functional theory calculations, with four different solvation models, either implicit or explicit. On average, after 50 fs, 30∼40 % of the bright excited state population has been transferred to GC-CT. This process is thus fast and effective, especially when transferring from the Guanine bright excited states, in line with the available experimental studies. Independent of the adopted solvation model, the population of GC-CT is however disfavoured in solution with respect to the gas phase. We show that dynamical solvation effects are responsible for this puzzling result and assess the different chemical-physical effects modulating the population of CT states on the ultrafast time-scale. We also propose some simple analyses to predict how solvent can affect the population transfer between bright and CT states, showing that the effect of the solute/solvent electrostatic interactions on the energy of the CT state can provide a rather reliable indication of its possible population.
Collapse
Affiliation(s)
- James A. Green
- Istituto di Biostrutture e Bioimmagini-CNRVia De Amicis 95I-80145Napoli
| | - Sandra Gómez
- Department of ChemistryUniversity College London20 Gordon StreetLondonWC1H 0AJUnited Kingdom,Departamento de Química FísicaUniversity of SalamancaSalamanca37008Spain
| | - Graham Worth
- Department of ChemistryUniversity College London20 Gordon StreetLondonWC1H 0AJUnited Kingdom
| | - Fabrizio Santoro
- Istituto di Chimica die Composti Organometallici (ICCOM-CNR), Area della Ricerca del CNRVia Moruzzi 1I-56124Pisa
| | - Roberto Improta
- Istituto di Biostrutture e Bioimmagini-CNRVia De Amicis 95I-80145Napoli
| |
Collapse
|
110
|
Lin K, Peng J, Xu C, Gu FL, Lan Z. Automatic Evolution of Machine-Learning-Based Quantum Dynamics with Uncertainty Analysis. J Chem Theory Comput 2022; 18:5837-5855. [DOI: 10.1021/acs.jctc.2c00702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Kunni Lin
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Jiawei Peng
- School of Chemistry, South China Normal University, Guangzhou510006, P. R. China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou510006, P. R. China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou510006, P. R. China
| |
Collapse
|
111
|
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
|
112
|
Alfonso Hernandez L, Freixas VM, Rodriguez-Hernandez B, Tretiak S, Fernandez-Alberti S, Oldani N. Exciton-vibrational dynamics induces efficient self-trapping in a substituted nanoring. Phys Chem Chem Phys 2022; 24:24095-24104. [PMID: 36178044 DOI: 10.1039/d2cp03162k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cycloparaphenylenes, being the smallest segments of carbon nanotubes, have emerged as prototypes of the simplest carbon nanohoops. Their unique structure-dynamics-optical properties relationships have motivated a wide variety of synthesis of new related nanohoop species. Studies of how chemical changes, introduced in these new materials, lead to systems with new structural, dynamics and optical properties, expand their functionalities for optoelectronics applications. Herein, we study the effect that conjugation extension of a cycloparaphenylene through the introduction of a satellite tetraphenyl substitution has on its structural and dynamical properties. Our non-adiabatic excited state molecular dynamics simulations suggest that this substitution accelerates the electronic relaxation from the high-energy band to the lowest excited state. This is partially due to efficient conjugation achieved between specific phenyl units as introduced by the tetraphenyl substitution. We observe a particular exciton redistribution during relaxation, in which the tetraphenyl substitution plays a significant role. As a result, an efficient inter-band energy transfer takes place. Besides, the observed phonon-exciton interplay induces a significant exciton self-trapping. Our results encourage and guide the future studies of new phenyl substitutions in carbon nanorings with desired optoelectronic properties.
Collapse
Affiliation(s)
- Laura Alfonso Hernandez
- Departamento de Ciencia Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina.
| | - Victor M Freixas
- Departamento de Ciencia Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina.
| | | | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | | | - Nicolas Oldani
- Departamento de Ciencia Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina.
| |
Collapse
|
113
|
Zhu Y, Peng J, Kang X, Xu C, Lan Z. The principal component analysis of the ring deformation in the nonadiabatic surface hopping dynamics. Phys Chem Chem Phys 2022; 24:24362-24382. [PMID: 36178471 DOI: 10.1039/d2cp03323b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The analysis of the leading active molecular motions in the on-the-fly trajectory surface hopping simulation provides the essential information to understand the geometric evolution in nonadiabatic dynamics. When the ring deformation is involved, the identification of the key active coordinates becomes challenging. A "hierarchical" protocol based on the dimensionality reduction and clustering approaches is proposed for the automatic analysis of the ring deformation in the nonadiabatic molecular dynamics. The representative system keto isocytosine is taken as the prototype to illustrate this protocol. The results indicate that the current hierarchical analysis protocol is a powerful way to clearly clarify both the major and minor active molecular motions of the ring distortion in nonadiabatic dynamics.
Collapse
Affiliation(s)
- Yifei Zhu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China. .,MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.,School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Xu Kang
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China. .,MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China. .,MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China. .,MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
114
|
Joubert-Doriol L. Variational Approach for Linearly Dependent Moving Bases in Quantum Dynamics: Application to Gaussian Functions. J Chem Theory Comput 2022; 18:5799-5809. [PMID: 36166838 DOI: 10.1021/acs.jctc.2c00461] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this paper, we present a variational treatment of the linear dependence for a non-orthogonal time-dependent basis set in solving the Schrödinger equation. The method is based on (i) the definition of a linearly independent working space and (ii) a variational construction of the propagator over finite time steps. The second point allows the method to properly account for changes in the dimensionality of the working space along the time evolution. In particular, the time evolution is represented by a semi-unitary transformation. Tests are carried out on a quartic double-well potential with Gaussian basis functions whose centers evolve according to classical equations of motion. We show that the resulting dynamics converges to the exact one and is unitary by construction.
Collapse
Affiliation(s)
- Loïc Joubert-Doriol
- Université Gustave Eiffel, Université Paris-Est Créteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| |
Collapse
|
115
|
Morzan UN, Díaz Mirón G, Grisanti L, González Lebrero MC, Kaminski Schierle GS, Hassanali A. Non-Aromatic Fluorescence in Biological Matter: The Exception or the Rule? J Phys Chem B 2022; 126:7203-7211. [PMID: 36128666 DOI: 10.1021/acs.jpcb.2c04280] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
While in the vast majority of cases fluorescence in biological matter has been attributed to aromatic or conjugated groups, peptides associated with neurodegenerative diseases, such as Alzheimer's, Parkinson's, or Huntington's, have been recently shown to display an intrinsic visible fluorescence even in the absence of aromatic residues. This has called the attention of researchers from many different fields, trying to understand the origin of this peculiar behavior and, at the same time, motivating the search for novel strategies to control the optical properties of new biophotonic materials. Today, after nearly 15 years of its discovery, there is a growing consensus about the mechanism underlying this phenomenon, namely, that electronic interactions between non-optically active molecules can result in supramolecular assemblies that are fluorescent. Despite this progress, many aspects of this phenomenon remain uncharted territory. In this Perspective, we lay down the state-of-the-art in the field highlighting the open questions from both experimental and theoretical fronts in this fascinating emerging area of non-aromatic fluorescence.
Collapse
Affiliation(s)
- Uriel N Morzan
- International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| | - Gonzalo Díaz Mirón
- Facultad de Ciencias Exactas y Naturales, Departamento de Química Inorgánica, Analítica y Química Física, Universidad de Buenos Aires, C1053 Buenos Aires, Argentina
| | - Luca Grisanti
- Division of Theoretical Physics, Ruđer Bos̆cković Institute, Bijenic̆ka cesta 54, 10000 Zagreb, Croatia
| | - Mariano C González Lebrero
- Facultad de Ciencias Exactas y Naturales, Departamento de Química Inorgánica, Analítica y Química Física, Universidad de Buenos Aires, C1053 Buenos Aires, Argentina
| | | | - Ali Hassanali
- International Centre for Theoretical Physics, Strada Costiera 11, 34151 Trieste, Italy
| |
Collapse
|
116
|
Freixas VM, Tretiak S, Fernandez-Alberti S. Infinitene: Computational Insights from Nonadiabatic Excited State Dynamics. J Phys Chem Lett 2022; 13:8495-8501. [PMID: 36066077 DOI: 10.1021/acs.jpclett.2c02296] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Progress in organic synthesis opens exploration of a rich diversity of molecules with interesting new structural topologies. This is the case of a recently synthesized helically twisted figure-eight molecule coined infinitene. The molecule belongs to a numerous family of looped polyarenes, where the degree of π-conjugation is controlled by high strain energies and steric hindrances. A particular balance of these ingredients leads to unusual optoelectronic properties potentially suitable for a range of applications in nanoelectronics and photonics. Due to its recent discovery, the photophysical properties of infinitene remain unexplored. In this Letter, atomistic nonadiabatic excited state molecular dynamics modeling unveils unique features of intramolecular electronic and vibrational energy relaxation and redistribution that take place after molecular photoexcitation. Our results detail relationships between optical and electronic properties providing useful knowledge for future molecular designs related to infinitene.
Collapse
Affiliation(s)
- Victor Manuel Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | | |
Collapse
|
117
|
Bondanza M, Demoulin B, Lipparini F, Barbatti M, Mennucci B. Trajectory Surface Hopping for a Polarizable Embedding QM/MM Formulation. J Phys Chem A 2022; 126:6780-6789. [PMID: 36107729 PMCID: PMC9527758 DOI: 10.1021/acs.jpca.2c04756] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We present the implementation of trajectory surface-hopping
nonadiabatic
dynamics for a polarizable embedding QM/MM formulation. Time-dependent
density functional theory was used at the quantum mechanical level
of theory, whereas the molecular mechanics description involved the
polarizable AMOEBA force field. This implementation has been obtained
by integrating the surface-hopping program Newton-X NS with an interface
between the Gaussian 16 and the Tinker suites of codes to calculate
QM/AMOEBA energies and forces. The implementation has been tested
on a photoinduced electron-driven proton-transfer reaction involving
pyrimidine and a hydrogen-bonded water surrounded by a small cluster
of water molecules and within a large water droplet.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | | | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, 13385 Marseille, France
- Institut Universitaire de France, 75231 Paris, France
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
118
|
Significant impact of deprotonated status on the photoisomerization dynamics of bacteriophytochrome chromophore. CHINESE CHEM LETT 2022. [DOI: 10.1016/j.cclet.2022.107850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
119
|
Bruce‐Chwatt T, Naidoo KJ. Molecular mechanisms from reaction coordinate graph enabled multidimensional free energies illustrated on water dimer hydrogen bonding. J Comput Chem 2022; 43:1802-1813. [PMID: 36054751 PMCID: PMC9543413 DOI: 10.1002/jcc.26982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 06/10/2022] [Accepted: 07/22/2022] [Indexed: 11/11/2022]
Abstract
Computing the free energies of molecular mechanisms in multidimensional space relies on combinations of geometrically complex reaction coordinates. We show how a graph theory implementation reduces complexity, and illustrate this on the arrangements of hydrogen bonding of a water dimer. The reaction coordinates and forces are computed using graphs that define the dependencies on the atoms in the Free Energy from Adaptive Reaction Coordinate Forces (FEARCF) library. The library can be interfaced with classical molecular dynamics as well as quantum molecular dynamics packages. Multidimensional interdependent reaction coordinates are constructed to produce complex free energy hypersurfaces. The reaction coordinates are graphed from atomic and molecular components to define points, distances, vectors, angles, planes and combinations thereof. The resultant free energy surfaces that are a function of the distance, angles, planes, and so on, can represent molecular mechanisms in reduced dimensions from the component atomic Cartesian coordinate degrees of freedom. The FEARCF library can be interfaced with any molecular package. Here, we demonstrate the link to NWChem to compute a hyperdimensional DFT (aug‐cc‐pVDZ basis set and X3LYP exchange correlation functionals) free energy space of a water dimer. Analysis of the water dimer free energy hypervolume reveals that while the chain and cyclic hydrogen bonding configurations are located in stable minimum energy wells, the bifurcated hydrogen bond configuration is a gateway to instability and dimer dissociation.
Collapse
Affiliation(s)
- Tomás Bruce‐Chwatt
- Scientific Computing Research Unit, Department of Chemistry University of Cape Town Cape Town South Africa
| | - Kevin J. Naidoo
- Scientific Computing Research Unit, Department of Chemistry University of Cape Town Cape Town South Africa
| |
Collapse
|
120
|
Stetzler J, Rassolov VA. Comparison of Born–Oppenheimer approximation and electron-nuclear correlation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2106321] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Affiliation(s)
- Julian Stetzler
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, SC, USA
| | - Vitaly A. Rassolov
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, SC, USA
| |
Collapse
|
121
|
Marsili E, Prlj A, Curchod BFE. A Theoretical Perspective on the Actinic Photochemistry of 2-Hydroperoxypropanal. J Phys Chem A 2022; 126:5420-5433. [PMID: 35900368 PMCID: PMC9393889 DOI: 10.1021/acs.jpca.2c03783] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The photochemical reactions triggered by the sunlight
absorption
of transient volatile organic compounds in the troposphere are notoriously
difficult to characterize experimentally due to the unstable and short-lived
nature of these organic molecules. Some members of this family of
compounds are likely to exhibit a rich photochemistry given the diversity
of functional groups they can bear. Even more interesting is the photochemical
fate of volatile organic compounds bearing more than one functional
group that can absorb light—this is the case, for example,
of α-hydroperoxycarbonyls, which are formed during the oxidation
of isoprene. Experimental observables characterizing the photochemistry
of these molecules like photoabsorption cross-sections or photolysis
quantum yields are currently missing, and we propose here to leverage
a recently developed computational protocol to predict in silico the
photochemical fate of 2-hydroperoxypropanal (2-HPP) in the actinic
region. We combine different levels of electronic structure methods—SCS-ADC(2)
and XMS-CASPT2—with the nuclear ensemble approach and trajectory
surface hopping to understand the mechanistic details of the possible
nonradiative processes of 2-HPP. In particular, we predict the photoabsorption
cross-section and the wavelength-dependent quantum yields for the
observed photolytic pathways and combine them to determine in silico
photolysis rate constants. The limitations of our protocol and possible
future improvements are discussed.
Collapse
Affiliation(s)
- Emanuele Marsili
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Antonio Prlj
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| |
Collapse
|
122
|
Avagliano D, Bonfanti M, Nenov A, Garavelli M. Automatized protocol and interface to simulate QM/MM time-resolved transient absorption at TD-DFT level with COBRAMM. J Comput Chem 2022; 43:1641-1655. [PMID: 35815854 PMCID: PMC9544370 DOI: 10.1002/jcc.26966] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2022] [Revised: 06/07/2022] [Accepted: 06/10/2022] [Indexed: 12/24/2022]
Abstract
We present a series of new implementations that we recently introduced in COBRAMM, the open-source academic software developed in our group. The goal of these implementations is to offer an automatized workflow and interface to simulate time-resolved transient absorption (TA) spectra of medium-to-big chromophore embedded in a complex environment. Therefore, the excited states absorption and the stimulated emission are simulated along nonadiabatic dynamics performed with trajectory surface hopping. The possibility of treating systems from medium to big size is given by the use of time-dependent density functional theory (TD-DFT) and the presence of the environment is taken into account employing a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme. The full implementation includes a series of auxiliary scripts to properly setup the QM/MM system, the calculation of the wavefunction overlap along the dynamics for the propagation, the evaluation of the transition dipole moment at linear response TD-DFT level, and scripts to setup, run and analyze the TA from an ensemble of trajectories. Altogether, we believe that our implementation will open the door to the easily simulate the time-resolved TA of systems so far computationally inaccessible.
Collapse
Affiliation(s)
- Davide Avagliano
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| | - Matteo Bonfanti
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy.,Fondazione Human Technopole - Viale Rita Levi-Montalcini, 1 - Area MIND - Cargo 6 - 20157, Milano, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| |
Collapse
|
123
|
Jones CM, List NH, Martínez TJ. Steric and Electronic Origins of Fluorescence in GFP and GFP-like Proteins. J Am Chem Soc 2022; 144:12732-12746. [PMID: 35786916 DOI: 10.1021/jacs.2c02946] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Fluorescent proteins have become routine tools for biological imaging. However, their nanosecond lifetimes on the excited state present computational hurdles to a full understanding of these photoactive proteins. In this work, we simulate approximately 0.5 nanoseconds of ab initio molecular dynamics to elucidate steric and electronic features responsible for fluorescent protein behavior. Using green fluorescent protein (GFP) and Dronpa2─widely used fluorescent proteins with contrasting functionality─as case studies, we leverage previous findings in the gas phase and solution to explore the deactivation mechanisms available to these proteins. Starting with ground-state analyses, we identify steric (the distribution of empty pockets near the chromophore) and electronic (electric fields exerted on chromophore moieties) factors that offer potential avenues for rational design. Picosecond timescale simulations on the excited state reveal that the chromophore can access twisted structures in Dronpa2, while the chromophore is largely confined to planarity in GFP. We couple ab initio multiple spawning (AIMS) and enhanced sampling simulations to discover and characterize conical intersection seams that facilitate internal conversion, which is a rare event in both systems. Our AIMS simulations correctly capture the relative fluorescence profiles of GFP and Dronpa2 within the first few picoseconds, and we attribute the diminished fluorescence intensity of Dronpa2, relative to GFP, to flexible chromophore intermediates on the excited state. Furthermore, we predict that twisted chromophore intermediates produce red-shifted intensities in the Dronpa2 fluorescence spectrum. If confirmed experimentally, this spectroscopic signature would provide valuable insights when screening and developing novel fluorescent proteins.
Collapse
Affiliation(s)
- Chey M Jones
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
124
|
Koessler ER, Mandal A, Huo P. Incorporating Lindblad Decay Dynamics into Mixed Quantum-Classical Simulations. J Chem Phys 2022; 157:064101. [DOI: 10.1063/5.0099922] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We derive the $\mathcal{L}$-MFE method to incorporate Lindblad jump operator dynamics into the mean-field Ehrenfest (MFE) approach. We map the density matrix evolution of Lindblad dynamics onto pure state coefficients using trajectory averages. We use simple assumptions to construct the $\mathcal{L}$-MFE method that satisfies this exact mapping. This establishes a method that uses independent trajectories which exactly reproduces Lindblad decay dynamics using a wavefunction description, with deterministic changes of the magnitudes of the quantum expansion coefficients, while only adding on a stochastic phase. We further demonstrate that when including nuclei in the Ehrenfest dynamics, the $\mathcal{L}$-MFE method gives semi-quantitatively accurate results, with the accuracy limited by the accuracy of the approximations present in the semiclassical MFE approach. This work provides a general framework to incorporate Lindblad dynamics into semiclassical or mixed quantum-classical simulations.
Collapse
Affiliation(s)
- Eric R Koessler
- Chemistry, University of Rochester, United States of America
| | | | - Pengfei Huo
- Department of Chemsitry, University of Rochester Department of Chemistry, United States of America
| |
Collapse
|
125
|
Schnack-Petersen AK, Koch H, Coriani S, Kjønstad EF. Efficient implementation of molecular CCSD gradients with Cholesky-decomposed electron repulsion integrals. J Chem Phys 2022; 156:244111. [DOI: 10.1063/5.0087261] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
We present an efficient implementation of ground and excited state coupled cluster singles and doubles (CCSD) gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition and density fitting are both inner projection methods, and, thus, similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3 O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. Furthermore, the computational time per optimization cycle is significantly lower compared to other implementations based on an inner projection method. We showcase the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.
Collapse
Affiliation(s)
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Eirik F. Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| |
Collapse
|
126
|
Liang R, Bakhtiiari A. Multiscale simulation unravels the light-regulated reversible inhibition of dihydrofolate reductase by phototrexate. J Chem Phys 2022; 156:245102. [DOI: 10.1063/5.0096349] [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
Molecular photoswitches are widely used in photopharmacology, where the biomolecular functions are photo-controlled reversibly with high spatiotemporal precision. Despite the success of this field, it remains elusive how the protein environment modulates the photochemical properties of photoswitches. Understanding this fundamental question is critical for designing more effective light-regulated drugs with mitigated side effects. In our recent work, we employed first-principles non-adiabatic dynamics simulations to probe the effects of protein on the trans to cis photoisomerization of phototrexate (PTX), a photochromic analog of the anticancer therapeutic methotrexate that inhibits the target enzyme dihydrofolate reductase (DHFR). Building upon this study, in this work, we employ multiscale simulations to unravel the full photocycle underlying the light-regulated reversible inhibition of DHFR by PTX, which remains elusive until now. First-principles non-adiabatic dynamics simulations reveal that the cis to trans photoisomerization quantum yield is hindered in the protein due to backward isomerization on the ground-state following non-adiabatic transition, which arises from the favorable binding of the cis isomer with the protein. However, free energy simulations indicate that cis to trans photoisomerization significantly decreases the binding affinity of the PTX. Thus, the cis to trans photoisomerization most likely precedes the ligand unbinding from the protein. We propose the most probable photocycle of the PTX-DHFR system. Our comprehensive simulations highlight the trade-offs among the binding affinity, photoisomerization quantum yield, and the thermal stability of the ligand's different isomeric forms. As such, our work reveals new design principles of light-regulated drugs in photopharmacology.
Collapse
Affiliation(s)
- Ruibin Liang
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Amirhossein Bakhtiiari
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| |
Collapse
|
127
|
Shakiba M, Stippell E, Li W, Akimov AV. Nonadiabatic Molecular Dynamics with Extended Density Functional Tight-Binding: Application to Nanocrystals and Periodic Solids. J Chem Theory Comput 2022; 18:5157-5180. [PMID: 35758936 DOI: 10.1021/acs.jctc.2c00297] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report a new methodology for nonadiabatic molecular dynamics calculations within the extended tight-binding (xTB) framework. We demonstrate the applicability of the developed approach to finite and periodic systems with thousands of atoms by modeling "hot" electron relaxation dynamics in silicon nanocrystals and electron-hole recombination in both a graphitic carbon nitride monolayer and a titanium-based metal-organic framework (MOF). This work reports the nonadiabatic dynamic simulations in the largest Si nanocrystals studied so far by the xTB framework, with diameters up to 3.5 nm. For silicon nanocrystals, we find a non-monotonic dependence of "hot" electron relaxation rates on the nanocrystal size, in agreement with available experimental reports. We rationalize this relationship by a combination of decreasing nonadiabatic couplings related to system size and the increase of available coherent transfer pathways in systems with higher densities of states. We emphasize the importance of proper treatment of coherences for obtaining such non-monotonic dependences. We characterize the electron-hole recombination dynamics in the graphitic carbon nitride monolayer and the Ti-containing MOF. We demonstrate the importance of spin-adaptation and proper sampling of surface hopping trajectories in modeling such processes. We also assess several trajectory surface hopping schemes and highlight their distinct qualitative behavior in modeling the excited-state dynamics in superexchange-like models depending on how they handle coherences between nearly parallel states.
Collapse
Affiliation(s)
- Mohammad Shakiba
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Elizabeth Stippell
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, China
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
128
|
Smith JET, Lee J, Sharma S. Near-Exact Nuclear Gradients of Complete Active Space Self-Consistent Field Wave Functions. J Chem Phys 2022; 157:094104. [DOI: 10.1063/5.0085515] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this paper, we study the nuclear gradients of heat bath configuration interaction self-consistent field (HCISCF) wave functions and use them to optimize molecular geometries for various molecules.We show that the HCISCF nuclear gradients are fairly insensitive to the size of the "selected" variational space, which allows us to reduce the computational cost without introducing significant error.The ability of HCISCF to treat larger active spaces combined with the flexibility for users to control the computational cost makes the method very attractive for studying strongly correlated systems which require a larger active space than possible with complete active space self-consistent field (CASSCF).Finally, we study the realistic catalyst, Fe(PDI), and highlight some of the challenges this system poses for density functional theory (DFT).We demonstrate how HCISCF can clarify the energetic stability of geometries obtained from DFT when the results are strongly dependent on the functional.We also use the HCISCF gradients to optimize geometries for this species and study the adiabatic singlet-triplet gap. During geometry optimization, we find that multiple near-degenerate local minima exist on the triplet potential energy surface.
Collapse
Affiliation(s)
- James E. T. Smith
- Center for Computational Quantum Phyics, Flatiron Institute, United States of America
| | - Joonho Lee
- Chemistry, Columbia University Department of Chemistry, United States of America
| | - Sandeep Sharma
- University of Colorado at Boulder, United States of America
| |
Collapse
|
129
|
Wang Y, Guo H, Yarkony DR. Internal conversion and intersystem crossing dynamics based on coupled potential energy surfaces with full geometry-dependent spin-orbit and derivative couplings. Nonadiabatic photodissociation dynamics of NH 3(A) leading to the NH(X 3Σ -, a 1Δ) + H 2 channel. Phys Chem Chem Phys 2022; 24:15060-15067. [PMID: 35696936 DOI: 10.1039/d2cp01271e] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We simulate the photodissociation of NH3 originating from its first excited singlet state S1 into the NH2 + H (radical) and NH + H2 (molecular) channels. The states considered are the ground singlet state S0, the first excited singlet state S1 and the lowest-lying triplet state T1, which permit for the first time a uniform treatment of the internal conversion and intersystem crossing. The simulations are based on a diabatic potential energy matrix (DPEM) of S0, S1 coupled by a conical intersection seam, as well as a potential energy surface (PES) for T1 coupled by spin-orbit coupling (SOC) to the two singlet states. The DPEM and PES are fitted to ab initio electronic structure data (ESD) including energies, energy gradients, and derivative couplings. The DPEM also defines an adiabatic to diabatic state (AtD) transformation, which is used to transform the singular adiabatic SOC into a smooth function of the nuclear coordinates in the diabatic representation, allowing the diabatic SOC to be fit to an analytical functional form. ESD and SOC data obtained from these surfaces can serve as input for either quantum or semi-classical characterization of the nonadiabatic dynamics. Using the SHARC suite of programs, nonadiabatic simulations based on over 40 000 semi-classical trajectories assess the convergence of our results. The production of NH + H2 is not direct, but is only achieved through a quasi-statistical dissociation mechanism after internal conversion to the ground electronic state. This leads to a much lower yield comparing with the main NH2 + H channel. The NH(X3Σ_) radical produced through the intersystem crossing from S0 to T1 is rare (∼0.2%) compared to NH(a1Δ) due to the process being spin forbidden.
Collapse
Affiliation(s)
- Yuchen Wang
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA.
| |
Collapse
|
130
|
Ten Brink M, Gräber S, Hopjan M, Jansen D, Stolpp J, Heidrich-Meisner F, Blöchl PE. Real-time non-adiabatic dynamics in the one-dimensional Holstein model: Trajectory-based vs exact methods. J Chem Phys 2022; 156:234109. [PMID: 35732530 DOI: 10.1063/5.0092063] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We benchmark a set of quantum-chemistry methods, including multitrajectory Ehrenfest, fewest-switches surface-hopping, and multiconfigurational-Ehrenfest dynamics, against exact quantum-many-body techniques by studying real-time dynamics in the Holstein model. This is a paradigmatic model in condensed matter theory incorporating a local coupling of electrons to Einstein phonons. For the two-site and three-site Holstein model, we discuss the exact and quantum-chemistry methods in terms of the Born-Huang formalism, covering different initial states, which either start on a single Born-Oppenheimer surface, or with the electron localized to a single site. For extended systems with up to 51 sites, we address both the physics of single Holstein polarons and the dynamics of charge-density waves at finite electron densities. For these extended systems, we compare the quantum-chemistry methods to exact dynamics obtained from time-dependent density matrix renormalization group calculations with local basis optimization (DMRG-LBO). We observe that the multitrajectory Ehrenfest method, in general, only captures the ultrashort time dynamics accurately. In contrast, the surface-hopping method with suitable corrections provides a much better description of the long-time behavior but struggles with the short-time description of coherences between different Born-Oppenheimer states. We show that the multiconfigurational Ehrenfest method yields a significant improvement over the multitrajectory Ehrenfest method and can be converged to the exact results in small systems with moderate computational efforts. We further observe that for extended systems, this convergence is slower with respect to the number of configurations. Our benchmark study demonstrates that DMRG-LBO is a useful tool for assessing the quality of the quantum-chemistry methods.
Collapse
Affiliation(s)
- M Ten Brink
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - S Gräber
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - M Hopjan
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - D Jansen
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - J Stolpp
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - F Heidrich-Meisner
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - P E Blöchl
- Institut für Theoretische Physik, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
131
|
Asaad M, Joubert-Doriol L, Izmaylov AF. Controlling energy conservation in quantum dynamics with independently moving basis functions: Application to multi-configuration Ehrenfest. J Chem Phys 2022; 156:204121. [PMID: 35649883 DOI: 10.1063/5.0087797] [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
Application of the time-dependent variational principle to a linear combination of frozen-width Gaussians describing the nuclear wavefunction provides a formalism where the total energy is conserved. The computational downside of this formalism is that trajectories of individual Gaussians are solutions of a coupled system of differential equations, limiting implementation to serial propagation algorithms. To allow for parallelization and acceleration of the computation, independent trajectories based on simplified equations of motion were suggested. Unfortunately, within practical realizations involving finite Gaussian bases, this simplification leads to breaking the energy conservation. We offer a solution for this problem by using Lagrange multipliers to ensure the energy and norm conservation regardless of basis function trajectories or basis completeness. We illustrate our approach within the multi-configurational Ehrenfest method considering a linear vibronic coupling model.
Collapse
Affiliation(s)
- Mina Asaad
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| | - Loïc Joubert-Doriol
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Artur F Izmaylov
- Department of Physical and Environmental Sciences, University of Toronto Scarborough, Toronto, Ontario M1C 1A4, Canada
| |
Collapse
|
132
|
Erickson M, Han Y, Rasulev B, Kilin D. Molecular Dynamics Study of the Photodegradation of Polymeric Chains. J Phys Chem Lett 2022; 13:4374-4380. [PMID: 35544382 DOI: 10.1021/acs.jpclett.2c00802] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The development of reusable polymeric materials inspires an attempt to combine renewable biomass with upcycling to form a biorenewable closed system. It has been reported that 2,5-furandicarboxylic acid (FDCA) can be recovered for recycling when incorporated as monomers into photodegradable polymeric systems. Here, we conduct density functional theory (DFT) studies with periodic boundary conditions on microscopic structures involved in the photodegradation of polymeric chains incorporating FDCA and 2-nitro-1,3-benzenedimethanol. The photodegradation process of polymeric chains is studied using time-dependent excited-state molecular dynamics (TDESMD) in vacuum and aqueous environments. Changes in the photophysical properties for reaction intermediates are characterized by ground-state observables. The distribution of reaction intermediates and products is obtained from TDESMD trajectories using cheminformatics techniques. Results show that a higher degree of polymeric chain degradation is achieved in the vacuum environment. Additionally, one finds that the FDCA molecule is recoverable in the aqueous environment, in qualitative agreement with experimental findings.
Collapse
Affiliation(s)
- Meade Erickson
- Department of Coatings and Polymeric Materials, North Dakota State University, Fargo, North Dakota 58108, United States
| | - Yulun Han
- Department of Chemistry and Biochemistry, North Dakota State University, Fargo, North Dakota 58108, United States
| | - Bakhtiyor Rasulev
- Department of Coatings and Polymeric Materials, North Dakota State University, Fargo, North Dakota 58108, United States
| | - Dmitri Kilin
- Department of Chemistry and Biochemistry, North Dakota State University, Fargo, North Dakota 58108, United States
| |
Collapse
|
133
|
Simulation of Diffusion-Controlled Growth of Interdependent Nuclei under Potentiostatic Conditions. MATERIALS 2022; 15:ma15103603. [PMID: 35629634 PMCID: PMC9147730 DOI: 10.3390/ma15103603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/23/2022] [Revised: 05/11/2022] [Accepted: 05/16/2022] [Indexed: 12/10/2022]
Abstract
The problem of diffusion-controlled growth following an instantaneous nucleation event was studied within the framework of a new numerical model, considering the spatial distribution of hemispherical nuclei on the electrode surface and the mutual influence of growing nuclei via the collision of 3D diffusion fields. The simulation of the diffusion-controlled growth of hexagonal and random ensembles was performed at the overpotential-dependent number density of nuclei. The diffusion flow to each nucleus within a random ensemble was simulated by the finite difference method using the derived analytical expressions for the surface areas and the volumes formed at the intersection of 3D diffusion fields with the side faces of a virtual right prism with a Voronoi polygon base. The implementation of this approach provides an accurate calculation of concentration profiles, time dependences of the size of nuclei, and current transients. The results, including total current density transients, growth exponents, and nucleus size distribution, were compared with models developed within the concept of planar diffusion zones, the mean-field approximation and the Brownian dynamics simulation method, as well as with experimental data from the literature. The prospects of the model for studying the initial stages of electrocrystallization were discussed.
Collapse
|
134
|
Yao N, Chen X, Fu ZH, Zhang Q. Applying Classical, Ab Initio, and Machine-Learning Molecular Dynamics Simulations to the Liquid Electrolyte for Rechargeable Batteries. Chem Rev 2022; 122:10970-11021. [PMID: 35576674 DOI: 10.1021/acs.chemrev.1c00904] [Citation(s) in RCA: 65] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Rechargeable batteries have become indispensable implements in our daily life and are considered a promising technology to construct sustainable energy systems in the future. The liquid electrolyte is one of the most important parts of a battery and is extremely critical in stabilizing the electrode-electrolyte interfaces and constructing safe and long-life-span batteries. Tremendous efforts have been devoted to developing new electrolyte solvents, salts, additives, and recipes, where molecular dynamics (MD) simulations play an increasingly important role in exploring electrolyte structures, physicochemical properties such as ionic conductivity, and interfacial reaction mechanisms. This review affords an overview of applying MD simulations in the study of liquid electrolytes for rechargeable batteries. First, the fundamentals and recent theoretical progress in three-class MD simulations are summarized, including classical, ab initio, and machine-learning MD simulations (section 2). Next, the application of MD simulations to the exploration of liquid electrolytes, including probing bulk and interfacial structures (section 3), deriving macroscopic properties such as ionic conductivity and dielectric constant of electrolytes (section 4), and revealing the electrode-electrolyte interfacial reaction mechanisms (section 5), are sequentially presented. Finally, a general conclusion and an insightful perspective on current challenges and future directions in applying MD simulations to liquid electrolytes are provided. Machine-learning technologies are highlighted to figure out these challenging issues facing MD simulations and electrolyte research and promote the rational design of advanced electrolytes for next-generation rechargeable batteries.
Collapse
Affiliation(s)
- Nan Yao
- Beijing Key Laboratory of Green Chemical Reaction Engineering and Technology, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Xiang Chen
- Beijing Key Laboratory of Green Chemical Reaction Engineering and Technology, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Zhong-Heng Fu
- Beijing Key Laboratory of Green Chemical Reaction Engineering and Technology, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Qiang Zhang
- Beijing Key Laboratory of Green Chemical Reaction Engineering and Technology, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| |
Collapse
|
135
|
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
|
136
|
Coonjobeeharry J, Spinlove KE, Sanz Sanz C, Sapunar M, Došlić N, Worth GA. Mixed-quantum-classical or fully-quantized dynamics? A unified code to compare methods. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200386. [PMID: 35341308 DOI: 10.1098/rsta.2020.0386] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2021] [Accepted: 09/06/2021] [Indexed: 06/14/2023]
Abstract
Three methods for non-adiabatic dynamics are compared to highlight their capabilities. Multi-configurational time-dependent Hartree is a full grid-based solution to the time-dependent Schrödinger equation, variational multi-configurational Gaussian (vMCG) uses a less flexible but unrestricted Gaussian wavepacket basis, and trajectory surface hopping (TSH) replaces the nuclear wavepacket with a swarm of classical trajectories. Calculations with all methods using a model Hamiltonian were performed. The vMCG and TSH were also then run in a direct dynamics mode, with the potential energy surfaces calculated on-the-fly using quantum chemistry calculations. All dynamics calculations used the Quantics package, with the TSH calculations using a new interface to a surface hopping code. A novel approach to calculate adiabatic populations from grid-based quantum dynamics using a time-dependent discrete variable representation is presented, allowing a proper comparison of methods. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
- J Coonjobeeharry
- Department of Chemistry, University College London, 20, Gordon St., London WC1H 0AJ, UK
| | - K E Spinlove
- Faculty of Science and Engineering, University of Groningen, Nijenborgh 4, Groningen 9747AG, Netherlands
| | - C Sanz Sanz
- Department of Applied Physical Chemistry, Faculty of Science, Autonoma University Madrid, Campus Cantoblanco, Madrid 28049, Spain
| | - M Sapunar
- Department of Physical Chemistry, Rudjer Bošković Institute, Bijenička 54, Zagreb 1000, Croatia
| | - N Došlić
- Department of Physical Chemistry, Rudjer Bošković Institute, Bijenička 54, Zagreb 1000, Croatia
| | - G A Worth
- Department of Chemistry, University College London, 20, Gordon St., London WC1H 0AJ, UK
| |
Collapse
|
137
|
Ashfold MNR, Kim SK. Non-Born-Oppenheimer effects in molecular photochemistry: an experimental perspective. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200376. [PMID: 35341307 DOI: 10.1098/rsta.2020.0376] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Accepted: 06/28/2021] [Indexed: 06/14/2023]
Abstract
Non-adiabatic couplings between Born-Oppenheimer (BO)-derived potential energy surfaces are now recognized as pivotal in describing the non-radiative decay of electronically excited molecules following photon absorption. This opinion piece illustrates how non-BO effects provide photostability to many biomolecules when exposed to ultraviolet radiation, yet in many other cases are key to facilitating 'reactive' outcomes like isomerization and bond fission. The examples are presented in order of decreasing molecular complexity, spanning studies of organic sunscreen molecules in solution, through two families of heteroatom containing aromatic molecules and culminating with studies of isolated gas phase H2O molecules that afford some of the most detailed insights yet available into the cascade of non-adiabatic couplings that enable the evolution from photoexcited molecule to eventual products. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
| | - Sang Kyu Kim
- Department of Chemistry, KAIST, Daejeon 34141, Republic of Korea
| |
Collapse
|
138
|
Marsili E, Agostini F, Nauts A, Lauvergnat D. Quantum dynamics with curvilinear coordinates: models and kinetic energy operator. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200388. [PMID: 35341305 DOI: 10.1098/rsta.2020.0388] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/29/2021] [Accepted: 11/01/2021] [Indexed: 06/14/2023]
Abstract
In order to simplify the numerical solution of the time-dependent or time-independent Schrödinger equations associated with atomic and molecular motions, the use of well-adapted coordinates is essential. Usually, this set of curvilinear coordinates leads to a Hamiltonian operator that is as separable as possible. Although their corresponding kinetic energy operator (KEO) expressions can be derived analytically for small systems or special kinds of coordinates, a numerical and exact approach allows one to compute them in terms of sophisticated curvilinear coordinates. Furthermore, the numerical approach enables one to easily define reduced-dimensionality or constrained models. We present here a recent implementation of this numerical approach that allows nested coordinate transformations, therefore leading to great flexibility in the definition of the curvilinear coordinates. Furthermore, this implementation has no limitations in terms of numbers of atoms or coordinate transformations. The quantum dynamics of the cis-trans photoisomerization of part of the retinal chromophore illustrates the construction of the coordinates and KEO part of a three-dimensional model. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
- Emanuele Marsili
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
- Department of Chemistry, Durham University, Durham DH1 3LE, UK
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| | - André Nauts
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
- Institute of Condensed Matter and Nanosciences (NAPS), Université Catholique de Louvain, BE-1348 Louvain-la-Neuve, Belgium
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| |
Collapse
|
139
|
Mukherjee S, Pinheiro M, Demoulin B, Barbatti M. Simulations of molecular photodynamics in long timescales. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2022; 380:20200382. [PMID: 35341303 PMCID: PMC8958277 DOI: 10.1098/rsta.2020.0382] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Nonadiabatic dynamics simulations in the long timescale (much longer than 10 ps) are the next challenge in computational photochemistry. This paper delimits the scope of what we expect from methods to run such simulations: they should work in full nuclear dimensionality, be general enough to tackle any type of molecule and not require unrealistic computational resources. We examine the main methodological challenges we should venture to advance the field, including the computational costs of the electronic structure calculations, stability of the integration methods, accuracy of the nonadiabatic dynamics algorithms and software optimization. Based on simulations designed to shed light on each of these issues, we show how machine learning may be a crucial element for long time-scale dynamics, either as a surrogate for electronic structure calculations or aiding the parameterization of model Hamiltonians. We show that conventional methods for integrating classical equations should be adequate to extended simulations up to 1 ns and that surface hopping agrees semiquantitatively with wave packet propagation in the weak-coupling regime. We also describe our optimization of the Newton-X program to reduce computational overheads in data processing and storage. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.
Collapse
Affiliation(s)
| | - Max Pinheiro
- Aix Marseille University, CNRS, ICR, Marseille, France
| | | | - Mario Barbatti
- Aix Marseille University, CNRS, ICR, Marseille, France
- Institut Universitaire de France, 75231 Paris, France
| |
Collapse
|
140
|
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
|
141
|
Liu W, Wang Z, Chen Z, Luo J, Li S, Wang L. Algorithm advances and applications of time‐dependent first‐principles simulations for ultrafast dynamics. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1577] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Wen‐Hao Liu
- State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors Chinese Academy of Sciences Beijing China
- Center of Materials Science and Optoelectronics Engineering University of Chinese Academy of Sciences Beijing China
| | - Zhi Wang
- State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors Chinese Academy of Sciences Beijing China
| | - Zhang‐Hui Chen
- Materials Science Division Lawrence Berkeley National Laboratory Berkeley California USA
| | - Jun‐Wei Luo
- State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors Chinese Academy of Sciences Beijing China
- Center of Materials Science and Optoelectronics Engineering University of Chinese Academy of Sciences Beijing China
- Beijing Academy of Quantum Information Sciences Beijing China
| | - Shu‐Shen Li
- State Key Laboratory of Superlattices and Microstructures, Institute of Semiconductors Chinese Academy of Sciences Beijing China
- Center of Materials Science and Optoelectronics Engineering University of Chinese Academy of Sciences Beijing China
- Beijing Academy of Quantum Information Sciences Beijing China
| | - Lin‐Wang Wang
- Materials Science Division Lawrence Berkeley National Laboratory Berkeley California USA
| |
Collapse
|
142
|
Yong H, Keefer D, Mukamel S. Imaging Purely Nuclear Quantum Dynamics in Molecules by Combined X-ray and Electron Diffraction. J Am Chem Soc 2022; 144:7796-7804. [DOI: 10.1021/jacs.2c01311] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Haiwang Yong
- Department of Chemistry, University of California, Irvine, Irvine, California 92697, United States
- Department of Physics and Astronomy, University of California, Irvine, Irvine, California 92697, United States
| | - Daniel Keefer
- Department of Chemistry, University of California, Irvine, Irvine, California 92697, United States
- Department of Physics and Astronomy, University of California, Irvine, Irvine, California 92697, United States
| | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine, Irvine, California 92697, United States
- Department of Physics and Astronomy, University of California, Irvine, Irvine, California 92697, United States
| |
Collapse
|
143
|
Talotta F, Lauvergnat D, Agostini F. Describing the photo-isomerization of a retinal chromophore model with coupled and quantum trajectories. J Chem Phys 2022; 156:184104. [DOI: 10.1063/5.0089415] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The exact factorization of the electron-nuclear wavefunction is applied to the study of the photo- isomerization of a retinal chromophore model. We describe such an ultrafast nonadiabatic process by analyzing the time-dependent potentials of the theory and by mimicking nuclear dynamics with quantum and coupled trajectories. The time-dependent vector and scalar potentials are the signature of the exact factorization, as they guide nuclear dynamics by encoding the complete electronic dynamics and including excited-state effects. Analysis of the potentials is, thus, essential - when possible - to predict the time-dependent behavior of the system of interest. In this work, we employ the exact time-dependent potentials, available for the numerically-exactly solvable model used here, to propagate quantum nuclear trajectories representing the isomerization reaction of the retinal chromophore. The quantum trajectories are the best possible trajectory-based description of the reaction when using the exact-factorization formalism, and thus allow us to assess the performance of the coupled-trajectory, fully approximate, schemes derived from the exact-factorization equations.
Collapse
Affiliation(s)
| | - David Lauvergnat
- Institut de Chimie Physique, UMR 8000, CNRS Délégation Ile-de-France Sud, France
| | | |
Collapse
|
144
|
Vandaele E, Mališ M, Luber S. The ΔSCF method for non-adiabatic dynamics of systems in the liquid phase. J Chem Phys 2022; 156:130901. [PMID: 35395890 DOI: 10.1063/5.0083340] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Computational studies of ultrafast photoinduced processes give valuable insights into the photochemical mechanisms of a broad range of compounds. In order to accurately reproduce, interpret, and predict experimental results, which are typically obtained in a condensed phase, it is indispensable to include the condensed phase environment in the computational model. However, most studies are still performed in vacuum due to the high computational cost of state-of-the-art non-adiabatic molecular dynamics (NAMD) simulations. The quantum mechanical/molecular mechanical (QM/MM) solvation method has been a popular model to perform photodynamics in the liquid phase. Nevertheless, the currently used QM/MM embedding techniques cannot sufficiently capture all solute-solvent interactions. In this Perspective, we will discuss the efficient ΔSCF electronic structure method and its applications with respect to the NAMD of solvated compounds, with a particular focus on explicit quantum mechanical solvation. As more research is required for this method to reach its full potential, some challenges and possible directions for future research are presented as well.
Collapse
Affiliation(s)
- Eva Vandaele
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
145
|
Xie BB, Jia PK, Wang KX, Chen WK, Liu XY, Cui G. Generalized Ab Initio Nonadiabatic Dynamics Simulation Methods from Molecular to Extended Systems. J Phys Chem A 2022; 126:1789-1804. [PMID: 35266391 DOI: 10.1021/acs.jpca.1c10195] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Nonadiabatic dynamics simulation has become a powerful tool to describe nonadiabatic effects involved in photophysical processes and photochemical reactions. In the past decade, our group has developed generalized trajectory-based ab initio surface-hopping (GTSH) dynamics simulation methods, which can be used to describe a series of nonadiabatic processes, such as internal conversion, intersystem crossing, excitation energy transfer and charge transfer of molecular systems, and photoinduced nonadiabatic carrier dynamics of extended systems with and without spin-orbit couplings. In this contribution, we will first give a brief introduction to our recently developed methods and related numerical implementations at different computational levels. Later, we will present some of our latest applications in realistic systems, which cover organic molecules, biological proteins, organometallic compounds, periodic organic and inorganic materials, etc. Final discussion is given to challenges and outlooks of ab initio nonadiabatic dynamics simulations.
Collapse
Affiliation(s)
- Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Pei-Ke Jia
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Ke-Xin Wang
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231, Zhejiang, P. R. China
| | - Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, Sichuan, P. R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| |
Collapse
|
146
|
Ren J, Li W, Jiang T, Wang Y, Shuai Z. Time‐dependent density matrix renormalization group method for quantum dynamics in complex systems. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1614] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Jiajun Ren
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Weitang Li
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Tong Jiang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Yuanheng Wang
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| | - Zhigang Shuai
- MOE Key Laboratory of Organic OptoElectronics and Molecular Engineering, Department of Chemistry Tsinghua University Beijing People's Republic of China
| |
Collapse
|
147
|
Ma XR, Zhang J, Xiong YC, Zhou W. Revising the performance of the Landau–Zener surface hopping on some typical one-dimensional nonadiabatic models. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2051761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Xiang-Rui Ma
- Department of Material Physics, School of Mathematics, Physics and Optoelectronic Engineering, Hubei University of Automotive Technology, Shiyan, People's Republic of China
- Collaborative Innovation Center for Optoelectronic Technology, Hubei University of Automotive Technology, Shiyan, People's Republic of China
| | - Jun Zhang
- Department of Material Physics, School of Mathematics, Physics and Optoelectronic Engineering, Hubei University of Automotive Technology, Shiyan, People's Republic of China
- Collaborative Innovation Center for Optoelectronic Technology, Hubei University of Automotive Technology, Shiyan, People's Republic of China
| | - Yong-Chen Xiong
- Department of Material Physics, School of Mathematics, Physics and Optoelectronic Engineering, Hubei University of Automotive Technology, Shiyan, People's Republic of China
- Collaborative Innovation Center for Optoelectronic Technology, Hubei University of Automotive Technology, Shiyan, People's Republic of China
| | - Wanghuai Zhou
- Department of Material Physics, School of Mathematics, Physics and Optoelectronic Engineering, Hubei University of Automotive Technology, Shiyan, People's Republic of China
- Collaborative Innovation Center for Optoelectronic Technology, Hubei University of Automotive Technology, Shiyan, People's Republic of China
| |
Collapse
|
148
|
Liang R, Bakhtiiari A. Effects of Enzyme-Ligand Interactions on the Photoisomerization of a Light-Regulated Chemotherapeutic Drug. J Phys Chem B 2022; 126:2382-2393. [PMID: 35297246 DOI: 10.1021/acs.jpcb.1c10819] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular photoswitches permit using light to control protein activity with high spatiotemporal resolutions, thereby alleviating the side effects of conventional chemotherapy. However, due to the challenges in probing ultrafast photoisomerization reactions in biological environments, it remains elusive how the protein influences the photochemistry of the photoswitches, which hampers the rational design of light-regulated therapeutics. To overcome this challenge, we employed first-principles nonadiabatic dynamics simulations to characterize the photodynamics of the phototrexate (PTX), a recently developed photoswitchable anticancer chemotherapeutic that reversibly inhibits its target enzyme dihydrofolate reductase (DHFR). Our simulations show that the protein environment impedes the trans to cis photoisomerization of the PTX. The confinement in the ligand-binding cavity slows down the isomerization kinetics and quantum yield of the photoswitch by reshaping its conical intersection, increasing its excited-state free-energy barrier and quenching its local density fluctuations. Also, the protein environment results in a suboptimal binding mode of the photoproduct that needs to undergo large structural rearrangement to effectively inhibit the enzyme. Therefore, we predict that the PTX's trans → cis photoisomerization in solution precedes its binding with the protein, despite the favorable binding energy of the trans isomer. Our findings highlight the importance of the protein environment on the photochemical reactions of the molecular photoswitches. As such, our work represents an important step toward the rational design of light-regulated drugs in photopharmacology.
Collapse
Affiliation(s)
- Ruibin Liang
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
| | - Amirhossein Bakhtiiari
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States
| |
Collapse
|
149
|
Cornetta LM, Martinez TJ, Varella MTDN. Dissociative electron attachment to 5-bromo-uracil: non-adiabatic dynamics on complex-valued potential energy surfaces. Phys Chem Chem Phys 2022; 24:6845-6855. [PMID: 35253036 DOI: 10.1039/d1cp05663h] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Electron induced dissociation reactions are relevant to many fields, ranging from prebiotic chemistry to cancer treatments. However, the simulation of dissociation electron attachment (DEA) dynamics is very challenging because the auto-ionization widths of the transient negative ions must be accounted for. We propose an adaptation of the ab initio multiple spawning (AIMS) method for complex-valued potential energy surfaces, along the lines of recent developments based on surface hopping dynamics. Our approach combines models for the energy dependence of the auto-ionization widths, obtained from scattering calculations, with survival probabilities computed for the trajectory basis functions employed in the AIMS dynamics. The method is applied to simulate the DEA dynamics of 5-bromo-uracil in full dimensionality, i.e., taking all the vibrational modes into consideration. The propagation starts on the resonance state and describes the formation of Br- anions mediated by non-adiabatic couplings. The potential energies, gradients and non-adiabatic couplings were computed with the fractional-occupancy molecular orbital complete-active-space configuration-interaction method, and the calculated DEA cross section are consistent with the observed DEA intensities.
Collapse
Affiliation(s)
- Lucas M Cornetta
- Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, Rua Sérgio Buarque de Holanda, 777 - Cidade Universitária, Campinas, São Paulo, Brazil.
| | - Todd J Martinez
- Department of Chemistry and PULSE Institute, Stanford University, Stanford, California 94305, USA.,SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA
| | - Marcio T do N Varella
- Instituto de Física, Universidade de São Paulo, Rua do Matão 1371, 05508-090, São Paulo, São Paulo, Brazil
| |
Collapse
|
150
|
Song C. State-averaged CASSCF with polarizable continuum model for studying photoreactions in solvents: Energies, analytical nuclear gradients, and non-adiabatic couplings. J Chem Phys 2022; 156:104102. [DOI: 10.1063/5.0085855] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This paper presents state-averaged complete active space self-consistent field in polarizable continuum model (PCM) for studies of photoreactions in solvents. The wavefunctions of the solute and the PCM surface charges of the solvent are optimized simultaneously such that the state-averaged free energy is variationally minimized. The method supports both fixed weights and dynamic weights where the weights are automatically adjusted based on the energy gaps. The corresponding analytical nuclear gradients and non-adiabatic couplings are also derived. Furthermore, we show how the new method can be entirely formulated in terms of seven basic operations, which allows the implementation to benefit from existing high-performance libraries on graphical processing units. Results demonstrating the accuracy and performance of the implementation are presented and discussed. We also apply the new method to the study of minimal conical intersection search and photoreaction energy pathways in solvents. Effects from the polarity of the solvents and different formulas of dynamic weights are compared and discussed.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry, University of California Davis, Davis, California 95616, USA
| |
Collapse
|