Charge migration engineered by localisation: electron-nuclear dynamics in polyenes and glycine

Coherent Excitation of the CH Stretching Vibrations in C2H4 +: The Role of the Derivative Coupling Studied by the Quantum Ehrenfest Method

We report nonadiabatic dynamics computations on C2H4 + initiated on a coherent superposition of the five lowest cationic states, employing the Quantum Ehrenfest method. In addition to the totally symmetric carbon-carbon double bond stretch and carbon-hydrogen stretches, we see that the three non-totally symmetric modes become stimulated; torsion and three different CH stretching patterns. Thus, a coherent superposition of states, of the type involved in an attochemistry experiment, leads to the stimulation of specific non-totally symmetric motions. The computations were also performed on the specific combination of the A and C states. In each case normal mode 9 (cis-asymmetric H2CCH2 stretch), out of the set of non-totally-symmetric normal modes, dominates. Thus, we can steer the nuclear motion along specific non-totally symmetric normal modes using a defined coherent superposition.

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.

Attochemistry Regulation of Charge Migration

Charge migration (CM) is a coherent attosecond process that involves the movement of localized holes across a molecule. To determine the relationship between a molecule's structure and the CM dynamics it exhibits, we perform systematic studies of para-functionalized bromobenzene molecules (X-C₆H₄-R) using real-time time-dependent density functional theory. We initiate valence-electron dynamics by emulating rapid strong-field ionization leading to a localized hole on the bromine atom. The resulting CM, which takes on the order of 1 fs, occurs via an X localized → C₆H₄ delocalized → R localized mechanism. Interestingly, the hole contrast on the acceptor functional group increases with increasing electron-donating strength. This trend is well-described by the Hammett σ value of the group, which is a commonly used metric for quantifying the effect of functionalization on the chemical reactivity of benzene derivatives. These results suggest that simple attochemistry principles and a density-based picture can be used to predict and understand CM.

Quantum coherence in molecular photoionization

The study of onset and decay, as well as control of ultrafast quantum coherence in many-electron systems is in the focus of interest of attosecond physics. Interpretation of attosecond experiments detecting the ultrafast quantum coherence requires application of advanced theoretical and computational tools combining many-electron theory, description of the electronic continuum, including in the strong laser field scenario, as well as nuclear dynamics theory. This perspective reviews the recent theoretical advances in understanding the attosecond dynamics of quantum coherence in photoionized molecular systems and outlines possible future directions of theoretical and experimental study of coherence and entanglement in the attosecond regime.

How electronic superpositions drive nuclear motion following the creation of a localized hole in the glycine radical cation

In this work, we have studied the nuclear and electron dynamics in the glycine cation starting from localized hole states using the quantum Ehrenfest method. The nuclear dynamics is controlled both by the initial gradient and by the instantaneous gradient that results from the oscillatory electron dynamics (charge migration). We have used the Fourier transform (FT) of the spin densities to identify the “normal modes” of the electron dynamics. We observe an isomorphic relationship between the electron dynamics normal modes and the nuclear dynamics, seen in the vibrational normal modes. The FT spectra obtained this way show bands that are characteristic of the energy differences between the adiabatic hole states. These bands contain individual peaks that are in one-to-one correspondence with atom pair (+·) ↔ (·+) resonances, which, in turn, stimulate nuclear motion involving the atom pair. With such understanding, we anticipate “designer” coherent superpositions that can drive nuclear motion in a particular direction.

Quantum electronic coherences by attosecond transient absorption spectroscopy: ab initio B-spline RCS-ADC study

Here I present a fully ab initio time-resolved study of X-ray attosecond transient absorption spectroscopy (ATAS) in a prototypical polyatomic molecule, pyrazine, and demonstrate the possibility of retrieving the many-electron quantum ionic coherences arising in attosecond molecular photoionisation and pre-determining the subsequent charge-directed photochemical reactivity. Advanced first-principles many-electron simulations are performed, within a hybrid XUV pump/X-ray probe setup, to describe the interaction of pyrazine with both XUV pump and X-ray probe pulses, and study the triggered correlated many-electron dynamics. The calculations are carried out by means of the recently-developed ab initio method for many-electron dynamics in polyatomic molecules, the time-dependent (TD) B-spline Restricted Correlation Space-Algebraic Diagrammatic Construction (RCS-ADC). RCS-ADC simulates molecular ionisation from first principles, combining the accurate description of electron correlation of quantum chemistry with the full account of the continuum dynamics of the photoelectron. Complete theoretical characterisation of the atto-ionised many-electron state and photo-induced attosecond charge dynamics is achieved by calculating the reduced ionic density matrix (R-IDM) of the bipartite ion-photoelectron system, with full inclusion of the correlated shakeup states. Deviations from the sudden approximation picture of photoionisation, in the low-photon-energy limit, are presented. The effect of the multi-channel interaction between the parent-ion and the emitted photoelectron on the onset of the quantum electronic coherences is analysed. Moreover, I show how the Schmidt decomposition of the R-IDM unravels the many-electron dynamics triggered by the pump, allowing for the identification of the key channels involved. Finally, I calculate the X-ray attosecond transient absorption spectra of XUV-ionised pyrazine. The results unveil the mapping of the ATAS measurement onto the quantum electronic coherences, and related non-zero R-IDM matrix elements, produced upon ionisation by the XUV pump laser pulse.

Real-time observation of a correlation-driven sub 3 fs charge migration in ionised adenine

Sudden ionisation of a relatively large molecule can initiate a correlation-driven process dubbed charge migration, where the electron density distribution is expected to rapidly move along the molecular backbone. Capturing this few-femtosecond or attosecond charge redistribution would represent the real-time observation of electron correlation in a molecule with the enticing prospect of following the energy flow from a single excited electron to the other coupled electrons in the system. Here, we report a time-resolved study of the correlation-driven charge migration process occurring in the nucleic-acid base adenine after ionisation with a 15-35 eV attosecond pulse. We find that the production of intact doubly charged adenine - via a shortly-delayed laser-induced second ionisation event - represents the signature of a charge inflation mechanism resulting from many-body excitation. This conclusion is supported by first-principles time-dependent simulations. These findings may contribute to the control of molecular reactivity at the electronic, few-femtosecond time scale.

First-Principles Simulations of Biological Molecules Subjected to Ionizing Radiation

Ionizing rays cause damage to genomes, proteins, and signaling pathways that normally regulate cell activity, with harmful consequences such as accelerated aging, tumors, and cancers but also with beneficial effects in the context of radiotherapies. While the great pace of research in the twentieth century led to the identification of the molecular mechanisms for chemical lesions on the building blocks of biomacromolecules, the last two decades have brought renewed questions, for example, regarding the formation of clustered damage or the rich chemistry involving the secondary electrons produced by radiolysis. Radiation chemistry is now meeting attosecond science, providing extraordinary opportunities to unravel the very first stages of biological matter radiolysis. This review provides an overview of the recent progress made in this direction, focusing mainly on the atto- to femto- to picosecond timescales. We review promising applications of time-dependent density functional theory in this context.

Competition between charge migration and charge transfer induced by nuclear motion following core ionization: Model systems and application to Li2. J Chem Phys 2019;151:124108. [PMID: 31575213 DOI: 10.1063/1.5117246]

Attosecond and femtosecond spectroscopies present opportunities for the control of chemical reaction dynamics and products, as well as for quantum information processing; we address the somewhat unique situation of core-ionization spectroscopy which, for dimeric chromophores, leads to strong valence charge localization and hence tightly paired potential-energy surfaces of very similar shape. Application is made to the quantum dynamics of core-ionized Li₂ ⁺. This system is chosen as Li₂ is the simplest stable molecule facilitating both core ionization and valence ionization. First, the quantum dynamics of some model surfaces are considered, with the surprising result that subtle differences in shape between core-ionization paired surfaces can lead to dramatic differences in the interplay between electronic charge migration and charge transfer induced by nuclear motion. Then, equation-of-motion coupled-cluster calculations are applied to determine potential-energy surfaces for 8 core-excited state pairs, calculations believed to be the first of their type for other than the lowest-energy core-ionized molecular pair. While known results for the lowest-energy pair suggest that Li₂ ⁺ is unsuitable for studying charge migration, higher-energy pairs are predicted to yield results showing competition between charge migration and charge transfer. Central is a focus on the application of Hush's 1975 theory for core-ionized X-ray photoelectron spectroscopy to understand the shapes of the potential-energy surfaces and hence predict key features of charge migration.

Onset of ionic coherence and ultrafast charge dynamics in attosecond molecular ionisation

Here is presented a fully ab initio theoretical framework for simulating the correlated many-electron dynamics occurring during and emerging from molecular ionisation by attosecond laser pulses. This is based on the time-dependent (TD) version of the B-spline restricted correlation space (RCS)-algebraic diagrammatic construction (ADC) method, with the full description of the photoelectron and inclusion of electron correlation effects, such as shakeup processes and inter-channel couplings. The nature of the ultrafast charge dynamics in the molecular ion is elucidated by quantitatively predicting the degree of electronic coherence and eigenstate content of the prepared molecular cationic state, beyond the commonly used sudden approximation. The results presented here for the acetylene and ethylene molecules show that even in the high photon energy regime the simulated hole dynamics is quantitatively different from the prediction of the sudden approximation. Moreover, for high-bandwidth ionising pulse, the residual interaction between the cation, in highly-excited shake-up states, and the emitted slow photoelectron gives rise to a loss of coherence in the ionic system which can persist for the first few femtoseconds after ionisation.

The Ehrenfest method with fully quantum nuclear motion (Qu-Eh): Application to charge migration in radical cations

An algorithm is described for quantum dynamics where an Ehrenfest potential is combined with fully quantum nuclear motion (Quantum-Ehrenfest, Qu-Eh). The method is related to the single-set variational multi-configuration Gaussian approach (vMCG) but has the advantage that only a single quantum chemistry computation is required at each time step since there is only a single time-dependent potential surface. Also shown is the close relationship to the "exact factorization method." The quantum Ehrenfest method is compared with vMCG for study of electron dynamics in a modified bismethylene-adamantane cation system. Illustrative examples of electron-nuclear dynamics are presented for a distorted allene system and for HCCI⁺ where one has a degenerate Π system.