Excited state dynamics initiated by an electromagnetic field within the Variational Multi-Configurational Gaussian (vMCG) method

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.

Taming Disulfide Bonds with Laser Fields. Nonadiabatic Surface-Hopping Simulations in a Ruthenium Complex

Laser control of chemical reactions is a challenging field of research. In particular, the theoretical description of coupled electronic and nuclear motion in the presence of laser fields is not a trivial task and simulations are mostly restricted to small systems or molecules treated within reduced dimensionality. Here, we demonstrate how the excited state dynamics of [Ru(^S-Sbpy)(bpy)₂]²⁺ can be controlled using explicit laser fields in the context of fewest-switches surface hopping. In particular, the transient properties along the excited state dynamics leading to population of the T₁ minimum energy structure are exploited to define simple laser fields capable of slowing and even completely stopping the onset of S-S bond dissociation. The use of a linear vibronic coupling model to parametrize the potential energy surfaces showcases the strength of the surface-hopping methodology to study systems including explicit laser fields using many nuclear degrees of freedom and a large amount of close-lying electronic excited states.

Validating fewest-switches surface hopping in the presence of laser fields

The capability of fewest-switches surface hopping (FSSH) to describe non-adiabatic dynamics under explicit excitation with external fields is evaluated. Different FSSH parameters are benchmarked against multi-configurational time dependent Hartree (MCTDH) reference calculations using SO₂ and 2-thiocytosine as model, yet realistic, molecular systems. Qualitatively, FSSH is able to reproduce the trends in the MCTDH dynamics with (also without) an explicit external field; however, no set of FSSH parameters is ideal. The adequate treatment of the overcoherence in FSSH is revealed as the driving factor to improve the description of the excitation process with respect to the MCTDH reference. Here, two corrections were tested: the augmented-FSSH (AFSSH) correction and the energy-based decoherence correction. A dependence on the employed basis is detected in AFSSH, performing better when spin-orbit and external laser field couplings are treated as off-diagonal elements instead of projecting them onto the diagonal of the Hamilton operator. In the presence of an electric field, the excited state dynamics was found to depend strongly on the vector used to rescale the kinetic energy along after a transition between surfaces. For SO₂, recurrence of the excited wave packet throughout the duration of the applied laser pulse is observed for laser pulses (>100 fs), resulting in additional interferences missed by FSSH and only visible in variational multi-configurational Gaussian when utilizing a large number of Gaussian basis functions. This feature vanishes when going toward larger molecules, such as 2-thiocytosine, where this effect is barely visible in a laser pulse 200 fs long.

Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies

We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.

The quantum-Ehrenfest method with the inclusion of an IR pulse: Application to electron dynamics of the allene radical cation

We describe the implementation of a laser control pulse in the quantum-Ehrenfest method, a molecular quantum dynamics method that solves the time-dependent Schrödinger equation for both electrons and nuclei. The oscillating electric field-dipole interaction is incorporated directly in the one-electron Hamiltonian of the electronic structure part of the algorithm. We then use the coupled electron-nuclear dynamics of the π-system in the allene radical cation (•CH₂=C=CH₂)⁺ as a simple model of a pump-control experiment. We start (pump) with a two-state superposition of two cationic states. The resulting electron dynamics corresponds to the rapid oscillation of the unpaired electron between the two terminal methylenes. This electron dynamics is, in turn, coupled to the torsional motion of the terminal methylenes. There is a conical intersection at 90° twist, where the electron dynamics collapses because the adiabatic states become degenerate. After passing the conical intersection, the electron dynamics revives. The IR pulse (control) in our simulations is timed to have its maximum at the conical intersection. Our simulations show that the effect of the (control) pulse is to change the electron dynamics at the conical intersection and, as a consequence, the concomitant nuclear dynamics, which is dominated by the change in the torsional angle.

Understanding and Designing Thermally Activated Delayed Fluorescence Emitters: Beyond the Energy Gap Approximation

In this article recent progress in the development of molecules exhibiting Thermally Activated Delayed Fluorescence (TADF) is discussed with a particular focus upon their application as emitters in highly efficient organic light emitting diodes (OLEDs). The key aspects controlling the desirable functional properties, e. g. fast intersystem crossing, high radiative rate and unity quantum yield, are introduced with a particular focus upon the competition between the key requirements needed to achieve high performance OLEDs. The design rules required for organic and metal organic materials are discussed, and the correlation between them outlined. Recent progress towards understanding the influence of the interaction between a molecule and its environment are explained as is the role of the mechanism for excited state formation in OLEDs. Finally, all of these aspects are combined to discuss the ability to implement high level design rules for achieving higher quality materials for commercial applications. This article highlights the significant progress that has been made in recent years, but also outlines the significant challenges which persist to achieve a full understanding of the TADF mechanism and improve the stability and performance of these materials.

Ultrafast nonadiabatic dynamics probed by nitrogen K-edge absorption spectroscopy

Quantum dynamics simulations are used to simulate the ultrafast X-ray Absorption Near-Edge Structure (XANES) spectra of photoexcited pyrazine including two strongly coupled electronically excited states and four normal mode degrees of freedom.

Time-resolved near-edge X-ray absorption fine structure of pyrazine from electronic structure and nuclear wave packet dynamics simulations

As a demonstration of the analysis of the electronic structure and the nuclear dynamics from time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS), we present the TR-NEXAFS spectra of pyrazine following the excitation to the ¹B_2u(ππ^*) state. The spectra are calculated combining the frozen-core/core-valence separated equation-of-motion coupled cluster singles and doubles approach for the spectral signatures and the multiconfiguration time-dependent Hartree method for the wave packet propagation. The population decay from the ¹B_2u(ππ^*) state to the ¹B_3u(nπ^*) and ¹A_u(nπ^*) states, followed by oscillatory flow of population between the ¹B_3u(nπ^*) and ¹A_u(nπ^*) states, is interpreted by means of visualization of the potential energy curves and the reduced nuclear densities. By examining the electronic structure of the three valence-excited states and the final core-excited states, we observe that the population dynamics is explicitly reflected in the TR-NEXAFS spectra, especially when the heteroatoms are selected as the X-ray absorption sites. This work illustrates the feasibility of extracting fine details of molecular photophysical processes from TR-NEXAFS spectra by using currently available theoretical methods.