Electron Dynamics In Molecules Research Articles

Cavity-modified molecular dipole switching dynamics.

Polaritonic states, which are formed by resonances between a molecular excitation and the photonic mode of a cavity, have a number of useful properties that offer new routes to control molecular photochemistry using electric fields. To provide a theoretical description of how polaritonic states affect the real-time electron dynamics in molecules, a new method is described where the effects of strong light-molecule coupling are implemented using real-time electronic structure theory. The coupling between the molecular electronic states and the cavity is described by the Pauli-Fierz Hamiltonian, and transitions between polaritonic states are induced via an external time-dependent electric field using time-dependent configuration interaction (TDCI) theory, producing quantum electrodynamics TDCI (QED-TDCI). This method is used to study laser-induced ultrafast charge transfer and dipole-switching dynamics of the LiCN molecule inside a cavity. The increase in cavity coupling strength is found to have a significant impact on the energies and transition dipole moments of the molecule-cavity system. The convergence of the polaritonic state energies as a function of the number of included electronic and photonic basis states is discussed.

Photoionization cross sections and photoelectron angular distributions of molecules with XCHEM-2.0

The XCHEM code was introduced in 2017 [1] to provide an accurate description of electron correlation and exchange in the electronic continuum of molecules at the same level as complete or restricted active-space self-consistent field (CASSCF or RASSCF) methods. This has allowed for an accurate description of molecular photoionization in the region of Feshbach resonances, shake up processes in which ionization is accompanied by excitation of one or several of the remaining electrons, and interchannel couplings. The success of XCHEM for small molecules has led us to improve its performance in several aspects, which now allows for the description of resonant molecular photoionization in larger systems. In addition, we have incorporated the possibility to calculate photoelectron angular distributions in the laboratory and molecular frames, which are essential to interpret angularly resolved photoionization experiments. Here we show its performance in the N2 and pyrazine molecules. The new version of the code, XCHEM-2.0, is freely available at https://gitlab.com/xchem/xchem_public. Program summaryProgram Title: XCHEM-2.0CPC Library link to program files:https://doi.org/10.17632/t8tbk9gdt2.1Developer's repository link:https://gitlab.com/xchem/xchem_public/-/archive/2.0.0/xchem_public-2.0.0.tar.gzLicensing provisions: LGPLProgramming language: Fortran, pythonNature of problem: Understanding photoionization of molecules is a recurrent goal in Atomic, Molecular and Optical physics, even more nowadays with the advent of XUV and X-ray high-harmonic generation sources and free electron lasers. Photoionization can be employed to induce electron dynamics in molecules, e.g., by using ultrashort XUV or X-ray pulses, or as a probe of these dynamics by recording time-resolved photoelectron spectra. At low photoelectron energies, photoionization may result from various competing processes, such as direct ionization, autoionization of Feshbach and shape resonances, shake up, etc., each of these leaving its trace in the measured spectra. So, interpretation of these spectra can be a difficult task and theoretical simulations are often required to separate the different contributions. However, a method able to do so must be able to describe electron correlation in the electronic continuum beyond the Hartree-Fock (HF) and configuration interaction singles (CIS) approximations. XCHEM-2.0 has been designed to account for all these processes by providing an accurate description of electron correlation in the ionization continuum of molecules.Solution method: XCHEM-2.0 employs a close coupling formalism to describe the electronic continuum of molecules at the level of multiference configuration interaction (MRCI) methods in combination with a hybrid Gaussian-Bspline basis set (the so-called GABs basis). In particular, it makes use of a restricted active space self-consistent field (RASSCF) approximation to evaluate the bound electronic states of the remaining molecular cation using localized Gaussian functions. These localized Gaussians are then supplemented with a single centered Gaussian expansion to allow the electron to leave the nuclear environment. Finally, the basis set includes a set of B-spline functions, centered at the same position as the single-centered Gaussians, which reach the asymptotic region. Thanks to the last B-spline function, which does not vanish at the box boundary, XCHEM-2.0 can determine continuum states fulfilling any required asymptotic behavior. The present version of the code improves on the original version by incorporating a more efficient augmentation procedure to build Ne-electron configurations from (Ne−1) ones, reducing the space to store data, providing a more efficient removal of linear dependencies resulting from the over-completeness of the polycentric+GABS basis and a more efficient solution of the scattering equations, and including new routines to calculate photoelectron angular distributions in the laboratory and molecular frames.Additional comments including restrictions and unusual features: The proposed method focuses on RASSCF methodology, so part of electron correlation is still not included. Also, one- and two electron integrals involving exclusively Gaussian functions are done using a modified version of the OpenMOLCAS software, so the installation of this code is also required and limitation associated with the latter is applicable. The code is limited to include a few ionization channels (around 20) with medium angular momentum (limited to ℓ = 15 as in OpenMOLCAS).

Hamiltonian formulation and symplectic split-operator schemes for time-dependent density-functional-theory equations of electron dynamics in molecules

We revisit Kohn–Sham time-dependent density-functional theory (TDDFT) equations and show that they derive from a canonical Hamiltonian formalism. We use this geometric description of the TDDFT dynamics to define families of symplectic split-operator schemes that accurately and efficiently simulate the time propagation for certain classes of DFT functionals. We illustrate these with numerical simulations of the far-from-equilibrium electronic dynamics of a one-dimensional carbon chain. In these examples, we find that an optimized 4th order scheme provides a good compromise between the numerical complexity of each time step and the accuracy of the scheme. We also discuss how the Hamiltonian structure changes when using a basis set to discretize TDDFT and the challenges this raises for using symplectic split-operator propagation schemes.

Approximation Schemes to Include Nuclear Motion in Laser-Driven Ab Initio Electron Dynamics: Application to High Harmonic Generation.

The quantum mechanical description of many-electron dynamics in molecules driven by short laser pulses is at the heart of theoretical attochemistry. In addition to the formidable time-dependent electronic structure problem, the field faces the challenge that nuclear motion, ideally also treated quantum mechanically, may not be negligible, but scales enormously in effort. As a consequence, most first-principles calculations on ultrafast electron dynamics in molecules are done within the fixed-nuclear approximation. For laser-pulse excitation in H2+, where an "exact" treatment of the coupled nuclear-electron dynamics is possible, it has been shown that nuclear motion can have a nonnegligible impact on high harmonic generation (HHG) spectra (Witzorky et al., J. Chem. Theor. Comput. 2021, 17, 7353-7365). It is not so clear, however, how to include (quantum) nuclear motion also for more complicated molecules, with more electrons and/or nuclei, in particular when the electronic structure is described by correlated, multistate wavefunction methods such as the time-dependent configuration interaction (TD-CI). In this work, we suggest a scheme in which the Born-Oppenheimer potential energy surfaces of a molecule are approximated by model potentials (harmonic and asymptotic, as an expansion in 1/R), obtained from only a few ab initio calculations, with the prospect to treat complex molecular systems. The method is tested successfully for HHG by few-cycle laser pulses for the "exact" H2+ reference. It is then applied for diatomic molecules with more electrons and for a two-dimensional model of the water molecule using TD-CIS (S = single) for the electronic structure part.

Iterative projection algorithm for retrieval of angle-resolved single-molecule dipoles from high-harmonic spectra

High-order harmonic spectroscopy (HHS), as an emerging technique, provides a way to detect molecular structures with \aa{}ngstr\"om spatial resolution and ultrafast electron dynamics in molecules on an attosecond timescale. However, the implementation of HHS still encounters great difficulties due to the coherent average of single-molecule emissions from different alignment angles. Here, we provide a machine learning algorithm based on the iterative projection method to fully retrieve the complex single-molecule dipole of a linear molecule from the harmonic spectra measured at different delays and polarization angles between the alignment and driving pulses. We demonstrated our algorithm both theoretically and experimentally using the ${\mathrm{N}}_{2}$ molecule. The results show that our algorithm can accurately retrieve the complex single-molecule dipole from the harmonic spectra with good noise stability and robustness. From the retrieved single-molecule dipole, the contributions of multiple orbitals in high harmonic generation are identified even with a low degree of alignment.

Attosecond imaging of photoinduced dynamics in molecules using time-resolved photoelectron momentum microscopy

We explore the novel capabilities offered by attosecond extreme ultraviolet and x-ray pulses that can be now generated by free-electron lasers and high-harmonics generation sources for probing photon-induced electron dynamics in molecules. We theoretically analyze how spatial and temporal dependence of charge migration in a pentacene molecule can be followed by means of time-resolved photoelectron microscopy on the attosecond time scale. Performing the analysis, we accurately take into account that an attosecond probe pulse leads to considerable spectral broadening. We demonstrate that the excited-state dynamics of a neutral pentacene molecule in the real space map onto unique features of photoelectron momentum maps.

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 N2 and CO2 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.

Isolating Attosecond Electron Dynamics in Molecules where Nuclei Move Fast.

Capturing electronic dynamics in real time has been the ultimate goal of attosecond science since its beginning. While for atomic targets the existing measurement techniques have been thoroughly validated, in molecules there are open questions due to the inevitable copresence of moving nuclei, which are not always mere spectators of the phototriggered electron dynamics. Previous work has shown that not only can nuclear motion affect the way electrons move in a molecule, but it can also lead to contradictory interpretations depending on the chosen experimental approach. In this Letter we investigate how nuclear motion affects and eventually distorts the electronic dynamics measured by using two of the most popular attosecond techniques, reconstruction of attosecond beating by interference of two-photon transitions and attosecond streaking. Both methods are employed, in combination with abinitio theoretical calculations, to retrieve photoionization delays in the dissociative ionization of H_{2}, H_{2}→H^{+}+H+e^{-}, in the region of the Q_{1} series of autoionizing states, where nuclear motion plays a prominent role. We find that the experimental reconstruction of attosecond beating by interference of two-photon transitions results are very sensitive to bond softening around the Q_{1} threshold (27.8eV), even at relatively low infrared (IR) intensity (I_{0}∼1.4×10^{11} W/cm^{2}), due to the long duration of the probe pulse that is inherent to this technique. Streaking, on the other hand, seems to be a better choice to isolate attosecond electron dynamics, since shorter pulses can be used, thus reducing the role of bond softening. This conclusion is supported by very good agreement between our streaking measurements and the results of accurate theoretical calculations. Additionally, the streaking technique offers the necessary energy resolution to accurately retrieve the fast-oscillating phase of the photoionization matrix elements, an essential requirement for extending this technique to even more complicated molecular targets.

Picometer-Resolved Photoemission Position within the Molecule by Strong-Field Photoelectron Holography.

Laser-induced tunneling ionization is one of the fundamental light-matter interaction processes. An accurate description of the tunnel-ionized electron wave packet is central to understanding and controlling subsequent electron dynamics. Because of the anisotropic molecular structure, tunneling ionization of molecules involves considerable challenges in accurately describing the tunneling electron wave packet. Up to now, some basic properties of the tunneling electron from molecules still remain unexplored. Here, we demonstrate that the tunneling electron from a molecule is not always emitted from the geometric center of the molecule along the tunnel direction. Rather, the photoemission position depends on the molecular orientation. Using a photoelectron holographic technique, we determine the photoemission position for a nitrogen molecule relative to the molecular geometric center to be 95±21 pm when the molecular axis is oriented along the tunnel direction. Our Letter poses, and answers experimentally, a fundamental question as to where the molecular photoionization actually begins, which has significant implications for time-resolved probing of valence electron dynamics in molecules.

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.

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.

Time-Resolved Excited-State Analysis of Molecular Electron Dynamics by TDDFT and Bethe-Salpeter Equation Formalisms.

In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe–Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.

Realization of a Continuously Phase-Locked Few-Cycle Deep-UV/XUV Pump-Probe Beamline with Attosecond Precision for Ultrafast Spectroscopy

Chemical and physical processes in molecules can be controlled through the manipulation of quantum interferences between rotational, vibrational, and electronic degrees of freedom. Most of the past efforts have been focused on the control of nuclear dynamics. Even though electronic coherence and its coupling to nuclear degrees of freedom may profoundly affect the outcome of these processes, electron dynamics have received less attention. Proper investigation of electron dynamics in materials demands ultrafast sources in the visible, ultraviolet (UV), and extreme ultraviolet (XUV) spectral region. For this purpose, a few-cycle deep-UV and XUV beamlines have been constructed for studying ultrafast electron dynamics in molecules. To ensure the required high temporal resolution on the attosecond time scale, vibration isolation from environmental mechanical noise and active stabilization have been implemented to achieve attosecond timing control between pump and probe pulses with excellent stability. This is achieved with an actively phase-stabilized double-layer Mach-Zehnder interferometer system capable of continuous time-delay scans over a range of 200 fs with a root-mean-square timing jitter of only 13 as over a few seconds and ~80 as of peak-to-peak drift over several hours.

First-Principles Simulations of X-ray Transient Absorption for Probing Attosecond Electron Dynamics.

X-ray transient absorption spectroscopy (XTAS) is a promising technique for measuring electron dynamics in molecules and solids with attosecond time resolutions. In XTAS, the elemental specificity and spatial locality of core-to-valence X-ray absorption is exploited to relate modulations in the time-resolved absorption spectra to local electron density variations around particular atoms. However, interpreting these absorption modulations and frequency shifts as a function of the time delay in terms of dynamics can be challenging. In this paper, we present a first-principles study of attosecond XTAS in a selection of simple molecules based on real-time time-dependent density functional theory (RT-TDDFT) with constrained DFT to emulate the state of the system following the interaction with a ultraviolet pump laser. In general, there is a decrease in the optical density and a blue shift in the frequency with increasing electron density around the absorbing atom. In carbon monoxide (CO), modulations in the O K-edge occur at the frequency of the valence electron dynamics, while for dioxygen (O2) they occur at twice the frequency, due to the indistinguishability of the oxygen atoms. In 4-aminophenol (H2NC6H4OH), likewise, there is a decrease in the optical density and a blue shift in the frequency for the oxygen and nitrogen K-edges with increasing charge density on the O and N, respectively. Similar effects are observed in the nitrogen K-edge for a long-range charge-transfer excitation in a benzene (C6H6)–tetracyanoethylene (C6N4; TCNE) dimer but with weaker modulations due to the delocalization of the charge across the entire TCNE molecule. Additionally, in all cases, there are pre-edge features corresponding to core transitions to depopulated orbitals. These potentially offer a background-free signal that only appears in pumped molecules.

Electron dynamics method using a locally projected group diabatic Fock matrix for molecules and aggregates

We propose a method using reduced size of Hilbert space to describe an electron dynamics in molecule and aggregate based on our previous theoretical scheme (Yonehara and Nakajima, 2017). The real-time time-dependent density functional theory is combined with newly introduced projected group diabatic Fock matrix. First, this projection method is applied to a test donor-acceptor dimer, namely, a naphthalene-tetracyanoethylene with and without initial local excitations and light fields. Secondly, we calculate an absorption spectrum of five-unit-polythiophene monomer. The importance of feedback of instantaneous density to Fock matrix is also clarified. In all cases, half of the orbitals were safely reduced without loss of accuracy in descriptions of properties. The present scheme provides one possible way to investigate and analyze a complex excited electron dynamics in molecular aggregates within a moderate computational cost.

Ultrafast ab Initio Quantum Chemistry Using Matrix Product States.

Ultrafast dynamics in chemical systems provide a unique access to fundamental processes at the molecular scale. A proper description of such systems is often very challenging because of the quantum nature of the problem. The concept of matrix product states (MPS), however, has proven its performance in describing such correlated quantum systems in recent years for a wide range of applications. In this work, we continue the development of the MPS approach to study ultrafast electron dynamics in quantum chemical systems. The method combines time evolution schemes, such as fourth-order Runge-Kutta and Krylov space time evolution, with MPS, in order to solve the time-dependent Schrödinger equation efficiently. This allows for describing electron dynamics in molecules on a full configurational interaction (CI) level for a few femtoseconds after excitation. As a benchmark, we compare MPS-based calculations to full CI calculations for a chain of hydrogen atoms and for the water molecule. Krylov space time evolution is in particular suited for the MPS approach, as it provides a wide range of opportunities to be adjusted to the reduced MPS dimension case. Finally, we apply the MPS approach to describe charge migration effects in iodoacetylene and find direct agreement between our results and experimental observations.

A new method to calculate ultrafast electron dynamics in molecules based on matrix product states

We propose a new method to describe electron dynamics in molecules on the scale of femtoseconds. It is based on factorizing the electronic wave function into a matrix product state and using this factorization to solve the time dependent Schrodinger equation.

SALMON: Scalable Ab-initio Light–Matter simulator for Optics and Nanoscience

SALMON (Scalable Ab-initio Light–Mattersimulator for Optics and Nanoscience, http://salmon-tddft.jp) is a software package for the simulation of electron dynamics and optical properties of molecules, nanostructures, and crystalline solids based on first-principles time-dependent density functional theory. The core part of the software is the real-time, real-space calculation of the electron dynamics induced in molecules and solids by an external electric field solving the time-dependent Kohn–Sham equation. Using a weak instantaneous perturbing field, linear response properties such as polarizabilities and photoabsorptions in isolated systems and dielectric functions in periodic systems are determined. Using an optical laser pulse, the ultrafast electronic response that may be highly nonlinear in the field strength is investigated in time domain. The propagation of the laser pulse in bulk solids and thin films can also be included in the simulation via coupling the electron dynamics in many microscopic unit cells using Maxwell’s equations describing the time evolution of the electromagnetic fields. The code is efficiently parallelized so that it may describe the electron dynamics in large systems including up to a few thousand atoms. The present paper provides an overview of the capabilities of the software package showing several sample calculations. Program summaryProgram Title: SALMON: Scalable Ab-initio Light–Matter simulator for Optics and NanoscienceProgram Files doi:http://dx.doi.org/10.17632/8pm5znxtsb.1Licensing provisions: Apache-2.0Programming language: Fortran 2003Nature of problem: Electron dynamics in molecules, nanostructures, and crystalline solids induced by an external electric field is calculated based on first-principles time-dependent density functional theory. Using a weak impulsive field, linear optical properties such as polarizabilities, photoabsorptions, and dielectric functions are extracted. Using an optical laser pulse, the ultrafast electronic response that may be highly nonlinear with respect to the exciting field strength is described as well. The propagation of the laser pulse in bulk solids and thin films is considered by coupling the electron dynamics in many microscopic unit cells using Maxwell’s equations describing the time evolution of the electromagnetic field.Solution method: Electron dynamics is calculated by solving the time-dependent Kohn–Sham equation in real time and real space. For this, the electronic orbitals are discretized on a uniform Cartesian grid in three dimensions. Norm-conserving pseudopotentials are used to account for the interactions between the valence electrons and the ionic cores. Grid spacings in real space and time, typically 0.02 nm and 1 as respectively, determine the spatial and temporal resolutions of the simulation results. In most calculations, the ground state is first calculated by solving the static Kohn–Sham equation, in order to prepare the initial conditions. The orbitals are evolved in time with an explicit integration algorithm such as a truncated Taylor expansion of the evolution operator, together with a predictor–corrector step when necessary. For the propagation of the laser pulse in a bulk solid, Maxwell’s equations are solved using a finite-difference scheme. By this, the electric field of the laser pulse and the electron dynamics in many microscopic unit cells of the crystalline solid are coupled in a multiscale framework.

The role of the CI expansion length in time-dependent studies.

With the recent advances in experimental attosecond science, theoretical predictions of electron dynamics can now be validated against experiment. Time-dependent studies of the electron motion in molecules can be used to obtain information about electronic transitions and the interaction of the electrons with electromagnetic fields. Often, these approaches rely on single-excited wave functions. Presented here is a first attempt to evaluate the accuracy of the time-dependent configuration interaction method so that the optimal representation of the electronic wave function for time-dependent studies can be assessed. A quantifier is determined that can aid in finding this optimal representation. The approach is demonstrated on a variety of molecules that include both localized and intramolecular charge transfer electron excitations. Observables including excitation energies, dipole moments, strengths, and static polarizabilities are obtained from time-independent and time-dependent calculations and are compared to experimental data. In this way, a rigorous routine is developed by which the reliability and accuracy of the CI wave function can be assessed and which represents a first step to a more quantitative description of electron dynamics in molecules.

Correlation-driven charge migration following double ionization and attosecond transient absorption spectroscopy

We theoretically investigate charge migration following prompt double ionization. Thereby, we extend the concept of correlation-driven charge migration, which was introduced by Cederbaum and coworkers for single ionization [Chem. Phys. Lett. 307, 205 (1999)], to doubly ionized molecules. This allows us to demonstrate that compared to singly ionized molecules, in multiply ionized molecules, electron dynamics originating from electronic relaxation and correlation are particularly prominent. In addition, we also discuss how these correlation-driven electron dynamics might be evidenced and traced experimentally using attosecond transient absorption spectroscopy. For this purpose, we determine the time-resolved absorption cross section and find that the correlated electron dynamics discussed are reflected in it with exceptionally great detail. Strikingly, we find that features in the cross section can be traced back to electron hole populations and time-dependent partial charges and hence, can be interpreted with surprising ease. By taking advantage of element-specific core-to-valence transitions even atomic spatial resolution can be achieved. Thus, with the theoretical considerations presented, not only do we predict particularly diverse and correlated electron dynamics in molecules to follow prompt multiple ionization but we also identify a promising route towards their experimental investigation.