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
|
Loreti A, Freixas VM, Avagliano D, Segatta F, Song H, Tretiak S, Mukamel S, Garavelli M, Govind N, Nenov A. WFOT: A Wave Function Overlap Tool between Single- and Multi-Reference Electronic Structure Methods for Spectroscopy Simulation. J Chem Theory Comput 2024; 20:4804-4819. [PMID: 38828948 DOI: 10.1021/acs.jctc.4c00310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2024]
Abstract
We report the development of a novel diagnostic tool, named wave function overlap tool (WFOT), designed to evaluate the overlap between wave functions computed at single-reference [i.e., time-dependent density functional theory or configuration interaction singles (CIS)] and multireference (i.e., CASSCF/CASPT2) electronic structure levels of theory. It relies on truncating the single- and multireference WFs to CIS-like expansions spanning the same configurational space and maximizing the molecular orbital overlap by means of a unitary transformation. To demonstrate the functionality of the tool, we calculate the transient spectrum of acetylacetone by evaluating excited state absorption signals with multireference quality on top of single-reference on-the-fly dynamics simulations. Semiautomatic spectra generation is facilitated by interfacing the tool with the COBRAMM package, which also allows one to use WFOT with several quantum chemistry codes such as Gaussian, NWChem, and OpenMolcas. Other exciting possibilities for the utilization of the code beyond the simulation of transient absorption spectroscopy are eventually discussed.
Collapse
Affiliation(s)
- Alessandro Loreti
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| | - Victor Manuel Freixas
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, United States
| | - Davide Avagliano
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| | - Francesco Segatta
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| | - Huajing Song
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Sergei Tretiak
- Physics and Chemistry of Materials, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Shaul Mukamel
- Department of Chemistry and Department of Physics and Astronomy, University of California, Irvine, California 92697, United States
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", University of Bologna, Via Piero Gobetti 85, Bologna 40129, Italy
| |
Collapse
|
3
|
Suchan J, Liang F, Durden AS, Levine BG. Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone. J Chem Phys 2024; 160:134310. [PMID: 38573851 DOI: 10.1063/5.0198333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 03/18/2024] [Indexed: 04/06/2024] Open
Abstract
Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum-classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?
Collapse
Affiliation(s)
- Jiří Suchan
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| | - Fangchun Liang
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Andrew S Durden
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| | - Benjamin G Levine
- Institute of Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
4
|
Palombo R, Barneschi L, Pedraza-González L, Yang X, Olivucci M. Picosecond quantum-classical dynamics reveals that the coexistence of light-induced microbial and animal chromophore rotary motion modulates the isomerization quantum yield of heliorhodopsin. Phys Chem Chem Phys 2024; 26:10343-10356. [PMID: 38501246 DOI: 10.1039/d4cp00193a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Rhodopsins are light-responsive proteins forming two vast and evolutionary distinct superfamilies whose functions are invariably triggered by the photoisomerization of a single retinal chromophore. In 2018 a third widespread superfamily of rhodopsins called heliorhodopsins was discovered using functional metagenomics. Heliorhodopsins, with their markedly different structural features with respect to the animal and microbial superfamilies, offer an opportunity to study how evolution has manipulated the chromophore photoisomerization to achieve adaptation. One question is related to the mechanism of such a reaction and how it differs from that of animal and microbial rhodopsins. To address this question, we use hundreds of quantum-classical trajectories to simulate the spectroscopically documented picosecond light-induced dynamics of a heliorhodopsin from the archaea thermoplasmatales archaeon (TaHeR). We show that, consistently with the observations, the trajectories reveal two excited state decay channels. However, inconsistently with previous hypotheses, only one channel is associated with the -C13C14- rotation of microbial rhodopsins while the second channel is characterized by the -C11C12- rotation typical of animal rhodopsins. The fact that such -C11C12- rotation is aborted upon decay and ground state relaxation, explains why illumination of TaHeR only produces the 13-cis isomer with a low quantum efficiency. We argue that the documented lack of regioselectivity in double-bond excited state twisting motion is the result of an "adaptation" that could be completely lost via specific residue substitutions modulating the steric hindrance experienced along the isomerization motion.
Collapse
Affiliation(s)
- Riccardo Palombo
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Siena, Italy.
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, USA.
| | - Leonardo Barneschi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Siena, Italy.
| | - Laura Pedraza-González
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Giuseppe Moruzzi, 13, I-56124 Pisa, Italy
| | - Xuchun Yang
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, USA.
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, via A. Moro 2, I-53100 Siena, Siena, Italy.
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, USA.
| |
Collapse
|
5
|
Chen X, Zhang X, Han J, Xia SH. Photochemical Mechanisms of Hydroxyquinoline Benzimidazole: Insights from Electronic Structure Calculations and Nonadiabatic Dynamics Simulations. J Phys Chem A 2024; 128:1984-1992. [PMID: 38446415 DOI: 10.1021/acs.jpca.3c07298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Excited-state intramolecular double proton transfer (ESIDPT) has received much attention because of its widespread existence in the life reactions of living organisms, and materials with this property are significant for their special luminescent properties. In this work, the complete active space self-consistent field (CASSCF) and OM2/multireference configuration interaction (OM2/MRCI) methods have been employed to study the static electronic structure calculations of the photochemistry and the possibility of ESIDPT process of hydroxyquinoline benzimidazole (HQB) molecule, along with the nonadiabatic dynamics simulations. The computational results show that the HQB molecule is relaxed to the S1-ENOL minimum after being excited to the Franck-Condon point in the S1 state. Subsequently, during the nonadiabatic deactivation process, the OH···N proton transfer and the twisting of benzimidazole occur before arriving at the single proton transfer conical intersection S1S0-KETO. Finally, the system can either return to the initial ground-state structure S0-ENOL or to the single proton transfer ground-state structure S0-KETO, both of which have almost the same probability. The dynamics simulations also show that no double proton transfer occurs. The excited-state lifetime of HQB is fitted to 1.1 ps, and only 64% of the dynamic trajectories return to the ground state within the 2.0 ps simulation time. We hope the detailed reaction mechanism of the HQB molecule will provide new insights into similar systems.
Collapse
Affiliation(s)
- Xiaohang Chen
- College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China
| | - XinYu Zhang
- College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China
| | - Juan Han
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Shu-Hua Xia
- College of Life and Environmental Sciences, Minzu University of China, Beijing 100081, China
| |
Collapse
|
6
|
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
|
7
|
Rose PA, Krich JJ. Interpretations of High-Order Transient Absorption Spectroscopies. J Phys Chem Lett 2023; 14:10849-10855. [PMID: 38032056 DOI: 10.1021/acs.jpclett.3c02491] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2023]
Abstract
Transient absorption (TA) spectroscopy is an invaluable tool for determining the energetics and dynamics of excited states. When pump intensities are sufficiently high, TA spectra include both the generally desired third-order response and responses that are higher in order in the field amplitudes. Recent work demonstrated that pump-intensity-dependent TA measurements allow separating the orders of response, but the information content in those higher orders has not been described. We give a general framework for understanding high-order TA spectra. We extend to higher order the fundamental processes of standard TA: ground-state bleach (GSB), stimulated emission (SE), and excited-state absorption (ESA). Each order introduces two new processes: SE and ESA from previously inaccessible highly excited states and negations of lower-order processes. We show the new spectral and dynamical information at each order and show how the relative signs of the signals in different orders can be used to identify which processes dominate.
Collapse
Affiliation(s)
- Peter A Rose
- Department of Physics, University of Ottawa, Ottawa ON K1N 6N5, Canada
| | - Jacob J Krich
- Department of Physics, University of Ottawa, Ottawa ON K1N 6N5, Canada
- Nexus for Quantum Technologies, University of Ottawa, Ottawa ON K1N 6N5, Canada
| |
Collapse
|
8
|
Kaczun T, Dempwolff AL, Huang X, Gelin MF, Domcke W, Dreuw A. Tuning UV Pump X-ray Probe Spectroscopy on the Nitrogen K Edge Reveals the Radiationless Relaxation of Pyrazine: Ab Initio Simulations Using the Quasiclassical Doorway-Window Approximation. J Phys Chem Lett 2023:5648-5656. [PMID: 37310800 DOI: 10.1021/acs.jpclett.3c01018] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Transient absorption UV pump X-ray probe spectroscopy has been established as a versatile technique for the exploration of ultrafast photoinduced dynamics in valence-excited states. In this work, an ab initio theoretical framework for the simulation of time-resolved UV pump X-ray probe spectra is presented. The method is based on the description of the radiation-matter interaction in the classical doorway-window approximation and a surface-hopping algorithm for the nonadiabatic nuclear excited-state dynamics. Using the second-order algebraic-diagrammatic construction scheme for excited states, UV pump X-ray probe signals were simulated for the carbon and nitrogen K edges of pyrazine, assuming a duration of 5 fs of the UV pump and X-ray probe pulses. It is predicted that spectra measured at the nitrogen K edge carry much richer information about the ultrafast nonadiabatic dynamics in the valence-excited states of pyrazine than those measured at the carbon K edge.
Collapse
Affiliation(s)
- Tobias Kaczun
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg D-69120, Germany
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg D-69120, Germany
| | - Xiang Huang
- Department of Chemistry, Technical University of Munich, Garching D-85747, Germany
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, Garching D-85747, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg D-69120, Germany
| |
Collapse
|
9
|
Zhan S, Gelin MF, Huang X, Sun K. Ab initio simulation of peak evolutions and beating maps for electronic two-dimensional signals of a polyatomic chromophore. J Chem Phys 2023; 158:2890773. [PMID: 37191214 DOI: 10.1063/5.0150387] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Accepted: 04/28/2023] [Indexed: 05/17/2023] Open
Abstract
By employing the doorway-window (DW) on-the-fly simulation protocol, we performed ab initio simulations of peak evolutions and beating maps of electronic two-dimensional (2D) spectra of a polyatomic molecule in the gas phase. As the system under study, we chose pyrazine, which is a paradigmatic example of photodynamics dominated by conical intersections (CIs). From the technical perspective, we demonstrate that the DW protocol is a numerically efficient methodology suitable for simulations of 2D spectra for a wide range of excitation/detection frequencies and population times. From the information content perspective, we show that peak evolutions and beating maps not only reveal timescales of transitions through CIs but also pinpoint the most relevant coupling and tuning modes active at these CIs.
Collapse
Affiliation(s)
- Siying Zhan
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - 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
| | - Kewei Sun
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| |
Collapse
|
10
|
Cerezo J, Gao S, Armaroli N, Ingrosso F, Prampolini G, Santoro F, Ventura B, Pastore M. Non-Phenomenological Description of the Time-Resolved Emission in Solution with Quantum-Classical Vibronic Approaches-Application to Coumarin C153 in Methanol. Molecules 2023; 28:molecules28093910. [PMID: 37175320 PMCID: PMC10180259 DOI: 10.3390/molecules28093910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2023] [Revised: 04/27/2023] [Accepted: 04/27/2023] [Indexed: 05/15/2023] Open
Abstract
We report a joint experimental and theoretical work on the steady-state spectroscopy and time-resolved emission of the coumarin C153 dye in methanol. The lowest energy excited state of this molecule is characterized by an intramolecular charge transfer thus leading to remarkable shifts of the time-resolved emission spectra, dictated by the methanol reorganization dynamics. We selected this system as a prototypical test case for the first application of a novel computational protocol aimed at the prediction of transient emission spectral shapes, including both vibronic and solvent effects, without applying any phenomenological broadening. It combines a recently developed quantum-classical approach, the adiabatic molecular dynamics generalized vertical Hessian method (Ad-MD|gVH), with nonequilibrium molecular dynamics simulations. For the steady-state spectra we show that the Ad-MD|gVH approach is able to reproduce quite accurately the spectral shapes and the Stokes shift, while a ∼0.15 eV error is found on the prediction of the solvent shift going from gas phase to methanol. The spectral shape of the time-resolved emission signals is, overall, well reproduced, although the simulated spectra are slightly too broad and asymmetric at low energies with respect to experiments. As far as the spectral shift is concerned, the calculated spectra from 4 ps to 100 ps are in excellent agreement with experiments, correctly predicting the end of the solvent reorganization after about 20 ps. On the other hand, before 4 ps solvent dynamics is predicted to be too fast in the simulations and, in the sub-ps timescale, the uncertainty due to the experimental time resolution (300 fs) makes the comparison less straightforward. Finally, analysis of the reorganization of the first solvation shell surrounding the excited solute, based on atomic radial distribution functions and orientational correlations, indicates a fast solvent response (≈100 fs) characterized by the strengthening of the carbonyl-methanol hydrogen bond interactions, followed by the solvent reorientation, occurring on the ps timescale, to maximize local dipolar interactions.
Collapse
Affiliation(s)
- Javier Cerezo
- Departamento de Química and Institute for Advanced Research in Chemical Sciences (IAdChem), Universidad Autónoma de Madrid, 28049 Madrid, Spain
- Institute of Chemistry of OrganoMetallic Compounds (ICCOM), National Research Council of Italy (CNR), Area di Ricerca di Pisa, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Sheng Gao
- Institute for Organic Synthesis and Photoreactivity (ISOF), National Research Council of Italy (CNR), Via P. Gobetti 101, I-40129 Bologna, Italy
| | - Nicola Armaroli
- Institute for Organic Synthesis and Photoreactivity (ISOF), National Research Council of Italy (CNR), Via P. Gobetti 101, I-40129 Bologna, Italy
| | - Francesca Ingrosso
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), F-54000 Nancy, France
| | - Giacomo Prampolini
- Institute of Chemistry of OrganoMetallic Compounds (ICCOM), National Research Council of Italy (CNR), Area di Ricerca di Pisa, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Fabrizio Santoro
- Institute of Chemistry of OrganoMetallic Compounds (ICCOM), National Research Council of Italy (CNR), Area di Ricerca di Pisa, Via Moruzzi 1, I-56124 Pisa, Italy
| | - Barbara Ventura
- Institute for Organic Synthesis and Photoreactivity (ISOF), National Research Council of Italy (CNR), Via P. Gobetti 101, I-40129 Bologna, Italy
| | - Mariachiara Pastore
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), F-54000 Nancy, France
| |
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
|
Xu Y, Liu C, Ma H. Hierarchical Mapping for Efficient Simulation of Strong System-Environment Interactions. J Chem Theory Comput 2023; 19:426-435. [PMID: 36626721 DOI: 10.1021/acs.jctc.2c00851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Quantum dynamics (QD) simulation is a powerful tool for interpreting ultrafast spectroscopy experiments and unraveling their microscopic mechanism in out-of-equilibrium excited state behaviors in various chemical, biological, and material systems. Although state-of-the-art numerical QD approaches such as the time-dependent density matrix renormalization group (TD-DMRG) already greatly extended the solvable system size of general linearly coupled exciton-phonon models with up to a few hundred phonon modes, the accurate simulation of larger system sizes or strong system-environment interactions is still computationally highly challenging. Based on quantum information theory (QIT), in this work, we realize that only a small number of effective phonon modes couple to the excitonic system directly regardless of a large or even infinite number of modes in the condensed phase environment. On top of the identified small number of direct effective modes, we propose a hierarchical mapping (HM) approach through performing block Lanczos transformations on the remaining indirect modes, which transforms the Hamiltonian matrix to a nearly block-tridiagonal form and eliminates the long-range interactions. Numerical tests on model spin-boson systems and realistic singlet fission models in a rubrene crystal environment with up to 7000 modes and strong system-environment interactions indicate HM can reduce the system size by 1-2 orders of magnitude and accelerate the calculation by ∼80% without losing accuracy.
Collapse
Affiliation(s)
- Yihe Xu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Chungen Liu
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, China
| |
Collapse
|
13
|
Gelin MF, Chen L, Domcke W. Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals. Chem Rev 2022; 122:17339-17396. [PMID: 36278801 DOI: 10.1021/acs.chemrev.2c00329] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Lipeng Chen
- Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Strasse 38, D-01187 Dresden, Germany
| | - Wolfgang Domcke
- Department of Chemistry, Technical University of Munich, D-85747 Garching,Germany
| |
Collapse
|
14
|
Avagliano D, Bonfanti M, Nenov A, Garavelli M. Automatized protocol and interface to simulate QM/MM time-resolved transient absorption at TD-DFT level with COBRAMM. J Comput Chem 2022; 43:1641-1655. [PMID: 35815854 PMCID: PMC9544370 DOI: 10.1002/jcc.26966] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2022] [Revised: 06/07/2022] [Accepted: 06/10/2022] [Indexed: 12/24/2022]
Abstract
We present a series of new implementations that we recently introduced in COBRAMM, the open-source academic software developed in our group. The goal of these implementations is to offer an automatized workflow and interface to simulate time-resolved transient absorption (TA) spectra of medium-to-big chromophore embedded in a complex environment. Therefore, the excited states absorption and the stimulated emission are simulated along nonadiabatic dynamics performed with trajectory surface hopping. The possibility of treating systems from medium to big size is given by the use of time-dependent density functional theory (TD-DFT) and the presence of the environment is taken into account employing a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme. The full implementation includes a series of auxiliary scripts to properly setup the QM/MM system, the calculation of the wavefunction overlap along the dynamics for the propagation, the evaluation of the transition dipole moment at linear response TD-DFT level, and scripts to setup, run and analyze the TA from an ensemble of trajectories. Altogether, we believe that our implementation will open the door to the easily simulate the time-resolved TA of systems so far computationally inaccessible.
Collapse
Affiliation(s)
- Davide Avagliano
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| | - Matteo Bonfanti
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy.,Fondazione Human Technopole - Viale Rita Levi-Montalcini, 1 - Area MIND - Cargo 6 - 20157, Milano, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Bologna, Italy
| |
Collapse
|
15
|
Lei Y, Zheng Z, Vasquez L, Zhao J, Ma J, Ma H. Enhanced Electron Transfer and Spin Flip through Spin-Orbital Couplings in Organic/Inorganic Heterojunctions: A Nonadiabatic Surface Hopping Simulation. J Phys Chem Lett 2022; 13:4840-4848. [PMID: 35616399 DOI: 10.1021/acs.jpclett.2c01177] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The circumstances of transferred electrons across organic/inorganic interfaces have attracted intensive interest because of the distinctive electronic structure properties of those two components. Leveraging ab initio nonadiabatic molecular dynamics methods in conjunction with spin dynamics induced by spin-orbital couplings (SOCs), this study reports two competitive channels during photoinduced dynamical processes in the prototypical ZnPc/monolayer MoS2 heterojunction. Interestingly, the electron-transfer and relaxation processes occur simultaneously because of the enhancement of electron-phonon couplings and expansion of dynamical pathways by SOCs, suggesting that the electron-transfer rate and relaxation processes can be tuned by SOCs, hence yielding the performance promotion of photovoltaic and photocatalytic devices. Additionally, approximately half of the transferred electrons flip their spin within 1.6 ps because of strong SOCs in MoS2, achieving great agreement with experimental measurements. This investigation provides instructive perspectives for designing novel devices and applications based on organic/inorganic heterojunctions, demonstrating the importance of spin dynamics simulations in exploring sophisticated photoinduced processes in materials.
Collapse
Affiliation(s)
- Yuli Lei
- Jiangsu Key Laboratory of Vehicle Emissions Control, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Chemistry and Biomedicine Innovation Center (ChemBIC), School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Zhenfa Zheng
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Luis Vasquez
- Jiangsu Key Laboratory of Vehicle Emissions Control, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Chemistry and Biomedicine Innovation Center (ChemBIC), School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Jin Zhao
- Department of Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jing Ma
- Jiangsu Key Laboratory of Vehicle Emissions Control, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Chemistry and Biomedicine Innovation Center (ChemBIC), School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Haibo Ma
- Jiangsu Key Laboratory of Vehicle Emissions Control, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Chemistry and Biomedicine Innovation Center (ChemBIC), School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| |
Collapse
|