1
|
Huang H, Peng J, Zhang Y, Gu FL, Lan Z, Xu C. The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications. J Chem Phys 2024; 160:234101. [PMID: 38884395 DOI: 10.1063/5.0215036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Accepted: 05/15/2024] [Indexed: 06/18/2024] Open
Abstract
Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.
Collapse
Affiliation(s)
- Haiyi Huang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
- MOE Key Laboratory of Theoretical and Computational Photochemistry, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yulin Zhang
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
2
|
Nulakani NVR, Ali MA. Unveiling the chemical kinetics of aminomethanol (NH 2CH 2OH): insights into O . H and O 2 photo-oxidation reactions and formamide dominance. Front Chem 2024; 12:1407355. [PMID: 38873406 PMCID: PMC11169873 DOI: 10.3389/fchem.2024.1407355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2024] [Accepted: 05/08/2024] [Indexed: 06/15/2024] Open
Abstract
Aminomethanol is released into the atmosphere through various sources, including biomass burning. In this study, we have expounded the chemical kinetics of aminomethanol in the reaction pathways initiated by the hydroxyl radical (O ˙ H) with the aid of ab initio//density functional theory (DFT) i.e., coupled-cluster theory (CCSD(T))//hybrid-DFT (M06-2X/6-311++G (3df, 3pd). We have explored various possible directions of theO ˙ H radical on aminomethanol, as well as the formation of distinct pre-reactive complexes. Our computational findings reveal that the H transfer necessitates activation energies ranging from 4.1 to 6.5 kcal/mol from the -CH2 group, 3.5-6.5 kcal/mol from the -NH2 group and 7-9.3 kcal/mol from the -OH group of three rotational conformers. The H transfer from -CH2, -NH2 and -OH exhibits an estimated total rate constant (k OH) of approximately 1.97 × 10-11 cm3 molecule-1 s-1 at 300 K. The branching fraction analysis indicates a pronounced dominance of C-centered NH2C ˙ HOH radicals with a favorability of 77%, surpassing the N-centeredN ˙ HCH2OH (20%) and O-centered NH2CH2O ˙ (3%) radicals. Moreover, our investigation delves into the oxidation of the prominently favored carbon-centered NH2C ˙ HOH radical through its interaction with atmospheric oxygen molecules. Intriguingly, our findings reveal that formamide (NH2CHO) emerges as the predominant product in the NH2C ˙ HOH + 3O2 reaction, eclipsing alternative outcomes such as amino formic acid (NH2COOH) and formimidic acid (HN = C(H)-OH). At atmospheric conditions pertinent to the troposphere, the branching fraction value for the formation of formamide is about 99%, coupled with a rate constant of 5.5 × 10-12 cm3 molecule-1 s-1. Finally, we have scrutinized the detrimental impact of formamide on the atmosphere. Interaction of formamide with atmospheric hydroxyl radicals could give rise to the production of potentially perilous compounds such as HNCO. Further, unreactedN ˙ HCH2OH radicals may initiate the formation of carcinogenic nitrosamines when reacting with trace N-oxides (namely, NO and NO2). This, in turn, escalates the environmental risk factors.
Collapse
Affiliation(s)
| | - Mohamad Akbar Ali
- Department of Chemistry, Khalifa University of Science and Technology, Abu Dhabi, United Arab Emirates
- Center for the Catalyst and Separations, Khalifa University of Science and Technology, Abu Dhabi, United Arab Emirates
| |
Collapse
|
3
|
Araujo L, Lasser C, Schmidt B. FSSH-2: Fewest Switches Surface Hopping with Robust Switching Probability. J Chem Theory Comput 2024; 20:3413-3419. [PMID: 38696709 DOI: 10.1021/acs.jctc.4c00089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2024]
Abstract
This study introduces the FSSH-2 scheme, a redefined and numerically stable adiabatic Fewest Switches Surface Hopping (FSSH) method for mixed quantum-classical dynamics. It reformulates the standard FSSH hopping probability without using nonadiabatic coupling vectors and allows for numerical time integration with larger step sizes. The advantages of FSSH-2 are demonstrated by numerical experiments for five different model systems in one and two spatial dimensions with up to three electronic states.
Collapse
Affiliation(s)
- Leonardo Araujo
- Department of Mathematics, TUM School of Computation, Information and Technology, Technische Universität München, Boltzmannstr. 3, 85748 Garching bei München, Germany
| | - Caroline Lasser
- Department of Mathematics, TUM School of Computation, Information and Technology, Technische Universität München, Boltzmannstr. 3, 85748 Garching bei München, Germany
| | - Burkhard Schmidt
- Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, 10117 Berlin, Germany
| |
Collapse
|
4
|
Martyka M, Jankowska J. New insights into the photocyclization reaction of a popular diarylethene switch: a nonadiabatic molecular dynamics study. Phys Chem Chem Phys 2024; 26:13383-13394. [PMID: 38646878 DOI: 10.1039/d3cp06256b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/23/2024]
Abstract
Diarylethene (DAE) molecular switches have continued to attract the attention of researchers for over 20 years. Their remarkable photophysical properties endow them with countless applications in photonics and molecular technologies. However, despite extensive experimental and theoretical research, the mechanism of DAE photoswitching is not yet fully rationalized. In this work, we investigate the ring closure dynamics of a popular DAE switch, 1,2-bis(3-methyl-5-phenyl-2 thienyl)perfluorocyclopentene (PT), using nonadiabatic molecular dynamics (NAMD) simulations. Employing the fewest switches surface hopping protocol, along with the semi-empirical multireference ODM2/MRCI-SD method, we investigate possible reaction pathways for this photoprocess, as well as their timescales and resulting photoproducts. Furthermore, using a dynamic configuration-space sampling procedure, we elucidate the role of triplet states in the photocyclization of PT, supporting available experimental data for the closely related DMPT molecule, which indicate an ultrafast intersystem crossing (ISC) transition competing with the singlet-driven photoswitching reaction. Our findings not only corroborate experimental studies on DAE switches, but also provide new mechanistic insights into the potential use in the rational design of DAE switches tailored for specific technological applications.
Collapse
Affiliation(s)
- Mikołaj Martyka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw, 02-093, Poland.
- Interdisciplinary Doctoral School, University of Warsaw, Dobra 56/66, Warsaw, 00-312, Poland
| | - Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw, 02-093, Poland.
| |
Collapse
|
5
|
Miao X, Diemer K, Mitrić R. A CASSCF/MRCI trajectory surface hopping simulation of the photochemical dynamics and the gas phase ultrafast electron diffraction patterns of cyclobutanone. J Chem Phys 2024; 160:124309. [PMID: 38526800 DOI: 10.1063/5.0197768] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2024] [Accepted: 02/29/2024] [Indexed: 03/27/2024] Open
Abstract
We present the simulation of the photochemical dynamics of cyclobutanone induced by the excitation of the 3 s Rydberg state. For this purpose, we apply the complete active space self-consistent field method together with the spin-orbit multireference configuration interaction singles treatment, combined with the trajectory surface hopping for the inclusion of nonadiabatic effects. The simulations were performed in the spin-adiabatic representation, including nine electronic states derived from three singlet and two triplet spin-diabatic states. Our simulations reproduce the two previously observed primary dissociation channels: the C2 pathway yielding C2H4 + CH2CO and the C3 pathway producing c-C3H6 + CO. In addition, two secondary products, CH2 + CO from the C2 pathway and C3H6 from the C3 pathway, both of them previously reported, are also observed in our simulation. We determine the ratio of the C3:C2 products to be about 2.8. Our findings show that most of the trajectories reach their electronic ground state within 200 fs, with dissociation events finished after 300 fs. We also identify the minimum energy conical intersections that are responsible for the relaxation and provide an analysis of the photochemical reaction mechanism based on multidimensional scaling. Furthermore, we demonstrate a minimal impact of triplet states on the photodissociation mechanism within the observed timescale. In order to provide a direct link to experiments, we simulate the gas phase ultrafast electron diffraction patterns and connect their features to the underlying structural dynamics.
Collapse
Affiliation(s)
- Xincheng Miao
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, 97074 Würzburg, Germany
| | - Kira Diemer
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, 97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, 97074 Würzburg, Germany
| |
Collapse
|
6
|
Pieri E, Weingart O, Huix-Rotllant M, Ledentu V, Garavelli M, Ferré N. Modeling pH-Dependent Biomolecular Photochemistry. J Chem Theory Comput 2024; 20:842-855. [PMID: 38198619 DOI: 10.1021/acs.jctc.3c00980] [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: 01/12/2024]
Abstract
The tuning mechanism of pH can be extremely challenging to model computationally in complex biological systems, especially with respect to the photochemical properties. This article reports a protocol aimed at modeling pH-dependent photodynamics using a combination of constant-pH molecular dynamics and semiclassical nonadiabatic molecular dynamics simulations. With retinal photoisomerization in Anabaena sensory rhodopsin (ASR) as a testbed, we show that our protocol produces pH-dependent photochemical properties, such as the isomerization quantum yield or decay rates. We decompose our results into single-titrated residue contributions, identifying some key tuning amino acids. Additionally, we assess the validity of the single protonation state picture to represent the system at a given pH and propose the most populated protein charge state as a compromise between cost and accuracy.
Collapse
Affiliation(s)
- Elisa Pieri
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, 13013 Marseille, France
| | - Oliver Weingart
- Faculty of Mathematics and Natural Sciences, Institute for Theoretical and Computational Chemistry, Heinrich Heine University, Universitätsstr. 1, 40225 Düsseldorf, Germany
| | - Miquel Huix-Rotllant
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, 13013 Marseille, France
| | - Vincent Ledentu
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, 13013 Marseille, France
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Nicolas Ferré
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, 13013 Marseille, France
| |
Collapse
|
7
|
Weight BM, Li X, Zhang Y. Theory and modeling of light-matter interactions in chemistry: current and future. Phys Chem Chem Phys 2023; 25:31554-31577. [PMID: 37842818 DOI: 10.1039/d3cp01415k] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2023]
Abstract
Light-matter interaction not only plays an instrumental role in characterizing materials' properties via various spectroscopic techniques but also provides a general strategy to manipulate material properties via the design of novel nanostructures. This perspective summarizes recent theoretical advances in modeling light-matter interactions in chemistry, mainly focusing on plasmon and polariton chemistry. The former utilizes the highly localized photon, plasmonic hot electrons, and local heat to drive chemical reactions. In contrast, polariton chemistry modifies the potential energy curvatures of bare electronic systems, and hence their chemistry, via forming light-matter hybrid states, so-called polaritons. The perspective starts with the basic background of light-matter interactions, molecular quantum electrodynamics theory, and the challenges of modeling light-matter interactions in chemistry. Then, the recent advances in modeling plasmon and polariton chemistry are described, and future directions toward multiscale simulations of light-matter interaction-mediated chemistry are discussed.
Collapse
Affiliation(s)
- Braden M Weight
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
- Department of Physics and Astronomy, University of Rochester, Rochester, NY, 14627, USA
| | - Xinyang Li
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| | - Yu Zhang
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| |
Collapse
|
8
|
Polonius S, Zhuravel O, Bachmair B, Mai S. LVC/MM: A Hybrid Linear Vibronic Coupling/Molecular Mechanics Model with Distributed Multipole-Based Electrostatic Embedding for Highly Efficient Surface Hopping Dynamics in Solution. J Chem Theory Comput 2023; 19:7171-7186. [PMID: 37788824 PMCID: PMC10601485 DOI: 10.1021/acs.jctc.3c00805] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Indexed: 10/05/2023]
Abstract
We present a theoretical framework for a hybrid linear vibronic coupling model electrostatically embedded into a molecular mechanics environment, termed the linear vibronic coupling/molecular mechanics (LVC/MM) method, for the surface hopping including arbitrary coupling (SHARC) molecular dynamics package. Electrostatic embedding is realized through the computation of interactions between environment point charges and distributed multipole expansions (DMEs, up to quadrupoles) that represent each electronic state and transition densities in the diabatic basis. The DME parameters are obtained through a restrained electrostatic potential (RESP) fit, which we extended to yield higher-order multipoles. We also implemented in SHARC a scheme for achieving roto-translational invariance of LVC models as well as a general quantum mechanics/molecular mechanics (QM/MM) interface, an OpenMM interface, and restraining potentials for simulating liquid droplets. Using thioformaldehyde in water as a test case, we demonstrate that LVC/MM can accurately reproduce the solvation structure and energetics of rigid solutes, with errors on the order of 1-2 kcal/mol compared to a BP86/MM reference. The implementation in SHARC is shown to be very efficient, enabling the simulation of trajectories on the nanosecond time scale in a matter of days.
Collapse
Affiliation(s)
- Severin Polonius
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
| | - Oleksandra Zhuravel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Vienna
Doctoral School in Chemistry (DoSChem), University of Vienna, Währinger Str. 42, 1090 Vienna, Austria
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| |
Collapse
|
9
|
Puentes-Milián L, García-Alfonso E, Márquez-Mijares M, Rubayo-Soneira J. NeI 2 Photofragmentation Dynamics Through Quasiclassical and Semiclassical Studies. Chemphyschem 2023; 24:e202300406. [PMID: 37503677 DOI: 10.1002/cphc.202300406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Revised: 07/26/2023] [Accepted: 07/26/2023] [Indexed: 07/29/2023]
Abstract
The triatomic system NeI2 is studied under the consideration that the diatom is found in an excited electronic state (B). The vibrational levels (v=13, …, 23) are considered within two well-known theoretical procedures: quasi-classical trajectories (QCT), where the classical equations of motion for nuclei are solved on a single potential energy surface (PES), and the trajectory surface hopping (TSH) method, where the same are solved in a bunch of crossed vibrational PES (diabatic representation). The trajectory surface hopping fewest switches (TSHFS) is implemented to minimize the number of hoppings, thus allowing the calculations of hopping probability between the different PES's, and the kinetic mechanism to track the dissociation path. From these calculations, several observables such as, the lifetimes, vibrational and rotational energies (I2 ), dissociation channels, are obtained. Our results are compared with previous experimental and theoretical work.
Collapse
Affiliation(s)
- Loidel Puentes-Milián
- Instituto Superior de Tecnologías y Ciencias Aplicadas (InSTEC), Universidad de La Habana (UH), Ave. Salvador Allende 1110, Plaza de la Revolución, 10400, La Habana, Cuba
| | - Ernesto García-Alfonso
- Université Toulouse 3 and CNRS, Laboratoire des Collisions, Agrégats et Réactivité, IRSAMC, 118 Route de Narbonne, 31062, Toulouse, Cedex 09, France
| | - Maykel Márquez-Mijares
- Instituto Superior de Tecnologías y Ciencias Aplicadas (InSTEC), Universidad de La Habana (UH), Ave. Salvador Allende 1110, Plaza de la Revolución, 10400, La Habana, Cuba
| | - Jesús Rubayo-Soneira
- Instituto Superior de Tecnologías y Ciencias Aplicadas (InSTEC), Universidad de La Habana (UH), Ave. Salvador Allende 1110, Plaza de la Revolución, 10400, La Habana, Cuba
| |
Collapse
|
10
|
Zhao X, Merritt ICD, Lei R, Shu Y, Jacquemin D, Zhang L, Xu X, Vacher M, Truhlar DG. Nonadiabatic Coupling in Trajectory Surface Hopping: Accurate Time Derivative Couplings by the Curvature-Driven Approximation. J Chem Theory Comput 2023; 19:6577-6588. [PMID: 37772732 DOI: 10.1021/acs.jctc.3c00813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/30/2023]
Abstract
Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.
Collapse
Affiliation(s)
- Xiaorui Zhao
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
- School of Aerospace Engineering, Tsinghua University, Beijing 100084, P. R. China
| | | | - Ruiqing Lei
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
| | - Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
- Institut Universitaire de France (IUF), Paris 75005, France
| | - Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, P. R. China
| | - Xuefei Xu
- Center for Combustion Energy, Tsinghua University, Beijing 100084, P. R. China
- Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, P. R. China
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, Nantes F-44000, France
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
11
|
Chen WK, Wang SR, Liu XY, Fang WH, Cui G. Nonadiabatic Derivative Couplings Calculated Using Information of Potential Energy Surfaces without Wavefunctions: Ab Initio and Machine Learning Implementations. Molecules 2023; 28:molecules28104222. [PMID: 37241962 DOI: 10.3390/molecules28104222] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Revised: 05/16/2023] [Accepted: 05/18/2023] [Indexed: 05/28/2023] Open
Abstract
In this work, we implemented an approximate algorithm for calculating nonadiabatic coupling matrix elements (NACMEs) of a polyatomic system with ab initio methods and machine learning (ML) models. Utilizing this algorithm, one can calculate NACMEs using only the information of potential energy surfaces (PESs), i.e., energies, and gradients as well as Hessian matrix elements. We used a realistic system, namely CH2NH, to compare NACMEs calculated by this approximate PES-based algorithm and the accurate wavefunction-based algorithm. Our results show that this approximate PES-based algorithm can give very accurate results comparable to the wavefunction-based algorithm except at energetically degenerate points, i.e., conical intersections. We also tested a machine learning (ML)-trained model with this approximate PES-based algorithm, which also supplied similarly accurate NACMEs but more efficiently. The advantage of this PES-based algorithm is its significant potential to combine with electronic structure methods that do not implement wavefunction-based algorithms, low-scaling energy-based fragment methods, etc., and in particular efficient ML models, to compute NACMEs. The present work could encourage further research on nonadiabatic processes of large systems simulated by ab initio nonadiabatic dynamics simulation methods in which NACMEs are always required.
Collapse
Affiliation(s)
- Wen-Kai Chen
- Hebei Key Laboratory of Inorganic Nano-Materials, College of Chemistry and Materials Science, Hebei Normal University, Shijiazhuang 050024, China
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Sheng-Rui Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Hefei National Laboratory, Hefei 230088, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
- Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
12
|
Weight BM, Sifain AE, Gifford BJ, Htoon H, Tretiak S. On-the-Fly Nonadiabatic Dynamics Simulations of Single-Walled Carbon Nanotubes with Covalent Defects. ACS NANO 2023; 17:6208-6219. [PMID: 36972076 DOI: 10.1021/acsnano.2c08579] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Single-walled carbon nanotubes (SWCNTs) with covalent surface defects have been explored recently due to their promise for use in single-photon telecommunication emission and in spintronic applications. The all-atom dynamic evolution of electrostatically bound excitons (the primary electronic excitations) in these systems has only been loosely explored from a theoretical perspective due to the size limitations of these large systems (>500 atoms). In this work, we present computational modeling of nonradiative relaxation in a variety of SWCNT chiralities with single-defect functionalizations. Our excited-state dynamics modeling uses a trajectory surface hopping algorithm accounting for excitonic effects with a configuration interaction approach. We find a strong chirality and defect-composition dependence on the population relaxation (varying over 50-500 fs) between the primary nanotube band gap excitation E11 and the defect-associated, single-photon-emitting E11* state. These simulations give direct insight into the relaxation between the band-edge states and the localized excitonic state, in competition with dynamic trapping/detrapping processes observed in experiment. Engineering fast population decay into the quasi-two-level subsystem with weak coupling to higher-energy states increases the effectiveness and controllability of these quantum light emitters.
Collapse
Affiliation(s)
- Braden M Weight
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, United States
- Center for Integrated Nanotechnologies, Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Andrew E Sifain
- Department of Chemistry, Princeton University, Princeton, New Jersey 08540 United States
| | - Brendan J Gifford
- Center for Integrated Nanotechnologies, Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Han Htoon
- Center for Integrated Nanotechnologies, Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Center for Integrated Nanotechnologies, Center for Nonlinear Studies, and Theoretical Division Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
13
|
Martyka M, Jankowska J. Nonadiabatic molecular dynamics study of a complete photoswitching cycle for a full-size diarylethene system. J Photochem Photobiol A Chem 2023. [DOI: 10.1016/j.jphotochem.2022.114513] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
14
|
Reiner M, Bachmair B, Tiefenbacher MX, Mai S, González L, Marquetand P, Dellago C. Nonadiabatic Forward Flux Sampling for Excited-State Rare Events. J Chem Theory Comput 2023; 19:1657-1671. [PMID: 36856706 PMCID: PMC10061683 DOI: 10.1021/acs.jctc.2c01088] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Indexed: 03/02/2023]
Abstract
We present a rare event sampling scheme applicable to coupled electronic excited states. In particular, we extend the forward flux sampling (FFS) method for rare event sampling to a nonadiabatic version (NAFFS) that uses the trajectory surface hopping (TSH) method for nonadiabatic dynamics. NAFFS is applied to two dynamically relevant excited-state models that feature an avoided crossing and a conical intersection with tunable parameters. We investigate how nonadiabatic couplings, temperature, and reaction barriers affect transition rate constants in regimes that cannot be otherwise obtained with plain, traditional TSH. The comparison with reference brute-force TSH simulations for limiting cases of rareness shows that NAFFS can be several orders of magnitude cheaper than conventional TSH and thus represents a conceptually novel tool to extend excited-state dynamics to time scales that are able to capture rare nonadiabatic events.
Collapse
Affiliation(s)
- Madlen
Maria Reiner
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Physics, University of
Vienna, Boltzmanngasse
5, 1090 Vienna, Austria
| | - Brigitta Bachmair
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry, University
of Vienna, Währinger
Strasse 42, 1090 Vienna, Austria
| | - Maximilian Xaver Tiefenbacher
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Doctoral School in Chemistry, University
of Vienna, Währinger
Strasse 42, 1090 Vienna, Austria
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Leticia González
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Christoph Dellago
- Research
Platform on Accelerating Photoreaction Discovery (ViRAPID), University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Faculty
of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| |
Collapse
|
15
|
Tichý O, Pederzoli M, Pittner J, Burda JV. Vertical Excitation Energies and Lifetimes of the Two Lowest Singlet Excited States of Cytosine, 5-Aza-cytosine, and the Triazine Family: Quantum Mechanics-Molecular Mechanics Studies. J Chem Theory Comput 2023; 19:1976-1985. [PMID: 36961980 PMCID: PMC10100535 DOI: 10.1021/acs.jctc.2c01262] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/26/2023]
Abstract
A swarm of semi-classical quantum mechanics/molecular mechanics molecular-dynamics simulations where OM2/MNDO is combined with the Gromacs program for consideration of explicit water is performed, solving the time-dependent Schrödinger equation in each step of the trajectories together with the Tully's fewest switches algorithm. Within this stochastic treatment, time dependent probabilities of the three lowest electronic states are determined. The fact that nucleobases are quickly deactivated is confirmed in the cytosine case where our best lifetime estimation is τ1=0.82 ps for the model with 100 water molecules with the SPCE force field and a time step of 0.1 fs. Lifetimes of the remaining molecules are visibly longer: 5-azacytosine, 2,4-diamino-1,3,5-triazine (DT), and 2,4,6-triamino-1,3,5-triazine (TT) molecules have an S1 → S0 de-excitation time of slightly above 10 ps. The lifetimes of the triazine family increases with the increasing number of exocyclic amino groups, that is, s-triazine < 2-amino-1,3,5-triazine < DT < TT. This can be explained by a higher mobility of the carbon-bonded hydrogen atoms in comparison with heavier amino groups since their movement is slowed down due to a substantially higher mass than hydrogen atoms, which can easier reach the out-of-plane positions required in the conical intersection structures. Moreover, bulkier NH2 ligands suffer due to greater friction caused by the surrounding water environment. These mechanical aspects caused a change in the explored lifetime dependences in comparison with our previous gas-phase study.
Collapse
Affiliation(s)
- Ondřej Tichý
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | - Marek Pederzoli
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences, Dolejškova 3, 182 23 Prague 8, Czech Republic
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences, Dolejškova 3, 182 23 Prague 8, Czech Republic
| | - Jaroslav V Burda
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| |
Collapse
|
16
|
Merritt ICD, Jacquemin D, Vacher M. Nonadiabatic Coupling in Trajectory Surface Hopping: How Approximations Impact Excited-State Reaction Dynamics. J Chem Theory Comput 2023; 19:1827-1842. [PMID: 36897995 DOI: 10.1021/acs.jctc.2c00968] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023]
Abstract
Photochemical reactions are widely modeled using the popular trajectory surface hopping (TSH) method, an affordable mixed quantum-classical approximation to the full quantum dynamics of the system. TSH is able to account for nonadiabatic effects using an ensemble of trajectories, which are propagated on a single potential energy surface at a time and which can hop from one electronic state to another. The occurrences and locations of these hops are typically determined using the nonadiabatic coupling between electronic states, which can be assessed in a number of ways. In this work, we benchmark the impact of some approximations to the coupling term on the TSH dynamics for several typical isomerization and ring-opening reactions. We have identified that two of the schemes tested, the popular local diabatization scheme and a scheme based on biorthonormal wave function overlap implemented in the OpenMOLCAS code as part of this work, reproduce at a much reduced cost the dynamics obtained using the explicitly calculated nonadiabatic coupling vectors. The other two schemes tested can give different results, and in some cases, even entirely incorrect dynamics. Of these two, the scheme based on configuration interaction vectors gives unpredictable failures, while the other scheme based on the Baeck-An approximation systematically overestimates hopping to the ground state as compared to the reference approaches.
Collapse
Affiliation(s)
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Morgane Vacher
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
17
|
Zhou Z, Wu Y, Bian X, Subotnik JE. Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping. J Chem Theory Comput 2023; 19:718-732. [PMID: 36655857 DOI: 10.1021/acs.jctc.2c00948] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Nonadiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (where the electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase-space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Thus, the F-PSSH algorithm sets the stage for future modeling of nonadiabatic dynamics under strong externally pumped circular polarization.
Collapse
Affiliation(s)
- Zeyu Zhou
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Yanze Wu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joseph Eli Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
18
|
Wu X, Liu B, Frauenheim T, Tretiak S, Yam C, Zhang Y. Investigation of plasmon relaxation mechanisms using nonadiabatic molecular dynamics. J Chem Phys 2022; 157:214201. [DOI: 10.1063/5.0127435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
Hot carriers generated from the decay of plasmon excitation can be harvested to drive a wide range of physical or chemical processes. However, their generation efficiency is limited by the concomitant phonon-induced relaxation processes by which the energy in excited carriers is transformed into heat. However, simulations of dynamics of nanoscale clusters are challenging due to the computational complexity involved. Here, we adopt our newly developed Trajectory Surface Hopping (TSH) nonadiabatic molecular dynamics algorithm to simulate plasmon relaxation in Au20 clusters, taking the atomistic details into account. The electronic properties are treated within the Linear Response Time-Dependent Tight-binding Density Functional Theory (LR-TDDFTB) framework. The relaxation of plasmon due to coupling to phonon modes in Au20 beyond the Born–Oppenheimer approximation is described by the TSH algorithm. The numerically efficient LR-TDDFTB method allows us to address a dense manifold of excited states to ensure the inclusion of plasmon excitation. Starting from the photoexcited plasmon states in Au20 cluster, we find that the time constant for relaxation from plasmon excited states to the lowest excited states is about 2.7 ps, mainly resulting from a stepwise decay process caused by low-frequency phonons of the Au20 cluster. Furthermore, our simulations show that the lifetime of the phonon-induced plasmon dephasing process is ∼10.4 fs and that such a swift process can be attributed to the strong nonadiabatic effect in small clusters. Our simulations demonstrate a detailed description of the dynamic processes in nanoclusters, including plasmon excitation, hot carrier generation from plasmon excitation dephasing, and the subsequent phonon-induced relaxation process.
Collapse
Affiliation(s)
- Xiaoyan Wu
- Shenzhen JL Computational Science and Applied Research Institute, Longhua District, Shenzhen 518110, China
| | - Baopi Liu
- Shenzhen JL Computational Science and Applied Research Institute, Longhua District, Shenzhen 518110, China
| | - Thomas Frauenheim
- Shenzhen JL Computational Science and Applied Research Institute, Longhua District, Shenzhen 518110, China
- Beijing Computational Science Research Center, Haidian District, Beijing 100193, China
- Bremen Center for Computational Materials Science, University of Bremen, Am Fallturm 1, 28359 Bremen, Germany
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
- Center of Integrated Nanotechnlogies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - ChiYung Yam
- Shenzhen Institute for Advanced Study, University of Electronic Science and Technology of China, Shenzhen 518000, China
- Hong Kong Quantum AI Lab Limited, Hong Kong, China
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
19
|
Yue L. Trajectory surface hopping molecular dynamics on Chemiluminescence of cyclic peroxides. J CHIN CHEM SOC-TAIP 2022. [DOI: 10.1002/jccs.202200329] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Ling Yue
- Key Laboratory for Non‐Equilibrium Synthesis and Modulation of Condensed Matter, Ministry of Education, School of Chemistry Xi'an Jiaotong University Xi'an China
| |
Collapse
|
20
|
Ma J, Zhao D, Jiang C, Lan Z, Li F. Effect of Temperature on Photoisomerization Dynamics of a Newly Designed Two-Stroke Light-Driven Molecular Rotary Motor. Int J Mol Sci 2022; 23:ijms23179694. [PMID: 36077091 PMCID: PMC9456002 DOI: 10.3390/ijms23179694] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2022] [Revised: 08/20/2022] [Accepted: 08/23/2022] [Indexed: 11/16/2022] Open
Abstract
The working mechanism of conventional light-driven molecular rotary motors, especially Feringa-type motors, contains two photoisomerization steps and two thermal helix inversion steps. Due to the existence of a thermal helix inversion step, both the ability to work at lower temperatures and the rotation speed are limited. In this work, a two-stroke light-driven molecular rotary motor, 2-(1,5-dimethyl-4,5-dihydrocyclopenta[b]pyrrol-6(1H)-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDPY), is proposed, which is capable of performing unidirectional and repetitive rotation by only two photoisomerization (EP→ZP and ZP→EP) steps. With trajectory surface-hopping simulation at the semi-empirical OM2/MRCI level, the EP→ZP and ZP→EP nonadiabatic dynamics of DDPY were systematically studied at different temperatures. Both EP→ZP and ZP→EP photoisomerizations are on an ultrafast timescale (ca. 200–300 fs). The decay mode of EP→ZP photoisomerization is approximately bi-exponential, while that of ZP→EP photoisomerization is found to be periodic. For EP and ZP isomers of DDPY, after the S0→S1 excitation, the dynamical processes of nonadiabatic decay are both followed by twisting about the central C=C double bond and the pyramidalization of the C atom at the stator-axle linkage. The effect of temperature on the nonadiabatic dynamics of EP→ZP and ZP→EP photoisomerizations of DDPY has been systematically investigated. The average lifetimes of the S1 excited state and quantum yields for both EP→ZP and ZP→EP photoisomerization are almost temperature-independent, while the corresponding unidirectionality of rotation is significantly increased (e.g., 74% for EP→ZP and 72% for ZP→EP at 300 K vs 100% for EP→ZP and 94% for ZP→EP at 50 K) with lowering the temperature.
Collapse
Affiliation(s)
- Jianzheng Ma
- Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi’an Jiaotong University, Xi’an 710049, China
| | - Di Zhao
- Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi’an Jiaotong University, Xi’an 710049, China
| | - Chenwei Jiang
- Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi’an Jiaotong University, Xi’an 710049, China
- Correspondence: (C.J.); (Z.L.)
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, China
- Correspondence: (C.J.); (Z.L.)
| | - Fuli Li
- Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi’an Jiaotong University, Xi’an 710049, China
| |
Collapse
|
21
|
He X, Wu B, Rivlin T, Liu J, Pollak E. Transition Path Flight Times and Nonadiabatic Electronic Transitions. J Phys Chem Lett 2022; 13:6966-6974. [PMID: 35877977 PMCID: PMC9358656 DOI: 10.1021/acs.jpclett.2c01425] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Accepted: 07/18/2022] [Indexed: 06/15/2023]
Abstract
Transition path flight times are studied for scattering on two electronic surfaces with a single crossing. These flight times reveal nontrivial quantum effects such as resonance lifetimes and nonclassical passage times and reveal that nonadiabatic effects often increase flight times. The flight times are computed using numerically exact time propagation and compared with results obtained from the Fewest Switches Surface Hopping (FSSH) method. Comparison of the two methods shows that the FSSH method is reliable for transition path times only when the scattering is classically allowed on the relevant adiabatic surfaces. However, where quantum effects such as tunneling and resonances dominate, the FSSH method is not adequate to accurately predict the correct times and transition probabilities. These results highlight limitations in methods which do not account for quantum interference effects, and suggest that measuring flight times is important for obtaining insights from the time-domain into quantum effects in nonadiabatic scattering.
Collapse
Affiliation(s)
- Xin He
- Beijing
National Laboratory for Molecular Sciences, Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Baihua Wu
- Beijing
National Laboratory for Molecular Sciences, Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Tom Rivlin
- Chemical
and Biological Physics Department, Weizmann
Institute of Science, 76100 Rehovot, Israel
| | - Jian Liu
- Beijing
National Laboratory for Molecular Sciences, Institute of Theoretical
and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Eli Pollak
- Chemical
and Biological Physics Department, Weizmann
Institute of Science, 76100 Rehovot, Israel
| |
Collapse
|
22
|
Wu X, Wen S, Song H, Frauenheim T, Tretiak S, Yam C, Zhang Y. Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method. J Chem Phys 2022; 157:084114. [DOI: 10.1063/5.0100339] [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
Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. The nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.
Collapse
Affiliation(s)
- Xiaoyan Wu
- Shenzhen JL Computational Science and Applied Research Institute, Shenzhen JL Computational Science and Applied Research Institute, China
| | | | - Huajing Song
- Los Alamos National Laboratory, United States of America
| | | | - Sergei Tretiak
- Theoretical Division, T-1, Los Alamos National Laboratory, United States of America
| | - ChiYung Yam
- Beijing Computational Science Research Center, Beijing Computational Science Research Center, China
| | - Yu Zhang
- Theoretical Division, Los Alamos National Laboratory, United States of America
| |
Collapse
|
23
|
Mukherjee S, Hazra S, Ghosh S, Mukherjee S, Adhikari S. Trajectory Surface Hopping vs. Quantum Scattering Calculations on D+ + H2 and H + H2+ Reactions using Ab Initio Surfaces and Couplings. Chem Phys 2022. [DOI: 10.1016/j.chemphys.2022.111588] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
24
|
Xie H, Xu X, Wang L, Zhuang W. Surface hopping dynamics in periodic solid-state materials with a linear vibronic coupling model. J Chem Phys 2022; 156:154116. [PMID: 35459287 DOI: 10.1063/5.0085759] [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 report a surface hopping approach in which the implemented linear vibronic coupling Hamiltonian is constructed and the electronic wavefunction is propagated in the reciprocal space. The parameters of the linear vibronic coupling model, including onsite energies, phonon frequencies, and electron-phonon couplings, are calculated with density-functional theory and density-functional perturbation theory and interpolated in fine sampling points of the Brillouin zone with maximally localized Wannier functions. Using this approach, we studied the relaxation dynamics of the photo-excited hot carrier in a one-dimensional periodic carbon chain. The results show that the completeness of the number of Hilbert space k points and the number of phonon q points plays an important role in the hot carrier relaxation processes. By calculating the relaxation times of hot carriers under different reciprocal space sampling and extrapolating with the stretched-compressed exponential function, the relaxation times of hot electrons and holes in the quasi-continuous energy band are obtained. By considering the feedback effect in the hopping processes and analyzing the time-dependent phonon energy in different normal modes, we found that the long-wave longitudinal optical phonons play a major role in the relaxation dynamics of hot electrons and holes. We, therefore, provided herein an efficient and accurate approach for modeling the photophysical processes in periodic solid-state material systems.
Collapse
Affiliation(s)
- Hua Xie
- Key Laboratory of Strongly-Coupled Quantum Matter Physics, Chinese Academy of Sciences, School of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiaoliang Xu
- Key Laboratory of Strongly-Coupled Quantum Matter Physics, Chinese Academy of Sciences, School of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Wei Zhuang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou, Fujian 350002, China
| |
Collapse
|
25
|
Design and Nonadiabatic Photoisomerization Dynamics Study of a Three-Stroke Light-Driven Molecular Rotary Motor. Int J Mol Sci 2022; 23:ijms23073908. [PMID: 35409268 PMCID: PMC8999534 DOI: 10.3390/ijms23073908] [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: 02/25/2022] [Revised: 03/29/2022] [Accepted: 03/30/2022] [Indexed: 01/01/2023] Open
Abstract
Working cycle of conventional light-driven molecular rotary motors (LDMRMs), especially Feringa-type motors, usually have four steps, two photoisomerization steps, and two thermal helix inversion (THI) steps. THI steps hinder the ability of the motor to operate at lower temperatures and limit the rotation speed of LDMRMs. A three-stroke LDMRM, 2-(2,7-dimethyl-2,3-dihydro-1H-inden-1-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDIY), is proposed, which is capable of completing an unidirectional rotation by two photoisomerization steps and one thermal helix inversion step at room temperature. On the basis of trajectory surface-hopping simulation at the semi-empirical OM2/MRCI level, the EP→ZP and ZP→EM nonadiabatic photoisomerization dynamics of DDIY were systematically analyzed. Quantum yields of EP→ZP and ZP→EM photoisomerization of DDIY are ca. 34% and 18%, respectively. Both EP→ZP and ZP→EM photoisomerization processes occur on an ultrafast time scale (ca. 100-300 fs). This three-stroke LDMRM may stimulate further research for the development of new families of more efficient LDMRMs.
Collapse
|
26
|
Tichý O, Burda JV. Estimation of electron absorption spectra and lifetime of the two lowest singlet excited states of pyrimidine nucleobases and their derivatives. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.131863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
27
|
Wang Z, Dong J, Qiu J, Wang L. All-Atom Nonadiabatic Dynamics Simulation of Hybrid Graphene Nanoribbons Based on Wannier Analysis and Machine Learning. ACS APPLIED MATERIALS & INTERFACES 2022; 14:22929-22940. [PMID: 35100503 DOI: 10.1021/acsami.1c22181] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Trajectory surface hopping combined with ab initio electronic structure calculations is a popular and powerful approach for on-the-fly nonadiabatic dynamics simulations. For large systems, however, this remains a significant challenge because of the unaffordable computational cost of large-scale electronic structure calculations. Here, we present an efficient divide-and-conquer approach to construct the system Hamiltonian based on Wannier analysis and machine learning. In detail, the large system under investigation is first decomposed into small building blocks, and then all possible segments formed by building blocks within a cutoff distance are found out. Ab initio molecular dynamics is carried out to generate a sequence of geometries for each equivalent segment with periodicity. The Hamiltonian matrices in the maximum localized Wannier function (MLWF) basis are obtained for all geometries and utilized to train artificial neural networks (ANNs) for the structure-dependent Hamiltonian elements. Taking advantage of the orthogonality and spatial locality of MLWFs, the one-electron Hamiltonian of a large system at arbitrary geometry can be directly constructed by the trained ANNs. As demonstrations, we study charge transport in a zigzag graphene nanoribbon (GNR), a coved GNR, and a series of hybrid GNRs with a state-of-the-art surface hopping method. The interplay between delocalized and localized states is found to determine the electron dynamics in hybrid GNRs. Our approach has successfully studied GNRs with >10 000 atoms, paving the way for efficient and reliable all-atom nonadiabatic dynamics simulation of general systems.
Collapse
Affiliation(s)
- Zedong Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jiawei Dong
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Jing Qiu
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| | - Linjun Wang
- Key Laboratory of Excited-State Materials of Zhejiang Province, Department of Chemistry, Zhejiang University, Hangzhou 310027, China
| |
Collapse
|
28
|
Liu J, Lan Z, Yang J. An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics. Phys Chem Chem Phys 2021; 23:22313-22323. [PMID: 34591049 DOI: 10.1039/d1cp03477d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.
Collapse
Affiliation(s)
- Jie Liu
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| |
Collapse
|
29
|
Lin K, Hu D, Peng J, Xu C, Gu FL, Lan Z. Prediction of the excited-state reaction channels in photo-induced processes of nitrofurantoin using first-principle calculations and dynamics simulations. CHEMOSPHERE 2021; 281:130831. [PMID: 34289597 DOI: 10.1016/j.chemosphere.2021.130831] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/21/2021] [Revised: 03/25/2021] [Accepted: 05/06/2021] [Indexed: 06/13/2023]
Abstract
The understanding of the photochemistry of antibiotic compounds is important because it gives the direct information on the possible environmental pollution caused by them. Due to their large size, the theoretical studies of their excited-state reactions are rather challenging. In current work, we combined the on-the-fly trajectory surface-hopping dynamics, conical-intersection optimizations and excited-state pathway calculations to study the photochemistry of the trans-isomer of nitrofurantoin, a widely-used drug to treat the urinary tract infections. The dynamics-then-pathway approach was taken. First the trajectory surface hopping dynamics at the state-averaged complete-active-space self-consistent-field (SA-CASSCF) level with small active space and small basis sets were run. Second, the minimum-energy conical-intersection optimizations were performed. Finally the excited pathways from the Frank-Condon region to different reaction channels were built at the multi-state multi-reference second-order perturbation (MS-CASPT2) level with large active space and large basis set. Several possible channels responsible for the photo-induced reaction mechanism of the trans-nitrofurantoin were obtained, including the cleavage of the NO bond of the NO2 moiety, the photoisomerization at the central CN bond, and other internal conversion channels. Our findings give some preliminary explanations on available experimental observations. It is also demonstrates that the current theoretical approach is a powerful tool to explore the excited-state reactions in the photochemistry of media-sized or large-sized drug compounds.
Collapse
Affiliation(s)
- Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou, 510006, PR China
| | - Deping Hu
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou, 510006, PR China
| | - Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou, 510006, PR China
| | - Chao Xu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou, 510006, PR China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou, 510006, PR China.
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou, 510006, PR China.
| |
Collapse
|
30
|
Freixas VM, Wilhelm P, Nelson T, Hinderer F, Höger S, Tretiak S, Lupton JM, Fernandez-Alberti S. Excitation Energy Transfer between bodipy Dyes in a Symmetric Molecular Excitonic Seesaw. J Phys Chem A 2021; 125:8404-8416. [PMID: 34542292 DOI: 10.1021/acs.jpca.1c06332] [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
We examine the redistribution of energy between electronic and vibrational degrees of freedom that takes place between a π-conjugated oligomer, a phenylene-butadiynylene, and two identical boron-dipyrromethene (bodipy) end-caps using femtosecond transient absorption spectroscopy, single-molecule spectroscopy, and nonadiabatic excited-state molecular dynamics (NEXMD) modeling techniques. The molecular structure represents an excitonic seesaw in that the excitation energy on the oligomer backbone can migrate to either one end-cap or the other, but not to both. The NEXMD simulations closely reproduce the characteristic time scale for redistribution of electronic and vibrational energy of 2.2 ps and uncover the vibrational modes contributing to the intramolecular relaxation. The calculations indicate that the dihedral angle between the bodipy dye and the oligomer change upon excitation of the oligomer. Single-molecule experiments reveal a difference in photoluminescence lifetime of the bodipy dyes depending on whether they are excited by direct absorption or by redistribution of energy from the backbone. This difference in lifetime may be attributed to the difference in dihedral angle. The simulations also suggest that a strong coupling can occur between the two end-caps, giving rise to a reversible shuttling of excitation energy between them. Strong coupling should lead to a pronounced loss in polarization memory of the fluorescence since the oligomer backbone tends to be slightly distorted and the two bodipy transition dipoles have different orientations. A sensitive single-molecule technique is presented to test for such coupling. However, although redistribution of electronic and vibrational energy between the end-caps can occur, it appears to be unidirectional and irreversible, suggesting that an additional localization mechanism is at play which is, as yet, not fully accounted for in the simulations.
Collapse
Affiliation(s)
- Victor M Freixas
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Philipp Wilhelm
- Institut für Angewandte und Experimentelle Physik, Universität Regensburg, Universitätsstrasse 31, 93053 Regensburg, Germany
| | - Tammie Nelson
- Theoretical Division and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Florian Hinderer
- Kekulé-Institut für Organische Chemie und Biochemie, Universität Bonn, Gerhard-Domagk-Straße 1, 53121 Bonn, Germany
| | - Sigurd Höger
- Kekulé-Institut für Organische Chemie und Biochemie, Universität Bonn, Gerhard-Domagk-Straße 1, 53121 Bonn, Germany
| | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - John M Lupton
- Institut für Angewandte und Experimentelle Physik, Universität Regensburg, Universitätsstrasse 31, 93053 Regensburg, Germany
| | | |
Collapse
|
31
|
Posenitskiy E, Spiegelman F, Lemoine D. On application of deep learning to simplified quantum-classical dynamics in electronically excited states. MACHINE LEARNING: SCIENCE AND TECHNOLOGY 2021. [DOI: 10.1088/2632-2153/abfe3f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Abstract
Deep learning (DL) is applied to simulate non-adiabatic molecular dynamics of phenanthrene, using the time-dependent density functional based tight binding (TD-DFTB) approach for excited states combined with mixed quantum–classical propagation. Reference calculations rely on Tully’s fewest-switches surface hopping (FSSH) algorithm coupled to TD-DFTB, which provides electronic relaxation dynamics in fair agreement with various available experimental results. Aiming at describing the coupled electron-nuclei dynamics in large molecular systems, we then examine the combination of DL for excited-state potential energy surfaces (PESs) with a simplified trajectory surface hopping propagation based on the Belyaev–Lebedev (BL) scheme. We start to assess the accuracy of the TD-DFTB approach upon comparison of the optical spectrum with experimental and higher-level theoretical results. Using the recently developed SchNetPack (Schütt et al 2019 J. Chem. Theory Comput.
15 448–55) for DL applications, we train several models and evaluate their performance in predicting excited-state energies and forces. Then, the main focus is given to the analysis of the electronic population of low-lying excited states computed with the aforementioned methods. We determine the relaxation timescales and compare them with experimental data. Our results show that DL demonstrates its ability to describe the excited-state PESs. When coupled to the simplified BL scheme considered in this study, it provides reliable description of the electronic relaxation in phenanthrene as compared with either the experimental data or the higher-level FSSH/TD-DFTB theoretical results. Furthermore, the DL performance allows high-throughput analysis at a negligible cost.
Collapse
|
32
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 154] [Impact Index Per Article: 51.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
33
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
34
|
Song H, Nam Y, Keefer D, Garavelli M, Mukamel S, Tretiak S. Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene. J Phys Chem Lett 2021; 12:5716-5722. [PMID: 34128675 DOI: 10.1021/acs.jpclett.1c01397] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S2/S3(D2d) → S1(D8h) → S0(D4h) → S0(D2d) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. The obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.
Collapse
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
| | - Yeonsig Nam
- Department of Chemistry, University of California, Irvine 92697, California, United States
| | - Daniel Keefer
- Department of Chemistry, University of California, Irvine 92697, California, United States
| | - Marco Garavelli
- Department of Industrial Chemistry, "T. Montanari", Università degli Studi di Bologna, Viale del Risorgimento, 4, Bologna 40136, Italy
| | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine 92697, California, United States
- Department of Physics and Astronomy, University of California, Irvine 92697, California, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos 87545, New Mexico, United States
| |
Collapse
|
35
|
Hu D, Xie Y, Peng J, Lan Z. On-the-Fly Symmetrical Quasi-Classical Dynamics with Meyer-Miller Mapping Hamiltonian for the Treatment of Nonadiabatic Dynamics at Conical Intersections. J Chem Theory Comput 2021; 17:3267-3279. [PMID: 34028268 DOI: 10.1021/acs.jctc.0c01249] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The on-the-fly version of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (SQC/MM) is implemented to study the nonadiabatic dynamics at conical intersections of polyatomic systems. The current on-the-fly implementation of the SQC/MM method is based on the adiabatic representation and the dressed momentum. To include the zero-point energy (ZPE) correction of the electronic mapping variables, we employ both the γ-adjusted and γ-fixed approaches. Nonadiabatic dynamics of the methaniminium cation (CH2NH2+) and azomethane are simulated using the on-the-fly SQC/MM method. For CH2NH2+, both ZPE correction approaches give reasonable and consistent results. However, for azomethane, the γ-adjusted version of the SQC/MM dynamics behaves much better than the γ-fixed version. Further analysis indicates that it is always recommended to use the γ-adjusted SQC/MM dynamics in the on-the-fly simulation of photoinduced dynamics of polyatomic systems, particularly when the excited state is well separated from the ground state in the Franck-Condon region. This work indicates that the on-the-fly SQC/MM method is a powerful simulation protocol to deal with the nonadiabatic dynamics of realistic polyatomic systems.
Collapse
Affiliation(s)
- Deping Hu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Jiawei Peng
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China.,School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
36
|
Barbatti M. Velocity Adjustment in Surface Hopping: Ethylene as a Case Study of the Maximum Error Caused by Direction Choice. J Chem Theory Comput 2021; 17:3010-3018. [PMID: 33844922 DOI: 10.1021/acs.jctc.1c00012] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The most common surface hopping dynamics algorithms require velocity adjustment after hopping to ensure total-energy conservation. Based on the semiclassical analysis, this adjustment must be made parallel to the nonadiabatic coupling vector's direction. Nevertheless, this direction is not always known, and the common practice has been to adjust the velocity in either the linear momentum or velocity directions. This paper benchmarks surface hopping dynamics of photoexcited ethylene with velocity adjustment in several directions, including those of the nonadiabatic coupling vector, the momentum, and the energy gradient difference. It is shown that differences in time constants and structural evolution fall within the statistical uncertainty of the method considering up to 500 trajectories in each dynamics set, rendering the three approaches statistically equivalent. For larger ensembles beyond 1000 trajectories, significant differences between the results arise, limiting the validity of adjustment in alternative directions. Other possible adjustment directions (velocity, single-state gradients, angular momentum) are evaluated as well. Given the small size of ethylene, the results reported in this paper should be considered an upper limit for the error caused by the choice of the velocity-adjustment direction on surface hopping dynamics.
Collapse
|
37
|
Jankowska J, Góra RW. Ultrafast nonradiative deactivation of photoexcited 8-oxo-hypoxanthine: a nonadiabatic molecular dynamics study. Phys Chem Chem Phys 2021; 23:1234-1241. [PMID: 33355573 DOI: 10.1039/d0cp05271j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
In the scientific endeavor to understand the chemical origins of life, the photochemistry of the smallest life building blocks, nucleobases, has been a constant object of focus and intense research. Here, we report the results of the first theoretical study on the photo-properties of an 8-oxo-hypoxanthine molecule, the chromophore of 8-oxo-inosine, which is relevant to the recently proposed, prebiotically plausible synthetic routes to the formation of purine- and pyrimidine-nucleotides. With ab initio and semi-empirical OM2/MRCI quantum-chemistry calculations, we predict a strong photostability of the 8-oxo-hypoxanthine system and see the origin of this effect in ultrafast nonradiative relaxation through puckering of the 6-membered heterocyclic ring.
Collapse
Affiliation(s)
- Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland.
| | - Robert W Góra
- Department of Physical and Quantum Chemistry, Wroclaw University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370, Wrocław, Poland.
| |
Collapse
|
38
|
Darghouth AAMHM, Casida ME, Zhu X, Natarajan B, Su H, Humeniuk A, Titov E, Miao X, Mitrić R. Effect of varying the TD-lc-DFTB range-separation parameter on charge and energy transfer in a model pentacene/buckminsterfullerene heterojunction. J Chem Phys 2021; 154:054102. [PMID: 33557554 DOI: 10.1063/5.0024559] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Atomistic modeling of energy and charge transfer at the heterojunction of organic solar cells is an active field with many remaining outstanding questions owing, in part, to the difficulties in performing reliable photodynamics calculations on very large systems. One approach to being able to overcome these difficulties is to design and apply an appropriate simplified method. Density-functional tight binding (DFTB) has become a popular form of approximate density-functional theory based on a minimal valence basis set and neglect of all but two center integrals. We report the results of our tests of a recent long-range correction (lc) [A. Humeniuk and R. Mitrić, J. Chem. Phys. 143, 134120 (2015)] for time-dependent (TD) lc-DFTB by carrying out TD-lc-DFTB fewest switches surface hopping calculations of energy and charge transfer times using the relatively new DFTBABY [A. Humeniuk and R. Mitrić, Comput. Phys. Commun. 221, 174 (2017)] program. An advantage of this method is the ability to run enough trajectories to get meaningful ensemble averages. Our interest in the present work is less in determining exact energy and charge transfer rates than in understanding how the results of these calculations vary with the value of the range-separation parameter (Rlc = 1/μ) for a model organic solar cell heterojunction consisting of a gas-phase van der Waals complex P/F made up of a single pentacene (P) molecule together with a single buckminsterfullerene (F) molecule. The default value of Rlc = 3.03 a0 is found to be much too small as neither energy nor charge transfer is observed until Rlc ≈ 10 a0. Tests at a single geometry show that the best agreement with high-quality ab initio spectra is obtained in the limit of no lc (i.e., very large Rlc). A plot of energy and charge transfer rates as a function of Rlc is provided, which suggests that a value of Rlc ≈ 15 a0 yields the typical literature (condensed-phase) charge transfer time of about 100 fs. However, energy and charge transfer times become as high as ∼300 fs for Rlc ≈ 25 a0. A closer examination of the charge transfer process P*/F → P+/F- shows that the initial electron transfer is accompanied by a partial delocalization of the P hole onto F, which then relocalizes back onto P, consistent with a polaron-like picture in which the nuclei relax to stabilize the resultant redistribution of charges.
Collapse
Affiliation(s)
| | - Mark E Casida
- Laboratoire de Spectrométrie, Interactions et Chimie Théorique (SITh), Département de Chimie Moléculaire (DCM), Institut de Chimie Moléculaire de Grenoble (ICMG), Université Grenoble-Alpes, 301 rue de la Chimie, CS 40700, 38058 Grenoble Cedex 9, France
| | - Xi Zhu
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Bhaarathi Natarajan
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Haibin Su
- Institute of Advanced Studies, Nanyang Technological University, 60 Nanyang View, 639673, Singapore
| | - Alexander Humeniuk
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Evgenii Titov
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Xincheng Miao
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| | - Roland Mitrić
- Institut für Physikalische und Theoretische Chemie, Julius-Maximilians-Universität Würzburg, Emil-Fischer-Straße 42, D-97074 Würzburg, Germany
| |
Collapse
|
39
|
Temen S, Akimov AV. A Simple Solution to Trivial Crossings: A Stochastic State Tracking Approach. J Phys Chem Lett 2021; 12:850-860. [PMID: 33427475 DOI: 10.1021/acs.jpclett.0c03428] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We present a new state tracking algorithm based on a stochastic state reassignment that reflects the quantum mechanical interpretation of the state time-overlaps. We assess the new method with a range of model Hamiltonians and demonstrate that it yields the results generally consistent with the deterministic min-cost algorithm. However, the stochastic state tracking algorithm reduces magnitudes of the state population fluctuations as the quantum system evolves toward its equilibrium. The new algorithm facilitates the thermalization of quantum state populations and suppresses the population revivals and oscillations near the equilibrium in many-state systems. The new stochastic algorithm has a favorable computational scaling, is easy to implement due to its conceptual transparency, and treats various types of state identity changes (trivial or avoided crossings and any intermediate cases) on equal footing.
Collapse
Affiliation(s)
- Story Temen
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| | - Alexey V Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260, United States
| |
Collapse
|
40
|
Westermayr J, Marquetand P. Machine learning and excited-state molecular dynamics. MACHINE LEARNING-SCIENCE AND TECHNOLOGY 2020. [DOI: 10.1088/2632-2153/ab9c3e] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
41
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
42
|
Song H, Fischer SA, Zhang Y, Cramer CJ, Mukamel S, Govind N, Tretiak S. First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem. J Chem Theory Comput 2020; 16:6418-6427. [DOI: 10.1021/acs.jctc.0c00295] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
| | - Sean A. Fischer
- Chemistry Division, U.S. Naval Research Laboratory, Washington, District of Columbia 20375, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
| | - Christopher J. Cramer
- Department of Chemistry, Supercomputing Institute and Chemical Theory Center, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Shaul Mukamel
- Departments of Chemistry, and physics and astronomy, University of California, Irvine, California 92697, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico, 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
43
|
Suchan J, Janoš J, Slavíček P. Pragmatic Approach to Photodynamics: Mixed Landau–Zener Surface Hopping with Intersystem Crossing. J Chem Theory Comput 2020; 16:5809-5820. [DOI: 10.1021/acs.jctc.0c00512] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Jiří Suchan
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Prague, Technická 5, 16628 Prague, Czech Republic
| |
Collapse
|
44
|
Shu Y, Zhang L, Sun S, Truhlar DG. Time-Derivative Couplings for Self-Consistent Electronically Nonadiabatic Dynamics. J Chem Theory Comput 2020; 16:4098-4106. [PMID: 32456433 DOI: 10.1021/acs.jctc.0c00409] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Electronically nonadiabatic dynamics methods based on a self-consistent potential, such as semiclassical Ehrenfest and coherent switching with decay of mixing, have a number of advantages but are computationally slower than approximations based on an unaveraged potential because they require evaluation of all components of the nonadiabatic coupling vector. Here we introduce a new approximation to the self-consistent potential that does not have this computational drawback. The new approximation uses time-derivative couplings evaluated by overlap integrals of electronic wave functions to approximate the nonadiabatic coupling terms in the equations of motion. We present a numerical test of the method for ethylene that shows there is little loss of accuracy in the ensemble-averaged results. This new approximation to the self-consistent potential makes direct dynamics calculations with self-consistent potentials more efficient for complex systems and makes them practically affordable for some cases where the cost was previously too high.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Linyao Zhang
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States.,School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People's Republic of China
| | - Donald G Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
45
|
Malone W, Nebgen B, White A, Zhang Y, Song H, Bjorgaard JA, Sifain AE, Rodriguez-Hernandez B, Freixas VM, Fernandez-Alberti S, Roitberg AE, Nelson TR, Tretiak S. NEXMD Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:5771-5783. [DOI: 10.1021/acs.jctc.0c00248] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Walter Malone
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Benjamin Nebgen
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Alexander White
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Yu Zhang
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Josiah A. Bjorgaard
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Andrew E. Sifain
- U.S. Army Research Laboratory, Aberdeen Proving Ground, Aberdeen, Maryland 21005, United States
| | | | - Victor M. Freixas
- Universidad Nacional de Quilmes/CONICET, Roque Saenz Peña 352, B1876BXD Bernal, Argentina
| | | | - Adrian E. Roitberg
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Tammie R. Nelson
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
46
|
Photoinduced Double Proton Transfer in the Glyoxal-Methanol Complex Revisited: The Role of the Excited States. J Chem Theory Comput 2020; 16:3273-3286. [PMID: 32275423 DOI: 10.1021/acs.jctc.0c00007] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Under irradiation in the visible range, the glyoxal-methanol complex in a cryogenic argon matrix undergoes a double proton transfer (DPT) reaction through which the glyoxal molecule isomerizes into hydroxyketene. In this work, we employ electronic structure simulations in order to shed more light on the underlying mechanism. Rewardingly, we find that the lowest singlet excited state (S1) of the complex acts as a gateway to two previously unknown isomerization pathways, of which one takes place entirely in the singlet manifold and the other also involves the lowest triplet state (T1). Both of these pathways are fully compatible with the available experimental data, implying that either or both are operative under experimental conditions. In either pathway, the methanol molecule acts as a proton shuttle between the proton-donating and proton-accepting sites of glyoxal, resulting in a dramatic lowering of the potential energy barrier to isomerization with respect to the case of isolated glyoxal. The occurrence of DPT in the singlet manifold is demonstrated directly with the use of nonadiabatic molecular dynamics simulations at the spin-flip time-dependent density functional theory level.
Collapse
|
47
|
Campetella M, Sanz García J. Following the evolution of excited states along photochemical reaction pathways. J Comput Chem 2020; 41:1156-1164. [DOI: 10.1002/jcc.26162] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2019] [Revised: 12/14/2019] [Accepted: 01/14/2020] [Indexed: 01/04/2023]
Affiliation(s)
- Marco Campetella
- Sorbonne Université, CNRS, Institut des Nanosciences de Paris, UMR7588 F‐75005 Paris France
| | - Juan Sanz García
- Sorbonne Université, Laboratoire de Chimie Théorique, UPMC Paris 06, UMR7616 F‐75005 Paris France
| |
Collapse
|
48
|
Lischka H, Shepard R, Müller T, Szalay PG, Pitzer RM, Aquino AJA, Araújo do Nascimento MM, Barbatti M, Belcher LT, Blaudeau JP, Borges I, Brozell SR, Carter EA, Das A, Gidofalvi G, González L, Hase WL, Kedziora G, Kertesz M, Kossoski F, Machado FBC, Matsika S, do Monte SA, Nachtigallová D, Nieman R, Oppel M, Parish CA, Plasser F, Spada RFK, Stahlberg EA, Ventura E, Yarkony DR, Zhang Z. The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry. J Chem Phys 2020; 152:134110. [PMID: 32268762 DOI: 10.1063/1.5144267] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.
Collapse
Affiliation(s)
- Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Ron Shepard
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Thomas Müller
- Institute for Advanced Simulation, Jülich Supercomputing Centre, Forschungszentrum Jülich, Jülich 52428, Germany
| | - Péter G Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Budapest, Hungary
| | - Russell M Pitzer
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Adelia J A Aquino
- School of Pharmaceutical Sciences and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | | | | | - Lachlan T Belcher
- Laser and Optics Research Center, Department of Physics, US Air Force Academy, Colorado 80840, USA
| | | | - Itamar Borges
- Departamento de Química, Instituto Militar de Engenharia, Rio de Janeiro, RJ 22290-270, Brazil
| | - Scott R Brozell
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Emily A Carter
- Office of the Chancellor and Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, Box 951405, Los Angeles, California 90095-1405, USA
| | - Anita Das
- Indian Institute of Engineering Science and Technology, Shibpur, Howrah, India
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, USA
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - William L Hase
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Gary Kedziora
- Air Force Research Laboratory, Wright-Patterson Air Force Base, Ohio 45433, USA
| | - Miklos Kertesz
- Department of Chemistry, Georgetown University, 37th and O Streets, NW, Washington, DC 20057-1227, USA
| | | | - Francisco B C Machado
- Departamento de Química, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, 1901 N. 13th St., Philadelphia, Pennsylvania 19122, USA
| | | | - Dana Nachtigallová
- Institute of Organic Chemistry and Biochemistry v.v.i., The Czech Academy of Sciences, Flemingovo nám. 2, 160610 Prague 6, Czech Republic
| | - Reed Nieman
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, USA
| | - Markus Oppel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, 1090 Vienna, Austria
| | - Carol A Parish
- Department of Chemistry, Gottwald Center for the Sciences, University of Richmond, Richmond, Virginia 23173, USA
| | - Felix Plasser
- Department of Chemistry, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Rene F K Spada
- Departamento de Física, Instituto Tecnológico de Aeronáutica, São José dos Campos 12228-900, São Paulo, Brazil
| | - Eric A Stahlberg
- Biomedical Informatics and Data Science, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, USA
| | - Elizete Ventura
- Universidade Federal da Paraíba, 58059-900 João Pessoa, PB, Brazil
| | - David R Yarkony
- Department of Chemistry, Johns Hopkins University, 3400 N. Charles Street, Baltimore, Maryland 21218, USA
| | - Zhiyong Zhang
- Stanford Research Computing Center, Stanford University, 255 Panama Street, Stanford, California 94305, USA
| |
Collapse
|
49
|
Park JW, Al-Saadon R, MacLeod MK, Shiozaki T, Vlaisavljevich B. Multireference Electron Correlation Methods: Journeys along Potential Energy Surfaces. Chem Rev 2020; 120:5878-5909. [PMID: 32239929 DOI: 10.1021/acs.chemrev.9b00496] [Citation(s) in RCA: 54] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Multireference electron correlation methods describe static and dynamical electron correlation in a balanced way and, therefore, can yield accurate and predictive results even when single-reference methods or multiconfigurational self-consistent field theory fails. One of their most prominent applications in quantum chemistry is the exploration of potential energy surfaces. This includes the optimization of molecular geometries, such as equilibrium geometries and conical intersections and on-the-fly photodynamics simulations, both of which depend heavily on the ability of the method to properly explore the potential energy surface. Because such applications require nuclear gradients and derivative couplings, the availability of analytical nuclear gradients greatly enhances the scope of quantum chemical methods. This review focuses on the developments and advances made in the past two decades. A detailed account of the analytical nuclear gradient and derivative coupling theories is presented. Emphasis is given to the software infrastructure that allows one to make use of these methods. Notable applications of multireference electron correlation methods to chemistry, including geometry optimizations and on-the-fly dynamics, are summarized at the end followed by a discussion of future prospects.
Collapse
Affiliation(s)
- Jae Woo Park
- Department of Chemistry, Chungbuk National University, Chungdae-ro 1, Cheongju 28644, Korea
| | - Rachael Al-Saadon
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Matthew K MacLeod
- Workday, 4900 Pearl Circle East, Suite 100, Boulder, Colorado 80301, United States
| | - Toru Shiozaki
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States.,Quantum Simulation Technologies, Inc., 625 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, 414 East Clark Street, Vermillion, South Dakota 57069, United States
| |
Collapse
|
50
|
Zhang YH, Sun XW, Zhang TS, Liu XY, Cui G. Nonadiabatic Dynamics Simulations on Early-Time Photochemistry of Spirobenzopyran. J Phys Chem A 2020; 124:2547-2559. [PMID: 32187492 DOI: 10.1021/acs.jpca.0c00791] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Photoinduced ring-opening, decay, and isomerization of spirobenzopyran have been explored by the OM2/MRCI nonadiabatic dynamics simulations based on Tully's fewest-switches surface hopping scheme. The efficient S1 to S0 internal conversion as observed in experiments is attributed to the existence of two efficient excited-state decay pathways. The first one is related to the C-N dissociation, and the second one is done to the C-O dissociation. The C-O dissociation pathway is dominant, and more than 90% trajectories decay to the S0 state via the C-O bond-fission related S1/S0 conical intersections. Near these regions in the S0 state, trajectories can either return to spirobenzopyran or proceed to various intermediates including merocyanine via a series of bond rotations. Our nonadiabatic dynamics simulations also demonstrate that the hydrogen-out-of-plane (HOOP) motion is important for efficient and ultrafast excited-state deactivation. On the other hand, we have also found that the replacement of methyl groups by hydrogen atoms in spirobenzopyran can artificially introduce different intramolecular hydrogen transfers leading to hydrogen-transferred intermediates. This finding is important for the community and demonstrates that such a kind of structural truncation, sometimes, could be problematic, leading to incorrect photodynamics. Our present work provides valuable insights into the photodynamics of spirobenzopyran, which could be helpful for the design of spiropyran-based photochromic materials.
Collapse
Affiliation(s)
- Ya-Hui Zhang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, 100875, China
| | - Xin-Wei Sun
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, 100875, China
| | - Teng-Shuo Zhang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, 100875, China
| | - Xiang-Yang Liu
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing, 100875, China
| |
Collapse
|