Multi-faceted spectroscopic mapping of ultrafast nonadiabatic dynamics near conical intersections: A computational study

Interconnection between Polarization-Detected and Population-Detected Signals: Theoretical Results and Ab Initio Simulations

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.

On-the-Fly Simulation of Two-Dimensional Fluorescence-Excitation Spectra

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.

Artificial-Intelligence-Enhanced On-the-Fly Simulation of Nonlinear Time-Resolved Spectra

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.

What can molecular assembly learn from catalysed assembly in living organisms?

Molecular assembly is the process of organizing individual molecules into larger structures and complex systems. The self-assembly approach is predominantly utilized in creating artificial molecular assemblies, and was believed to be the primary mode of molecular assembly in living organisms as well. However, it has been shown that the assembly of many biological complexes is "catalysed" by other molecules, rather than relying solely on self-assembly. In this review, we summarize these catalysed-assembly (catassembly) phenomena in living organisms and systematically analyse their mechanisms. We then expand on these phenomena and discuss related concepts, including catalysed-disassembly and catalysed-reassembly. Catassembly proves to be an efficient and highly selective strategy for synergistically controlling and manipulating various noncovalent interactions, especially in hierarchical molecular assemblies. Overreliance on self-assembly may, to some extent, hinder the advancement of artificial molecular assembly with powerful features. Furthermore, inspired by the biological catassembly phenomena, we propose guidelines for designing artificial catassembly systems and developing characterization and theoretical methods, and review pioneering works along this new direction. Overall, this approach may broaden and deepen our understanding of molecular assembly, enabling the construction and control of intelligent assembly systems with advanced functionality.

Probing avoided crossings and conical intersections by two-pulse femtosecond stimulated Raman spectroscopy: Theoretical study

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.

Finite-Temperature Hole-Magnon Dynamics in an Antiferromagnet

Employing the numerically accurate multiple Davydov Ansatz in combination with the thermo-field dynamics approach, we delve into the interplay of the finite-temperature dynamics of holes and magnons in an antiferromagnet, which allows for scrutinizing previous predictions from the self-consistent Born approximation while offering, for the first time, accurate finite-temperature computation of detailed magnon dynamics as a response and a facilitator to the hole motion. The study also uncovers a pronounced temperature dependence of the magnon and hole populations, pointing to the feasibility of potential thermal manipulation and control of hole dynamics. Our methodology can be applied not only to the calculation of steady-state angular-resolved photoemission spectra but also to the simulation of femtosecond terahertz pump-probe and other nonlinear signals for the characterization of antiferromagnetic materials.

Ab initio simulation of peak evolutions and beating maps for electronic two-dimensional signals of a polyatomic chromophore

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.

The hierarchy of Davydov's Ansätze: From guesswork to numerically "exact" many-body wave functions

This Perspective presents an overview of the development of the hierarchy of Davydov's Ansätze and a few of their applications in many-body problems in computational chemical physics. Davydov's solitons originated in the investigation of vibrational energy transport in proteins in the 1970s. Momentum-space projection of these solitary waves turned up to be accurate variational ground-state wave functions for the extended Holstein molecular crystal model, lending unambiguous evidence to the absence of formal quantum phase transitions in Holstein systems. The multiple Davydov Ansätze have been proposed, with increasing Ansatz multiplicity, as incremental improvements of their single-Ansatz parents. For a given Hamiltonian, the time-dependent variational formalism is utilized to extract accurate dynamic and spectroscopic properties using Davydov's Ansätze as its trial states. A quantity proven to disappear for large multiplicities, the Ansatz relative deviation is introduced to quantify how closely the Schrödinger equation is obeyed. Three finite-temperature extensions to the time-dependent variation scheme are elaborated, i.e., the Monte Carlo importance sampling, the method of thermofield dynamics, and the method of displaced number states. To demonstrate the versatility of the methodology, this Perspective provides applications of Davydov's Ansätze to the generalized Holstein Hamiltonian, variants of the spin-boson model, and systems of cavity-assisted singlet fission, where accurate dynamic and spectroscopic properties of the many-body systems are given by the Davydov trial states.

Exciton Dynamics and Time-Resolved Fluorescence in Nanocavity-Integrated Monolayers of Transition-Metal Dichalcogenides

We have developed an ab initio-based, fully quantum, numerically accurate methodology for the simulation of the exciton dynamics and time- and frequency-resolved fluorescence spectra of the cavity-controlled two-dimensional materials at finite temperatures and applied this methodology to the single-layer WSe₂ system. Specifically, the multiple Davydov D₂ Ansatz has been employed in combination with the method of thermofield dynamics for the finite-temperature extension of accurate time-dependent variation. This allowed us to establish dynamical and spectroscopic signatures of the polaronic and polaritonic effects as well as uncover their characteristic time scales in the relevant range of temperatures. Our study reveals the pivotal role of multidimensional conical intersections in controlling the many-body dynamics of highly intertwined excitonic, phononic, and photonic modes.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

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.

Accurate Simulation of Spectroscopic Signatures of Cavity-Assisted, Conical-Intersection-Controlled Singlet Fission Processes

A numerically accurate, fully quantum methodology has been developed for the simulation of the dynamics and nonlinear spectroscopic signals of cavity-assisted, conical-intersection-controlled singlet fission systems. The methodology is capable of handling several molecular systems strongly coupled to the photonic mode of the cavity and treats the intrinsic conical intersection and cavity-induced polaritonic conical intersections in a numerically exact manner. Contributions of higher-lying molecular electronic states are accounted for comprehensively. The intriguing process of cavity-modified fission dynamics, including all of its electronic, vibrational, and photonic degrees of freedom, together with its two-dimensional spectroscopic manifestation, is simulated for two rubrene dimers strongly coupled to the cavity mode.

Engineering Cavity Singlet Fission in Rubrene

By employing the numerically exact multiple Davydov D₂ ansatz, we study cavity-manipulated singlet fission that is mediated by polaritonic conical intersections for both one- and two-molecule systems. The population evolution of the TT state and the cavity photons is carefully examined in search for a high fission efficiency via cavity engineering. Several interesting mechanisms have been uncovered, such as photon-assisted singlet fission, system localization via a displaced photon state, and collective enhancement of the fission efficiency for the two-molecule system. It is also found that the system localization process in the two-molecule system differs substantially from that in the one-molecule system because of the appearance of a novel central polaritonic conical intersection in the two-molecule system. It has been demonstrated that the cavity-controlled singlet fission process can be switched on and off by controlling the average pumping photon number.

Probing nonadiabatic dynamics with attosecond pulse trains and soft x-ray Raman spectroscopy

Linear off-resonant x-ray Raman techniques are capable of detecting the ultrafast electronic coherences generated when a photoexcited wave packet passes through a conical intersection. A hybrid femtosecond or attosecond probe pulse is employed to excite the system and stimulate the emission of the signal photon, where both fields are components of a hybrid pulse scheme. In this paper, we investigate how attosecond pulse trains, as provided by high-harmonic generation processes, perform as probe pulses in the framework of this spectroscopic technique, instead of single Gaussian pulses. We explore different combination schemes for the probe pulse as well as the impact of parameters of the pulse trains on the signals. Furthermore, we show how Raman selection rules and symmetry consideration affect the spectroscopic signal, and we discuss the importance of vibrational contributions to the overall signal. We use two different model systems, representing molecules of different symmetries, and quantum dynamics simulations to study the difference in the spectra. The results suggest that such pulse trains are well suited to capture the key features associated with the electronic coherence.

Dynamics of disordered Tavis-Cummings and Holstein-Tavis-Cummings models

By employing the time-dependent variational principle and the versatile multi-D₂ Davydov trial states, in combination with the Green's function method, we study the dynamics of the Tavis-Cummings model and the Holstein-Tavis-Cummings model in the presence of diagonal disorder and cavity-qubit coupling disorder. For the Tavis-Cummings model, time evolution of the photon population, the optical absorption spectra, and the hetero-entanglement between the qubits and the cavity mode are calculated by using the Green's function method to corroborate numerically exact results of Davydov's Ansätze. For the Holstein-Tavis-Cummings model, only the latter is utilized to simulate effects of disorder on the photon population dynamics and the absorption spectra. We have demonstrated that the complementary employment of analytical and numerical methods permits uncovering a fairly comprehensive picture of a variety of complex behaviors in disordered multidimensional polaritonic cavity quantum electrodynamics systems.

Ab Initio Quasiclassical Simulation of Femtosecond Time-Resolved Two-Dimensional Electronic Spectra of Pyrazine

Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.

Efficient quantum dynamics simulations of complex molecular systems: A unified treatment of dynamic and static disorder

We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.

Time-Resolved Photoelectron Spectroscopy of Conical Intersections with Attosecond Pulse Trains

Conical Intersections (CIs), which are believed to be ubiquitous in molecular and biological systems, open up ultrafast nonradiative decay channels. A superposition of electronic states is created when a molecule passes through a CI and the nuclear wave packet branches. The resulting electronic coherence can be considered a unique signature of the CI. The involved electronic states can be resolved in the energy domain with photoelectron spectroscopy using a femtosecond pulse as a probe. However, the observation of the created electronic coherence in the time domain requires probe pulses with several electron volts of bandwidth. Attosecond pulses can probe the electronic coherence but are unable to resolve the involved electronic states. In this Letter, we propose to address this restriction by using time-resolved photoelectron spectroscopy with an attosecond pulse train as a probe. We theoretically demonstrate that the resulting photoelectron spectrum may yield energy resolution as well as the information on the created coherences in the time domain.

X-ray transient absorption reveals the 1Au (nπ*) state of pyrazine in electronic relaxation

Electronic relaxation in organic chromophores often proceeds via states not directly accessible by photoexcitation. We report on the photoinduced dynamics of pyrazine that involves such states, excited by a 267 nm laser and probed with X-ray transient absorption spectroscopy in a table-top setup. In addition to the previously characterized ¹B_2u (ππ*) (S₂) and ¹B_3u (nπ*) (S₁) states, the participation of the optically dark ¹A_u (nπ*) state is assigned by a combination of experimental X-ray core-to-valence spectroscopy, electronic structure calculations, nonadiabatic dynamics simulations, and X-ray spectral computations. Despite ¹A_u (nπ*) and ¹B_3u (nπ*) states having similar energies at relaxed geometry, their X-ray absorption spectra differ largely in transition energy and oscillator strength. The ¹A_u (nπ*) state is populated in 200 ± 50 femtoseconds after electronic excitation and plays a key role in the relaxation of pyrazine to the ground state.

Combined Surface-Hopping, Dyson Orbital, and B-Spline Approach for the Computation of Time-Resolved Photoelectron Spectroscopy Signals: The Internal Conversion in Pyrazine

A computational protocol for simulating time-resolved photoelectron signals of medium-sized molecules is presented. The procedure is based on a trajectory surface-hopping description of the excited-state dynamics and a combined Dyson orbital and multicenter B-spline approach for the computation of cross sections and asymmetry parameters. The accuracy of the procedure has been illustrated for the case of ultrafast internal conversion of gas-phase pyrazine excited to the ¹B_2u(ππ*) state. The simulated spectra and the asymmetry map are compared to the experimental data, and a very good agreement was obtained without applying any energy-dependent rescaling or broadening. An interesting side result of this work is the finding that the signature of the ¹A_u(nπ*) state is indistinguishable from that of the ¹B_3u(nπ*) state in the time-resolved photoelectron spectrum. By locating four symmetrically equivalent minima on the lowest-excited (S₁) adiabatic potential energy surface of pyrazine, we revealed the strong vibronic coupling of the ¹A_u(nπ*) and ¹B_3u(nπ*) states near the S₁ ← S₀ band origin.

Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex

Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.

Simulation of Time- and Frequency-Resolved Four-Wave-Mixing Signals at Finite Temperatures: A Thermo-Field Dynamics Approach

We propose a new approach to simulate four-wave-mixing signals of molecular systems at finite temperatures by combining the multiconfigurational Ehrenfest method with the thermo-field dynamics theory. In our approach, the four-time correlation functions at finite temperatures are mapped onto those at zero temperature in an enlarged Hilbert space with twice the vibrational degrees of freedom. As an illustration, we have simulated three multidimensional spectroscopic signals, time- and frequency-resolved fluorescence spectra, transient-absorption pump-probe spectra, and electronic two-dimensional (2D) spectra at finite temperatures, for a conical intersection-mediated singlet fission model of a rubrene crystal. It is shown that a detailed dynamical picture of the singlet fission process can be extracted from the three spectroscopic signals. An increasing temperature leads to lower intensities of the signals and broadened vibrational peaks, which can be attributed to faster singlet-triplet population transfer and stronger bath-induced electronic dephasing at higher temperatures.

Identification of an ultrafast internal conversion pathway of pyrazine by time-resolved vacuum ultraviolet photoelectron spectrum simulations

The internal conversion from the optically bright S₂ (¹B_2u, ππ*) state to the dark S₁ (¹B_3u, nπ*) state in pyrazine is a standard benchmark for experimental and theoretical studies on ultrafast radiationless decay. Since 2008, a few theoretical groups have suggested significant contributions of other dark states S₃ (¹A_u, nπ*) and S₄ (¹B_2g, nπ*) to the decay of S₂. We have previously reported the results of nuclear wave packet simulations [Kanno et al., Phys. Chem. Chem. Phys. 17, 2012 (2015)] and photoelectron spectrum calculations [Mignolet et al., Chem. Phys. 515, 704 (2018)] that support the conventional two-state picture. In this article, the two different approaches, i.e., wave packet simulation and photoelectron spectrum calculation, are combined: We computed the time-resolved vacuum ultraviolet photoelectron spectrum and photoelectron angular distribution for the ionization of the wave packet transferred from S₂ to S₁. The present results reproduce almost all the characteristic features of the corresponding experimental time-resolved spectrum [Horio et al., J. Chem. Phys. 145, 044306 (2016)], such as a rapid change from a three-band to two-band structure. This further supports the existence and character of the widely accepted pathway (S₂ → S₁) of ultrafast internal conversion in pyrazine.

Spiers Memorial Lecture: Introduction to ultrafast spectroscopy and imaging of photochemical reactions

A brief overview is presented on ultrafast spectroscopy and imaging of photochemical reactions by highlighting several experimental studies reported in the last five years. A particular focus is placed on new experiments performed using high-order harmonic generation, X-ray free electron lasers, and relativistic electron beams. Exploration of fundamental chemical reaction dynamics using these advanced experimental methodologies is in an early stage, and exciting new research opportunities await in this rapidly expanding and advancing research field. At the same time, there is no experimental methodology that provides all aspects of the electronic and structural dynamics in a single experiment, and investigations using different methodologies with various perspectives need to be considered in a comprehensive manner.

Capturing fingerprints of conical intersection: Complementary information of non-adiabatic dynamics from linear x-ray probes

Many recent experimental ultrafast spectroscopy studies have hinted at non-adiabatic dynamics indicating the existence of conical intersections, but their direct observation remains a challenge. The rapid change of the energy gap between the electronic states complicated their observation by requiring bandwidths of several electron volts. In this manuscript, we propose to use the combined information of different x-ray pump-probe techniques to identify the conical intersection. We theoretically study the conical intersection in pyrrole using transient x-ray absorption, time-resolved x-ray spontaneous emission, and linear off-resonant Raman spectroscopy to gather evidence of the curve crossing.

Ab Initio Surface-Hopping Simulation of Femtosecond Transient-Absorption Pump-Probe Signals of Nonadiabatic Excited-State Dynamics Using the Doorway-Window Representation

An ab initio theoretical framework for the simulation of femtosecond time-resolved transient absorption (TA) pump-probe (PP) spectra with quasi-classical trajectories is presented. The simulations are based on the classical approximation to the doorway-window (DW) representation of third-order four-wave-mixing signals. The DW formula accounts for the finite duration and spectral shape of the pump and probe pulses. In the classical DW formalism, classical trajectories are stochastically sampled from a positive definite doorway distribution, and the signals are evaluated by averaging over a positive definite window distribution. Nonadiabatic excited-state dynamics is described by a stochastic surface-hopping algorithm. The method has been implemented for the pyrazine molecule with the second-order algebraic-diagrammatic construction (ADC(2)) ab initio electronic-structure method. The methodology is illustrated by ab initio simulations of the ground-state bleach, stimulated emission, and excited-state absorption contributions to the TA PP spectrum of gas-phase pyrazine.

Efficient simulation of time- and frequency-resolved four-wave-mixing signals with a multiconfigurational Ehrenfest approach

We have extended the multiconfigurational Ehrenfest approach to the simulation of four-wave-mixing signals of systems involving multiple electronic and vibrational degrees of freedom. As an illustration, we calculate signals of three widely used spectroscopic techniques, time- and frequency-resolved fluorescence spectroscopy, transient absorption spectroscopy, and two-dimensional (2D) electronic spectroscopy, for a two-electronic-state, twenty-four vibrational-mode conical intersection model. It has been shown that all these three spectroscopic signals characterize fast population transfer from the higher excited electronic state to the lower excited electronic state. While the time- and frequency-resolved spectrum maps the wave packet propagation exclusively on the electronically excited states, the transient absorption and 2D electronic spectra reflect the wave packet dynamics on both electronically excited states and the electronic ground state. Combining trajectory-guided Gaussian basis functions and the nonlinear response function formalism, the present approach provides a promising general technique for the applications of various Gaussian basis methods to the calculations of four-wave-mixing spectra of polyatomic molecules.