Imaging Electron Dynamics with Ultrashort Light Pulses: A Theory Perspective

A wide range of ultrafast phenomena in various atomic, molecular and condense matter systems is governed by electron dynamics. Therefore, the ability to image electronic motion in real space and real time would provide a deeper understanding of such processes and guide developments of tools to control them. Ultrashort light pulses, which can provide unprecedented time resolution approaching subfemtosecond time scale, are perspective to achieve real-time imaging of electron dynamics. This task is challenging not only from an experimental view, but also from a theory perspective, since standard theories describing light-matter interaction in a stationary regime can provide erroneous results in an ultrafast case as demonstrated by several theoretical studies. We review the theoretical framework based on quantum electrodynamics, which has been shown to be necessary for an accurate description of time-resolved imaging of electron dynamics with ultrashort light pulses. We compare the results of theoretical studies of time-resolved nonresonant and resonant X-ray scattering, and time- and angle-resolved photoelectron spectroscopy and show that the corresponding time-resolved signals encode analogous information about electron dynamics. Thereby, the information about an electronic system provided by these time-resolved techniques is different from the information provided by their time-independent analogues.


Introduction
Real-time imaging of electron dynamics is one of the most important and challenging tasks for modern ultrafast science [1][2][3][4][5][6][7][8]. Valence electron dynamics occurs on a subfemtosecond to few-femtosecond timescale (1 fs = 10 −15 s) and its real-time measurement requires attosecond temporal resolution (1 as = 10 −18 s). At the same time, Ångstrom spatial resolution is required to access inter-atomic distances in molecular structures and solids.
Alternatively to using photons for imaging electron dynamics with X-ray scattering, one can consider employing photoelectrons dislodged by ultrashort light pulses for this goal. Timeand angle-resolved photoelectron spectroscopy (TRARPES), i.e., time-and energy-resolved molecular-frame photoelectron angular distributions, using ultrashort extreme ultraviolet (XUV) probe pulses inducing single-photon ionization has been proposed for imaging coherent electron dynamics in molecules [29][30][31][32]. High-energy photoelectrons generated by XUV pulses, on the one hand, provide Ångstrom spatial resolution and, on the other hand, allow for a less sophisticated interpretation than that in the case of low-energy photoelectrons, where multiple scattering effects can be substatial [33].
The goal of this article is to show that imaging a nonstationary electronic system by ultrashort light pulses is nontrivial not only from an experimental but also from a theoretical perspective. We aim to demonstrate that a simple extrapolation of an observable of a stationary measurement to a time-dependent quantity to obtain a time-resolved signal does not work and can lead to erroneous results. Therefore, a careful analysis of time-dependent processes governing a nonstationary measurement is necessary for a correct interpretation of its outcome.
For a measurement of a stationary electronic system by means of X-ray scattering, its electron density is the central quantity that determines a scattering pattern. Thus, one may naively expect that time-resolved scattering patterns from a nonstationary electronic system depend on its time-dependent electron density. However, it has been demonstrated that this assumption is not correct and time-resolved scattering patterns are not connected to a time-dependent electron density due to a considerable contribution of inelastic (Compton-type) processes that cannot be physically avoided [24,[26][27][28]34]. Only taking into account all transitions that can be induced by a broadband, ultrashort probe pulse, one can obtain correct time-resolved scattering patterns. Such an analysis can be accurately performed within the quantum electrodynamics (QED) framework.
Analogously, one may assume that electron dynamics can be measured by means of time-dependent chemical shifts using time-resolved photoelectron spectroscopy. Since electron binding energies, which determine photoelectron spectra, depend on the local chemical environment, one would think that photoelectron peaks would shift following a time-dependent electron density. However, an analysis of a time-resolved photoelectron probability within the QED framework has shown that this expectation is also incorrect [32]. Thus, an intuitive approach to interpret an outcome of a time-resolved measurement can lead to erroneous results.
In this article, we review the theoretical framework based on the QED, which has been developed to describe the interaction of an ultrashort light pulse with a coherently evolving electronic system. In particular, we will concentrate on three techniques, which have been suggested to image electron dynamics in real space and real time, and analyzed within the QED framework, namely, time-resolved NRXS, time-resolved RXS and TRARPES. We will show that results of the QED analyses of the corresponding time-resolved signals can be represented in a common way although these techniques relay on different processes. Time-resolved signals obtained by means of these techniques have a similar temporal dependence on the evolution of the electronic state of a system being probed, which does not coincide with the temporal dependence of its electron density. We will discuss that accounting for certain features characteristic to the corresponding processes, analogous time-and space-dependent quantities of a nonstationary electronic system are encoded in corresponding timeand momentum-resolved signals. A common procedure based on a Fourier analysis of these signals can be applied to disentangle these quantities.
This review article is organized as follows. We describe a coherently evolving electronic system from a quantum-mechanical prospective in Section 2. In Section 3, we introduce the basics of the QED treatment of the interaction of an ultrashort light pulse with such an electronic system. We review the QED analysis of time-resolved NRXS from a nonstationary electronic system and its results in Section 4. The theory and suggested applications of time-resolved RXS employed for imaging of nonstationary electronic systems are reviewed in Section 5. The TRARPES technique is considered in Section 6. We discuss a limitation on a probe-pulse duration for an appropriate time-resolved measurement in Section 7. We summarize the review in Section 8.

Coherently Evolving Electronic System
Nowadays, it is possible to trigger and observe coherent electron dynamics in atoms [3,6], molecules [2,4,5,35] and crystals [36][37][38][39][40], and even control the outcome of a simple chemical reaction [41]. Sub-femtosecond timing synchronization between pump and probe pulses required for observation of such dynamics has been achieved [3,42]. It has been demonstrated that valence-electron wave packets can evolve with a high degree of coherence for much longer than 10 fs [3].
Let us consider such coherently evolving electronic system that had been excited by a pump pulse into a coherent superposition of the electronic eigenstates at time t 0 . Its time evolution is given by Here, |Φ I are the eigenstates and E I are the corresponding eigenenergies of the many-body Hamiltonian of the electronic system in the absence of an X-ray field,Ĥ m . We use atomic units for this and the following expressions. Since the purpose of this article is to discuss techniques that would allow probing such electron dynamics in real space and real time, we do not concentrate on a specific pump process giving rise to |Ψ(t) . We assume that the pump and probe pulses do not overlap in time that makes it possible to describe the probe and the pump steps separately [43].
The density matrix of this nonstationary electronic system given bŷ contains time-dependent complex off-diagonal elements. The time-dependent electron density of the electronic system is given by the relation ρ(

Interaction with a Probe Pulse within the QED Framework
When the electronic system is probed by an electromagnetic pulse, the total Hamiltonian of the whole system, matter and light, can be written as [44] whereâ † k,p andâ k,p are creation and annihilation operators of a photon in the k, p mode of the radiation field with energy ω k = |k|c, where c is the speed of light.Ĥ int is the minimal coupling interaction Hamiltonian between the matter and the electromagnetic field in Coulomb gaugê whereÂ is the vector potential operator of the electromagnetic field, p is the canonical momentum of an electron,ψ † andψ are electron creation and annihilation field operators. Depending on the character and conditions of the interaction of the electronic system with the electromagnetic pulse, one of two terms in Equation (5) becomes relevant for the interaction process and the other one turns to be insignificant. TRARPES, which is governed by absorption, is determined by HamiltonianĤ (1) int . Scattering driven by the second term in Equation (5) is dominant in the case of high-energy X-rays and, thus,Ĥ (2) int determines NRXS. However, the scattering cross section given byĤ (1) int becomes much larger than that given byĤ (2) int in the case of a resonant X-ray pulse and, thus,Ĥ (1) int determines RXS. In the case of X-ray scattering, the signal is determined by the probability P XS (k s ) to observe a scattered photon with momentum k s , which differs from the incoming photon momenta. Since we consider energy-unresolved scattering patterns, the signal is given by the energy average of P XS (k s ), i.e., by the differential scattering probability (DSP), where V is the quantization volume.
Within the density-matrix formalism, the expectation value of some observableÔ in a state, represented by a density matrixρ, is given by the relation Ô = Tr[ρÔ] [45]. Consequently, the probability of X-ray scattering is connected to the operatorÔ k s , which describes the observation of a photon in the scattering mode k s independently on its polarization, whereρ f (t f ) is the total density matrix of the electronic system and the electromagnetic field at time t f after the action of the probe pulse [24,44,45]. The operatorÔ k s is given bŷ Here, the sum is over two possible polarizations of a scattered photon. It is assumed that a photon detector has some acceptance range represented by the function W(ω k s ).
In the case of TRARPES, the signal is directly connected to the probability to observe an electron with momentum Q el given by where the operatorÔ describes the observation of a photoelectron with momentum Q el .ĉ † Q el ,σ andĉ Q el ,σ are creation and annihilation operators of an electron with momentum Q el and spin σ [32,44,45]. Since we consider spin-unresolved photoelectron spectra, the sum is over spin σ.
Let us now consider the total density matrixρ f (t f ) of the matter and the electromagnetic field, which is obtained by the propagation of the initial total density matrixρ 0 with the time evolution . The initial density matrix is given bŷ is the initial density operator of the electromagnetic field with {n} and { n} being complete sets that specify the number of photons in all initially occupied field modes with a distribution ρ X {n},{ n} [44,45]. Thus, the total density matrix can be represented asρ where an appropriate wave function |Ψ {n} , t f , which is an entangled state of the electronic and photonic states, is substituted dependent on the process considered as shown in the next Sections.

QED Description of Time-Resolved Nonresonant X-ray Scattering
Nonresonant X-ray scattering (usually, the specification 'nonresonant' is omitted) is an established technique that is employed to reveal structural information about a sample. In the scattering process, an electronic system in an initial state I interacting with an X-ray pulse absorbs and emits a photon, which leaves a system either in the same state I or brings it to a different final state F. The former case means that the photon has been scattered elastically and its energy is equal to the energy of the incoming X-ray beam ω in . In the latter case, inelastic scattering has taken place and a photon with energy ω in − (E F − E I ) has been scattered. Thus, elastic and inelastic scattering events from a stationary system can be distinguished by the spectroscopy of a scattered photon. Contributions due to elastic scattering to a scattering pattern dominate over that due to inelastic scattering, since transition amplitudes of elastic scattering events sum up coherently. Thus, scattering patterns from a stationary object are determined by elastic scattering and, as a result, encode the electron density ρ(r) of the object via the relation where Q is the scattering vector [46]. Solving the phase problem, one can reconstruct the electron density with the spatial resolution determined by the wave length of the X-ray pulse.
Let us now consider the interaction of an X-ray pulse with a nonstationary electronic system in the state |Ψ(t) = ∑ I C I e −iE I (t−t 0 ) |Φ I . In References [47,48], where this problem has been studied, it was assumed that the time-dependent electron density of a nonstationary electron system at the time of measurement t p , ρ(r, t p ), should be substituted for the electron density ρ(r) in Equation (12) to obtain time-resolved scattering patterns However, the connection of scattering patterns to the electron density stem from elastic scattering. If one would try to consider elastic scattering from an electronic system in the state |Ψ(t) = ∑ I C I e −iE I (t−t 0 ) |Φ I , one would have to assume that the final state of a system F is exactly the same superposition |Ψ(t) as before the interaction with the X-ray pulse. In addition to the fact that the probability of such a situation is extremely weak, it is impossible to spectroscopically distinguish this scattering event from other transitions with final states within the eigenstates of the wave packet |Ψ(t) . In view of this, the concept of 'elastic scattering' in the case of the interaction of an X-ray pulse with a nonstationary electronic system is ambiguous, and the connection of scattering patterns to the electron density must be carefully reconsidered for the time-dependent case [24].
In order to accurately describe the interaction of an X-ray pulse with an electronic system evolving coherently, a theoretical analysis employing the QED framework has been performed by Dixit et al. in Reference [24]. Let us briefly review the result of this work. As discussed in the previous Section, nonresonant X-ray scattering is governed by the HamiltonianĤ (2) int in Equation (5). In order to obtain the probability of nonresonant X-ray scattering with Equation (7), the total density matrixρ f (t f ) must be evaluated within the first-order time-dependent perturbation theory usingĤ (2) int as the perturbation. The resulting density matrix is connected to the first-order wave function which must be substituted for |Ψ {n} , t f in Equation (11). Thus, according to Equations (6) and (7), the DSP of NRXS from an electronic system evolving coherently is with the function Here, −Ê+ ] is the first-order radiation field correlation function, wherê p p e −iω k t+ik·r with p being the unit vector corresponding to the polarization p [44,49]. ω in is the photon energy and in is the mean polarization vector of the incoming X-ray beam. The sum is over all possible final states F with the corresponding energies E F . The role of the function W(ω k s ) is to represent the acceptance range of the photon detector, which can, for example, resolve only certain scattered energies ω k s . This expression considers a general case, which is independent on a probe-pulse duration and coherence properties.
Equation (15) can be simplified under certain conditions. Dixit et al. considered a perfectly coherent Gaussian-shaped probe pulse with the amplitude of the electric field where r 0 is the position of the object, t p is the time of the measurement, τ p is the pulse duration (FWHM of the pulse intensity) and I 0 (r 0 ) = cE 2 (r 0 , t = 0)/(8π) is the peak intensity. It was further assumed that the probe-pulse duration is much shorter than the characteristic time scale of changes in the electron density. Hereafter, we will refer to this assumption as the frozen-density approximation. As we discuss in Section 7, this is a necessary condition for an appropriate time-resolved measurement. The resulting time-dependent DSP can be represented as with the operatorĜ Here, Q = k in − k s is the scattering vector with k in being the mean wave vector of the incoming X-ray beam, P NRXS Let us now compare the resulting relation for the time-resolved DSP in Equation (18) with the relation in Equation (13). Equation (13) that can be represented as dP NRXS (Q)/dΩ ∝ d 3 r 1 d 3 r 2 e iQ·(r 1 −r 2 ) Ψ(t p )|ψ † (r 1 )ψ(r 1 )|Ψ(t p ) Ψ(t p )|ψ † (r 2 )ψ(r 2 )|Ψ(t p ) assumes that elastic scattering has taken place and the final state coincides with the initial electronic state |Ψ(t p ) . Equation (18) contains a sum over final states F, which are eigenstates of the Hamiltonian of the electronic systemĤ m , but not a superposition of electronic states. Contributions due to scattering to final states F are summed incoherently and are weighted by the function e −Ω 2 F τ 2 p /(4 ln 2) originating from the probe-pulse spectral density. Even if one could spectroscopically resolve scattering exclusively to final states, which are involved in the wave packet Ψ(t p ), a contribution due to scattering to each eigenstate would have to be added incoherently, and it would not be possible to substitute Ψ(t p ) for Φ F in Equation (19). As a result, the state Ψ(t p ) enters Equation (18) twice instead of four times as in Equation (13). Thus, the assumption that time-resolved X-ray scattering is given the Fourier transform of the time-dependent electron density is incorrect and even fails to reflect the correct time-dependence of scattering patterns on the evolution of the electronic system.
In order to illustrate the difference between Equations (13) and (18), Dixit et al. considered an electronic wave packet prepared by a pump pulse as a coherent superposition |4 f of the 3d and 4 f eigenstates of atomic hydrogen with the projection of orbital angular momentum equal to zero, the polarization direction along the z-axis and the corresponding energies E 3d and E 4 f . The probe X-ray pulse was assumed to have duration of 1 fs, 4 keV photon energy and propagate along the y direction. They calculated scattering patterns from the electronic wave packet with Equation (13) and with Equation (18) and compared them in Figure 1.  Their structure is much more diverse than the structure of the semiclassical scattering patterns calculated with Equation (13) shown in Figure 1d due to inelastic processes that have to be taken into account in Equation (18). Namely, the calculation of the scattering patterns in Figure 1a considered all possible final states matching the assumed photon detection width of 0.5 eV. However, the most important difference appears at delay times T/4 and 3T/4, when the electronic charge distributions are identical, but the electronic wave packet carry a different phase. The correct scattering patterns in Figure 1a are different from each other at delay times T/4 and 3T/4 unlike the scattering patterns at delay times T/4 and 3T/4 in Figure 1d, which are equal to each other. Thus, the correct scattering patterns depend on the phase of the wave packet, but not on the electron density at the time of measurement as the scattering patterns in Figure 1d do. The further principal difference is that the scattering patterns in Figure 1a are not centrosymmetric (i.e., they are not equal for Q and −Q) at delay times T/4 and 3T/4 breaking the Friedel's law [50]. For the sake of comparison, Dixit et al. have shown scattering patterns calculated with Equation (18) within the QED framework, but restricting the calculation to the two electronic states involved in the dynamics [see Figure 1b]. The structure of these patterns is similar to the structure of the semiclassical patterns in Figure 1d. However, analogously to the patterns in Figure 1a, they depend on the phase of the electronic wave packet and are not centrosymmetric at delay times T/4 and 3T/4. Thus, not only inelastic transitions to a manifold of possible final states do lead to a disagreement between the QED and semiclassical scattering patterns, but also an incorrect description of the time dependence of the DSP in Equation (13) does.
Obviously, the connection to the electronic state encoded in time-resolved scattering pattern is much less straightforward as suggested by Equation (13). In order to interpret information encoded in time-resolved scattering patterns, a Fourier analysis has been proposed for time-resolved resonant scattering patterns [28]. Let us apply such an analysis to the nonresonant case. Equation (18) can be represented as The function entering the expression above is connected to the real part of the elements of the electron density matrixρ m (t p ) [cf. Equation (2)]. Therefore, its temporal dependence correlates with the electron density [cf. Equation (3)] and related quantities such as charge distributions at the time t p . Please notice that this temporal dependence is different to that given by the semiclassical Equation (13), which correlates with the time-dependent electron density squared. The corresponding term in Equation (20) is centrosymmetric with respect to Q. The function is connected to the imaginary parts of the electron density matrix elements and its temporal dependence correlates with the electron current density, which is given by The corresponding term in Equation (20) is an odd function with respect to Q and is responsible for the time-resolved scattering patterns being noncentrosymmetric. Thus, the fact that the scattering patterns in Figure 1 are not centrosymmetric reflects that the electron current is nonzero, i.e., that electrons are moving.
The two contributions to the scattering patterns can be decomposed performing the Fourier transform of the DSP F NRXS (r, t p ) = 1/(2π) 3 d 3 Qe −iQ·r dP NRXS (Q)/dΩ. The real part of the Fourier transform is connected to the function R NRXS via the relation and the imaginary part is connected to the function J NRXS according to Thus, performing the Fourier transform of the scattering patterns from Q space to real space obtained at different time delays t p − t 0 , one obtains a complex function, which depends on the time of measurement t p and the space coordinates r. The real part of this function correlates with charge distributions at the time of measurement t p , and its imaginary part correlates with electron currents at t p . A more specific connection of Re F NRXS (r, t p ) and Im F NRXS (r, t p ) to certain time-dependent quantities has to be identified for a given electronic system. Such an analysis has been performed for the case of RXS [27,28] as described in the next Section.
Several alternative approaches to image electron dynamics by means of time-resolved NRXS have been developed through the application of the QED framework. Dixit et al. has proposed to image the instantaneous electron density of an electronic wave packet via ultrafast X-ray phase contrast imaging [51]. They have shown that ultrafast phase contrast imaging provides the Laplacian of the electron density, which reveals complex bonding and topology of the charge distributions in an electronic system. However, this technique is quite experimentally challenging, since it requires detector pixels that are small enough to resolve the image of small objects such as molecules.
Another approach has been suggested by Grosser et al. in Reference [52], who have shown that the Compton-scattering cross section, in the impulse approximation, depends solely on the electron momentum distribution. Thus, time-resolved Compton scattering can be used to obtain momentum-space images of the sample to be probed. Kowalewski et al. have suggested time-resolved X-ray diffraction as a probe of molecules in the gas phase undergoing nonadiabatic avoided-crossing dynamics involving strongly coupled electrons and nuclei [53]. They have shown that it provides signatures of a created electronic coherence on top of dominant ground-and excited-state wavepacket motions. A photon-coincidence measurement based on the time-and wavevector-resolved detection of photons generated by the scattering of multiple X-ray pulses with variable delays has been proposed by Biggs et al. in Reference [54]. They have shown that it directly measures multipoint correlation functions of the charge density through superpositions of valence excitations which are created impulsively by the scattering process.

QED Description of Time-Resolved Resonant X-ray Scattering
RXS is an element specific technique that provides insight into charge, orbital and spin degrees of freedom [55][56][57][58]. Due to the resonant nature of this process, the scattering cross section can be considerably enhanced in comparison to the nonresonant case discussed in the previous Section. In contrast to NRXS, which probes simultaneously electrons involved in the dynamics and electrons that are essentially stationary [59], RXS can directly probe (quasi-)particles involved in the dynamics. Therefore, RXS particularly suits for measurement of heavy elements, where the vast majority of electrons are stationary, and resonant energies approach several keV corresponding to Ångstrom spatial resolution.
RXS is a two-step process consisting of an absorption and an emission process, which involves an intermediate state in distinction to NRXS. In the first step, the absorption of a photon from an X-ray pulse induces a resonant transition of an electron from a core shell to a valence shell of an atomic system (e.g., a molecule or a crystal) being measured. Thus, after the first step, the electronic system is brought from its initial state I to an intermediate state J C with an electron hole in a core shell of a scattering atom C and an additional electron in the valence shell. In the second step, an electron from either the same valence shell or some other shell fills the electron hole in the core shell leading to a spontaneous emission of a photon, which is detected. The electronic system is in a final state F afterwards. If elastic scattering has taken place, then the final state F coincides with the initial state I, and the energy of the scattered photon is equal to ω in . An inelastic scattering event results in the final state F being different from the initial state I and leads to scattering of a photon with the energy ω in − (E F − E I ). As in the case of NRXS, elastic and inelastic scattering events from a stationary system can be distinguished by the spectroscopy of a scattered photon, and the contribution of elastic scattering events to RXS scattering patterns are dominating. Through elastic scattering, resonant scattering patterns are determined by the electron density of the object being measured, but are indirectly connected to it via the relation [55] where f C is a scattering amplitude of an atom C located at a position R C . Since RXS involves transitions of electrons from core shells, which are highly localized in comparison to X-ray wavelengths, the spatial distribution of the X-ray electric field is treated within the dipole approximation for each absorbing atom. A scattering amplitude f C is proportional to the product of the dipole matrix elements describing the absorption process and the emission process: , let us consider the interaction of a resonant X-ray pulse with an atomic system with coherent electron dynamics in the valence shell described by the state |Ψ(t) = ∑ I C I e −iE I (t−t 0 ) |Φ I . In the absorption step, an electron from a core shell is excited into the valence shell, whereby it destroys the coherent electron dynamics. The system is brought to an intermediate state J C , and the following step, emission, is determined by this intermediate state and is not connected to the state |Ψ(t) . Thus, as in the case of NRXS, it is, first, very unlikely that the final state after the emission process would coincide with the superposition |Ψ(t) and, second, this situation would be spectroscopically indistinguishable. This again demonstrates that the notion of 'elastic scattering' for an interaction of a light pulse with a nonstationary electronic system is unclear. Since the relation in Equation (26) relays on the assumption that elastic scattering provides a dominating contribution to a scattering pattern, it must be reconsidered for time-resolved RXS.
The interaction of an ultrashort resonant X-ray pulse with an electronic system evolving coherently has been described in Reference [27] as follows. Since RXS is a two-step process governed by the HamiltonianĤ (1) int in Equation (5), the second-order wave function must be substituted for |Ψ {n} , t f in Equation (11). Thus, according to Equations (6) and (7), the DSP of RXS from an electronic system evolving coherently is whereT s(in) (r) =ψ † (r)( s(in) · p)ψ(r) and Γ J C is the decay width of the intermediate state J C . This expression is general for a resonant X-ray pulse with arbitrary coherence properties, and assumes neither the frozen-density approximation nor the dipole approximation. The relation in Equation (28) has much in common with the general expression for time-resolved NRXS from a coherently evolving electronic system in Equation (15). Analogously to Equation (15), the DSP for time-resolved RXS involves the sum over final states F with corresponding transition amplitudes summed incoherently and also depends on the first-order correlation function G (1) (r 1 , t 1 , r 2 , t 2 ). Applying the dipole approximation for each absorbing atom, and assuming a coherent resonant X-ray pulse with the amplitude of the electric field defined in Equation (17) and duration that is much shorter than the characteristic time scale of the changes in the electron density, the DSP of RXS can be represented as with the operatorĜ where P RXS 0 = τ 2 p I 0 /(4 ln 2c 4 ) andr = d 3 rψ † (r) r ψ(r). The expression in Equation (30) is also very similar to the expression for the DSP of NRXS in Equation (18); but, instead of the space integrals in Equation (18), it involves the sum over scattering atoms C 1 and C 2 , and their positions R C 1 and R C 2 due to the dipole approximation. In particular, it has a similar temporal dependence on the electronic state at the time of measurement Ψ(t p ).
Let us compare Equation (30) with the expression for the DSP from a stationary system in Equation (26), which can be represented as dP stRXS /dΩ ∝ ∑ C 1 ,C 2 e iQ·(R C 1 −R C 2 ) f C 1 f * C 2 . The DSP according to Equation (30) could be represented as dP RXS /dΩ ∝ . Obviously, it is not possible to obtain the correct expression for the time-resolved DSP by simply substituting the scattering amplitudes This assumption suggests that the absorption process first brings the coherently evolving electronic system into an intermediate state J C , and then the emission process brings the system from the intermediate state J C into the same coherent superposition Ψ(t p ) as before the absorption. Equation (30) indeed takes into account that the absorption process bringing the system into the intermediate state J C destroys electron dynamics, and final states after the emission process are eigenstates of the HamiltonianĤ m . Therefore, the information about electron dynamics is contained in the absorption term of f CF (t p ), but not in the emission term.
The difference between the semiclassical approach to obtain the time-resolved DSP according to Equation (32) and the correct expression in Equation (30) has been illustrated in Reference [27] by the calculation of the DSP from an ionized Br 2 molecule with coherent electron dynamics in the valence shell. The two highest occupied molecular orbitals of Br 2 are the πg and πu orbitals of Br 4p character. They assumed that a photoionizing pump pulse created an electron hole initially localized at the 4pz orbital of one of Br 2 atoms (at site A in Figure 2) at time t 0 = 0, which is actually a coherent superposition of the electron hole being in the 4pπg orbital and the electron hole being in the 4pπu orbital. Such electron-hole localization in a molecule by a photoionizing pump pulse is possible as has been demonstrated in Reference [4]. Then, the electron hole, which was initially localized on one of atoms, starts oscillating coherently between the two atoms of the Br 2 molecule moving from site A to site B and back (see Figure 2). The electron hole dynamics is probed by an X-ray pulse with the photon energy of 13.5 keV tuned to the K edge of Br.  (32) (right) and the corresponding schematic representation of the hole density (middle) at different t p . The Br 2 molecule is aligned along the x direction parameters, the interatomic distance is R x = 2.3 Å and the oscillation period T is 1.7 fs [60]. Probe pulse is polarized along the z direction and propagates along the y-axis with corresponding parameters τ p = 200 as, ω in = 13.5 keV, I 0 = 10 18 W/cm 2 . N is a positive integer. The dependence on polarization is not shown. The ranges are limited by Q 2 x + Q 2 z ≤ ω 2 in /c 2 . It is assumed that the detector acceptance range W(ω k s ) allows detecting only 1s − 4pπg and 1s − 4pπu scattering events. Figure is reproduced with permission from Reference [27]. Copyright American Physical Society, 2015.
The DSP from the Br 2 molecule has been calculated according to the semiclassical expression in Equation (32) and the correct expression in Equation (30), and compared in Figure 2. In both cases, the signals are constant when the hole is localized on one of the atoms (at times 0 and T/2), since the scattering channel for the other atom is blocked at these moments, and there is no interference. In contrast to the semiclassical scattering patterns, the scattering patterns according to Equation (30) are different from each other at times T/4 and 3T/4, when the electron hole densities are equal, and depend on whether the hole is moving from site A to site B or vice versa. These scattering patterns are not centrosymmetric with respect to Q as in the case of NRXS (see Section 4). The QED patterns at times T/4 and 3T/4 are not just phase-shifted with respect to the semiclassical patterns, but also have different amplitudes. For example, their intensity is nonzero at any Q in contrast to the semiclassical patterns, which are zero at Q x R x = π/2. This comparison again demonstrates that the semiclassical approach to describe the interaction of a light pulse with a nonstationary electronic system can lead to erroneous time-and momentum-dependence of a time-resolved signal. Correct results demand a thorough analysis within the QED framework.
The analysis of a feasibility to realize the experiment as described above has been performed in Reference [27]. It has been shown that if one would prepare a beam of aligned molecules as shown in References [61,62] with the density of 10 10 cm −3 and size 0.4 cm, and excite a coherent electron dynamics in 10% of the molecules, then approximately four molecules per shot would contribute to a scattering pattern assuming the interaction area of 10 −8 cm 2 . Thus, summation of signals from roughly 10 5 scattering patterns provides a signal of 0.01 photons per pixel, which are enough to reconstruct a structure of a single molecule [63]. This means that such an experiment is in principle feasible with the forthcoming European X-ray free electron laser facility, which will provide 27,000 X-ray pulses per second [64] and allow collecting the necessary number of scattering patterns within hours.
The scattering patterns in Figure 2 depend on the direction in which the electron hole is moving. This means that the scattering patterns encode the interatomic electron current and it must be possible to reconstruct it. It has been shown in References [27,28] that this is a general property of scattering patterns obtained by means of ultrafast RXS and developed the following method to reconstruct the interatomic electron current. We applied the analogous analysis to that employed by this method for the case of NRXS in the previous Section 4 and obtained very similar results.
Equation (30) can be represented as where the first term is centrosymmetric with respect to Q and is determined by the fuction which depends on the real parts of the density matrix elements and, therefore, evolves in time similarly to the electron density [cf. Equation (3)] and charge distributions. The second term is an odd function with respect to Q and is determined by the function connected to the imaginary parts of the density matrix elements and correlates with the probability current density j(r, t p ) [cf. Equation (23)]. j(r, t p ) can be decomposed into intra-atomic and inter-atomic contributions as where the sums are over all atoms in the system. The volume-integrated probability current between scattering atoms C 1 and C 2 is given by where n C 1 C 2 is the unit vector pointing from site C 1 to site C 2 , andĜ C 1 C 2 = d 3 rξ † C 2 (r)(∇ · n C 1 C 2 )ξ C 1 (r) with the operatorξ C annihilating a particle at site R C . It has been shown in Reference [27] that the time evolution of the interatomic electron current j C 1 C 2 (t p ) is substantially reproduced by the function J RXS C 1 ,C 2 (t p ) and can be reconstructed from the Fourier transform of the DSP as shown below.
Performing the Fourier transform of the DSP in Equation (33), F RXS (r, t p ) = 1/(2π) 3 d 3 Qe −iQ·r dP RXS (Q)/dΩ, its centrosymmetric and noncetrosymmetric contributions can be decomposed. The real part of the Fourier transform is connected to the functions R RXS C 1 ,C 2 via the relation and its imaginary part is connected to J RXS C 1 ,C 2 via Thus, the Fourier transform of a time-resolved resonant scattering pattern is a sum of delta peaks at the interatomic distances R C 1 − R C 2 weighted by complex amplitudes, which depend on the time of measurement t p . In particular, the imaginary part of the amplitude of a peak at R C 1 − R C 2 is given by J RXS C 1 ,C 2 (t p ) and provides the electron current between atoms C 1 and C 2 . If there are several atomic pairs at the same interatomic distance, when the imaginary part corresponds to the sum of the interatomic currents between the corresponding pairs.
The procedure to reconstruct interatomic electron currents in crystals by means of ultrafast RXS has been illustrated by the calculation of time-resolved scattering patterns from KBr and Ge crystals with coherent electron dynamics in the valence bands by D. Popova-Gorelova and R. Santra in Reference [28]. Here, we will review their calculation of coherent electron dynamics in KBr. It is an ionic crystal, where the 4s electrons of K atoms are transferred to Br atoms [65]. The p-character electrons centered on the Br atoms form the outermost valence band of KBr. It was assumed that a photoionizing pump pulse induced coherent electron dynamics in this band by creating p-type electron holes centered on Br atoms at time t 0 . A probe pulse propagating along the x direction of the duration of 200 as and the photon energy of ≈ 13 keV tuned to the K edge of Br polarized along the y-axis parallel to one of the vectors connecting two Br atoms (see Figure 3) was assumed. Figure 3. Snapshot of the electron-hole density at time t p in a fragment of (KBr) 108 cluster visualized using VESTA software [66]. The orange arrows represent the electron-hole currents j C 1 C 2 (t p ) between the Br atoms. The blue-framed arrow represents a current parallel to the probe-pulse polarization. Figure is reproduced with permission from Reference [28]. Copyright American Physical Society, 2015.
Each coherently evolving electron hole is delocalized and distributed over many Br atoms in some region (see Figure 3). It was assumed that the concentration of the electron holes is sufficiently low to consider these regions isolated and the holes noninteracting. A region, where a single electron hole is distributed, was simulated by a (KBr) 108 cluster. Since the delocalized electron hole is coherently evolving, the charge is redistributing and flowing from one atom to another resulting in nonzero interatomic electron current.
A scattering pattern at a time t p is shown in Figure 4. Since scattered polarizations are different at each k s vector, an additional Q-dependence appears in the DSP and leads to a nonperiodic scattering pattern in Figure 4a. This dependence is eliminated in the scattering pattern in Figure 4b, which is determined solely by the trigonometric functions in Equation (33). Figure 5a shows the imaginary part of the Fourier transform from Q space to real space of the scattering pattern in Figure 4b extrapolated to a region of infinite Q y and Q z . In accordance with Equation (38), it consists of delta peaks at points corresponding to vectors connecting pairs of Br atoms (see Figure 3). The amplitudes of peaks at points (R y , R z ) are opposite to the amplitudes of peaks at (−R y , −R z ). Performing an analysis of the electronic structure of KBr, the authors found a way to resolve interatomic currents between Br atoms in a direction parallel to the probe-pulse polarization. They suggested to measure scattering patterns at different time delays t p − t 0 , perform the Fourier transform of these scattering patterns from Q space to real space and follow the amplitude of the peak outlined in the circle in Figure 5a in the imaginary part of the Fourier transform. The evolution of the amplitude of this peak depending on the time of measurement t p is shown in Figure 5b with the solid violet line. It precisely reproduces the time evolution of the computed sum of the currents between all pairs of nearest-neighbor Br atoms connected by the vector (0, R Br-Br , 0) shown in the orange dashed line in Figure 5b, where R Br-Br is the distance between the nearest-neighbor Br atoms. Thus, it has been shown in Reference [28] that if a proper polarization of the incoming probe pulse has been chosen for measurement of scattering patterns, the time evolution of the amplitude of a certain delta peak in their Fourier transform follows the time evolution of the sum of the interatomic currents between atoms connected by the vector corresponding to this peak. The required polarization of the probe pulse can be determined by an analysis of the electron structure of the crystal. To sum up, time-resolved RXS is a robust technique that is in many ways advantageous for probing electron dynamics in molecules and crystals. On top of the usual structural information about an atomic system, it can provide interatomic electron current in it. Feasibility analyses of the experiments to measure electron dynamics in Br 2 [27], and KBr and Ge crystals [28] by means of ultrafast RXS showed that they are achievable at X-ray free electron laser facilities. However, the main obstacle of such experiments is the need to produce hard X-ray pulses of subfemtosecond duration for an interatomic spatial resolution and sufficient temporal resolution. Although generation of such pulses has been recently achieved [23], these constraints on the probe-pulse parameters are still quite demanding. These restrictions can be partially overcome with TRARPES technique [32], which we review in the next Section.

QED Description of Time-Resolved Photoelectron Probability
TRARPES, i.e., time-and energy-resolved molecular-frame photoelectron angular distributions, has been proposed for imaging coherent electron dynamics in molecules [29][30][31][32]. An advantage of this technique over X-ray scattering is that it allows achieving Ångstrom spatial resolution with light pulses of a much lower photon energy. An analysis of time-and angle-resolved photoelectron distributions obtained by ultrashort XUV probe pulses inducing single-photon ionization of a coherently evolving electron system within the QED framework has been performed in Reference [32]. This study provided a correct way to take into account the consequence of the broad probe-pulse bandwidth on a signal, which has not been taken into account in earlier studies of time-resolved photoelectron probability [29], but is critical for a correct interpretation of time-resolved photoelectron spectra. This work is reviewed in this Section.
In order to describe the absorption process, the total density matrix must be evaluated within the first-order time-dependent perturbation theory using the HamiltonianĤ (1) int . Thus, in order to obtain the photoelectron probability, the first-order wave function must be substituted for |Ψ {n} , t f in Equation (11). Thus, according to Equation (9), the photoelectron probability of a coherently evolving electronic system consisting of N el electrons is Here, φ σe (Q el , r) is the wave function of a photoelectron with spin σ and momentum Q el , ε e = |Q el | 2 /2 is the photoelectron energy, |Φ is a final state of the electronic system with N el − 1 electrons and energy E N el −1 F , which by assumption does not interact with the emitted photoelectron. Equation (40) does not include any further assumptions concerning probe-pulse parameters, such as coherence properties or duration.
Applying the dipole and the frozen-density approximation, assuming a coherent XUV pulse with the amplitude of the electric field defined in Equation (17), and applying the plane-wave approximation to the photoelectron wave function, the photoelectron probability can be represented as with the operatorĜ  (18) and (30). The general Equations (15), (28) and (40) depend on the first-order radiation field correlation function G (1) (r 1 , t 1 , r 2 , t 2 ) multiplied by a function dependent on Ψ(t 1 ) and its conjugate dependent on Ψ(t 2 ). Equations (18), (30) and (42) describe analogous dependence of the corresponding time-resolved signals on the evolution of the electronic state Ψ(t p ) and the spectral density of an ultrashort probe pulse. All these relations involve an incoherent sum over final states. The only prominent difference of the relations for the photoelectron probability from the relations for X-ray scattering is that the former comprise the photoelectron wave function φ σe (Q el , r) instead of the wave function of a scattered photon e ik s ·r .
In addition, analogously to time-resolved NRXS and time-resolved RXS, there are two contributions to the photoelectron probability in Equation (42), which can be represented as ε e , r 1 , r 2 ).
The first contribution is centrosymmetric with respect to Q el and is determined by the function which is connected to the real parts of electron density matrix elements. The second contribution is an odd function with respect to Q el and is determined by which is dependent on the imaginary parts of electron density matrix elements. The time-and angle-resolved photoelectron probability is a sum of an even function with respect to Q el determined by R PE (t p , r 1 , r 2 ) and an odd function with respect to Q el determined by J PE (t p , r 1 , r 2 ) not only in the case of the photoelectron wave function being a plane wave, but in the more general case of φ σe (Q el , r) being a Hermitian function with respect to both Q el and r.
By an analogy to the concept of a chemical shift in stationary photoelectron spectroscopy, one may assume that electron dynamics can be measured by means of time-dependent chemical shifts, i.e., that photoelectron peaks would shift following the time-dependent electron density in a time-resolved measurement. However, according to Equation (42), time-dependent photoelectron spectra at each emission angle consist of a series of photoelectron peaks centered at energies Ω F corresponding to a transition to a final state F, whereby the position of the peaks are time independent. Therefore, quite counter-intuitively, electron dynamics cannot be followed by means of time-dependent chemical shifts. The quantity that depends on the electronic state at the time of measurement is the strength of the photoelectron peaks.
Equation (42) was applied for the calculation of time-and angle-resolved photoelectron spectra of an indole molecule with coherent electron dynamics in valence orbitals. It was assumed that a broadband photoionizing pump pulse launched coherent electron dynamics in indole by creating an electron hole in a superposition of the HOMO (highest occupied molecular orbital) and HOMO-1 orbitals at time t 0 (see Figure 6). Then, the electronic state of the indole molecular cation after the interaction with the pump pulse evolves in time as is an electronic state with an electron hole in the HOMO, |Φ ion H−1 is an electronic state with an electron hole in the HOMO-1, and C 1 and C 2 are time-independent complex coefficients, which are determined by the pump process. After the pump process, the electron density starts oscillating with the period T = 2π/(E 2 − E 1 ) ≈ 10.2 fs. Since there are two electronic states involved in the dynamics, there are two time points during the density oscillation period, at which the time-dependent electron densities are identical. The zero time point was adjusted such that the time-dependent electron densities coincide at times t = T/4 and t = 3T/4.
It was suggested to probe the electron dynamics by an XUV pulse, which creates a second electron hole in the indole molecule by a single-photon ionization at time t p , and analyze angular distributions of photoemission probabilities at a fixed photoelectron energy depending on the time delay t p − t 0 . A y-polarized probe pulse of 1 fs duration, 100 eV photon energy, and 10 12 W/cm 2 intensity was considered. An application of the XUV probe pulse provided photoelectron angular distributions at photoelectron energies of ∼ 80 eV-90 eV, which allowed for Ångstrom spatial resolution.
In addition, an analysis of such distributions was simplified, since the orthogonalization correction to the plane-wave approximation was suppressed at these energies in the case of the indole molecule.   Each q x − q z plane at the corresponding photoelectron energy ε e is a projection of the semisphere with radius |q| = √ 2ε e at q y > 0, where the color of the q point on the semisphere corresponds to the probability of detecting an electron with momentum q. Figure is reproduced with permission from Reference [32]. Copyright American Physical Society, 2016.
It was suggested to analyze the angular distributions at a fixed photoelectron energy ε e by performing the Fourier transform according to the following equation where δ(|Q el | − Q 0 ) is the Dirac delta function. This equation describes the Fourier transform of the spherical surfaces in Q el space of a fixed radius Q 0 = √ 2ε e , the projections of which are shown in Figure 7. It was found that the Fourier transform is connected to the time-dependent Dyson orbitals where ] with the function s(r) = sinc(Q 0 |r|). Thus, the Fourier transform of the time-and angle-resolved photoelectron probability at a fixed photoelectron energy ε e is determined by a linear combination of autocorrelation functions of Dyson orbitals A[φ D F (t p )] with the coefficients given by the exponential factor e −(Ω F −ε e ) 2 τ 2 p /(4 ln 2) . Alternatively to the representation of F PE (r, ε e , t p ) via the Dyson orbitals in Equation (48) provided in Reference [32], the Fourier transform can be expressed via the functions R PE and J PE defined in Equations (45) and (46) as Thus, analogously to the results of the Fourier analysis of scattering patterns, the real part of the Fourier transform of the photoelectron angular distributions is determined by the real part of electron density matrix elements, and correlates with the electron density and related quantities such as charge distributions at the time of the probe-pulse arrival. The imaginary part of the Fourier transform is determined by the imaginary part of electron density matrix elements and correlates with the instantaneous electron current.

Applicability of the Frozen-Density Approximation
Let us now consider the treatment of electron dynamics during the interaction of an electronic system with the probe pulse. In all examples discussed in the previous Sections, duration of a probe pulse was chosen such short that changes in the electron density during the probe were negligible. As a result, the frozen-density approximation, which neglects electron dynamics during the interaction with the probe pulse, could be applied for the QED derivation of Equation (18)  In order to check the applicability of the frozen-density approximation, a situation when the probe-pulse duration is comparable to the characteristic time scale of electron dynamics was considered in Reference [32]. It was assumed that the pump pulse launched coherent electron dynamics by creating an electron hole in a superposition of HOMO, HOMO-1, and HOMO-2 of indole, which was measured by the probe pulse of the same parameters as in the previous section. The shortest beating period of the time-dependent electron density in this case was 1.5 fs, which was still larger than the probe-pulse duration of 1 fs.
In order to accurately describe the situation, when the probe pulse duration is comparable with or longer than the characteristic time scale of electron dynamics, one has to take into account the evolution of the electronic system during the action of the probe pulse. The corresponding relation for TRARPES can be derived from the general expression in Equation (40), resulting in demonstrated that the frozen-density approximation broke down and provided completely incorrect results. Thus, they demonstrated that if the probe-pulse duration is comparable with or longer than the shortest beating period of the electron density of a given system, the generated photoelectron spectra are not connected to its instantaneous electronic state as assumed in Equation (42). At the same time, they verified that a probe-pulse duration being ten times shorter than the oscillation period of the electron density is enough to satisfy the frozen-density approximation. where Ω FI(FK) = ω in − ω k s − E F + E I(K) . In addition, such a relation for RXS results from the general relation in Equation (28), leading to FI τ 2 p /(8 ln 2) . (53) In both cases, the relations are not connected to the electronic state at the time of measurement Ψ(t p ) that cannot be singled out. Instead of that, all three relations in Equations (42), (52) and (53) include the electronic state Ψ(t p ) convoluted with the amplitude of the electric field [cf. Equation (17)]. Thus, an assumption that a time-resolved signal obtained by a probe pulse not satisfying the frozen-density approximation depends on the electronic state simply averaged over the pulse duration, t p +τ p /2 Actually, Equations (18), (30) and (42) assuming the frozen-density approximation can be obtained from corresponding Equations (51)-(53) by substituting the mean energy E of the superposition Ψ(t p ) for E I(K) in factors Ω el FI(FK) and Ω FI(FK) , correspondingly. This means that the frozen-density approximation is valid as long as the maximum energy difference of the eigenstates involved in the dynamics is negligible compared to the probe-pulse bandwidth.
Since the dependence of Equations (52) and (53) (51), it is possible to transfer the conclusion of Reference [32] to the case of time-resolved X-ray scattering. Thus, if the duration of an X-ray probe-pulse is about ten times shorter than the shortest beating period of the electron density, then time-resolved scattering patterns are connected to the electronic state Ψ(t p ) and can be described by corresponding Equations (18) and (30). Whereas, if the duration of the X-ray probe pulse is comparable with the shortest beating period, then this connection is lost, and scattering patterns are determined by Ψ(t p ) convoluted with the amplitude of the probe-pulse electric field as described by corresponding Equations (52) and (53).

Discussion
We reviewed the theoretical framework based on QED developed for an accurate description of an interaction between a coherently evolving electronic system and an ultrashort light probe pulse. We concentrated on its application to three techniques, namely, time-resolved NRXS, time-resolved RXS and TRARPES that can be employed for a measurement of electron dynamics in real space and real time [24,[26][27][28]32,[51][52][53][54]68]. It turned out that the corresponding theories and their results have much in common.
The outcomes of the theoretical analyses appeared to be quite counter-intuitive in all considered cases. Time-and momentum-resolved signals from a nonstationary electronic system obtained by means of the considered time-resolved techniques encode different information from their stationary analogues. X-ray scattering, which is determined by the electron density of a stationary object being measured, is not determined by the time-dependent electron density in a time-resolved measurement. The connection of scattering patterns to the electron density relay on elastic scattering dominating over inelastic scattering. In the case of an interaction of a light pulse with a nonstationary electronic system, the concept of elastic scattering is ambiguous, since it is extremely improbable that a final state after the interaction would be the same nonstationary electronic state. Moreover, such a situation would not be spectroscopically distinguishable from other possible inelastic scattering events. As a result, time-resolved scattering patterns encode spatial and temporal correlations substantially deviating from quantities encoded in stationary scattering patterns.
A similar consequence of a time-resolved measurement on photoelectron spectroscopy has been demonstrated. Contrary to an intuitive assumption that electron dynamics can be measured by means of time-dependent chemical shifts (i.e., by following temporal changes of electronic binding energies), time-resolved photoelectron spectra consist of a series of photoelectron peaks centered on time-independent positions. Electron dynamics is indeed encoded in time-dependent amplitudes of these peaks, which exhibit prominent temporal dependence in photoelectron angular distributions at fixed photoelectron energies.
In all cases considered, time-and momentum-resolved signals, i.e., scattering patterns in the case of NRXS and RXS, and photoelectron angular distributions in the case of TRARPES, do not follow instantaneous electron density. They depend on the phase of the electronic wave packet being measured and are not centrosymmetric with respect to the momentum (except for certain time points). It turns out that they have a quite similar dependence on the electronic state at the time of measurement Ψ(t p ) and the spectral density of a probe pulse. In all cases, they are connected to a function, which can be represented as Ψ(t p )|Ĥ(r 1 , r 2 )|Ψ(t p ) , where an operatorĤ(r 1 , r 2 ) includes a sum over final states of corresponding transition amplitudes weighted by a function determined by the spectral density of a probe pulse. The exact form of the operatorĤ(r 1 , r 2 ) depends on the technique by which the signal was obtained.
There are two contributions to the time-and momentum-resolved signals. The first contribution is centrosymmetric with respect to the momentum and is determined by a function connected to the real parts of electron density matrix elements. This means that the temporal evolution of this function correlates with the time-dependent electron density and charge distributions. The second contribution is an odd function with respect to the momentum and is determined by a function connected to the imaginary parts of electron density matrix elements. Its temporal dependence correlates with electron currents. Thus, the fact that the signals are not centrosymmetric reflects that electrons are moving. The two contributions can be disentangled by performing the Fourier transform from momentum space to real space. Specific connections of these functions to certain time-dependent quantities have to be identified for a given electronic system and technique. For instance, such an analysis has been performed for the case of ultrafast RXS in References [27,28], where a method to measure interatomic electron currents was introduced.
The connections of the time-resolved signals to the electronic state Ψ(t p ) hold as long as the bandwidth of a probe pulse is much larger than the maximum energy difference between electronic states involved in the dynamics. This condition is satisfied for a probe-pulse duration being ten times shorter than the shortest beating period of the electron density. If the probe-pulse duration is comparable with this period, when the time-resolved signals are determined by the convolution of the electronic state with the electric-field amplitude of the probe pulse, but not by its temporal average over the duration of the probe-pulse.

Conflicts of Interest:
The author declares no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: QED quantum electrodynamics NRXS nonresonant X-ray scattering RXS resonant X-ray scattering TRARPES time-and angle-resolved photoelectron spectroscopy XUV extreme ultraviolet DSP differential scattering probability HOMO highest occupied molecular orbital

List of symbols
The following symbols are used in this manuscript: state of a coherently evolving electronic system Φ I , Φ K electronic states comprising the coherent superposition C I , C K complex time-independent coefficients of the coherent superposition E I , E K energies of the electronic states comprising the coherent superposition E mean energy of the coherent superposition Φ F final state E F final-state energy ρ electron densitŷ ρ m density matrix of an electronic system ρ f total density matrix of the matter and the electromagnetic field j electron current density ψ electron annihilation field operator T oscillation period of an electronic wavepacket with two eigenstates t 0 time, when the coherent superposition was created t p time of measurement τ p probe-pulse duration t f time after the action of the probe pulse k photon momentum p index referring to a photon polarization a k,p annihilation operator of the photon in the k, p mode ω k photon energy c speed of light G (1) first-order radiation field correlation function E electric-field amplitude of a probe pulse I 0 probe-pulse peak intensity interaction Hamiltonian determined by the A 2 term P XS probability of X-ray scattering dP/dΩ differential scattering probability dP st /dΩ differential scattering probability according to the quasistationary treatment R C position of an atom C f C scattering amplitude of an atom C C, C 1 , C 2 atomic indices Φ J C intermediate state with an electron hole in a core shell of an atom C in the RXS process Γ J C decay width of the intermediate state J C E J C energy of the intermediate state J C j C 1 C 2 interatomic electron current between atoms C 1 and C 2 W(ω k s ) function representing photon-detector acceptance range P PE photoelectron probability Q el photoelectron momentum σ photoelectron spin φ σe photoelectron wave function ε e photoelectron energŷ c Q el ,σ annihilation operator of an electron with momentum Q el and spin σ N el number of electrons in a system before the action of a photoionizing probe pulse F Fourier transform from momentum space to real spacê U time-evolution operator |Ψ {n} , t f wave function being an entangled state of the electronic and photonic states Ω F factor equal to ω in − ω k s − E F + E Ω el