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
|
Peng J, Liu H, Lan Z. The photodissociation dynamics and ultrafast electron diffraction image of cyclobutanone from the surface hopping dynamics simulation. J Chem Phys 2024; 160:224305. [PMID: 38856062 DOI: 10.1063/5.0203462] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Accepted: 05/23/2024] [Indexed: 06/11/2024] Open
Abstract
The comprehension of nonadiabatic dynamics in polyatomic systems relies heavily on the simultaneous advancements in theoretical and experimental domains. The gas-phase ultrafast electron diffraction (UED) technique has attracted significant attention as a unique tool for monitoring photochemical and photophysical processes at the all-atomic level with high temporal and spatial resolutions. In this work, we simulate the UED spectra of cyclobutanone using the trajectory surface hopping method at the extended multi-state complete active space second order perturbation theory (XMS-CASPT2) level and thereby predict the results of the upcoming UED experiments in the Stanford Linear Accelerator Laboratory. The simulated results demonstrate that a few pathways, including the C2 and C3 dissociation channels, as well as the ring opening channel, play important roles in the nonadiabatic reactions of cyclobutanone. We demonstrate that the simulated UED signal can be directly interpreted in terms of atomic motions, which provides a unique way of monitoring the evolution of the molecular structure in real time. Our work not only provides numerical data that help to determine the accuracy of the well-known surface hopping dynamics at the high XMS-CASPT2 electronic-structure level but also facilitates the understanding of the microscopic mechanisms of the photoinduced reactions in cyclobutanone.
Collapse
Affiliation(s)
- Jiawei Peng
- School of Chemistry, South China Normal University, Guangzhou 510006, China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Hong Liu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
3
|
Mannouch JR, Kelly A. Quantum Quality with Classical Cost: Ab Initio Nonadiabatic Dynamics Simulations Using the Mapping Approach to Surface Hopping. J Phys Chem Lett 2024; 15:5814-5823. [PMID: 38781480 PMCID: PMC11163471 DOI: 10.1021/acs.jpclett.4c00535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2024] [Revised: 05/07/2024] [Accepted: 05/09/2024] [Indexed: 05/25/2024]
Abstract
Nonadiabatic dynamics methods are an essential tool for investigating photochemical processes. In the context of employing first-principles electronic structure techniques, such simulations can be carried out in a practical manner using semiclassical trajectory-based methods or wave packet approaches. While all approaches applicable to first-principles simulations are necessarily approximate, it is commonly thought that wave packet approaches offer inherent advantages over their semiclassical counterparts in terms of accuracy and that this trait simply comes at a higher computational cost. Here we demonstrate that the mapping approach to surface hopping (MASH), a recently introduced trajectory-based nonadiabatic dynamics method, can be efficiently applied in tandem with ab initio electronic structure. Our results even suggest that MASH may provide more accurate results than on-the-fly wave packet techniques, all at a much lower computational cost.
Collapse
Affiliation(s)
- Jonathan R. Mannouch
- Hamburg Center for Ultrafast
Imaging, Universität Hamburg and
the Max Planck Institute
for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| | - Aaron Kelly
- Hamburg Center for Ultrafast
Imaging, Universität Hamburg and
the Max Planck Institute
for the Structure and Dynamics of Matter, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
4
|
Shi W, Ma J, Jiang C, Taketsugu T. Advanced theoretical design of light-driven molecular rotary motors: enhancing thermal helix inversion and visible-light activation. Phys Chem Chem Phys 2024; 26:15672-15680. [PMID: 38766713 DOI: 10.1039/d4cp00037d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
In this study, we have advanced the field of light-driven molecular rotary motors (LDMRMs) by achieving two pivotal goals: lowering the thermal helix inversion (THI) barrier and extending the absorption wavelength into the visible spectrum. This study involves the structural reengineering of a second-generation visible LDMRM, resulting in the synthesis of a novel class, specifically, 2-((2S)-5-methoxy-2-methyl-2,3-dihydro-1H-cyclopenta[a]naphthalen-1-yl)-3-oxo-2,3-dihydro-1H-dibenzo[e,g]indole-6,9-dicarbonitrile. This redesigned motor stands out with its two photoisomerization stages and two thermal helix inversions, featuring exceptionally low THI barriers (4.00 and 2.05 kcal mol-1 at the OM2/MRCI level for the EM → EP and ZM → ZP processes, respectively). Moreover, it displays absorption wavelengths in the visible light range (482.98 and 465.76 nm for the EP and ZP isomers, respectively, at the TD-PBE0-D3/6-31G(d,p) level), surpassing its predecessors in efficiency, as indicated by the narrow HOMO-LUMO energy gap. Ultrafast photoisomerization kinetics (approximately 0.8-1.6 ps) and high quantum yields (around 0.3-0.6) were observed through trajectory surface hopping simulations. Additionally, the simulated time-resolved fluorescence emission spectrum indicates a significantly reduced "dark state" duration (0.09-0.26 ps) in these newly designed LDMRMs compared to the original ones, marking a substantial leap forward in the design and efficiency of LDMRMs.
Collapse
Affiliation(s)
- Weiliang Shi
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo 060-0810, Japan
| | - 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.
| | - 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.
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo 060-0810, Japan.
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo 001-0021, Japan
| |
Collapse
|
5
|
Xu C, Lin C, Peng J, Zhang J, Lin S, Gu FL, Gelin MF, Lan Z. On-the-fly simulation of time-resolved fluorescence spectra and anisotropy. J Chem Phys 2024; 160:104109. [PMID: 38477337 DOI: 10.1063/5.0201204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2024] [Accepted: 02/23/2024] [Indexed: 03/14/2024] Open
Abstract
We combine on-the-fly trajectory surface hopping simulations and the doorway-window representation of nonlinear optical response functions to create an efficient protocol for the evaluation of time- and frequency-resolved fluorescence (TFRF) spectra and anisotropies of the realistic polyatomic systems. This approach gives the effective description of the proper (e.g., experimental) pulse envelopes, laser field polarizations, and the proper orientational averaging of TFRF signals directly from the well-established on-the-fly nonadiabatic dynamic simulations without extra computational cost. To discuss the implementation details of the developed protocol, we chose cis-azobenzene as a prototype to simulate the time evolution of the TFRF spectra governed by its nonadiabatic dynamics. The results show that the TFRF is determined by the interplay of several key factors, i.e., decays of excited-state populations, evolution of the transition dipole moments along with the dynamic propagation, and scaling factor of the TFRF signals associated with the cube of emission frequency. This work not only provides an efficient and effective approach to simulate the TFRF and anisotropies of realistic polyatomic systems but also discusses the important relationship between the TFRF signals and the underlining nonadiabatic dynamics.
Collapse
Affiliation(s)
- Chao Xu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety; School of Environment, South China Normal University, Guangzhou 510006, People's Republic of China
| | - Congru Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety; School of Environment, South China Normal University, Guangzhou 510006, People's Republic of China
| | - Jiawei Peng
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety; School of Environment, South China Normal University, Guangzhou 510006, People's Republic of China
| | - Juanjuan Zhang
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety; School of Environment, South China Normal University, Guangzhou 510006, People's Republic of China
| | - Shichen Lin
- Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety; School of Environment, South China Normal University, Guangzhou 510006, People's Republic of China
| | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, People's Republic of China
| | - Zhenggang Lan
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education and Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety; School of Environment, South China Normal University, Guangzhou 510006, People's Republic of China
| |
Collapse
|
6
|
Kim SS, Rhee YM. Potential energy interpolation with target-customized weighting coordinates: application to excited-state dynamics of photoactive yellow protein chromophore in water. Phys Chem Chem Phys 2024; 26:9021-9036. [PMID: 38440829 DOI: 10.1039/d3cp05643k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2024]
Abstract
Interpolation of potential energy surfaces (PESs) can provide a practical route to performing molecular dynamics simulations with a reliability matching a high-level quantum chemical calculation. An obstacle to its widespread use is perhaps the lack of general and optimal interpolation settings that can be applied in a black-box manner for any given molecular system. How to set up the weights for interpolation is one such task, and we still need to diversify the approaches in order to treat various systems. Here, we develop a new interpolation weighting scheme, which allows us to choose the weighting coordinates in a system-specific manner, by amplifying the contribution from specific internal coordinates. The new weighting scheme with an appropriate selection of coordinates is proved to be effective in reducing the interpolation error along the reaction pathway. As a demonstration, we consider the photoactive yellow protein chromophore system, as it constitutes itself as an interesting target that bears long-standing questions related to excited-state dynamics inside protein environments. We build its two-state diabatic interpolated PES with the new weighting scheme. We indeed see the utility of our scheme by conducting nonadiabatic molecular dynamics simulations with the required semi-global PES based on a limited number of data points.
Collapse
Affiliation(s)
- Seung Soo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
7
|
Yao Y, Li Q, Chu W, Ding YM, Yan L, Gao Y, Neogi A, Govorov A, Zhou L, Wang Z. Exploration of the origin of the excellent charge-carrier dynamics in Ruddlesden-Popper oxysulfide perovskite Y 2Ti 2O 5S 2. Phys Chem Chem Phys 2023. [PMID: 38051151 DOI: 10.1039/d3cp02860g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2023]
Abstract
Although the efficient separation of electron-hole (e-h) pairs is one of the most sought-after electronic characteristics of materials, due to thermally induced atomic motion and other factors, they do not remain separated during the carrier transport process, potentially leading to rapid carrier recombination. Here, we utilized real-time time-dependent density functional theory in combination with nonadiabatic molecular dynamics (NAMD) to explore the separated dynamic transport path within Ruddlesden-Popper oxysulfide perovskite Y2Ti2O5S2 caused by the dielectric layer and phonon frequency difference. The underlying origin of the efficient overall water splitting in Y2Ti2O5S2 is systematically explored. We report the existence of the bi-directional e-h separate-path transport, in which, the electrons transport in the Ti2O5 layer and the holes diffuse in the rock-salt layer. This is in contrast to the conventional e-h separated distribution with a crowded transport channel, as observed in SrTiO3 and hybrid perovskites. Such a unique feature finally results in a long carrier lifetime of 321 ns, larger than that in the SrTiO3 perovskite (160 ns) with only one carrier transport channel. This work provides insights into the carrier transport in lead-free perovskites and yields a novel design strategy for next-generation functionalized optoelectronic devices.
Collapse
Affiliation(s)
- Yisen Yao
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China.
- Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Qiaoqiao Li
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Weibin Chu
- Key Laboratory of Computational Physical Sciences (Ministry of Education), Institute of Computational Physical Sciences, Fudan University, Shanghai 200433, China
| | - Yi-Min Ding
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China.
| | - Luo Yan
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Yang Gao
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China.
- Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Arup Neogi
- Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Alexander Govorov
- Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, Chengdu 610054, China
- Department of Physics and Astronomy, Ohio University, Athens, Ohio 45701, USA.
| | - Liujiang Zhou
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China.
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Zhiming Wang
- Yangtze Delta Region Institute (Huzhou), University of Electronic Science and Technology of China, Huzhou 313001, China.
- Institute of Fundamental and Frontier Sciences, University of Electronic Science and Technology of China, Chengdu 610054, China
| |
Collapse
|
8
|
Issler K, Mitrić R, Petersen J. HORTENSIA, a program package for the simulation of nonadiabatic autoionization dynamics in molecules. J Chem Phys 2023; 159:134801. [PMID: 37787145 DOI: 10.1063/5.0167412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2023] [Accepted: 09/11/2023] [Indexed: 10/04/2023] Open
Abstract
We present a program package for the simulation of ultrafast vibration-induced autoionization dynamics in molecular anions in the manifold of the adiabatic anionic states and the discretized ionization continuum. This program, called HORTENSIA (Hopping Real-time Trajectories for Electron-ejection by Nonadiabatic Self-Ionization in Anions), is based on the nonadiabatic surface-hopping methodology, wherein nuclei are propagated as an ensemble along classical trajectories in the quantum-mechanical potential created by the electronic density of the molecular system. The electronic Schrödinger equation is numerically integrated along the trajectory, providing the time evolution of electronic state coefficients, from which switching probabilities into discrete electronic states are determined. In the case of a discretized continuum state, this hopping event is interpreted as the ejection on an electron. The derived diabatic and nonadiabatic couplings in the time-dependent electronic Schrödinger equation are calculated from anionic and neutral wavefunctions obtained from quantum-chemical calculations with commercially available program packages interfaced with our program. Based on this methodology, we demonstrate the simulation of autoionization electron kinetic energy spectra that are both time- and angle-resolved. In addition, the program yields data that can be interpreted easily with respect to geometric characteristics, such as bonding distances and angles, which facilitate the detection of molecular configurations important for the autoionization process. Furthermore, several useful extensions are included, namely, tools for the generation of initial conditions and input files as well as for the evaluation of output files, all of this both through console commands and a graphical user interface.
Collapse
Affiliation(s)
- Kevin Issler
- Julius-Maximilians-Universität Würzburg, Institut für Physikalische und Theoretische Chemie, Emil-Fischer-Str. 42, 97074 Würzburg, Germany
| | - Roland Mitrić
- Julius-Maximilians-Universität Würzburg, Institut für Physikalische und Theoretische Chemie, Emil-Fischer-Str. 42, 97074 Würzburg, Germany
| | - Jens Petersen
- Julius-Maximilians-Universität Würzburg, Institut für Physikalische und Theoretische Chemie, Emil-Fischer-Str. 42, 97074 Würzburg, Germany
| |
Collapse
|
9
|
Freixas VM, Malone W, Li X, Song H, Negrin-Yuvero H, Pérez-Castillo R, White A, Gibson TR, Makhov DV, Shalashilin DV, Zhang Y, Fedik N, Kulichenko M, Messerly R, Mohanam LN, Sharifzadeh S, Bastida A, Mukamel S, Fernandez-Alberti S, Tretiak S. NEXMD v2.0 Software Package for Nonadiabatic Excited State Molecular Dynamics Simulations. J Chem Theory Comput 2023; 19:5356-5368. [PMID: 37506288 DOI: 10.1021/acs.jctc.3c00583] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/30/2023]
Abstract
We present NEXMD version 2.0, the second release of the NEXMD (Nonadiabatic EXcited-state Molecular Dynamics) software package. Across a variety of new features, NEXMD v2.0 incorporates new implementations of two hybrid quantum-classical dynamics methods, namely, Ehrenfest dynamics (EHR) and the Ab-Initio Multiple Cloning sampling technique for Multiconfigurational Ehrenfest quantum dynamics (MCE-AIMC or simply AIMC), which are alternative options to the previously implemented trajectory surface hopping (TSH) method. To illustrate these methodologies, we outline a direct comparison of these three hybrid quantum-classical dynamics methods as implemented in the same NEXMD framework, discussing their weaknesses and strengths, using the modeled photodynamics of a polyphenylene ethylene dendrimer building block as a representative example. We also describe the expanded normal-mode analysis and constraints for both the ground and excited states, newly implemented in the NEXMD v2.0 framework, which allow for a deeper analysis of the main vibrational motions involved in vibronic dynamics. Overall, NEXMD v2.0 expands the range of applications of NEXMD to a larger variety of multichromophore organic molecules and photophysical processes involving quantum coherences and persistent couplings between electronic excited states and nuclear velocity.
Collapse
Affiliation(s)
- Victor M Freixas
- Departments of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697-2025, United States
| | - Walter Malone
- Department of Physics, Tuskegee University, Tuskegee, Alabama 36088, United States
| | - Xinyang Li
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Huajing Song
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Hassiel Negrin-Yuvero
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Royle Pérez-Castillo
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET, B1876BXD Bernal, Argentina
| | - Alexander White
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tammie R Gibson
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Dmitry V Makhov
- School of Chemistry, University of Leeds, Leeds LS2 9JT, United Kingdom
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | | | - Yu Zhang
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Nikita Fedik
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Maksim Kulichenko
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Richard Messerly
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Luke Nambi Mohanam
- Department of Electrical and Computer Engineering, College of Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Sahar Sharifzadeh
- Department of Electrical and Computer Engineering, College of Engineering, Boston University, Boston, Massachusetts 02215, United States
| | - Adolfo Bastida
- Departamento de Química Física, Universidad de Murcia, Murcia 30100, Spain
| | - Shaul Mukamel
- Departments of Chemistry and Physics and Astronomy, University of California, Irvine, California 92697-2025, United States
| | | | - Sergei Tretiak
- Theoretical Division, Center for Nonlinear Studies (CNLS), and Center for Integrated Nanotechnologies (CINT), Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
10
|
Zhang J, Peng J, Zhu Y, Hu D, Lan Z. Influence of Mode-Specific Excitation on the Nonadiabatic Dynamics of Methyl Nitrate (CH 3ONO 2). J Phys Chem Lett 2023:6542-6549. [PMID: 37450883 DOI: 10.1021/acs.jpclett.3c00664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2023]
Abstract
The impact of mode-specific vibrational excitations on initial-preparation conditions was studied by examining the excited-state population decay rates in the nonadiabatic dynamics of methyl nitrate (CH3ONO2). In particular, exciting a few specific modes by adding a single quantum of energy clearly decelerated the nonadiabatic dynamics population decay rates. The underlying reason for this slower population decay was explained by analyzing the profiles of the excited-state potential energy surfaces in the Franck-Condon regions and the topology of the S1/S0 conical intersection. This study not only provides physical insights into the key mechanisms controlling nonadiabatic dynamics but also shows the possibility of controlling nonadiabatic dynamics via mode-specific vibrational excitations.
Collapse
Affiliation(s)
- Juanjuan Zhang
- 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 Chemistry, South China Normal University, Guangzhou 510006, China
| | - Yifei Zhu
- 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
| | - Deping Hu
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, 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
|
11
|
Hu D, Huo P. Ab Initio Molecular Cavity Quantum Electrodynamics Simulations Using Machine Learning Models. J Chem Theory Comput 2023; 19:2353-2368. [PMID: 37000936 PMCID: PMC10134431 DOI: 10.1021/acs.jctc.3c00137] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
Abstract
We present a mixed quantum-classical simulation of polariton dynamics for molecule-cavity hybrid systems. In particular, we treat the coupled electronic-photonic degrees of freedom (DOFs) as the quantum subsystem and the nuclear DOFs as the classical subsystem and use the trajectory surface hopping approach to simulate non-adiabatic dynamics among the polariton states due to the coupled motion of nuclei. We use the accurate nuclear gradient expression derived from the Pauli-Fierz quantum electrodynamics Hamiltonian without making further approximations. The energies, gradients, and derivative couplings of the molecular systems are obtained from the on-the-fly simulations at the level of complete active space self-consistent field (CASSCF), which are used to compute the polariton energies and nuclear gradients. The derivatives of dipoles are also necessary ingredients in the polariton nuclear gradient expression but are often not readily available in electronic structure methods. To address this challenge, we use a machine learning model with the Kernel ridge regression method to construct the dipoles and further obtain their derivatives, at the same level as the CASSCF theory. The cavity loss process is modeled with the Lindblad jump superoperator on the reduced density of the electronic-photonic quantum subsystem. We investigate the azomethane molecule and its photoinduced isomerization dynamics inside the cavity. Our results show the accuracy of the machine-learned dipoles and their usage in simulating polariton dynamics. Our polariton dynamics results also demonstrate the isomerization reaction of azomethane can be effectively tuned by coupling to an optical cavity and by changing the light-matter coupling strength and the cavity loss rate.
Collapse
|
12
|
Taddei M, Garavelli M, Amirjalayer S, Conti I, Nenov A. Modus Operandi of a Pedalo-Type Molecular Switch: Insight from Dynamics and Theoretical Spectroscopy. MOLECULES (BASEL, SWITZERLAND) 2023; 28:molecules28020816. [PMID: 36677872 PMCID: PMC9863296 DOI: 10.3390/molecules28020816] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Revised: 01/09/2023] [Accepted: 01/11/2023] [Indexed: 01/15/2023]
Abstract
Molecular switches which can be triggered by light to interconvert between two or more well-defined conformation differing in their chemical or physical properties are fundamental for the development of materials with on-demand functionalities. Recently, a novel molecular switch based on a the azodicarboxamide core has been reported. It exhibits a volume-conserving conformational change upon excitation, making it a promising candidate for embedding in confined environments. In order to rationally implement and efficiently utilize the azodicarboxamide molecular switch, detailed insight into the coordinates governing the excited-state dynamics is needed. Here, we report a detailed comparative picture of the molecular motion at the atomic level in the presence and absence of explicit solvent. Our hybrid quantum mechanics/molecular mechanics (QM/MM) excited state simulations reveal that, although the energy landscape is slightly modulated by the solvation, the light-induced motion is dominated by a bending-assisted pedalo-type motion independent of the solvation. To support the predicted mechanism, we simulate time-resolved IR spectroscopy from first principles, thereby resolving fingerprints of the light-induced switching process. Our calculated time-resolved data are in good agreement with previously reported measured spectra.
Collapse
Affiliation(s)
- Mario Taddei
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, 40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, 40136 Bologna, Italy
- Correspondence: (M.G.); (I.C.); (A.N.)
| | - Saeed Amirjalayer
- Center for Nanotechnology, Center for Multiscale Theory and Computation, Physikalisches Institut, Westfälische Wilhelms-Universität Münster, 48149 Münster, Germany
| | - Irene Conti
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, 40136 Bologna, Italy
- Correspondence: (M.G.); (I.C.); (A.N.)
| | - Artur Nenov
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, 40136 Bologna, Italy
- Correspondence: (M.G.); (I.C.); (A.N.)
| |
Collapse
|
13
|
Jain A, Sindhu A. Pedagogical Overview of the Fewest Switches Surface Hopping Method. ACS OMEGA 2022; 7:45810-45824. [PMID: 36570264 PMCID: PMC9773185 DOI: 10.1021/acsomega.2c04843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/31/2022] [Accepted: 11/18/2022] [Indexed: 06/17/2023]
Abstract
The fewest switches surface hopping method continues to grow in popularity to capture electronic nonadiabaticity and quantum nuclear effects due to its simplicity and accuracy. Knowing the basics of the method is essential for the correct implementation and interpretation of results. This review covers the fundamentals of the fewest switches surface hopping method with a detailed discussion of the nuances such as decoherence schemes and frustrated hops and the correct approach to calculating populations. The consequences of incorrect implementation are further discussed toward calculating kinetic and thermodynamic properties. Some tips for practitioners and a step-by-step algorithm for developers are provided. Finally, some of the finer technicalities of the fewest switches surface hopping method that are buried deep in the literature are pointed out to help graduate students better appreciate this method.
Collapse
|
14
|
Fang Y, Huang H, Lin K, Xu C, Gu FL, Lan Z. The impact of different geometrical restrictions on the nonadiabatic photoisomerization of biliverdin chromophores. Phys Chem Chem Phys 2022; 24:26190-26199. [PMID: 36278817 DOI: 10.1039/d2cp02941c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The photoisomerization mechanism of the chromophore of bacterial biliverdin (BV) phytochromes is explored via nonadiabatic dynamics simulation by using the on-the-fly trajectory surface-hopping method at the semi-empirical OM2/MRCI level. Particularly, the current study focuses on the influence of geometrical constrains on the nonadiabatic photoisomerization dynamics of the BV chromophore. Here a rather simplified approach is employed in the nonadiabatic dynamics to capture the features of geometrical constrains, which adds mechanical restrictions to the specific moieties of the BV chromophore. This simplified method provides a rather quick approach to examine the influence of geometrical restrictions on photoisomerization. As expected, different constrains bring distinctive influences on the photoisomerization mechanism of the BV chromophore, giving either strong or minor modification of both involved reaction channels and excited-state lifetimes after the constrains are added in different ring moieties. These observations not only contribute to the primary understanding of the role of the spatial restriction caused by biological environments in photoinduced dynamics of the BV chromophore, but also provide useful ideas for the artificial regulation of the photoisomerization reaction channels of phytochrome proteins.
Collapse
Affiliation(s)
- Yuan Fang
- School of Chemistry, South China Normal University, Guangzhou 510006, China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
| | - Haiyi Huang
- School of Chemistry, South China Normal University, Guangzhou 510006, China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
| | - Kunni Lin
- School of Chemistry, South China Normal University, Guangzhou 510006, China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
15
|
Liu J, Hu W, Yang J. Accelerating Linear-Response Time-Dependent Hybrid Density Functional Theory with Low-Rank Decomposition Techniques in the Plane-Wave Basis. J Chem Theory Comput 2022; 18:6713-6721. [PMID: 36242561 DOI: 10.1021/acs.jctc.2c00763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present an efficient low-rank implementation of linear-response time-dependent density functional theory for hybrid functionals (hybrid-LR-TDDFT) within the plane-wave pseudopotential framework. The adaptively compressed exchange (ACE) operator and the natural transition orbitals (NTOs) are introduced to build the low-rank representation of the nonlocal exchange operator in the hybrid-LR-TDDFT Hamiltonian. Numerical tests demonstrate that the ACE approximation significantly reduces the computational cost of applying the nonlocal exchange operator without loss of accuracy, and the NTO approximation can further accelerate the hybrid-LR-TDDFT calculations by introducing an NTO cutoff parameter. This new method enables us to effectively study the excitonic properties of two-dimensional MoS2 consisting of 216 atoms and ∼1900 electrons with range-separated hybrid functionals on a single graphics processing unit.
Collapse
Affiliation(s)
- Jie Liu
- Hefei National Laboratory, University of Science and Technology of China, Hefei230088, China
| | - Wei Hu
- Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui230026, China
| | - Jinlong Yang
- Hefei National Laboratory, University of Science and Technology of China, Hefei230088, China.,Hefei National Research Center for Physical Sciences at the Microscale, University of Science and Technology of China, Hefei, Anhui230026, China
| |
Collapse
|
16
|
Barbatti M, Bondanza M, Crespo-Otero R, Demoulin B, Dral PO, Granucci G, Kossoski F, Lischka H, Mennucci B, Mukherjee S, Pederzoli M, Persico M, Pinheiro M, Pittner J, Plasser F, Sangiogo Gil E, Stojanovic L. Newton-X Platform: New Software Developments for Surface Hopping and Nuclear Ensembles. J Chem Theory Comput 2022; 18:6851-6865. [PMID: 36194696 PMCID: PMC9648185 DOI: 10.1021/acs.jctc.2c00804] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Abstract
![]()
Newton-X is an open-source computational platform to
perform nonadiabatic
molecular dynamics based on surface hopping and spectrum simulations
using the nuclear ensemble approach. Both are among the most common
methodologies in computational chemistry for photophysical and photochemical
investigations. This paper describes the main features of these methods
and how they are implemented in Newton-X. It emphasizes the newest
developments, including zero-point-energy leakage correction, dynamics
on complex-valued potential energy surfaces, dynamics induced by incoherent
light, dynamics based on machine-learning potentials, exciton dynamics
of multiple chromophores, and supervised and unsupervised machine
learning techniques. Newton-X is interfaced with several third-party
quantum-chemistry programs, spanning a broad spectrum of electronic
structure methods.
Collapse
Affiliation(s)
- Mario Barbatti
- Aix Marseille University, CNRS, ICR, 13013Marseille, France.,Institut Universitaire de France, 75231Paris, France
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | - Rachel Crespo-Otero
- Department of Chemistry, Queen Mary University of London, Mile End Road, E1 4NSLondon, U.K
| | | | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, 361005Xiamen, China
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000Toulouse, France
| | - Hans Lischka
- Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas79409, United States
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | | | - Marek Pederzoli
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223Prague 8, Czech Republic
| | - Maurizio Persico
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | - Max Pinheiro
- Aix Marseille University, CNRS, ICR, 13013Marseille, France
| | - Jiří Pittner
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223Prague 8, Czech Republic
| | - Felix Plasser
- Department of Chemistry, Loughborough University, LE11 3TULoughborough, U.K
| | - Eduarda Sangiogo Gil
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via Moruzzi 13, 56124Pisa, Italy
| | - Ljiljana Stojanovic
- Department of Physics and Astronomy, University College London, Gower Street, WC1E 6BTLondon, U.K
| |
Collapse
|
17
|
Zhu Y, Peng J, Kang X, Xu C, Lan Z. The principal component analysis of the ring deformation in the nonadiabatic surface hopping dynamics. Phys Chem Chem Phys 2022; 24:24362-24382. [PMID: 36178471 DOI: 10.1039/d2cp03323b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The analysis of the leading active molecular motions in the on-the-fly trajectory surface hopping simulation provides the essential information to understand the geometric evolution in nonadiabatic dynamics. When the ring deformation is involved, the identification of the key active coordinates becomes challenging. A "hierarchical" protocol based on the dimensionality reduction and clustering approaches is proposed for the automatic analysis of the ring deformation in the nonadiabatic molecular dynamics. The representative system keto isocytosine is taken as the prototype to illustrate this protocol. The results indicate that the current hierarchical analysis protocol is a powerful way to clearly clarify both the major and minor active molecular motions of the ring distortion in nonadiabatic dynamics.
Collapse
Affiliation(s)
- Yifei Zhu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China. .,MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.,School of Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Xu Kang
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China. .,MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China. .,MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China. .,MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
18
|
Chakraborty P, Liu Y, McClung S, Weinacht T, Matsika S. Nonadiabatic Excited State Dynamics of Organic Chromophores: Take-Home Messages. J Phys Chem A 2022; 126:6021-6031. [PMID: 36069531 DOI: 10.1021/acs.jpca.2c04671] [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
Nonadiabatic excited state dynamics are important in a variety of processes. Theoretical and experimental developments have allowed for a great progress in this area, while combining the two is often necessary and the best approach to obtain insight into the photophysical behavior of molecules. In this Feature Article we use examples of our recent work combining time-resolved photoelectron spectroscopy with theoretical nonadiabatic dynamics to highlight important lessons we learned. We compare the nonadiabatic excited state dynamics of three different organic molecules with the aim of elucidating connections between structure and dynamics. Calculations and measurements are compared for uracil, 1,3-cyclooctadiene, and 1,3-cyclohexadiene. The comparison highlights the role of rigidity in influencing the dynamics and the difficulty of capturing the dynamics accurately with calculations.
Collapse
Affiliation(s)
- Pratip Chakraborty
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States.,Division of Theoretical Chemistry and Biology, Department of Chemistry, School of Engineering Sciences in Chemistry, Biotechnology and Health (CBH), KTH Royal Institute of Technology, SE-10044 Stockholm, Sweden
| | - Yusong Liu
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, United States.,Stanford PULSE Institute, SLAC National Laboratory, Menlo Park, California 94025, United States
| | - Samuel McClung
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, United States
| | - Thomas Weinacht
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, United States
| | - Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
19
|
Significant impact of deprotonated status on the photoisomerization dynamics of bacteriophytochrome chromophore. CHINESE CHEM LETT 2022. [DOI: 10.1016/j.cclet.2022.107850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
20
|
T. do Casal M, Toldo JM, Pinheiro Jr M, Barbatti M. Fewest switches surface hopping with Baeck-An couplings. OPEN RESEARCH EUROPE 2022; 1:49. [PMID: 37645211 PMCID: PMC10446015 DOI: 10.12688/openreseurope.13624.2] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/01/2022] [Indexed: 08/31/2023]
Abstract
In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.
Collapse
|
21
|
T. do Casal M, Toldo JM, Pinheiro Jr M, Barbatti M. Fewest switches surface hopping with Baeck-An couplings. OPEN RESEARCH EUROPE 2022; 1:49. [PMID: 37645211 PMCID: PMC10446015 DOI: 10.12688/openreseurope.13624.1] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 03/01/2022] [Indexed: 03/12/2024]
Abstract
In the Baeck-An (BA) approximation, first-order nonadiabatic coupling vectors are given in terms of adiabatic energy gaps and the second derivative of the gaps with respect to the coupling coordinate. In this paper, a time-dependent (TD) BA approximation is derived, where the couplings are computed from the energy gaps and their second time-derivatives. TD-BA couplings can be directly used in fewest switches surface hopping, enabling nonadiabatic dynamics with any electronic structure methods able to provide excitation energies and energy gradients. Test results of surface hopping with TD-BA couplings for ethylene and fulvene show that the TD-BA approximation delivers a qualitatively correct picture of the dynamics and a semiquantitative agreement with reference data computed with exact couplings. Nevertheless, TD-BA does not perform well in situations conjugating strong couplings and small velocities. Considered the uncertainties in the method, TD-BA couplings could be a competitive approach for inexpensive, exploratory dynamics with a small trajectories ensemble. We also assessed the potential use of TD-BA couplings for surface hopping dynamics with time-dependent density functional theory (TDDFT), but the results are not encouraging due to singlet instabilities near the crossing seam with the ground state.
Collapse
|
22
|
Xu C, Lin K, Hu D, Gu FL, Gelin MF, Lan Z. Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods. J Phys Chem Lett 2022; 13:661-668. [PMID: 35023755 DOI: 10.1021/acs.jpclett.1c03373] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.
Collapse
Affiliation(s)
- Chao Xu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Chemistry, South China Normal University, Guangzhou, 510006, P. R. China
| | - Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Chemistry, South China Normal University, Guangzhou, 510006, P. R. 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, P. R. China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education; School of Chemistry, South China Normal University, Guangzhou, 510006, P. R. China
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou, 310018, P. R. 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
| |
Collapse
|
23
|
Understanding the Photolysis of CH 3ONO 2 with the On-the-fly Nonadiabatic Dynamics Simulation at the ADC(2) Level. CHINESE J CHEM PHYS 2022. [DOI: 10.1063/1674-0068/cjcp2201013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
|
24
|
Borrego-Varillas R, Nenov A, Kabaciński P, Conti I, Ganzer L, Oriana A, Jaiswal VK, Delfino I, Weingart O, Manzoni C, Rivalta I, Garavelli M, Cerullo G. Tracking excited state decay mechanisms of pyrimidine nucleosides in real time. Nat Commun 2021; 12:7285. [PMID: 34907186 PMCID: PMC8671501 DOI: 10.1038/s41467-021-27535-7] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Accepted: 11/24/2021] [Indexed: 11/18/2022] Open
Abstract
DNA owes its remarkable photostability to its building blocks-the nucleosides-that efficiently dissipate the energy acquired upon ultraviolet light absorption. The mechanism occurring on a sub-picosecond time scale has been a matter of intense debate. Here we combine sub-30-fs transient absorption spectroscopy experiments with broad spectral coverage and state-of-the-art mixed quantum-classical dynamics with spectral signal simulations to resolve the early steps of the deactivation mechanisms of uridine (Urd) and 5-methyluridine (5mUrd) in aqueous solution. We track the wave packet motion from the Franck-Condon region to the conical intersections (CIs) with the ground state and observe spectral signatures of excited-state vibrational modes. 5mUrd exhibits an order of magnitude longer lifetime with respect to Urd due to the solvent reorganization needed to facilitate bulky methyl group motions leading to the CI. This activates potentially lesion-inducing dynamics such as ring opening. Involvement of the 1nπ* state is found to be negligible.
Collapse
Affiliation(s)
| | - Artur Nenov
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, Bologna, Italy
| | - Piotr Kabaciński
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - Irene Conti
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, Bologna, Italy
| | - Lucia Ganzer
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - Aurelio Oriana
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy
| | - Vishal Kumar Jaiswal
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, Bologna, Italy
| | - Ines Delfino
- Dipartimento di Scienze Ecologiche e Biologiche, Università della Tuscia, Via San Camillo de Lellis, snc, Viterbo, Italy
| | - Oliver Weingart
- Institut für Theoretische Chemie und Computerchemie, Heinrich Heine Universität Düsseldorf, Universitätsstrasse 1, Düsseldorf, Germany
| | | | - Ivan Rivalta
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, Bologna, Italy
- Université de Lyon, École Normale Supérieure de Lyon, CNRS UMR 5182, Laboratoire de Chimie, 46 allée d'Italie, F69364, Lyon, France
| | - Marco Garavelli
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, Bologna, Italy.
| | - Giulio Cerullo
- IFN-CNR, Piazza Leonardo da Vinci 32, Milano, Italy.
- Dipartimento di Fisica, Politecnico di Milano, Piazza Leonardo da Vinci 32, Milano, Italy.
| |
Collapse
|
25
|
Chen WK, Sun XW, Fang Q, Liu XY, Cui GL. GW/BSE nonadiabatic dynamics simulations on excited-state relaxation processes of zinc phthalocyanine-fullerene dyads: Roles of bridging chemical bonds. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2109162] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Wen-kai Chen
- 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
| | - Qiu Fang
- 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
| | - Gang-long Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
26
|
Gómez S, Soysal EN, Worth GA. Micro-Solvated DMABN: Excited State Quantum Dynamics and Dual Fluorescence Spectra. Molecules 2021; 26:7247. [PMID: 34885829 PMCID: PMC8658867 DOI: 10.3390/molecules26237247] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Revised: 11/23/2021] [Accepted: 11/24/2021] [Indexed: 12/04/2022] Open
Abstract
In this work, we report a complete analysis by theoretical and spectroscopic methods of the short-time behaviour of 4-(dimethylamino)benzonitrile (DMABN) in the gas phase as well as in cyclohexane, tetrahydrofuran, acetonitrile, and water solution, after excitation to the La state. The spectroscopic properties of DMABN were investigated experimentally using UV absorption and fluorescence emission spectroscopy. The computational study was developed at different electronic structure levels and using the Polarisable Continuum Model (PCM) and explicit solvent molecules to reproduce the solvent environment. Additionally, excited state quantum dynamics simulations in the diabatic picture using the direct dynamics variational multiconfigurational Gaussian (DD-vMCG) method were performed, the largest quantum dynamics "on-the-fly" simulations performed with this method until now. The comparison with fully converged multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) dynamics on parametrised linear vibronic coupling (LVC) potentials show very similar population decays and evolution of the nuclear wavepacket. The ring C=C stretching and three methyl tilting modes are identified as the responsible motions for the internal conversion from the La to the Lb states. No major differences are observed in the ultrafast initial decay in different solvents, but we show that this effect depends strongly on the level of electronic structure used.
Collapse
Affiliation(s)
- Sandra Gómez
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK;
| | | | - Graham A. Worth
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK;
| |
Collapse
|
27
|
Lin K, Peng J, Gu FL, Lan Z. Simulation of Open Quantum Dynamics with Bootstrap-Based Long Short-Term Memory Recurrent Neural Network. J Phys Chem Lett 2021; 12:10225-10234. [PMID: 34647736 DOI: 10.1021/acs.jpclett.1c02672] [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/13/2023]
Abstract
The recurrent neural network with the long short-term memory cell (LSTM-NN) is employed to simulate the long-time dynamics of open quantum systems. The bootstrap method is applied in the LSTM-NN construction and prediction, which provides a Monte Carlo estimation of a forecasting confidence interval. Within this approach, a large number of LSTM-NNs are constructed by resampling time-series sequences that were obtained from the early stage quantum evolution given by numerically exact multilayer multiconfigurational time-dependent Hartree method. The built LSTM-NN ensemble is used for the reliable propagation of the long-time quantum dynamics, and the simulated result is highly consistent with the exact evolution. The forecasting uncertainty that partially reflects the reliability of the LSTM-NN prediction is also given. This demonstrates the bootstrap-based LSTM-NN approach is a practical and powerful tool to propagate the long-time quantum dynamics of open systems with high accuracy and low computational cost.
Collapse
Affiliation(s)
- Kunni Lin
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. 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, P. R. China
| | - Feng Long Gu
- Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry, South China Normal University, Guangzhou 510006, P. R. 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
| |
Collapse
|
28
|
Hu D, Peng J, Chen L, Gelin MF, Lan Z. Spectral Fingerprint of Excited-State Energy Transfer in Dendrimers through Polarization-Sensitive Transient-Absorption Pump-Probe Signals: On-the-Fly Nonadiabatic Dynamics Simulations. J Phys Chem Lett 2021; 12:9710-9719. [PMID: 34590858 DOI: 10.1021/acs.jpclett.1c02640] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The time-resolved polarization-sensitive transient-absorption (TA) pump-probe (PP) spectra are simulated using on-the-fly surface-hopping nonadiabatic dynamics and the doorway-window representation of nonlinear spectroscopy. A dendrimer model system composed of two linear phenylene ethynylene units (2-ring and 3-ring) is taken as an example. The ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA) contributions as well as the total TA PP signals are obtained and carefully analyzed. It is shown that intramolecular excited-state energy transfer from the 2-ring unit to the 3-ring unit can be conveniently identified by employing pump and probe pulses with different polarizations. Our results demonstrate that time-resolved polarization-sensitive TA PP signals provide a powerful tool for the elucidation of excited-state energy-transfer pathways, notably in molecular systems possessing several optically bright nonadiabatically coupled electronic states with different orientations of transition dipole moments.
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
| | - 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
| | - Lipeng Chen
- Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, 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
|
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
|
Nag P, Anand N, Vennapusa SR. Ultrafast nonadiabatic excited-state intramolecular proton transfer in 3-hydroxychromone: A surface hopping approach. J Chem Phys 2021; 155:094301. [PMID: 34496583 DOI: 10.1063/5.0060934] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We employ the ab initio molecular dynamics within the surface hopping method to explore the excited-state intramolecular proton transfer taking place on the coupled "bright" S1 (ππ*) and "dark" S2 (nπ*) states of 3-hydroxychromone. The nonadiabatic population transfer between these states via an accessible conical intersection would open up multiple proton transfer pathways. Our findings reveal the keto tautomer formation via S1 on a timescale similar to the O-H in-plane vibrational period (<100 fs). Structural analysis indicates that a few parameters of the five-membered proton transfer geometry that constitute the donor (hydroxyl) and acceptor (carbonyl) groups would be adequate to drive the enol to keto transformation. We also investigate the role of O-H in-plane and out-of-plane vibrational motions in the excited-state dynamics of 3-hydroxychromone.
Collapse
Affiliation(s)
- Probal Nag
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram 695551, India
| | - Neethu Anand
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram 695551, India
| | - Sivaranjana Reddy Vennapusa
- School of Chemistry, Indian Institute of Science Education and Research Thiruvananthapuram, Maruthamala PO, Vithura, Thiruvananthapuram 695551, India
| |
Collapse
|
31
|
Jankowska J, Sobolewski AL. Modern Theoretical Approaches to Modeling the Excited-State Intramolecular Proton Transfer: An Overview. Molecules 2021; 26:molecules26175140. [PMID: 34500574 PMCID: PMC8434569 DOI: 10.3390/molecules26175140] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2021] [Revised: 08/20/2021] [Accepted: 08/23/2021] [Indexed: 02/02/2023] Open
Abstract
The excited-state intramolecular proton transfer (ESIPT) phenomenon is nowadays widely acknowledged to play a crucial role in many photobiological and photochemical processes. It is an extremely fast transformation, often taking place at sub-100 fs timescales. While its experimental characterization can be highly challenging, a rich manifold of theoretical approaches at different levels is nowadays available to support and guide experimental investigations. In this perspective, we summarize the state-of-the-art quantum-chemical methods, as well as molecular- and quantum-dynamics tools successfully applied in ESIPT process studies, focusing on a critical comparison of their specific properties.
Collapse
Affiliation(s)
- Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
- Correspondence:
| | | |
Collapse
|
32
|
Matsika S. Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections. Chem Rev 2021; 121:9407-9449. [PMID: 34156838 DOI: 10.1021/acs.chemrev.1c00074] [Citation(s) in RCA: 50] [Impact Index Per Article: 16.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.
Collapse
Affiliation(s)
- Spiridoula Matsika
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
33
|
Wang Y, Simine L. Solving the Wigner equation with signed particle Monte Carlo for chemically relevant potentials. J Chem Phys 2021; 155:034109. [PMID: 34293875 DOI: 10.1063/5.0055603] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Expanding the set of stable, accurate, and scalable methods for simulating molecular quantum dynamics is important for accelerating the computational exploration of molecular processes. In this paper, we adapt the signed particles Monte Carlo algorithm for solving the transient Wigner equation to scenarios of chemical interest. This approach was used in the past to study electronic processes in semi-conductors, but to the best of our knowledge, it had never been applied to molecular modeling. We present the algorithm and demonstrate its excellent performance on harmonic and double well potentials for electronic and nuclear systems. We explore the stability of the algorithm, discuss the choice of hyper-parameters, and cautiously speculate that it may be used in quantum molecular dynamics simulations.
Collapse
Affiliation(s)
- Yu Wang
- Department of Chemistry, McGill University, 801 Sherbrooke Street W, Montreal, Quebec H3A 0B8, Canada
| | - Lena Simine
- Department of Chemistry, McGill University, 801 Sherbrooke Street W, Montreal, Quebec H3A 0B8, Canada
| |
Collapse
|
34
|
Lee IS, Ha JK, Han D, Kim TI, Moon SW, Min SK. PyUNIxMD: A Python-based excited state molecular dynamics package. J Comput Chem 2021; 42:1755-1766. [PMID: 34197646 PMCID: PMC8362049 DOI: 10.1002/jcc.26711] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 06/10/2021] [Accepted: 06/16/2021] [Indexed: 01/17/2023]
Abstract
Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.
Collapse
Affiliation(s)
- In Seong Lee
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Jong-Kwon Ha
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Daeho Han
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Tae In Kim
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Sung Wook Moon
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| | - Seung Kyu Min
- Department of Chemistry, Ulsan National Institute of Science and Technology (UNIST), Ulsan, South Korea
| |
Collapse
|
35
|
Li W, She Y, Vasenko AS, Prezhdo OV. Ab initio nonadiabatic molecular dynamics of charge carriers in metal halide perovskites. NANOSCALE 2021; 13:10239-10265. [PMID: 34031683 DOI: 10.1039/d1nr01990b] [Citation(s) in RCA: 51] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron-phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge-charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.
Collapse
Affiliation(s)
- Wei Li
- School of Chemistry and Materials Science, Hunan Agricultural University, Changsha 410128, People's Republic of China.
| | | | | | | |
Collapse
|
36
|
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
|
37
|
Chu W, Prezhdo OV. Concentric Approximation for Fast and Accurate Numerical Evaluation of Nonadiabatic Coupling with Projector Augmented-Wave Pseudopotentials. J Phys Chem Lett 2021; 12:3082-3089. [PMID: 33750138 DOI: 10.1021/acs.jpclett.0c03853] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We develop an efficient and accurate method for numerical evaluation of nonadiabatic (NA) coupling in the Kohn-Sham representation with projector augmented-wave (PAW) pseudopotentials that are commonly used in electronic structure calculations on nanoscale, condensed matter, and molecular systems. Without additional cost, the method provides an order of magnitude improvement in accuracy compared to the current technique, while it is 3-4 orders of magnitude faster than the exact evaluation. Atomic displacements over typical time steps in molecular dynamics (MD) simulations are much smaller than the size of the PAW core region, and therefore, evaluation of the NA in the core is simplified. The accuracy is demonstrated with three condensed matter systems. The method is robust to variation in the MD time step. The accurate NA coupling evaluation also helps in maintaining phase-consistency of the NA coupling and identifying trivial crossings of adiabatic states. The approach stimulates NAMD applications to modeling of modern materials and processes.
Collapse
Affiliation(s)
- Weibin Chu
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| | - Oleg V Prezhdo
- Departments of Chemistry, and Physics and Astronomy, University of Southern California, Los Angeles, California 90089, United States
| |
Collapse
|
38
|
Gelin MF, Huang X, Xie W, Chen L, Došlić NA, Domcke W. Ab Initio Surface-Hopping Simulation of Femtosecond Transient-Absorption Pump-Probe Signals of Nonadiabatic Excited-State Dynamics Using the Doorway-Window Representation. J Chem Theory Comput 2021; 17:2394-2408. [PMID: 33755464 DOI: 10.1021/acs.jctc.1c00109] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
An ab initio theoretical framework for the simulation of femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectra with quasi-classical trajectories is presented. The simulations are based on the classical approximation to the doorway-window (DW) representation of third-order four-wave-mixing signals. The DW formula accounts for the finite duration and spectral shape of the pump and probe pulses. In the classical DW formalism, classical trajectories are stochastically sampled from a positive definite doorway distribution, and the signals are evaluated by averaging over a positive definite window distribution. Nonadiabatic excited-state dynamics is described by a stochastic surface-hopping algorithm. The method has been implemented for the pyrazine molecule with the second-order algebraic-diagrammatic construction (ADC(2)) ab initio electronic-structure method. The methodology is illustrated by ab initio simulations of the ground-state bleach, stimulated emission, and excited-state absorption contributions to the TA PP spectrum of gas-phase pyrazine.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Xiang Huang
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| | - Weiwei Xie
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Lipeng Chen
- Max Planck Institute for the Physics of Complex Systems, D-01187 Dresden, Germany
| | - Nad A Došlić
- Department of Physical Chemistry, Ruder Boscovic Institute, HR-10000 Zagreb, Croatia
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| |
Collapse
|
39
|
Liu J, Hu W, Yang J. Two-level iterative solver for linear response time-dependent density functional theory with plane wave basis set. J Chem Phys 2021; 154:064101. [PMID: 33588554 DOI: 10.1063/5.0032464] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a two-level iterative solver for linear response time-dependent density functional theory by combining two forms of the Casida equation in the Kohn-Sham orbital representation and in Hutter's formulation. This two-level iterative solver has been implemented with the plane wave pseudopotential method for excited-state simulations of molecular and low-dimensional solid materials. Numerical studies with the Davidson algorithm demonstrate that this two-level iterative solver yields excited-state properties for molecules (benzene C6H6 and fullerene C60) and low-dimensional semiconductors [two-dimensional molybdenum disulfide MoS2 monolayer and rutile titanium dioxide TiO2(110) surface] with significantly reduced computational cost and storage requirement compared with standard iterative algorithms. We apply our approach to investigate the photoinduced charge separation of methanol molecules adsorption on the rutile TiO2(110) surface from the exciton perspective and validate that the photogenerated hole can be captured by methanol molecules.
Collapse
Affiliation(s)
- Jie Liu
- Hefei National Laboratory for Physical Sciences at the Microscale, Department of Chemical Physics, and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Wei Hu
- Hefei National Laboratory for Physical Sciences at the Microscale, Department of Chemical Physics, and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at the Microscale, Department of Chemical Physics, and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
40
|
Zhang Q, Li Y, Cao Z, Zhu C. Aggregation-induced emission spectra of triphenylamine salicylaldehyde derivatives via excited-state intramolecular proton transfer revealed by molecular spectral and dynamics simulations. RSC Adv 2021; 11:37171-37180. [PMID: 35496419 PMCID: PMC9043594 DOI: 10.1039/d1ra07388e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Accepted: 11/03/2021] [Indexed: 11/21/2022] Open
Abstract
Aggregation-induced emission (AIE) spectra accompanied by excited state intramolecular proton transfer (ESIPT) for two triphenylamine salicylaldehyde derivatives (namely, TS and TS-OMe) are investigated by performing molecular spectral and dynamics simulations associated with the hybrid quantum mechanics/molecular mechanics (QM/MM) at the quantum level of the time-dependent density functional theory. The simulated emission spectral peaks and Stokes' shifts are in good agreement with the experimental results for both TS and TS-OMe. Furthermore, the AIE spectral mechanisms are well explained to be associated with the ESIPT processes for both TS and TS-OMe monomers in the aggregated crystal state, while the AIE spectra mechanism for the TS-OMe (TS) dimer is accompanied by intermolecular charge-transfer excitation process. Besides, the TS dimers also contributed to the AIE mechanisms in the crystal with the intermolecular charge-transfer from one monomer to another. In addition, the TS dimers are contributed to the AIE mechanisms in the crystal with the intermolecular charge-transfer from one monomer to another. On the other hand, simulated emission spectra for both the TS and TS-OMe monomers in acetonitrile solution are involved in mixed emission with and without the ESIPT process, as interpreted by nonadiabatic molecular dynamics simulation. It is also briefly addressed that the emission spectra in the solution are weak and enhanced in the crystal. The present study provides a great physical insight into the design of highly efficient AIE compounds. The emission spectra of triphenylamine salicylaldehyde derivatives are weak in solution and strong in the crystal, as interpreted by molecular spectral and dynamics simulations.![]()
Collapse
Affiliation(s)
- Qing Zhang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, China
| | - Yuanyuan Li
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, China
| | - Zexing Cao
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 360015, China
| | - Chaoyuan Zhu
- Department of Applied Chemistry, Institute of Molecular Science, National Chiao-Tung University, Hsinchu 30010, Taiwan
- Department of Applied Chemistry, Center for Emergent Functional Matter Science, National Yang Ming Chiao Tung University, Hsinchu 30010, Taiwan
| |
Collapse
|
41
|
Kochman MA, Durbeej B. Simulating the Nonadiabatic Relaxation Dynamics of 4-( N, N-Dimethylamino)benzonitrile (DMABN) in Polar Solution. J Phys Chem A 2020; 124:2193-2206. [PMID: 32083861 DOI: 10.1021/acs.jpca.9b10588] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The compound 4-(N,N-dimethylamino)benzonitrile (DMABN) represents the archetypal system for dual fluorescence, a rare photophysical phenomenon in which a given fluorophore shows two distinct emission bands. Despite extensive studies, the underlying mechanism remains the subject of debate. In the present contribution, we address this issue by simulating the excited-state relaxation process of DMABN as it occurs in polar solution. The potential energy surfaces for the system are constructed with the use of the additive quantum mechanics/molecular mechanics (QM/MM) method, and the coupled dynamics of the electronic wave function and the nuclei is propagated with the semiclassical fewest switches surface hopping method. The DMABN molecule, which comprises the QM subsystem, is treated with the use of the second-order algebraic diagrammatic construction (ADC(2)) method with the imposition of spin-opposite scaling (SOS). It is verified that this level of theory achieves a realistic description of the excited-state potential energy surfaces of DMABN. The simulation results qualitatively reproduce the main features of the experimentally observed fluorescence spectrum, thus allowing the unambiguous assignment of the two fluorescence bands: the normal band is due to the near-planar locally excited (LE) structure of DMABN, while the so-called "anomalous" second band arises from the twisted intramolecular charge transfer (TICT) structure. The transformation of the LE structure into the TICT structure takes place directly via intramolecular rotation, and is not mediated by another excited-state structure. In particular, the oft-discussed rehybridized intramolecular charge transfer (RICT) structure, which is characterized by a bent nitrile group, does not play a role in the relaxation process.
Collapse
Affiliation(s)
- Michał Andrzej Kochman
- Division of Theoretical Chemistry, Department of Physics, Chemistry and Biology (IFM), Linköping University, 581 83 Linköping, Sweden.,Department of Chemistry, University College London (UCL), 20 Gordon Street, London WC1H 0A, United Kingdom
| | - Bo Durbeej
- Division of Theoretical Chemistry, Department of Physics, Chemistry and Biology (IFM), Linköping University, 581 83 Linköping, Sweden
| |
Collapse
|
42
|
Smith B, Akimov AV. Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:073001. [PMID: 31661681 DOI: 10.1088/1361-648x/ab5246] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States of America
| | | |
Collapse
|
43
|
Nelson TR, White AJ, Bjorgaard JA, Sifain AE, Zhang Y, Nebgen B, Fernandez-Alberti S, Mozyrsky D, Roitberg AE, Tretiak S. Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials. Chem Rev 2020; 120:2215-2287. [PMID: 32040312 DOI: 10.1021/acs.chemrev.9b00447] [Citation(s) in RCA: 198] [Impact Index Per Article: 49.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.
Collapse
Affiliation(s)
- Tammie R Nelson
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Alexander J White
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Josiah A Bjorgaard
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Andrew E Sifain
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States.,U.S. Army Research Laboratory , Aberdeen Proving Ground , Maryland 21005 , United States
| | - Yu Zhang
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Benjamin Nebgen
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | | | - Dmitry Mozyrsky
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Adrian E Roitberg
- Department of Chemistry , University of Florida , Gainesville , Florida 32611 , United States
| | - Sergei Tretiak
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| |
Collapse
|
44
|
Ibele LM, Curchod BFE. A molecular perspective on Tully models for nonadiabatic dynamics. Phys Chem Chem Phys 2020; 22:15183-15196. [DOI: 10.1039/d0cp01353f] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
We present a series of standardized molecular tests for nonadiabatic dynamics, reminiscent of the one-dimensional Tully models proposed in 1990.
Collapse
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry
- Durham University
- Durham DH1 3LE
- UK
| | | |
Collapse
|
45
|
Ibele LM, Nicolson A, Curchod BFE. Excited-state dynamics of molecules with classically driven trajectories and Gaussians. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1665199] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Lea M. Ibele
- Department of Chemistry, Durham University, Durham, UK
| | | | | |
Collapse
|
46
|
Wang L, Qiu J, Bai X, Xu J. Surface hopping methods for nonadiabatic dynamics in extended systems. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1435] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Linjun Wang
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Jing Qiu
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Xin Bai
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| | - Jiabo Xu
- Center for Chemistry of Novel & High‐Performance Materials, Department of Chemistry Zhejiang University Hangzhou China
| |
Collapse
|
47
|
Minezawa N, Nakajima T. Trajectory surface hopping molecular dynamics simulation by spin-flip time-dependent density functional theory. J Chem Phys 2019; 150:204120. [DOI: 10.1063/1.5096217] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Affiliation(s)
- Noriyuki Minezawa
- Computational Molecular Science Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Takahito Nakajima
- Computational Molecular Science Research Team, RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
48
|
Peng J, Xie Y, Hu D, Lan Z. Performance of trajectory surface hopping method in the treatment of ultrafast intersystem crossing dynamics. J Chem Phys 2019; 150:164126. [PMID: 31042919 DOI: 10.1063/1.5079426] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We carried out extensive studies to examine the performance of the fewest-switches surface hopping method in the description of the ultrafast intersystem crossing dynamic of various singlet-triplet (S-T) models by comparison with the results of the exact full quantum dynamics. Different implementation details and some derivative approaches were examined. As expected, it is better to perform the trajectory surface hopping calculations in the spin-adiabatic representation or by the local diabatization approach, instead of in the spin-diabatic representation. The surface hopping method provides reasonable results for the short-time dynamics in the S-T model with weak spin-orbital coupling (diabatic coupling), although it does not perform well in the models with strong spin-orbital coupling (diabatic coupling). When the system accesses the S-T potential energy crossing with rather high kinetic energy, the trajectory surface hopping method tends to produce a good description of the nonadiabatic intersystem crossing dynamics. The impact of the decoherence correction on the performance of the trajectory surface hopping is system dependent. It improves the result accuracy in many cases, while its influence may also be minor for other cases.
Collapse
Affiliation(s)
- Jiawei Peng
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, South China Normal University, Guangzhou 510006, China
| | - Yu Xie
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, South China Normal University, Guangzhou 510006, China
| | - Deping Hu
- MOE Key Laboratory of Organic Optoelectronics and Molecular Engineering, Department of Chemistry, Tsinghua University, Beijing 100084, 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, South China Normal University, Guangzhou 510006, China
| |
Collapse
|
49
|
Zheng Q, Chu W, Zhao C, Zhang L, Guo H, Wang Y, Jiang X, Zhao J. Ab initio nonadiabatic molecular dynamics investigations on the excited carriers in condensed matter systems. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1411] [Citation(s) in RCA: 87] [Impact Index Per Article: 17.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Qijing Zheng
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly‐Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics University of Science and Technology of China Hefei China
| | - Weibin Chu
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly‐Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics University of Science and Technology of China Hefei China
| | - Chuanyu Zhao
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly‐Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics University of Science and Technology of China Hefei China
| | - Lili Zhang
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly‐Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics University of Science and Technology of China Hefei China
| | - Hongli Guo
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly‐Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics University of Science and Technology of China Hefei China
- School of Physics and Technology, Center for Nanoscience and Nanotechnology, and Key Laboratory of Artificial Micro‐ and Nano‐structures of Ministry of Education Wuhan University Wuhan China
| | - Yanan Wang
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly‐Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics University of Science and Technology of China Hefei China
| | - Xiang Jiang
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly‐Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics University of Science and Technology of China Hefei China
| | - Jin Zhao
- ICQD/Hefei National Laboratory for Physical Sciences at Microscale, and Key Laboratory of Strongly‐Coupled Quantum Matter Physics, Chinese Academy of Sciences, and Department of Physics University of Science and Technology of China Hefei China
- Department of Physics and Astronomy University of Pittsburgh Pittsburgh Pennsylvania
- Synergetic Innovation Center of Quantum Information & Quantum Physics University of Science and Technology of China Hefei China
| |
Collapse
|
50
|
Wang J, Oruganti B, Durbeej B. A Straightforward Route to Aromatic Excited States in Molecular Motors that Improves Photochemical Efficiency. CHEMPHOTOCHEM 2019. [DOI: 10.1002/cptc.201800268] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Jun Wang
- Division of Theoretical Chemistry, IFMLinköping University 581 83 Linköping Sweden
| | - Baswanth Oruganti
- Division of Theoretical Chemistry, IFMLinköping University 581 83 Linköping Sweden
| | - Bo Durbeej
- Division of Theoretical Chemistry, IFMLinköping University 581 83 Linköping Sweden
| |
Collapse
|