1
|
Zhu Y, Peng J, Xu C, Lan Z. Unsupervised Machine Learning in the Analysis of Nonadiabatic Molecular Dynamics Simulation. J Phys Chem Lett 2024; 15:9601-9619. [PMID: 39270134 DOI: 10.1021/acs.jpclett.4c01751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2024]
Abstract
The all-atomic full-dimensional-level simulations of nonadiabatic molecular dynamics (NAMD) in large realistic systems has received high research interest in recent years. However, such NAMD simulations normally generate an enormous amount of time-dependent high-dimensional data, leading to a significant challenge in result analyses. Based on unsupervised machine learning (ML) methods, considerable efforts were devoted to developing novel and easy-to-use analysis tools for the identification of photoinduced reaction channels and the comprehensive understanding of complicated molecular motions in NAMD simulations. Here, we tried to survey recent advances in this field, particularly to focus on how to use unsupervised ML methods to analyze the trajectory-based NAMD simulation results. Our purpose is to offer a comprehensive discussion on several essential components of this analysis protocol, including the selection of ML methods, the construction of molecular descriptors, the establishment of analytical frameworks, their advantages and limitations, and persistent challenges.
Collapse
Affiliation(s)
- Yifei Zhu
- MOE Key Laboratory of Environmental Theoretical Chemistry, 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
| | - Jiawei Peng
- MOE Key Laboratory of Environmental Theoretical Chemistry, 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
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, 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, 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
|
2
|
Sun K, Vasquez L, Borrelli R, Chen L, Zhao Y, Gelin MF. Interconnection between Polarization-Detected and Population-Detected Signals: Theoretical Results and Ab Initio Simulations. J Chem Theory Comput 2024; 20:7560-7573. [PMID: 39185737 DOI: 10.1021/acs.jctc.4c00592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/27/2024]
Abstract
Most of spectroscopic signals are specified by the nonlinear laser-induced polarization. In recent years, population-detection of signals becomes a trend in femtosecond spectroscopy. Polarization-detected (PD) and population-detected signals are fundamentally different, because they are determined by photoinduced processes acting on disparate time scales. In this work, we consider the fluorescence-detected (FD) N-wave-mixing (NWM) signal as a representative example of population-detected signals, derive a rigorous expression for this signal, and discuss its approximate variants suitable for numerical simulations. This leads us to the definition of the phenomenological FD (PFD) signal, which contains as a special case all definitions of FD signals available in the literature. Then we formulate and prove the population-polarization equivalence (PPE) theorem, which states that PFD NWM signals produced by (possibly strong) laser pulses can be evaluated as conventional PD signals in which the effective polarization is determined by the PFD transition dipole moment operator. We use the PPE theorem for the construction of the ab initio protocol for the simulation of PFD 4WM signals. As an example, we calculate electronic two-dimensional (2D) PFD spectra of the gas-phase pyrazine and compare them with the corresponding PD 2D spectra.
Collapse
Affiliation(s)
- Kewei Sun
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Luis Vasquez
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | | | | | - Yang Zhao
- School of Materials Science and Engineering, Nanyang Technological University, Singapore 639798, Singapore
| | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| |
Collapse
|
3
|
Pios SV, Gelin MF, Luis Vasquez, Hauer J, Chen L. On-the-Fly Simulation of Two-Dimensional Fluorescence-Excitation Spectra. J Phys Chem Lett 2024; 15:8728-8735. [PMID: 39162319 DOI: 10.1021/acs.jpclett.4c01842] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/21/2024]
Abstract
Two-dimensional (2D) fluorescence-excitation (2D-FLEX) spectroscopy is a recently proposed nonlinear femtosecond technique for the detection of photoinduced dynamics. The method records a time-resolved fluorescence signal in its excitation- and detection-frequency dependence and hence combines the exclusive detection of excited state dynamics (fluorescence) with signals resolved in both excitation and emission frequencies (2D electronic spectroscopy). In this work, we develop an on-the-fly protocol for the simulation of 2D-FLEX spectra of molecular systems, which is based on interfacing the classical doorway-window representation of spectroscopic responses with trajectory surface hopping simulations. Applying this methodology to gas-phase pyrazine, we show that femtosecond 2D-FLEX spectra can deliver detailed information that is otherwise obtainable via attosecond spectroscopy.
Collapse
Affiliation(s)
| | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Luis Vasquez
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Jürgen Hauer
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| | | |
Collapse
|
4
|
Perez-Castillo R, Freixas VM, Mukamel S, Martinez-Mesa A, Uranga-Piña L, Tretiak S, Gelin MF, Fernandez-Alberti S. Transient-absorption spectroscopy of dendrimers via nonadiabatic excited-state dynamics simulations. Chem Sci 2024; 15:13250-13261. [PMID: 39183915 PMCID: PMC11339953 DOI: 10.1039/d4sc01019a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Accepted: 07/10/2024] [Indexed: 08/27/2024] Open
Abstract
The efficiency of light-harvesting and energy transfer in multi-chromophore ensembles underpins natural photosynthesis. Dendrimers are highly branched synthetic multi-chromophoric conjugated supra-molecules that mimic these natural processes. After photoexcitation, their repeated units participate in a number of intramolecular electronic energy relaxation and redistribution pathways that ultimately funnel to a sink. Here, a model four-branched dendrimer with a pyrene core is theoretically studied using nonadiabatic molecular dynamics simulations. We evaluate excited-state photoinduced dynamics of the dendrimer, and demonstrate on-the-fly simulations of its transient absorption pump-probe (TA-PP) spectra. We show how the evolutions of the simulated TA-PP spectra monitor in real time photoinduced energy relaxation and redistribution, and provide a detailed microscopic picture of the relevant energy-transfer pathways. To the best of our knowledge, this is the first of this kind of on-the-fly atomistic simulation of TA-PP signals reported for a large molecular system.
Collapse
Affiliation(s)
- Royle Perez-Castillo
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET B1876BXD Bernal Argentina
| | - Victor M Freixas
- Department of Chemistry and Physics and Astronomy, University of California Irvine California 92697-2025 USA
| | - Shaul Mukamel
- Department of Chemistry and Physics and Astronomy, University of California Irvine California 92697-2025 USA
| | - Aliezer Martinez-Mesa
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET B1876BXD Bernal Argentina
- DynAMoS (Dynamical Processes in Atomic and Molecular Systems), Facultad de Física, Universidad de La Habana San Lázaro y L La Habana 10400 Cuba
| | - Llinersy Uranga-Piña
- Departamento de Ciencia y Tecnologia, Universidad Nacional de Quilmes/CONICET B1876BXD Bernal Argentina
- DynAMoS (Dynamical Processes in Atomic and Molecular Systems), Facultad de Física, Universidad de La Habana San Lázaro y L La Habana 10400 Cuba
| | - Sergei Tretiak
- Theoretical Division and Center for Integrated Nanotechnologies, Los Alamos National Laboratory Los Alamos New Mexico 87545 USA
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University Hangzhou 310018 China
| | | |
Collapse
|
5
|
Green D, Bressan G, Heisler IA, Meech SR, Jones GA. Vibrational coherences in half-broadband 2D electronic spectroscopy: Spectral filtering to identify excited state displacements. J Chem Phys 2024; 160:234104. [PMID: 38884412 DOI: 10.1063/5.0214023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Accepted: 06/03/2024] [Indexed: 06/18/2024] Open
Abstract
Vibrational coherences in ultrafast pump-probe (PP) and 2D electronic spectroscopy (2DES) provide insights into the excited state dynamics of molecules. Femtosecond coherence spectra and 2D beat maps yield information about displacements of excited state surfaces for key vibrational modes. Half-broadband 2DES uses a PP configuration with a white light continuum probe to extend the detection range and resolve vibrational coherences in the excited state absorption (ESA). However, the interpretation of these spectra is difficult as they are strongly dependent on the spectrum of the pump laser and the relative displacement of the excited states along the vibrational coordinates. We demonstrate the impact of these convoluting factors for a model based upon cresyl violet. A careful consideration of the position of the pump spectrum can be a powerful tool in resolving the ESA coherences to gain insights into excited state displacements. This paper also highlights the need for caution in considering the spectral window of the pulse when interpreting these spectra.
Collapse
Affiliation(s)
- Dale Green
- Physics, Faculty of Science, University of East Anglia, Norwich NR4 7TJ, United Kingdom
| | - Giovanni Bressan
- School of Chemistry, University of East Anglia, Norwich NR4 7TJ, United Kingdom
| | - Ismael A Heisler
- Instituto de Física, Universidade Federal do Rio Grande do Sul, 91509-900 Porto Alegre, RS, Brazil
| | - Stephen R Meech
- School of Chemistry, University of East Anglia, Norwich NR4 7TJ, United Kingdom
| | - Garth A Jones
- School of Chemistry, University of East Anglia, Norwich NR4 7TJ, United Kingdom
| |
Collapse
|
6
|
Kang M, Nuomin H, Chowdhury SN, Yuly JL, Sun K, Whitlow J, Valdiviezo J, Zhang Z, Zhang P, Beratan DN, Brown KR. Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics. Nat Rev Chem 2024; 8:340-358. [PMID: 38641733 DOI: 10.1038/s41570-024-00595-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/07/2024] [Indexed: 04/21/2024]
Abstract
Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.
Collapse
Affiliation(s)
- Mingyu Kang
- Duke Quantum Center, Duke University, Durham, NC, USA.
- Department of Physics, Duke University, Durham, NC, USA.
| | - Hanggai Nuomin
- Department of Chemistry, Duke University, Durham, NC, USA
| | | | - Jonathon L Yuly
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, USA
| | - Ke Sun
- Duke Quantum Center, Duke University, Durham, NC, USA
- Department of Physics, Duke University, Durham, NC, USA
| | - Jacob Whitlow
- Duke Quantum Center, Duke University, Durham, NC, USA
- Department of Electrical and Computer Engineering, Duke University, Durham, NC, USA
| | - Jesús Valdiviezo
- Kenneth S. Pitzer Theory Center, University of California, Berkeley, CA, USA
- Department of Chemistry, University of California, Berkeley, CA, USA
- Departamento de Ciencias, Sección Química, Pontificia Universidad Católica del Perú, Lima, Peru
| | - Zhendian Zhang
- Department of Chemistry, Duke University, Durham, NC, USA
| | - Peng Zhang
- Department of Chemistry, Duke University, Durham, NC, USA
| | - David N Beratan
- Department of Physics, Duke University, Durham, NC, USA.
- Department of Chemistry, Duke University, Durham, NC, USA.
- Department of Biochemistry, Duke University, Durham, NC, USA.
| | - Kenneth R Brown
- Duke Quantum Center, Duke University, Durham, NC, USA.
- Department of Physics, Duke University, Durham, NC, USA.
- Department of Chemistry, Duke University, Durham, NC, USA.
- Department of Electrical and Computer Engineering, Duke University, Durham, NC, USA.
| |
Collapse
|
7
|
Takahashi H, Borrelli R, Gelin MF, Chen L. Finite temperature dynamics in a polarized sub-Ohmic heat bath: A hierarchical equations of motion-tensor train study. J Chem Phys 2024; 160:164106. [PMID: 38656440 DOI: 10.1063/5.0202312] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2024] [Accepted: 04/08/2024] [Indexed: 04/26/2024] Open
Abstract
The dynamics of the sub-Ohmic spin-boson model under polarized initial conditions at finite temperatures is investigated by employing both analytical tools and the numerically accurate hierarchical equations of motion-tensor train method. By analyzing the features of nonequilibrium dynamics, we discovered a bifurcation phenomenon, which separates two regimes of the dynamics. It is found that before the bifurcation time, increasing temperature slows down the population dynamics, while the opposite effect occurs after the bifurcation time. The dynamics is highly sensitive to both initial preparation of the bath and thermal effects.
Collapse
Affiliation(s)
| | | | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | | |
Collapse
|
8
|
Zhang J, Benavides-Riveros CL, Chen L. Artificial-Intelligence-Based Surrogate Solution of Dissipative Quantum Dynamics: Physics-Informed Reconstruction of the Universal Propagator. J Phys Chem Lett 2024; 15:3603-3610. [PMID: 38527271 DOI: 10.1021/acs.jpclett.4c00598] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/27/2024]
Abstract
The accurate (or even approximate) solution of the equations that govern the dynamics of dissipative quantum systems remains a challenging task in quantum science. While several algorithms have been designed to solve those equations with different degrees of flexibility, they rely mainly on highly expensive iterative schemes. Most recently, deep neural networks have been used for quantum dynamics, but current architectures are highly dependent on the physics of the particular system and usually limited to population dynamics. Here we introduce an artificial-intelligence-based surrogate model that solves dissipative quantum dynamics by parametrizing quantum propagators as Fourier neural operators, which we train using both data set and physics-informed loss functions. Compared with conventional algorithms, our quantum neural propagator avoids time-consuming iterations and provides a universal superoperator that can be used to evolve any initial quantum state for arbitrarily long times. To illustrate the wide applicability of the approach, we employ our quantum neural propagator to compute the population dynamics and time-correlation functions of the Fenna-Matthews-Olson complex.
Collapse
|
9
|
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
|
10
|
Pios SV, Gelin MF, Ullah A, Dral PO, Chen L. Artificial-Intelligence-Enhanced On-the-Fly Simulation of Nonlinear Time-Resolved Spectra. J Phys Chem Lett 2024; 15:2325-2331. [PMID: 38386692 DOI: 10.1021/acs.jpclett.4c00107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2024]
Abstract
Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes in molecules of biological and technological significance. The nonlinear femtosecond signals detected for such systems must be interpreted, but it is a challenging task for which theoretical simulations are often indispensable. Accurate simulations of transient absorption or two-dimensional electronic spectra are, however, computationally very expensive, prohibiting the wider adoption of existing first-principles methods. Here, we report an artificial-intelligence-enhanced protocol to drastically reduce the computational cost of simulating nonlinear time-resolved electronic spectra, which makes such simulations affordable for polyatomic molecules of increasing size. The protocol is based on the doorway-window approach for the on-the-fly surface-hopping simulations. We show its applicability for the prototypical molecule of pyrazine for which it produces spectra with high precision with respect to ab initio reference while cutting the computational cost by at least 95% compared to pure first-principles simulations.
Collapse
Affiliation(s)
- Sebastian V Pios
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou, Zhejiang 310018, People's Republic of China
| | - Arif Ullah
- School of Physics and Optoelectronic Engineering, Anhui University, Hefei, Anhui 230601, People's Republic of China
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, People's Republic of China
| | - Lipeng Chen
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, People's Republic of China
| |
Collapse
|
11
|
Hou E, Sun K, Gelin MF, Zhao Y. Finite temperature dynamics of the Holstein-Tavis-Cummings model. J Chem Phys 2024; 160:084116. [PMID: 38421073 DOI: 10.1063/5.0193471] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 02/05/2024] [Indexed: 03/02/2024] Open
Abstract
By employing the numerically accurate multiple Davydov Ansatz (mDA) formalism in combination with the thermo-field dynamics (TFD) representation of quantum mechanics, we systematically explore the influence of three parameters-temperature, photonic-mode detuning, and qubit-phonon coupling-on population dynamics and absorption spectra of the Holstein-Tavis-Cummings (HTC) model. It is found that elevated qubit-phonon couplings and/or temperatures have a similar impact on all dynamic observables: they suppress the amplitudes of Rabi oscillations in photonic populations as well as broaden the peaks and decrease their intensities in the absorption spectra. Our results unequivocally demonstrate that the HTC dynamics is very sensitive to the concerted variation of the three aforementioned parameters, and this finding can be used for fine-tuning polaritonic transport. The developed mDA-TFD methodology can be efficiently applied for modeling, predicting, optimizing, and comprehensively understanding dynamic and spectroscopic responses of actual molecular systems in microcavities.
Collapse
Affiliation(s)
- Erqin Hou
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Kewei Sun
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
- School of Materials Science and Engineering, Nanyang Technological University, Singapore 639798, Singapore
| | - Yang Zhao
- School of Materials Science and Engineering, Nanyang Technological University, Singapore 639798, Singapore
| |
Collapse
|
12
|
Bruschi M, Gallina F, Fresch B. A Quantum Algorithm from Response Theory: Digital Quantum Simulation of Two-Dimensional Electronic Spectroscopy. J Phys Chem Lett 2024; 15:1484-1492. [PMID: 38295347 DOI: 10.1021/acs.jpclett.3c03499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
Multidimensional optical spectroscopies are powerful techniques to investigate energy transfer pathways in natural and artificial systems. Because of the high information content of the spectra, numerical simulations of the optical response are of primary importance to assist the interpretation of spectral features. However, the increasing complexity of the investigated systems and their quantum dynamics call for the development of novel simulation strategies. In this work, we consider using digital quantum computers. By combining quantum dynamical simulation and nonlinear response theory, we present a quantum algorithm for computing the optical response of molecular systems. The quantum advantage stems from the efficient quantum simulation of the dynamics governed by the molecular Hamiltonian, and it is demonstrated by explicitly considering exciton-vibrational coupling. The protocol is tested on a near-term quantum device, providing the digital quantum simulation of the linear and nonlinear response of simple molecular models.
Collapse
Affiliation(s)
- Matteo Bruschi
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, via Marzolo 1, Padua 35131, Italy
| | - Federico Gallina
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, via Marzolo 1, Padua 35131, Italy
| | - Barbara Fresch
- Dipartimento di Scienze Chimiche, Università degli Studi di Padova, via Marzolo 1, Padua 35131, Italy
- Padua Quantum Technologies Research Center, Università degli Studi di Padova, via Gradenigo 6/A, Padua 35131, Italy
| |
Collapse
|
13
|
Qiang Y, Sun K, Palacino-González E, Shen K, Rao BJ, Gelin MF, Zhao Y. Probing avoided crossings and conical intersections by two-pulse femtosecond stimulated Raman spectroscopy: Theoretical study. J Chem Phys 2024; 160:054107. [PMID: 38341700 DOI: 10.1063/5.0186583] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 01/15/2024] [Indexed: 02/13/2024] Open
Abstract
This study leverages two-pulse femtosecond stimulated Raman spectroscopy (2FSRS) to characterize molecular systems with avoided crossings (ACs) and conical intersections (CIs) in their low-lying excited electronic states. By simulating 2FSRS spectra of microscopically inspired ACs and CIs models, we demonstrate that 2FSRS not only delivers valuable information on the molecular parameters characterizing ACs and CIs but also helps distinguish between these two systems.
Collapse
Affiliation(s)
- Yijia Qiang
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Kewei Sun
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Elisa Palacino-González
- Zernike Institute for Advanced Materials, University of Groningen, 9747 AG Groningen, The Netherlands
| | - Kaijun Shen
- School of Materials Science and Engineering, Nanyang Technological University, Singapore 639798, Singapore
| | - B Jayachander Rao
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Maxim F Gelin
- School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Yang Zhao
- School of Materials Science and Engineering, Nanyang Technological University, Singapore 639798, Singapore
| |
Collapse
|
14
|
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
|
15
|
Gelin MF, Borrelli R. Thermo-Field Dynamics Approach to Photo-induced Electronic Transitions Driven by Incoherent Thermal Radiation. J Chem Theory Comput 2023; 19:6402-6413. [PMID: 37656914 DOI: 10.1021/acs.jctc.3c00590] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/03/2023]
Abstract
The effects of thermal light-matter interaction on the dynamics of photo-induced electronic transitions in molecules are investigated using a novel first principles approach based on the thermo-field dynamics description of both the molecular vibrational modes and of the radiation field. The developed approach permits numerically accurate simulations of quantum dynamics of electronic/excitonic systems coupled to nuclear and photonic baths kept at different temperatures. The baths can be described by arbitrary spectral densities and can have any system-bath coupling strengths. In agreement with the results obtained previously by less rigorous methods, we show that the excitation process obtained by the continuous interaction with the suddenly turned-on thermal radiation field creates a mixed ensemble having a nonnegligible component consisting of a superposition of vibronic eigenstates which can sustain coherent oscillations for relatively long times. The results become especially relevant for the dynamics of electronic transitions upon sunlight excitation. Analytical results based on time-dependent perturbation theory support the numerical simulations and provide a simple interpretation of the time evolution of quantum observables.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Raffaele Borrelli
- DISAFA, University of Torino, Largo Paolo Braccini 2, Grugliasco I-10095, Italy
| |
Collapse
|
16
|
Yang J, Gelin MF, Chen L, Šanda F, Thyrhaug E, Hauer J. Two-dimensional fluorescence excitation spectroscopy: A novel technique for monitoring excited-state photophysics of molecular species with high time and frequency resolution. J Chem Phys 2023; 159:074201. [PMID: 37581414 DOI: 10.1063/5.0156297] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2023] [Accepted: 07/27/2023] [Indexed: 08/16/2023] Open
Abstract
We propose a novel UV/Vis femtosecond spectroscopic technique, two-dimensional fluorescence-excitation (2D-FLEX) spectroscopy, which combines spectral resolution during the excitation process with exclusive monitoring of the excited-state system dynamics at high time and frequency resolution. We discuss the experimental feasibility and realizability of 2D-FLEX, develop the necessary theoretical framework, and demonstrate the high information content of this technique by simulating the 2D-FLEX spectra of a model four-level system and the Fenna-Matthews-Olson antenna complex. We show that the evolution of 2D-FLEX spectra with population time directly monitors energy transfer dynamics and can thus yield direct qualitative insight into the investigated system. This makes 2D-FLEX a highly efficient instrument for real-time monitoring of photophysical processes in polyatomic molecules and molecular aggregates.
Collapse
Affiliation(s)
- Jianmin Yang
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | | | - František Šanda
- Institute of Physics, Faculty of Mathematics and Physics, Charles University, 12116 Prague, Czech Republic
| | - Erling Thyrhaug
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| | - Jürgen Hauer
- Department of Chemistry, Technical University of Munich, D-85747 Garching, Germany
| |
Collapse
|
17
|
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: 6] [Impact Index Per Article: 6.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
|
18
|
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
|
19
|
Gera T, Chen L, Eisfeld A, Reimers JR, Taffet EJ, Raccah DIGB. Simulating optical linear absorption for mesoscale molecular aggregates: An adaptive hierarchy of pure states approach. J Chem Phys 2023; 158:2887556. [PMID: 37125709 DOI: 10.1063/5.0141882] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2023] [Accepted: 03/07/2023] [Indexed: 05/02/2023] Open
Abstract
In this paper, we present dyadic adaptive HOPS (DadHOPS), a new method for calculating linear absorption spectra for large molecular aggregates. This method combines the adaptive HOPS (adHOPS) framework, which uses locality to improve computational scaling, with the dyadic HOPS method previously developed to calculate linear and nonlinear spectroscopic signals. To construct a local representation of dyadic HOPS, we introduce an initial state decomposition that reconstructs the linear absorption spectra from a sum over locally excited initial conditions. We demonstrate the sum over initial conditions can be efficiently Monte Carlo sampled and that the corresponding calculations achieve size-invariant [i.e., O(1)] scaling for sufficiently large aggregates while trivially incorporating static disorder in the Hamiltonian. We present calculations on the photosystem I core complex to explore the behavior of the initial state decomposition in complex molecular aggregates as well as proof-of-concept DadHOPS calculations on an artificial molecular aggregate inspired by perylene bis-imide to demonstrate the size-invariance of the method.
Collapse
Affiliation(s)
- Tarun Gera
- Department of Chemistry, Southern Methodist University, P.O. Box, Dallas, Texas 750314, USA
| | - Lipeng Chen
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, Dresden, Germany
| | - Alexander Eisfeld
- Max Planck Institute for the Physics of Complex Systems, Nöthnitzer Str. 38, Dresden, Germany
| | - Jeffrey R Reimers
- International Centre for Quantum and Molecular Structures and the School of Physics, Shanghai University, 200444 Shanghai, China
- School of Mathematical and Physical Sciences, University of Technology Sydney, Sydney NSW 2007, Australia
| | - Elliot J Taffet
- Department of Chemistry, Southern Methodist University, P.O. Box, Dallas, Texas 750314, USA
| | - Doran I G B Raccah
- Department of Chemistry, Southern Methodist University, P.O. Box, Dallas, Texas 750314, USA
| |
Collapse
|
20
|
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
|
21
|
Leng X, Yan Y, Zhu R, Zou J, Zhang W, Shi Q. Revealing Intermolecular Electronic and Vibronic Coherence with Polarization-Dependent Two-Dimensional Beating Maps. J Phys Chem Lett 2023; 14:838-845. [PMID: 36656105 DOI: 10.1021/acs.jpclett.2c03413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Two-dimensional electronic spectroscopy (2DES) has been widely employed as an efficient tool to reveal the impact of intermolecular electronic and/or vibronic quantum coherence on excitation energy transfer in light-harvesting complexes. However, intramolecular vibrational coherence would also contribute to oscillating signals in 2D spectra, along with the intermolecular coherence signals that are directly related to energy transfer. In this work, the possibility of screening the vibrational coherence signals is explored through polarization-dependent 2DES. The all-parallel (AP) and double-crossed (DC) polarization-dependent two-dimensional rephasing spectra (2DRS) are simulated for a minimalist heterodimer model with vibrational coupling. By combining the DC-2DRS and the 2D beating maps, we demonstrate that the population and vibrational coherence signals can be largely suppressed, resulting in highlighted intermolecular electronic and vibronic coherence signals. Moreover, the AP- and DC-2DBMs show rather different patterns at the vibrational frequency, indicating a possible way to identify pure vibrational coherence.
Collapse
Affiliation(s)
- Xuan Leng
- Department of Physics, School of Physical Science and Technology, Ningbo University, Ningbo 315211, China
| | - Yaming Yan
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China
| | - Ruidan Zhu
- Beijing National Laboratory for Condensed Matter Physics, Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| | - Jiading Zou
- Beijing National Laboratory for Condensed Matter Physics, Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Wenzhao Zhang
- Department of Physics, School of Physical Science and Technology, Ningbo University, Ningbo 315211, China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|