Real-Time CASSCF (Ehrenfest) Modeling of Electron Dynamics in Organic Semiconductors. Dynamics Reaction Paths Driven by Quantum Coherences. Application to a Radical Organic Semiconductor

We present a strategy for the modeling of charge carrier dynamics in organic semiconductors using conventional quantum chemistry methods, including the analytic gradient for nuclear motion. The theoretical approach uses real-time CASSCF (Ehrenfest) all-electron dynamics coupled to classical nuclear dynamics for the special case of a small number (4-8) of molecular units. The objective is to obtain mechanistic/atomistic insight at the electronic structure level, relating to spin density dynamics, to the effect of crystal structure (e.g., slippage between spin/charge carriers), and to ferromagnetic and antiferromagnetic effects. The initial conditions for our simulations use the equilibrium structures of all the molecular units. At this geometry, a localized hole on one of the units corresponds to a coherent superposition of adiabatic states. We thus generate a dynamics reaction path driven by quantum coherences. Our aim is to inform experiment and to compare with parametrized theoretical models. The methodology is demonstrated for a perfectly π-stacked ethylene model (up to 8 eclipsed molecular units) for both hole transfer and localized exciton transfer. An application for hole transfer is presented for bisdithiazolyl (S,S) and bisdiselenazolyl (Se,Se) radicals for the special case of ferromagnetic coupling. For these examples, the embedded pyridine radical model organic chromophore (up to 6 eclipsed π-stacked molecular units) has been studied on its own as well as the target bisdithiazolyl (S,S) and bisdiselenazolyl (Se,Se) systems. A significant difference between these systems and the ethylene and pyridine stacks is that the (S,S) and (Se,Se) systems exhibit molecular slippage rather than being perfectly eclipsed. This slippage may result from crystal defects or intermolecular vibrations. For the model systems, the electron dynamics is dominated by the initial and final molecular units, irrespective of the length of the chain. The intervening units act as a "superexchange bridge". Our simulations reveal that, in the presence of slippage, charge migration cannot propagate across the entire system; instead, the coherence length is limited to 3 molecular units. The results also suggest that the mechanism of charge transport is different for bisdiselenazolyl (Se,Se) (superexchange-like A -[B]→ C) and bisdithiazolyl (S,S) (direct A → C). An analysis of the spin density suggests that, in the charge carrier dynamics, the additional charge carried by the Se versus S in the "scaffold" is small. Since we use a small number of molecular units, the coupled nuclear dynamics is seen to be complementary to the electron dynamics (i.e., creating a hole causes bond length contraction while filling a hole with an electron lengthens the bond). In all the cases studied, the mechanism of charge mobility is wave-like, rather than hopping, because we use the time dependent Schrödinger equation to propagate the electronic wave function.

Simulating Attochemistry: Which Dynamics Method to Use?

Attochemistry aims to exploit the properties of coherent electronic wavepackets excited via attosecond pulses to control the formation of photoproducts. Such molecular processes can, in principle, be simulated with various nonadiabatic dynamics methods, yet the impact of the approximations underlying the methods is rarely assessed. The performances of widely used mixed quantum-classical approaches, Tully surface hopping, and classical Ehrenfest methods are evaluated against the high-accuracy DD-vMCG quantum dynamics. This comparison is conducted for the valence ionization of fluorobenzene. Analyzing the nuclear motion induced in the branching space of the nearby conical intersection, the results show that the mixed quantum-classical methods reproduce quantitatively the average motion of a quantum wavepacket when initiated on a single electronic state. However, they fail to properly capture the nuclear motion induced by an electronic wavepacket along the derivative coupling, the latter originating from the quantum electronic coherence property, key to attochemistry.

Control of nuclear dynamics in the benzene cation by electronic wavepacket composition

The study of coupled electron-nuclear dynamics driven by coherent superpositions of electronic states is now possible in attosecond science experiments. The objective is to understand the electronic control of chemical reactivity. In this work we report coherent 8-state non-adiabatic electron-nuclear dynamics simulations of the benzene radical cation. The computations were inspired by the extreme ultraviolet (XUV) experimental results in which all 8 electronic states were prepared with significant population. Our objective was to study the nuclear dynamics using various bespoke coherent electronic state superpositions as initial conditions in the Quantum-Ehrenfest method. The original XUV measurements were supported by Multi-configuration time-dependent Hartree (MCTDH) simulations, which suggested a model of successive passage through conical intersections. The present computations support a complementary model where non-adiabatic events are seen far from a conical intersection and are controlled by electron dynamics involving non-adjacent adiabatic states. It proves to be possible to identify two superpositions that can be linked with two possible fragmentation paths.

Sub-Femtosecond Stark Control of Molecular Photoexcitation with Near Single-Cycle Pulses

Electric fields can tailor molecular potential energy surfaces by interaction with the electronic state-dependent molecular dipole moment. Recent developments in optics have enabled the creation of ultrashort few-cycle optical pulses with precise control of the carrier envelope phase (CEP) that determines the offset of the maxima in the field and the pulse envelope. This opens news ways of controlling ultrafast molecular dynamics by exploiting the CEP. In this work, we show that the photoabsorption efficiency of oriented H₂CSO (sulfine) can be controlled by tuning the CEP. We further show that this control emanates from a resonance condition related to Stark shifting of the electronic energy levels.

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.

The importance of Rydberg orbitals in dissociative ionization of small hydrocarbon molecules in intense laser fields

Much of our intuition about strong-field processes is built upon studies of diatomic molecules, which typically have electronic states that are relatively well separated in energy. In polyatomic molecules, however, the electronic states are closer together, leading to more complex interactions. A combined experimental and theoretical investigation of strong-field ionization followed by hydrogen elimination in the hydrocarbon series C₂D₂, C₂D₄ and C₂D₆ reveals that the photofragment angular distributions can only be understood when the field-dressed orbitals rather than the field-free orbitals are considered. Our measured angular distributions and intensity dependence show that these field-dressed orbitals can have strong Rydberg character for certain orientations of the molecule relative to the laser polarization and that they may contribute significantly to the hydrogen elimination dissociative ionization yield. These findings suggest that Rydberg contributions to field-dressed orbitals should be routinely considered when studying polyatomic molecules in intense laser fields.

Communication: XFAIMS-eXternal Field Ab Initio Multiple Spawning for electron-nuclear dynamics triggered by short laser pulses

Attoscience is an emerging field where attosecond pulses or few cycle IR pulses are used to pump and probe the correlated electron-nuclear motion of molecules. We present the trajectory-guided eXternal Field Ab Initio Multiple Spawning (XFAIMS) method that models such experiments "on-the-fly," from laser pulse excitation to fragmentation or nonadiabatic relaxation to the ground electronic state. For the photoexcitation of the LiH molecule, we show that XFAIMS gives results in close agreement with numerically exact quantum dynamics simulations, both for atto- and femtosecond laser pulses. We then show the ability of XFAIMS to model the dynamics in polyatomic molecules by studying the effect of nuclear motion on the photoexcitation of a sulfine (H₂CSO).

Probing time-dependent molecular dipoles on the attosecond time scale

Photoinduced molecular processes start with the interaction of the instantaneous electric field of the incident light with the electronic degrees of freedom. This early attosecond electronic motion impacts the fate of the photoinduced reactions. We report the first observation of attosecond time scale electron dynamics in a series of small- and medium-sized neutral molecules (N(2), CO(2), and C(2)H(4)), monitoring time-dependent variations of the parent molecular ion yield in the ionization by an attosecond pulse, and thereby probing the time-dependent dipole induced by a moderately strong near-infrared laser field. This approach can be generalized to other molecular species and may be regarded as a first example of molecular attosecond Stark spectroscopy.

Capturing Ultrafast Quantum Dynamics with Femtosecond and Attosecond X-ray Core-Level Absorption Spectroscopy

Recent technical advances in ultrafast laser sources enable the generation of femtosecond and attosecond soft X-ray pulses in tabletop laser setups as well as accelerator-based synchrotron and free-electron laser sources. These new light sources can be harnessed via pump-probe spectroscopy to elucidate ultrafast quantum dynamics in atoms, molecules, and condensed matter with unprecedented time resolution and chemical sensitivity. Employing such ultrashort pulses in transient X-ray absorption spectroscopy combines the unique advantages of core-level absorption probing of chemical environments and oxidation states with the ability to obtain ultimately freeze-frame snapshots of electronic and nuclear dynamics. In this Perspective, we provide an overview of the progress in applying the recently developed technique of femtosecond to attosecond time-resolved soft X-ray transient absorption spectroscopy to the study of ultrafast phenomena, including some of our own efforts to elucidate the interaction of intense laser pulses with atoms and molecules in the strong-field, nonperturbative limit. Possible avenues for future work are outlined.