Ultrafast charge migration by electron correlation

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 Real-Time Molecular Electron Dynamics Efficiently Using the Time-Dependent Density Matrix Renormalization Group

Compared to ground-state electronic structure optimizations, accurate simulations of molecular real-time electron dynamics are usually much more difficult to perform. To simulate electron dynamics, the time-dependent density matrix renormalization group (TDDMRG) has been shown to offer an attractive compromise between accuracy and cost. However, many simulation parameters significantly affect the quality and efficiency of a TDDMRG simulation. So far, it is unclear whether common wisdom from ground-state DMRG carries over to the TDDMRG, and a guideline on how to choose these parameters is missing. Here, in order to establish such a guideline, we investigate the convergence behavior of the main TDDMRG simulation parameters, such as time integrator, the choice of orbitals, and the choice of matrix-product-state representation for complex-valued nonsinglet states. In addition, we propose a method to select orbitals that are tailored to optimize the dynamics. Lastly, we showcase the TDDMRG by applying it to charge migration ionization dynamics in furfural, where we reveal a rapid conversion from an ionized state with a σ character to one with a π character within less than a femtosecond.

Ab-initio variational wave functions for the time-dependent many-electron Schrödinger equation

Understanding the real-time evolution of many-electron quantum systems is essential for studying dynamical properties in condensed matter, quantum chemistry, and complex materials, yet it poses a significant theoretical and computational challenge. Our work introduces a variational approach for fermionic time-dependent wave functions, surpassing mean-field approximations by accurately capturing many-body correlations. We employ time-dependent Jastrow factors and backflow transformations, enhanced through neural networks parameterizations. To compute the optimal time-dependent parameters, we employ the time-dependent variational Monte Carlo technique and introduce a new method based on Taylor-root expansions of the propagator, enhancing the accuracy of our simulations. The approach is demonstrated in three distinct systems. In all cases, we show clear signatures of many-body correlations in the dynamics. The results showcase the ability of our variational approach to accurately describe the time evolution, providing insight into quantum dynamical effects in interacting electronic systems, beyond the capabilities of mean-field.

Cavity Manipulation of Attosecond Charge Migration in Conjugated Dendrimers

Dendrimers are branched polymers with wide applications to photosensitization, photocatalysis, photodynamic therapy, photovoltaic conversion, and light sensor amplification. The primary step of numerous photophysical and photochemical processes in many molecules involves ultrafast coherent electronic dynamics and charge oscillations triggered by photoexcitation. This electronic wavepacket motion at short times where the nuclei are frozen is known as attosecond charge migration. We show how charge migration in a dendrimer can be manipulated by placing it in an optical cavity and monitored by time-resolved X-ray diffraction. Our simulations demonstrate that the dendrimer charge migration modes and the character of photoexcited wave function can be significantly influenced by the strong light-matter interaction in the cavity. This presents a new avenue for modulating initial ultrafast charge dynamics and subsequently controlling coherent energy transfer in dendritic nanostructures.

Self-Heterodyne Diffractive Imaging of Ultrafast Electron Dynamics Monitored by Single-Electron Pulses

The direct imaging of time-evolving molecular charge densities on atomistic scale and at femtosecond resolution has long been an elusive task. In this theoretical study, we propose a self-heterodyne electron diffraction technique based on single electron pulses. The electron is split into two beams, one passes through the sample and its interference with the second beam produces a heterodyne diffraction signal that images the charge density. Application to probing the ultrafast electronic dynamics in Mg-phthalocyanine demonstrates its potential for imaging chemical dynamics.

Core-Hole Coherent Spectroscopy in Molecules

We study the ultrafast dynamics initiated by a coherent superposition of core-excited states of nitrous oxide molecule. Using high-level ab initio methods, we show that the decoherence caused by the electronic decay and the nuclear dynamics is substantially slower than the induced ultrafast quantum beatings, allowing the system to undergo several oscillations before it dephases. We propose a proof-of-concept experiment using the harmonic up-conversion scheme available at x-ray free-electron laser facilities to trace the evolution of the created core-excited-state coherence through a time-resolved x-ray photoelectron spectroscopy.

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.

Attosecond Probing of Nuclear Vibrations in the D2+ and HeH+ Molecular Ions

We study the ultrafast photodissociation of small diatomic molecules using attosecond laser pulses of moderate intensity in the (extreme) ultraviolet regime. The simultaneous application of subfemtosecond laser pulses with different photon energies─resonant in the region of the molecular motion─allows one to monitor the vibrational dynamics of simple diatomics, like the D2+ and HeH+ molecular ions. In our real-time wave packet simulations, the nuclear dynamics is initiated either by sudden ionization (D2+) or by explicit pump pulses (HeH+) via distortion of the potential energy of the molecule. The application of time-delayed attosecond pulses leads to the breakup of the molecules, and the information on the underlying bound-state dynamics is imprinted in the kinetic energy release (KER) spectra of the outgoing fragments. We show that the KER-delay spectrograms generated in our ultrafast pump-probe schemes are able to reconstruct the most important features of the molecular motion within a given electronic state, such as the time period or amplitude of oscillations, interference patterns, or the revival and splitting of the nuclear wave packet. The impact of probe pulse duration, which is key to the applicability of the presented mapping scheme, is investigated in detail.

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump-probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron-nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. This work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

Tracking Charge Migration with Frequency-Matched Strobo-Spectroscopy

We present frequency-matched strobo-spectroscopy (FMSS) of charge migration (CM) in bromobutadiyne, simulated with time-dependent density functional theory. CM + FMSS is a pump-probe scheme that uses a frequency-matched high harmonic generation (HHG)-driving laser as an independent probe step, following the creation of a localized hole on the bromine atom that induces CM dynamics. We show that the delay-dependent harmonic yield tracks the phase of the CM dynamics through its sensitivity to the amount of electron density on the bromine end of the molecule. FMSS takes advantage of the intrinsic attosecond time resolution of the HHG process in which different harmonics are emitted at different times and thus probe different locations of the electron hole. Finally, we show that the CM-induced modulation of the HHG signal is dominated by the recombination step of the HHG process, with a negligible contribution from the ionization step.

Size Effect in Correlation-Driven Charge Migration in Correlation Bands of Alkyne Chains

Correlation-driven charge migration initiated by inner-valence ionization leading to the population of the correlation bands of alkyne chains containing between 4 and 12 carbon atoms is explored through ab initio simulations. Scaling laws are observed, both for the time scale of the charge migration and for the slope of the density of states of the correlation bands. These can be used for predicting the relaxation time scale in much larger systems from the same molecular family and for finding promising candidates for the development of an attochemistry scheme taking advantages of the specificity of the dynamics in the correlation bands of molecules.

Approximate Exponential Integrators for Time-Dependent Equation-of-Motion Coupled Cluster Theory

With a growing demand for time-domain simulations of correlated many-body systems, the development of efficient and stable integration schemes for the time-dependent Schrödinger equation is of keen interest in modern electronic structure theory. In this work, we present two approaches for the formation of the quantum propagator for time-dependent equation-of-motion coupled cluster theory based on the Chebyshev and Arnoldi expansions of the complex, nonhermitian matrix exponential, respectively. The proposed algorithms are compared with the short-iterative Lanczos method of Cooper et al. [J. Phys. Chem. A 2021 125, 5438-5447], the fourth-order Runge-Kutta method, and exact dynamics for a set of small but challenging test problems. For each of the cases studied, both of the proposed integration schemes demonstrate superior accuracy and efficiency relative to the reference simulations.

Direct temporal characterization of sub-3-fs deep UV pulses generated by resonant dispersive wave emission

We report on the complete temporal characterization of ultrashort pulses, generated by resonant dispersive wave emission in gas-filled hollow-capillary fibers, with energy in the microjoule range and continuously tunable from the deep-ultraviolet to the ultraviolet. Temporal characterization of such ultrabroad pulses, particularly challenging in this spectral region, was performed using an all-in-vacuum setup for self-diffraction frequency resolved optical gating (SD-FROG). Sub-3-fs pulses were measured, tunable from 250 nm to 350 nm, with a minimum pulse duration of 2.4 ± 0.1 fs.

Symmetry breaking in core-valence double ionisation of allene

Conventional electron spectroscopy is an established one-electron-at-the-time method for revealing the electronic structure and dynamics of either valence or inner shell ionized systems. By combining an electron-electron coincidence technique with the use of soft X-radiation we have measured a double ionisation spectrum of the allene molecule in which one electron is removed from a C1s core orbital and one from a valence orbital, well beyond Siegbahns Electron-Spectroscopy-for-Chemical-Analysis method. This core-valence double ionisation spectrum shows the effect of symmetry breaking in an extraordinary way, when the core electron is ejected from one of the two outer carbon atoms. To explain the spectrum we present a new theoretical approach combining the benefits of a full self-consistent field approach with those of perturbation methods and multi-configurational techniques, thus establishing a powerful tool to reveal molecular orbital symmetry breaking on such an organic molecule, going beyond Löwdins standard definition of electron correlation.

Manipulating Attosecond Charge Migration in Molecules by Optical Cavities

The ultrafast electronic charge dynamics in molecules upon photoionization while the nuclear motions are frozen is known as charge migration. In a theoretical study of the quantum dynamics of photoionized 5-bromo-1-pentene, we show that the charge migration process can be induced and enhanced by placing the molecule in an optical cavity, and can be monitored by time-resolved photoelectron spectroscopy. The collective nature of the polaritonic charge migration process is investigated. We find that, unlike spectroscopy, molecular charge dynamics in a cavity is local and does not show many-molecule collective effects. The same conclusion applies to cavity polaritonic chemistry.

Coherent Mixing of Singlet and Triplet States in Acrolein and Ketene: A Computational Strategy for Simulating the Electron-Nuclear Dynamics of Intersystem Crossing

We present a theoretical study of intersystem crossing (ISC) in acrolein and ketene with the Ehrenfest method that can describe a superposition of singlet and triplet states. Our simulations illustrate a new mechanistic effect of ISC, namely, that a superposition of singlets and triplets yields nonadiabatic dynamics characteristic of that superposition rather than the constituent state potential energy surfaces. This effect is particularly significant in ketene, where mixing of singlet and triplet states along the approach to a singlet/singlet conical intersection occurs, with the spin-orbit coupling (SOC) remaining small throughout. In both cases, the effects require many recrossings of the singlet/triplet state crossing seam, consistent with the textbook treatment of ISC.

The role of symmetric vibrational modes in the decoherence of correlation-driven charge migration

Due to the electron correlation, the fast removal of an electron from a molecule may create a coherent superposition of cationic states and in this way initiate pure electronic dynamics in which the hole-charge left by the ionization migrates throughout the system on an ultrashort time scale. The coupling to the nuclear motion introduces a decoherence that eventually traps the charge, and crucial questions in the field of attochemistry include how long the electronic coherence lasts and which nuclear degrees of freedom are mostly responsible for the decoherence. Here, we report full-dimensional quantum calculations of the concerted electron-nuclear dynamics following outer-valence ionization of propynamide, which reveal that the pure electronic coherences last only 2-3 fs before being destroyed by the nuclear motion. Our analysis shows that the normal modes that are mostly responsible for the fast electronic decoherence are the symmetric in-plane modes. All other modes have little or no effect on the charge migration. This information can be useful to guide the development of reduced dimensionality models for larger systems or the search for molecules with long coherence times.

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.

Photo-ionization initiated differential ultrafast charge migration: impacts of molecular symmetries and tautomeric forms

Photo-ionization induced ultrafast electron dynamics is considered as a precursor for the slower nuclear dynamics associated with molecular dissociation. Here, using the ab initio multielectron wave-packet propagation method, we study the overall many-electron dynamics, triggered by ionizing the outer-valence orbitals of different tautomers for a prototype molecule with more than one symmetry element. From the time evolution of the initially created averaged hole density of each system, we identify distinctly different charge dynamics responses in the tautomers. We observe that the keto form shows a charge migration direction away from the nitrogen bonded with hydrogen, while in enol-U - away from oxygen bonded to hydrogen. Additionally, the dynamics following the ionization of molecular orbitals with different symmetries reveals that a' orbitals show a fast and highly delocalized charge density in comparison to a'' symmetry. These observations indicate why different tautomers respond differently to an XUV ionization, and might explain the subsequent different fragmentation pathways. An experimental schematics allowing the detection and reconstruction of such charge dynamics is also proposed. Although the present study uses a simple, prototypical bio-relevant molecule, it reveals the explicit role of molecular symmetry and tautomerism in the ionization-triggered charge migration that controls many ultrafast physical, chemical, and biological processes, making tautomeric forms a promising tool of molecular design for desired charge migration.

Probing Electronic Symmetry Reduction during Charge Migration via Time-Resolved X-Ray Diffraction

The present work focuses on probing ultrafast charge migration after symmetry-breaking excitation using ultrashort laser pulses. LiCN is chosen as prototypical system because it can be oriented in the laboratory frame and it possesses optically-accessible charge transfer states at low energies. The charge migration is simulated within the hybrid time-dependent density functional theory/configuration interaction framework. Time-resolved electronic current densities and simulated time-resolved x-ray diffraction signals are used to unravel the mechanism of charge migration. Our simulations demonstrate that specific choices of laser polarization lead to a control over the symmetry of the induced charge migration. Moreover, time-resolved x-ray diffraction signals are shown to encode transient symmetry reduction at intermediate times.

Ultrafast chirality: the road to efficient chiral measurements

Today we are witnessing the electric-dipole revolution in chiral measurements. Here we reflect on its lessons and outcomes, such as the perspective on chiral measurements using the complementary principles of "chiral reagent" and "chiral observer", the hierarchy of scalar, vectorial and tensorial enantio-sensitive observables, the new properties of the chiro-optical response in the ultrafast and non-linear domains, and the geometrical magnetism associated with the chiral response in photoionization. The electric-dipole revolution is a landmark event. It has opened routes to extremely efficient enantio-discrimination with a family of new methods. These methods are governed by the same principles but work in vastly different regimes - from microwaves to optical light; they address all molecular degrees of freedom - electronic, vibrational and rotational, and use flexible detection schemes, i.e. detecting photons or electrons, making them applicable to different chiral phases, from gases to liquids to amorphous solids. The electric-dipole revolution has also enabled enantio-sensitive manipulation of chiral molecules with light. This manipulation includes exciting and controlling ultrafast helical currents in vibronic states of chiral molecules, enantio-sensitive control of populations in electronic, vibronic and rotational molecular states, and opens the way to efficient enantio-separation and enantio-sensitive trapping of chiral molecules. The word "perspective" has two meanings: an "outlook" and a "point of view". In this perspective article, we have tried to cover both meanings.

Charge-induced chemical dynamics in glycine probed with time-resolved Auger electron spectroscopy

In the present contribution, we use x-rays to monitor charge-induced chemical dynamics in the photoionized amino acid glycine with femtosecond time resolution. The outgoing photoelectron leaves behind the cation in a coherent superposition of quantum mechanical eigenstates. Delayed x-ray pulses track the induced coherence through resonant x-ray absorption that induces Auger decay. Temporal modulation of the Auger electron signal correlated with specific ions is observed, which is governed by the initial electronic coherence and subsequent vibronic coupling to nuclear degrees of freedom. In the time-resolved x-ray absorption measurement, we monitor the time-frequency spectra of the resulting many-body quantum wave packets for a period of 175 fs along different reaction coordinates. Our experiment proves that by measuring specific fragments associated with the glycine dication as a function of the pump-probe delay, one can selectively probe electronic coherences at early times associated with a few distinguishable components of the broad electronic wave packet created initially by the pump pulse in the cation. The corresponding coherent superpositions formed by subsets of electronic eigenstates and evolving along parallel dynamical pathways show different phases and time periods in the range of ( - 0.3 ± 0.1 ) π ≤ ϕ ≤ ( 0.1 ± 0.2 ) π and 18.2 - 1.4 + 1.7 ≤ T ≤ 23.9 - 1.1 + 1.2 fs. Furthermore, for long delays, the data allow us to pinpoint the driving vibrational modes of chemical dynamics mediating charge-induced bond cleavage along different reaction coordinates.

Probing Resonant Photoionization Time Delay by Self-Referenced Molecular Attoclock

Attosecond time-resolved electron tunneling dynamics have been investigated by using attosecond angular streaking spectroscopy, where a clock reference to the laser field vector is required in atomic strong-field ionization and the situation becomes complicated in molecules. Here we reveal a resonant ionization process via a transient state by developing an electron-tunneling-site-resolved molecular attoclock in Ar-Kr^{+}. Two distinct deflection angles are observed in the photoelectron angular distribution in the molecular frame, corresponding to the direct and resonant ionization pathways. We find the electron is temporally trapped in the Coulomb potential wells of the Ar-Kr^{+} before finally releasing into the continuum when the electron tunnels through the internal barrier. By utilizing the direct tunneling ionization as a self-referenced arm of the attoclock, the time delay of the electron trapped in the resonant state is revealed to be 3.50±0.04  fs. Our results give an impetus to exploring the ultrafast electron dynamics in complex systems and also endow a semiclassical presentation of the electron trapping dynamics in a quantum resonant state.

Quantum Interference Paves the Way for Long-Lived Electronic Coherences

The creation and dynamical fate of a coherent superposition of electronic states generated in a polyatomic molecule by broadband ionization with extreme ultraviolet pulses is studied using the multiconfiguration time-dependent Hartree method together with an ionization continuum model Hamiltonian. The electronic coherence between the hole states usually lasts until the nuclear dynamics leads to decoherence. A key goal of attosecond science is to control the electronic motion and design laser control schemes to retain this coherence for longer timescales. Here, we investigate this possibility using time-delayed pulses and show how this opens up the prospect of coherent control of charge migration phenomenon.

Decoherence and Revival in Attosecond Charge Migration Driven by Non-adiabatic Dynamics

Attosecond charge migration is a periodic evolution of the charge density at specific sites of a molecule on a time scale defined by the energy intervals between the electronic states involved. Here, we report the observation of charge migration in neutral silane (SiH₄) in 690 as, its decoherence within 15 fs, and its revival after 40-50 fs, using X-ray attosecond transient absorption spectroscopy. We observe the migration of charge as pairs of quantum beats with a characteristic spectral phase in the transient spectrum, in agreement with theory. The decay and revival of the degree of electronic coherence is found to be a result of both adiabatic and non-adiabatic dynamics in the populated Rydberg and valence states. The experimental results are supported by fully quantum-mechanical ab-initio calculations that include both electronic and nuclear dynamics, which additionally support the experimental evidence that conical intersections can mediate the transfer of electronic coherence from an initial superposition state to another one involving a different lower-lying state.

Filming movies of attosecond charge migration in single molecules with high harmonic spectroscopy

Electron migration in molecules is the progenitor of chemical reactions and biological functions after light-matter interaction. Following this ultrafast dynamics, however, has been an enduring endeavor. Here we demonstrate that, by using machine learning algorithm to analyze high-order harmonics generated by two-color laser pulses, we are able to retrieve the complex amplitudes and phases of harmonics of single fixed-in-space molecules. These complex dipoles enable us to construct movies of laser-driven electron migration after tunnel ionization of N₂ and CO₂ molecules at time steps of 50 attoseconds. Moreover, the angular dependence of the migration dynamics is fully resolved. By examining the movies, we observe that electron holes do not just migrate along the laser polarization direction, but may swirl around the atom centers. Our result establishes a general scheme for studying ultrafast electron dynamics in molecules, paving a way for further advance in tracing and controlling photochemical reactions by femtosecond lasers.

Light-Induced Ultrafast Molecular Dynamics: From Photochemistry to Optochemistry

By precisely controlling the waveform of ultrashort laser fields, electronic and nuclear motions in molecules can be steered on extremely short time scales, even in the attosecond regime. This new research field, termed "optochemistry", presents the light field in the time-frequency domain and opens new avenues for tailoring molecular reactions beyond photochemistry. This Perspective summarizes the ultrafast laser techniques employed in recent years for manipulating the molecular reactions based on waveform control of intense ultrashort laser pulses, where the chemical reactions can take place in isolated molecules, clusters, and various nanosystems. The underlying mechanisms for the coherent control of molecular dynamics are explicitly explored. Challenges and opportunities coexist in the field of optochemistry. Advanced technologies and theoretical modeling are still being pursued, with great prospects for controlling chemical reactions with unprecedented spatiotemporal precision.

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.

Electronic quantum coherence in glycine molecules probed with ultrashort x-ray pulses in real time

Here, we use x-rays to create and probe quantum coherence in the photoionized amino acid glycine. The outgoing photoelectron leaves behind the cation in a coherent superposition of quantum mechanical eigenstates. Delayed x-ray pulses track the induced coherence through resonant x-ray absorption that induces Auger decay and by photoelectron emission from sequential double photoionization. Sinusoidal temporal modulation of the detected signal at early times (0 to 25 fs) is observed in both measurements. Advanced ab initio many-electron simulations allow us to explain the first 25 fs of the detected coherent quantum evolution in terms of the electronic coherence. In the kinematically complete x-ray absorption measurement, we monitor its dynamics for a period of 175 fs and observe an evolving modulation that may implicate the coupling of electronic to vibronic coherence at longer time scales. Our experiment provides a direct support for the existence of long-lived electronic coherence in photoionized biomolecules.

High harmonic spectra computed using time-dependent Kohn-Sham theory with Gaussian orbitals and a complex absorbing potential

High harmonic spectra for H₂ and H₂ ⁺ are simulated by solving the time-dependent Kohn-Sham equation in the presence of a strong laser field using an atom-centered Gaussian representation of the density and a complex absorbing potential. The latter serves to mitigate artifacts associated with the finite extent of the basis functions, including spurious reflection of the outgoing electronic wave packet. Interference between the outgoing and reflected waves manifests as peak broadening in the spectrum as well as the appearance of spurious high-energy peaks after the harmonic progression has terminated. We demonstrate that well-resolved spectra can be obtained through the use of an atom-centered absorbing potential. As compared to grid-based algorithms, the present approach is more readily extensible to larger molecules.

Attosecond spectroscopy for the investigation of ultrafast dynamics in atomic, molecular and solid-state physics

Since the first demonstration of the generation of attosecond pulses (1 as = 10^-18s) in the extreme-ultraviolet spectral region, several measurement techniques have been introduced, at the beginning for the temporal characterization of the pulses, and immediately after for the investigation of electronic and nuclear ultrafast dynamics in atoms, molecules and solids with unprecedented temporal resolution. The attosecond spectroscopic tools established in the last two decades, together with the development of sophisticated theoretical methods for the interpretation of the experimental outcomes, allowed to unravel and investigate physical processes never observed before, such as the delay in photoemission from atoms and solids, the motion of electrons in molecules after prompt ionization which precede any notable nuclear motion, the temporal evolution of the tunneling process in dielectrics, and many others. This review focused on applications of attosecond techniques to the investigation of ultrafast processes in atoms, molecules and solids. Thanks to the introduction and ongoing developments of new spectroscopic techniques, the attosecond science is rapidly moving towards the investigation, understanding and control of coupled electron-nuclear dynamics in increasingly complex systems, with ever more accurate and complete investigation techniques. Here we will review the most common techniques presenting the latest results in atoms, molecules and solids.

Systematic Investigation of the Reliability of the Frozen Nuclei Approximation for Short-Pulse Excitation: The Example of HCCI+

In this work we quantitatively study the reliability of the frozen nuclei approximation for ultrafast dynamics. Specifically we study laser excitation of HCCI⁺ from its ground state to the first electronically excited state. The population of the first excited state is obtained by both the frozen nuclei approximation and by multidimensional nuclear dynamics. Detailed comparison of the results by the two methods are performed to provide quantitative criteria for the reliability of the frozen nuclei approximation for this system.

Attosecond coherent electron motion in Auger-Meitner decay

In quantum systems, coherent superpositions of electronic states evolve on ultrafast time scales (few femtoseconds to attoseconds; 1 attosecond = 0.001 femtoseconds = 10^-18 seconds), leading to a time-dependent charge density. Here we performed time-resolved measurements using attosecond soft x-ray pulses produced by a free-electron laser, to track the evolution of a coherent core-hole excitation in nitric oxide. Using an additional circularly polarized infrared laser pulse, we created a clock to time-resolve the electron dynamics and demonstrated control of the coherent electron motion by tuning the photon energy of the x-ray pulse. Core-excited states offer a fundamental test bed for studying coherent electron dynamics in highly excited and strongly correlated matter.

Coherent ring-current migration of Mg-phthalocyanine probed by time-resolved X-ray circular dichroism

The coherent ring current of Mg-phthalocyanine created by a broad band UV-visible pump pulse shows variation with time, where the ring currents at the corner benzene rings, around the Mg cation and on the outer ring oscillate with different time periods and the current density migrates among these regions. The 7 pairs of E_u degenerate excited states populated upon photoexcitation, generate 21 distinct coherent ring currents. We further calculate the time-resolved X-ray circular dichroism (TRXCD) spectrum of the coherences contributing to the ring current obtained by an attosecond X-ray probe pulse resonant with the nitrogen K-edge. A frequency domain TRXCD signal obtained by a Fourier transform of the signal with respect to the pump-probe delay time clearly separates the currents induced by different state pairs.

The coherent ring current of Mg-phthalocyanine are created by a broad band UV-visible pump pulse and migrate into different regions within the molecule. This coherent ring current dynamics is probed by time-resolved X-ray circular dichroism.

Attosecond charge migration following oxygen K-shell ionization in DNA bases and base pairs

Core ionization of DNA begins a cascade of events which could lead to cellular inactivation or death. The created core-hole following an impulse inner-shell ionization of molecules naturally decays in the auger timescale. We simulated charge migration (CM) phenomena following an impulsive core ionization of individual DNA bases at the oxygen K-edge which occurs before Auger decay of the oxygen. Our approach is based on real-time time dependent density functional theory (RT-TDDFT). It is shown that the pronounced hole fluctuation observed around bonds of the initial core-hole results in various valence orbital migrations. Also, the same photo-core-ionized dynamics is studied for the related base pairs. We investigate the role of base pairing and H-bonding interactions in the attosecond CM dynamics. In particular, the creation of a core-hole in the oxygen involved in H-bonding leads to an enhancement of charge migration relative to the respective single bases. Importantly, the hole oscillation of the adenine-thymine base pair upon creation of a core-hole at the oxygen, which does not contribute to the donor-acceptor interactions (not H-bonded), decreases compared to the single thymine base. Understanding the detailed dynamics of the localized core-hole initiating CM process would open the way for chemically controlling DNA damage/repair in the future.

Ultrafast Valence-Electron Dynamics in Oxazole Monitored by X-ray Diffraction Following a Stimulated X-ray Raman Excitation

Direct imaging of the ultrafast quantum motion of valence electrons in molecules is essential for understanding many elementary chemical and physical processes. We present a simulation study of valence-electron dynamics of oxazole. A valence-state electronic wavepacket is prepared with an attosecond soft X-ray pulse through a stimulated resonant X-ray Raman process and then probed with time-resolved off-resonant single-molecule X-ray diffraction. We find that the time dependent diffraction signal originates solely from the electronic coherences and can be detected by existing experimental techniques. We thus provide a feasible way of imaging electron dynamics in molecules. Moreover, the created electronic coherences and subsequent electron dynamics can be manipulated by the resonant X-ray Raman excitation tuned to different core-excited states.

The mystery of sub-picosecond charge transfer following irradiation of hydrated uridine monophosphate

The early mechanisms by which ionizing rays damage biological structures by so-called direct effects are largely elusive. In a recent picosecond pulse radiolysis study of concentrated uridine monophosphate solutions [J. Ma, S. A. Denisov, J.-L. Marignier, P. Pernot, A. Adhikary, S. Seki and M. Mostafavi, J. Phys. Chem. Lett., 2018, 9, 5105], unexpected results were found regarding the oxidation of the nucleobase. The signature of the oxidized nucleobase could not be detected 5 ps after the electron pulse, but only the oxidized phosphate, raising intriguing questions about the identity of charge-transfer mechanisms that could explain the absence of U⁺. We address here this question by means of advanced first-principles atomistic simulations of solvated uridine monophosphate, combining Density Functional Theory (DFT) with polarizable embedding schemes. We contrast three very distinct mechanisms of charge transfer covering the atto-, femto- and pico-second timescales. We first investigate the ionization mechanism and subsequent hole/charge migrations on a timescale of attoseconds to a few femtoseconds under the frozen nuclei approximation. We then consider a nuclear-driven phosphate-to-oxidized-nucleobase electron transfer, showing that it is an uncompetitive reaction channel on the sub-picosecond timescale, despite its high exothermicity and significant electronic coupling. Finally, we show that non-adiabatic charge transfer is enabled by femtosecond nuclear relaxation after ionization. We show that electronic decoherence and the electronic coupling strength are the key parameters that determine the hopping probabilities. Our results provide important insight into the interplay between electronics and nuclear motions in the early stages of the multiscale responses of biological matter subjected to ionizing radiation.

Ultrafast attosecond-magnetic-field generation of the charge migration process based on HeH2+ and H2 + electronically excited by circularly polarized laser pulses

The ultrafast process by the electron in molecular ions from one site or region to another that has come to be known as charge migration (CM), which is of fundamental importance to photon induced chemical or physical reactions. In this work, we study the electron current and ultrafast magnetic-field generation based on CM process of oriented asymmetric (HeH²⁺) and symmetric (H₂ ⁺) molecular ions. Calculated results show that they are ascribed to quantum interference of electronic states for these molecular ions under intense circularly polarized (CP) laser pulses. The two scenarios of (i) resonance excitation and (ii) direct ionization are considered through appropriately utilizing designed laser pulses. By comparison, the magnetic field induced by the scenario (i) is stronger than that of scenario (ii) for molecular ions. However, the scheme (ii) is very sensitive to the helicity of CP field, which is opposite to the scenario (i). Moreover, the magnetic field generated by H₂ ⁺ is stronger than that by HeH²⁺ through scenario (i). Our findings provide a guiding principle for producing ultrafast magnetic fields in molecular systems for future research in ultrafast magneto-optics.

Core-Valence Attosecond Transient Absorption Spectroscopy of Polyatomic Molecules

Tracing ultrafast processes induced by interaction of light with matter is often very challenging. In molecular systems, the initially created electronic coherence becomes damped by the slow nuclear rearrangement on a femtosecond timescale which makes real-time observations of electron dynamics in molecules particularly difficult. In this work, we report an extension of the theory underlying the attosecond transient absorption spectroscopy (ATAS) for the case of molecules, including a full account for the coupled electron-nuclear dynamics in the initially created wave packet, and apply it to probe the oscillations of the positive charge created after outer-valence ionization of the propiolic acid molecule. By taking advantage of element-specific core-to-valence transitions induced by x-ray radiation, we show that the resolution of ATAS makes it possible to trace the dynamics of electron density with atomic spatial resolution.

Attochemistry: Is Controlling Electrons the Future of Photochemistry?

Controlling matter with light has always been a great challenge, leading to the ever-expanding field of photochemistry. In addition, since the first generation of light pulses of attosecond (1 as = 10^-18 s) duration, a great deal of effort has been devoted to observing and controlling electrons on their intrinsic time scale. Because of their short duration, attosecond pulses have a large spectral bandwidth populating several electronically excited states in a coherent manner, i.e., an electronic wavepacket. Because of interference, such a wavepacket has a new electronic distribution implying a potentially different and totally new reactivity as compared to traditional photochemistry, leading to the novel concept of "attochemistry". This nascent field requires the support of theory right from the start. In this Perspective, we discuss the opportunities offered by attochemistry, the related challenges, and the current and future state-of-the-art developments in theoretical chemistry needed to model it accurately.

Short Iterative Lanczos Integration in Time-Dependent Equation-of-Motion Coupled-Cluster Theory

A time-dependent (TD) formulation of equation-of-motion coupled-cluster (EOM-CC) theory can provide excited-state information over an arbitrarily wide energy window with a reduced memory footprint relative to conventional, frequency-domain EOM-CC theory. However, the floating-point costs of the time-integration required by TD-EOM-CC are generally far larger than those of the frequency-domain form of the approach. This work considers the potential of the short iterative Lanczos (SIL) integration scheme [J. Chem. Phys. 1986, 85, 5870-5876] to reduce the floating-point costs of TD-EOM-CC simulations. Low-energy and K-edge absorption features for small molecules are evaluated using TD-EOM-CC with single and double excitations, with the time-integrations carried out via SIL and fourth-order Runge-Kutta (RK4) schemes. Spectra derived from SIL- and RK4-driven simulations are nearly indistinguishable, and with an appropriately chosen subspace dimension, the SIL requires far fewer floating-point operations than are required by RK4. For K-edge spectra, SIL is the more efficient scheme by an average factor of 7.2.

Molecular fragmentation as a way to reveal early electron dynamics induced by attosecond pulses

We present a theoretical study of the electron and nuclear dynamics that would arise in an attosecond two-color XUV-pump/XUV-probe experiment in glycine. In this scheme, the broadband pump pulse suddenly ionizes the molecule and creates an electronic wave packet that subsequently evolves under the influence of the nuclear motion until it is finally probed by the second XUV pulse. To describe the different steps of such an experiment, we have combined a multi-reference static-exchange scattering method with a trajectory surface hopping approach. We show that by changing the central frequency of the pump pulse, i.e., by engineering the initial electronic wave packet with the pump pulse, one can drive the cation dynamics into a specific fragmentation pathway. Reminiscence of this early electron dynamics can be observed in specific fragmentation channels (not all of them) as a function of the pump-probe delay and in time-resolved photoelectron spectra at specific photoelectron energies. The optimum conditions to visualize the initial electronic coherence in photoelectron and photo-ion spectra depend very much on the characteristics of the pump pulse as well as on the electronic structure of the molecule under study.