-
Tracking the Electron Density Changes in Excited States -- A Computational Study on Pyrazine
Authors:
Sebastian V. Pios,
Jiaji Zhang,
Maxim F. Gelin,
Hong-Guang Duan,
Lipeng Chen
Abstract:
The development of X-ray free-electron lasers (XFELs) has enabled ultrafast X-ray diffraction (XRD) experiments, which are capable of resolving electronic/vibrational transitions and structural changes in molecules, or capturing molecular movies. While time-resolved XRD has received increasing attention, the extraction of information content from signals is challenging and requires theoretical sup…
▽ More
The development of X-ray free-electron lasers (XFELs) has enabled ultrafast X-ray diffraction (XRD) experiments, which are capable of resolving electronic/vibrational transitions and structural changes in molecules, or capturing molecular movies. While time-resolved XRD has received increasing attention, the extraction of information content from signals is challenging and requires theoretical support. In this work, we combined X-ray scattering theory and trajectory surface hopping approach to resolve dynamical changes in the electronic structure of photo-excited molecules by studying time evolution of electron density changes between electronic excited states and ground state. Using pyrazine molecule as an example, we show that key features of reaction pathways can be identified, enabling the capture of structural changes associated with electronic transitions for a photo-excited molecule.
△ Less
Submitted 27 August, 2024;
originally announced August 2024.
-
On-the-Fly Simulation of Two-Dimensional Fluorescence-Excitation Spectra
Authors:
Sebastian V. Pios,
Maxim F. Gelin,
Luis Vasquez,
Jürgen Hauer,
Lipeng Chen
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 excitatio…
▽ More
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 the gas-phase pyrazine, we show that femtosecond 2D-FLEX spectra can deliver detailed information otherwise obtainable via attosecond spectroscopy.
△ Less
Submitted 21 June, 2024;
originally announced June 2024.
-
Accelerating molecular vibrational spectra simulations with a physically informed deep learning model
Authors:
Yuzhuo Chen,
Sebastian V. Pios,
Maxim F. Gelin,
Lipeng Chen
Abstract:
In recent years, machine learning (ML) surrogate models have emerged as an indispensable tool to accelerate simulations of physical and chemical processes. However, there is still a lack of ML models that can accurately predict molecular vibrational spectra. Here, we present a highly efficient high-dimensional neural network potentials (HD-NNP) architecture to accurately calculate infrared (IR) an…
▽ More
In recent years, machine learning (ML) surrogate models have emerged as an indispensable tool to accelerate simulations of physical and chemical processes. However, there is still a lack of ML models that can accurately predict molecular vibrational spectra. Here, we present a highly efficient high-dimensional neural network potentials (HD-NNP) architecture to accurately calculate infrared (IR) and Raman spectra based on dipole moments and polarizabilities obtained on-the-fly via ML-molecular dynamics (MD) simulations. The methodology is applied to pyrazine, a prototypical polyatomic chromophore. The HD-NNP predicted energies are well within the chemical accuracy (1 kcal/mol), and the errors for HD-NNP predicted forces are only one-half of those obtained from a popular high-performance ML model. Compared to the ab initio reference, the HD-NNP predicted frequencies of IR and Raman spectra differ only by less than 8.3 cm^(-1), and the intensities of IR spectra and the depolarizaiton ratios of Raman spectra are well reproduced. The HD-NNP architecture developed in this work highlights importance of constructing highly accurate NNPs for predicting molecular vibrational spectra.
△ Less
Submitted 10 February, 2024;
originally announced February 2024.
-
Finite temperature dynamics in a polarized sub-Ohmic heat bath: a hierarchical equations of motion-tensor train study
Authors:
Hideaki Takahashi,
Raffaele Borrelli,
Maxim F. Gelin,
Lipeng Chen
Abstract:
Dynamics of the sub-Ohmic spin-boson model under polarized initial conditions at finite temperature 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…
▽ More
Dynamics of the sub-Ohmic spin-boson model under polarized initial conditions at finite temperature 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.
△ Less
Submitted 3 February, 2024;
originally announced February 2024.
-
AI-enhanced on-the-fly simulation of nonlinear time-resolved spectra
Authors:
Sebastian V. Pios,
Maxim F. Gelin,
Arif Ullah,
Pavlo O. Dral,
Lipeng Chen
Abstract:
Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes of 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 elec…
▽ More
Time-resolved spectroscopy is an important tool for unraveling the minute details of structural changes of 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 AI-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 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.
△ Less
Submitted 14 January, 2024;
originally announced January 2024.
-
Dynamics of the spin-boson model: the effect of bath initial conditions
Authors:
Lipeng Chen,
Yiying Yan,
Maxim F. Gelin,
Zhiguo Lü
Abstract:
Dynamics of the (sub-)Ohmic spin-boson model under various bath initial conditions is investigated by employing the Dirac-Frenkel time-dependent variational approach with the multiple Davydov $\mathrm{D_1}$ ansatz in the interaction picture. The validity of our approach is carefully checked by comparing results with those of the hierarchy equations of motion method. By analyzing the features of no…
▽ More
Dynamics of the (sub-)Ohmic spin-boson model under various bath initial conditions is investigated by employing the Dirac-Frenkel time-dependent variational approach with the multiple Davydov $\mathrm{D_1}$ ansatz in the interaction picture. The validity of our approach is carefully checked by comparing results with those of the hierarchy equations of motion method. By analyzing the features of nonequilibrium dynamics, we identify the phase diagrams for different bath initial conditions. We find that for spectral exponent $s<s_c$, there exists a transition from coherent to quasicoherent dynamics with increasing the coupling strengths. For $s_c<s\leq{1}$, the coherent to incoherent crossover occurs at a certain coupling strength, and the quasicoherent dynamics emerges at much larger couplings. The initial preparation of the bath has considerable influence on the dynamics.
△ Less
Submitted 12 December, 2022;
originally announced December 2022.
-
Ultrafast Internal Conversion Dynamics Through the on-the-fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods
Authors:
Chao Xu,
Kunni Lin,
Deping Hu,
Feng Long Gu,
Maxim F. Gelin,
Zhenggang Lan
Abstract:
The ultrafast nonadiabatic internal conversion in azomethane is explored by the on-the-fly trajectory surface-hopping simulations of photoinduced dynamics and femtosecond transient absorption (TA) pump-probe (PP) spectra at three electronic-structure theory levels, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. All these dynamics simulations predict ultrafast internal conversion. On the one hand, the OM2/MR…
▽ More
The ultrafast nonadiabatic internal conversion in azomethane is explored by the on-the-fly trajectory surface-hopping simulations of photoinduced dynamics and femtosecond transient absorption (TA) pump-probe (PP) spectra at three electronic-structure theory levels, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. All these dynamics simulations predict ultrafast internal conversion. On the one hand, the OM2/MRCI and SA-CASSCF methods yield similar excited-state dynamics, while the XMS-CASPT2 method predicts a much slower population decay. On the other hand, the TA PP signals simulated at the SA-CASSCF and XMS-CASPT2 levels show the similar spectral features, particularly for the similar stimulated emission contributions, while the OM2/MRCI signals are quite different. This demonstrates that the nonadiabatic population dynamics and time-resolved stimulated emission signals may reflect different aspects of photoinduced processes. The combination of the dynamical and spectral simulations definitely provides more accurate and detailed information which sheds light on the microscopic mechanisms of photophysical and photochemical processes.
△ Less
Submitted 28 October, 2021;
originally announced October 2021.
-
First-passage time theory of activated rate chemical processes in electronic molecular junctions
Authors:
Riley J. Preston,
Maxim F. Gelin,
Daniel S. Kosov
Abstract:
Confined nanoscale spaces, electric fields and tunneling currents make the molecular electronic junction an experimental device for the discovery of new, out-of-equilibrium chemical reactions. Reaction-rate theory for current-activated chemical reactions is developed by combining a Keldysh nonequilibrium Green's functions treatment of electrons, Fokker-Planck description of the reaction coordinate…
▽ More
Confined nanoscale spaces, electric fields and tunneling currents make the molecular electronic junction an experimental device for the discovery of new, out-of-equilibrium chemical reactions. Reaction-rate theory for current-activated chemical reactions is developed by combining a Keldysh nonequilibrium Green's functions treatment of electrons, Fokker-Planck description of the reaction coordinate, and Kramers' first-passage time calculations. The NEGF provide an adiabatic potential as well as a diffusion coefficient and temperature with local dependence on the reaction coordinate. Van Kampen's Fokker-Planck equation, which describes a Brownian particle moving in an external potential in an inhomogeneous medium with a position-dependent friction and diffusion coefficient, is used to obtain an analytic expression for the first-passage time. The theory is applied to several transport scenarios: a molecular junction with a single, reaction coordinate dependent molecular orbital, and a model diatomic molecular junction. We demonstrate the natural emergence of Landauer's blowtorch effect as a result of the interplay between the configuration dependent viscosity and diffusion coefficients. The resultant localized heating in conjunction with the bond-deformation due to current-induced forces are shown to be the determining factors when considering chemical reaction rates; each of which result from highly tunable parameters within the system.
△ Less
Submitted 28 January, 2021;
originally announced January 2021.
-
A model for dynamical solvent control of molecular junction electronic properties
Authors:
Maxim F. Gelin,
Daniel S. Kosov
Abstract:
Experimental measurements of electron transport properties of molecular junctions are often performed in solvents. Solvent-molecule coupling and physical properties of the solvent can be used as the external stimulus to control electric current through a molecule. In this paper, we propose a model, which includes dynamical effects of solvent-molecule interaction in the non-equilibrium Green's func…
▽ More
Experimental measurements of electron transport properties of molecular junctions are often performed in solvents. Solvent-molecule coupling and physical properties of the solvent can be used as the external stimulus to control electric current through a molecule. In this paper, we propose a model, which includes dynamical effects of solvent-molecule interaction in the non-equilibrium Green's function calculations of electric current. The solvent is considered as a macroscopic dipole moment that reorients stochastically and interacts with the electrons tunnelling through the molecular junction. The Keldysh-Kadanoff-Baym equations for electronic Green's functions are solved in time-domain with subsequent averaging over random realisations of rotational variables using Furutsu-Novikov method for exact closure of infinite hierarchy of stochastic correlation functions. The developed theory requires the use of wide-band approximation as well as classical treatment of solvent degrees of freedom. The theory is applied to a model molecular junction. It is demonstrated that not only electrostatic interaction between molecular junction and solvent but also solvent viscosity can be used to control electrical properties of the junction. Aligning of the rotating dipole moment breaks particle-hole symmetry of the transmission favouring either hole or electron transport channels depending upon the aligning potential.
△ Less
Submitted 28 January, 2021; v1 submitted 2 December, 2020;
originally announced December 2020.
-
Effects of high pulse intensity and chirp in two-dimensional electronic spectroscopy of an atomic vapor
Authors:
Marcel Binz,
Lukas Bruder,
Lipeng Chen,
Maxim F. Gelin,
Wolfgang Domcke,
Frank Stienkemeier
Abstract:
The effects of high pulse intensity and chirp on two-dimensional electronic spectroscopy signals are experimentally investigated in the highly non-perturbative regime using atomic rubidium vapor as clean model system. Data analysis is performed based on higher-order Feynman diagrams and non-perturbative numerical simulations of the system response. It is shown that higher-order contributions may l…
▽ More
The effects of high pulse intensity and chirp on two-dimensional electronic spectroscopy signals are experimentally investigated in the highly non-perturbative regime using atomic rubidium vapor as clean model system. Data analysis is performed based on higher-order Feynman diagrams and non-perturbative numerical simulations of the system response. It is shown that higher-order contributions may lead to a fundamental change of the static appearance and beating-maps of the 2D spectra and that chirped pulses enhance or suppress distinct higher-order pathways. We further give an estimate of the threshold intensity beyond which the high-intensity effects become visible for the system under consideration.
△ Less
Submitted 9 September, 2021; v1 submitted 28 April, 2020;
originally announced April 2020.
-
What can be learned about molecular reorientation from single molecule polarization microscopy?
Authors:
M. F. Gelin,
D. S. Kosov
Abstract:
We have developed a general approach for the calculation of the single molecule polarization correlation function C(t), which delivers a correlation of the emission dichroisms at time 0 and t. The approach is model independent and valid for general asymmetric top molecules. The key dynamic quantities of our analysis are the even-rank orientational correlation functions, the weighted sum of which…
▽ More
We have developed a general approach for the calculation of the single molecule polarization correlation function C(t), which delivers a correlation of the emission dichroisms at time 0 and t. The approach is model independent and valid for general asymmetric top molecules. The key dynamic quantities of our analysis are the even-rank orientational correlation functions, the weighted sum of which yields C(t). We have demonstrated that the use of non-orthogonal schemes for the detection of the single molecule polarization responses makes it possible to manipulate the weighting coefficients in the expansion of C(t). Thus valuable information about the orientational correlation functions of the rank higher than second can be extracted from C(t).
△ Less
Submitted 11 November, 2006;
originally announced November 2006.