Relation of exact Gaussian basis methods to the dephasing representation: theory and application to time-resolved electronic spectra

Simulating ultrafast transient absorption spectra from first principles using a time-dependent configuration interaction probe

Transient absorption spectroscopy (TAS) is among the most common ultrafast photochemical experiments, but its interpretation remains challenging. In this work, we present an efficient and robust method for simulating TAS signals from first principles. Excited-state absorption and stimulated emission (SE) signals are computed using time-dependent complete active space configuration interaction (TD-CASCI) simulations, leveraging the robustness of time-domain simulation to minimize electronic structure failure. We demonstrate our approach by simulating the TAS signal of 1'-hydroxy-2'-acetonapthone (HAN) from ab initio multiple spawning nonadiabatic molecular dynamics simulations. Our results are compared to gas-phase TAS data recorded from both jet-cooled (T ∼ 40 K) and hot (∼403 K) molecules via cavity-enhanced TAS (CE-TAS). Decomposition of the computed spectrum allows us to assign a rise in the SE signal to excited-state proton transfer and the ultimate decay of the signal to relaxation through a twisted conical intersection. The total cost of computing the observable signal (∼1700 graphics processing unit hours for ∼4 ns of electron dynamics) was markedly less than that of performing the ab initio multiple spawning calculations used to compute the underlying nonadiabatic dynamics.

High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics

Gaussian wavepacket dynamics has proven to be a useful semiclassical approximation for quantum simulations of high-dimensional systems with low anharmonicity. Compared to Heller's original local harmonic method, the variational Gaussian wavepacket dynamics is more accurate, but much more difficult to apply in practice because it requires evaluating the expectation values of the potential energy, gradient, and Hessian. If the variational approach is applied to the local cubic approximation of the potential, these expectation values can be evaluated analytically, but they still require the costly third derivative of the potential. To reduce the cost of the resulting local cubic variational Gaussian wavepacket dynamics, we describe efficient high-order geometric integrators, which are symplectic, time-reversible, and norm-conserving. For small time steps, they also conserve the effective energy. We demonstrate the efficiency and geometric properties of these integrators numerically on a multidimensional, nonseparable coupled Morse potential.

Probing the sensitivity of ab initio multiple spawning to its parameters

Full multiple spawning (FMS) offers a strategy to simulate the nonadiabatic dynamics of molecular systems by describing their nuclear wavefunctions by a linear combination of coupled trajectory basis functions (TBFs). Applying a series of controlled approximations to the full multiple spawning (FMS) equations leads to the ab initio multiple spawning (AIMS), which is compatible with an on-the-fly propagation of the TBFs and an accurate description of nonadiabatic processes. The AIMS strategy and its numerical implementations, however, rely on a series of user-defined parameters. Herein, we investigate the influence of these parameters on the electronic-state population of two molecular systems- trans-azomethane and a two-dimensional model of the butatriene cation. This work highlights the stability of AIMS with respect to most of its parameters, underlines the specific parameters that require particular attention from the user of the method, and offers prescriptions for an informed selection of their value.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00214-023-03004-w.

Automatized protocol and interface to simulate QM/MM time-resolved transient absorption at TD-DFT level with COBRAMM

We present a series of new implementations that we recently introduced in COBRAMM, the open-source academic software developed in our group. The goal of these implementations is to offer an automatized workflow and interface to simulate time-resolved transient absorption (TA) spectra of medium-to-big chromophore embedded in a complex environment. Therefore, the excited states absorption and the stimulated emission are simulated along nonadiabatic dynamics performed with trajectory surface hopping. The possibility of treating systems from medium to big size is given by the use of time-dependent density functional theory (TD-DFT) and the presence of the environment is taken into account employing a hybrid quantum mechanics/molecular mechanics (QM/MM) scheme. The full implementation includes a series of auxiliary scripts to properly setup the QM/MM system, the calculation of the wavefunction overlap along the dynamics for the propagation, the evaluation of the transition dipole moment at linear response TD-DFT level, and scripts to setup, run and analyze the TA from an ensemble of trajectories. Altogether, we believe that our implementation will open the door to the easily simulate the time-resolved TA of systems so far computationally inaccessible.

Dynamics near a conical intersection-A diabolical compromise for the approximations of ab initio multiple spawning

Full multiple spawning (FMS) offers an exciting framework for the development of strategies to simulate the excited-state dynamics of molecular systems. FMS proposes to depict the dynamics of nuclear wavepackets by using a growing set of traveling multidimensional Gaussian functions called trajectory basis functions (TBFs). Perhaps the most recognized method emanating from FMS is the so-called ab initio multiple spawning (AIMS). In AIMS, the couplings between TBFs-in principle exact in FMS-are approximated to allow for the on-the-fly evaluation of required electronic-structure quantities. In addition, AIMS proposes to neglect the so-called second-order nonadiabatic couplings and the diagonal Born-Oppenheimer corrections. While AIMS has been applied successfully to simulate the nonadiabatic dynamics of numerous complex molecules, the direct influence of these missing or approximated terms on the nonadiabatic dynamics when approaching and crossing a conical intersection remains unknown to date. It is also unclear how AIMS could incorporate geometric-phase effects in the vicinity of a conical intersection. In this work, we assess the performance of AIMS in describing the nonadiabatic dynamics through a conical intersection for three two-dimensional, two-state systems that mimic the excited-state dynamics of bis(methylene)adamantyl, butatriene cation, and pyrazine. The population traces and nuclear density dynamics are compared with numerically exact quantum dynamics and trajectory surface hopping results. We find that AIMS offers a qualitatively correct description of the dynamics through a conical intersection for the three model systems. However, any attempt at improving the AIMS results by accounting for the originally neglected second-order nonadiabatic contributions appears to be stymied by the hermiticity requirement of the AIMS Hamiltonian and the independent first-generation approximation.

Finite-Temperature, Anharmonicity, and Duschinsky Effects on the Two-Dimensional Electronic Spectra from Ab Initio Thermo-Field Gaussian Wavepacket Dynamics

Accurate description of finite-temperature vibrational dynamics is indispensable in the computation of two-dimensional electronic spectra. Such simulations are often based on the density matrix evolution, statistical averaging of initial vibrational states, or approximate classical or semiclassical limits. While many practical approaches exist, they are often of limited accuracy and difficult to interpret. Here, we use the concept of thermo-field dynamics to derive an exact finite-temperature expression that lends itself to an intuitive wavepacket-based interpretation. Furthermore, an efficient method for computing finite-temperature two-dimensional spectra is obtained by combining the exact thermo-field dynamics approach with the thawed Gaussian approximation for the wavepacket dynamics, which is exact for any displaced, distorted, and Duschinsky-rotated harmonic potential but also accounts partially for anharmonicity effects in general potentials. Using this new method, we directly relate a symmetry breaking of the two-dimensional signal to the deviation from the conventional Brownian oscillator picture.

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.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

On-the-Fly Ab Initio Semiclassical Evaluation of Electronic Coherences in Polyatomic Molecules Reveals a Simple Mechanism of Decoherence

Irradiation of a molecular system by an intense laser field can trigger dynamics of both electronic and nuclear subsystems. The lighter electrons usually move on much faster, attosecond timescale but the slow nuclear rearrangement damps ultrafast electronic oscillations, leading to the decoherence of the electronic dynamics within a few femtoseconds. We show that a simple, single-trajectory semiclassical scheme can evaluate the electronic coherence time in polyatomic molecules accurately by demonstrating an excellent agreement with full-dimensional quantum calculations. In contrast to numerical quantum methods, the semiclassical one reveals the physical mechanism of decoherence beyond the general blame on nuclear motion. In the propiolic acid, the rate of decoherence and the large deviation from the static frequency of electronic oscillations are quantitatively described with just two semiclassical parameters-the phase space distance and signed area between the trajectories moving on two electronic surfaces. Because it evaluates the electronic structure on the fly, the semiclassical technique avoids the "curse of dimensionality" and should be useful for preselecting molecules for experimental studies.

SSAIMS-Stochastic-Selection Ab Initio Multiple Spawning for Efficient Nonadiabatic Molecular Dynamics

Ab initio multiple spawning provides a powerful and accurate way of describing the excited-state dynamics of molecular systems, whose strength resides in the proper description of coherence effects during nonadiabatic processes thanks to the coupling of trajectory basis functions. However, the simultaneous propagation of a large number of trajectory basis functions can be numerically inconvenient. We propose here an elegant and simple solution to this issue, which consists of (i) detecting uncoupled groups of coupled trajectory basis functions and (ii) selecting stochastically one of these groups to continue the ab initio multiple spawning dynamics. We show that this procedure can reproduce the results of full ab initio multiple spawning dynamics in cases where the uncoupled groups of trajectory basis functions stay uncoupled throughout the dynamics (which is often the case in high-dimensional problems). We present and discuss the aforementioned idea in detail and provide simple numerical applications on indole, ethylene, and protonated formaldimine, highlighting the potential of stochastic-selection ab initio multiple spawning.

A walk through the approximations of ab initio multiple spawning

Full multiple spawning offers an in principle exact framework for excited-state dynamics, where nuclear wavefunctions in different electronic states are represented by a set of coupled trajectory basis functions that follow classical trajectories. The couplings between trajectory basis functions can be approximated to treat molecular systems, leading to the ab initio multiple spawning method which has been successfully employed to study the photochemistry and photophysics of several molecules. However, a detailed investigation of its approximations and their consequences is currently missing in the literature. In this work, we simulate the explicit photoexcitation and subsequent excited-state dynamics of a simple system, LiH, and we analyze (i) the effect of the ab initio multiple spawning approximations on different observables and (ii) the convergence of the ab initio multiple spawning results towards numerically exact quantum dynamics upon a progressive relaxation of these approximations. We show that, despite the crude character of the approximations underlying ab initio multiple spawning for this low-dimensional system, the qualitative excited-state dynamics is adequately captured, and affordable corrections can further be applied to ameliorate the coupling between trajectory basis functions.

Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches

Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.

Path integral approach to the quantum fidelity amplitude

The Loschmidt echo is a measure of quantum irreversibility and is determined by the fidelity amplitude of an imperfect time-reversal protocol. Fidelity amplitude plays an important role both in the foundations of quantum mechanics and in its applications, such as time-resolved electronic spectroscopy. We derive an exact path integral formula for the fidelity amplitude and use it to obtain a series of increasingly accurate semiclassical approximations by truncating an exact expansion of the path integral exponent. While the zeroth-order expansion results in a remarkably simple, yet non-trivial approximation for the fidelity amplitude, the first-order expansion yields an alternative derivation of the so-called 'dephasing representation,' circumventing the use of a semiclassical propagator as in the original derivation. We also obtain an approximate expression for fidelity based on the second-order expansion, which resolves several shortcomings of the dephasing representation. The rigorous derivation from the path integral permits the identification of sufficient conditions under which various approximations obtained become exact.

Spectral lineshapes in nonlinear electronic spectroscopy

We outline a computational approach for nonlinear electronic spectra, which accounts for the electronic energy fluctuations due to nuclear degrees of freedom and explicitly incorporates the fluctuations of higher excited states, induced by the dynamics in the photoactive state(s). This approach is based on mixed quantum-classical dynamics simulations. Tedious averaging over multiple trajectories is avoided by employing the linearly displaced Brownian harmonic oscillator to model the correlation functions. The present strategy couples accurate computations of the high-lying excited state manifold with dynamics simulations. The application is made to the two-dimensional electronic spectra of pyrene, a polycyclic aromatic hydrocarbon characterized by an ultrafast (few tens of femtoseconds) decay from the bright S2 state to the dark S1 state. The spectra for waiting times t2 = 0 and t2 = 1 ps demonstrate the ability of this approach to model electronic state fluctuations and realistic lineshapes. Comparison with experimental spectra [Krebs et al., New Journal of Physics, 2013, 15, 085016] shows excellent agreement and allows us to unambiguously assign the excited state absorption features.

Efficient on-the-fly ab initio semiclassical method for computing time-resolved nonadiabatic electronic spectra with surface hopping or Ehrenfest dynamics

We derive a somewhat crude, yet very efficient semiclassical approximation for computing nonadiabatic spectra. The resulting method, which is a generalization of the multiple-surface dephasing representation, includes quantum effects through interference of mixed quantum-classical trajectories and through quantum treatment of the collective electronic degree of freedom. The method requires very little computational effort beyond the fewest-switches surface hopping or Ehrenfest locally mean-field dynamics and is very easy to implement. The proposed approximation is tested by computing the absorption and time-resolved stimulated emission spectra of pyrazine using the four-dimensional three-surface model which allows for comparison with the numerically exact quantum spectra. As expected, the multiple-surface dephasing representation is not suitable for high-resolution linear spectra, yet it seems to capture all the important features of pump-probe spectra. Finally, the method is combined with on-the-fly ab initio evaluation of the electronic structure (i.e., energies, forces, electric-dipole, and nonadiabatic couplings) in order to compute fully dimensional nonadiabatic spectra of pyrazine without approximations inherent to analytical, including vibronic-coupling models. The Appendix provides derivations of perturbative expressions for linear and pump-probe spectra of arbitrary mixed states and for arbitrary laser pulse shapes.

A new trajectory branching approximation to propagate the mixed quantum-classical Liouville equation

Starting from the mixed quantum-classical Liouville (MQCL) equation, we derive a new trajectory branching method as a modification to the conventional mean field approximation. In the new method, the mean field approximation is used to propagate the mixed quantum-classical dynamics for short times. When the mean field description becomes invalid, new trajectories are added in the simulation by branching the single trajectory into multiple ones. To achieve this, a new set of variables are defined to monitor the deviations of the dynamics on different potential energy surfaces from the reference mean field trajectory, and their equations of motion are derived from the MQCL equation based on the method of first moment expansion. The new method is tested on several one-dimensional two surface problems and is shown to correctly solve the problem of the mean field approximation in several cases.