Attosecond-Resolved Electron Dynamics in Many-Electron Atoms: Quantitative Theory and Comparison with Measurements

: A variety of processes originating from the interaction of atomic or molecular N-electron states with strong and/or hypershort radiation pulses can be understood quantitatively only by ﬁrst determining with good accuracy the solutions of the many-electron time-dependent Schrödinger equation (METDSE) that describe the corresponding physics. The METDSE is solvable nonperturbatively via the state-speciﬁc expansion approach (SSEA). SSEA solutions have been used, or can be used, for quantitative explanation and numerically reliable predictions of quantities that have been measured or are measurable in modern laser-driven experiments that can track, with hypershort (attosecond) time resolution, the effects of electron rearrangements in atoms and molecules. The calculations take into account in a transparent way the interplay between the phenomena and the electronic structures of the physically signiﬁcant states in discrete and multichannel continuous spectra, including multiply- and inner-hole–excited resonance states. The discussion focuses on novel topics of time-resolved many-electron physics and includes a comparison of our predictions to recent quantitative measurements of attosecond-resolved generation of the proﬁle of the ( 2 s 2 p ) 1 P o doubly excited resonance state of helium during photoionization and of the relative time delay in photoemission of the (2 s ,2 p ) electrons of neon. on


Prologue
The paper was written in response to an invitation by the editor (M. Schultze) to contribute a brief account of theoretical work and corresponding results that my colleagues Komninos and Mercouris and I have produced since the 1990s on topics that were announced for this special issue. The announcement states that the contribution may be a review "on studies of ultrafast phenomena and tools and methods to perform them". The following discussion aims at doing just that, by commenting briefly on certain of our proposals, computational tools, and numerical results on themes concerning electron dynamics in many-electron atoms and molecules that can be studied, theoretically and experimentally, with attosecond resolution.
Solving such past and future prototypical problems involves calculating and using the many-electron, time-dependent wavefunction,Ψ(t), describing the physics of interest. This is obtainable by nonperturbatively solving the many-electron time-dependent Schrödinger equation (METDSE) from first principles, in terms of the state-specific expansion approach (SSEA), which was introduced in 1993-1994 [1][2][3].
The label "many-electron" is meant to include atoms and molecules with more than two electrons. The initial states can have arbitrary electronic structures, where the couplings may create zero-order

Introduction
In modern studies of the interaction of hypershort and/or strong radiation pulses for a large range of wavelengths with atoms and molecules, a desideratum of crucial importance is the possibility of extracting quantitative (i.e., not just descriptive) information and understanding about nonlinear processes for different wavelengths, (e.g., high-harmonics generation) or time-resolved phenomena of hyperfast electron dynamics, where both the discrete spectrum and the multichannel continua, including resonance (autoionizing) electronic states, are involved.
As regards theory, such information and understanding can become accessible via formal constructions and computational methods that can go beyond the level of phenomenology or one-electron models and solve a variety of TDMEPs on the time axis. This endeavor depends directly on the possibility of determining and using the time-dependent solution,Ψ(t), of the METDSE. The matter-field interaction contains the pulse characteristics of duration, shape, intensity, and frequency.
If the many-electron Ψ(t) for a particular problem can be calculated accurately, its use within an appropriate theoretical construction allows the reliable description of the physics of interest, including, of course, quantitative predictions, the hallmark of any generally applicable theory. The same physics is obtainable via the appropriate measurement.
For ordinary ground states of neutral atoms and molecules interacting with pulsed radiation, if the field is strong, say 10 12 to 10 13 W/cm 2 and above, and the pulse is hypershort, say in the range of attoseconds (as) to a few decades of femtoseconds (fs) (normally not too many cycles of the field), Ψ(t) must, in general, be obtained nonperturbatively, taking into account, correctly and systematically, the significance of each problem's contributions to the states of both the discrete and the continuous spectra. This criterion is satisfied by the SSEA, which is conceptually simple, yet general and transparent with respect to the interplay between dynamics and electronic structures and spectra. A brief description of this computational tool in the theoretical study of time-resolved electron dynamics is given in Section 3. For details, the reader is referred to [1][2][3] and their references.
Implementation of the SSEA has allowed different types of TDMEPs to be solved, with quantitative results that can be compared with a variety of experiments. For example, see [4,5] for results from the calculation and quantitative understanding of observables in the multiphoton ionization of helium by UV and EUV femtosecond pulses.
Lest it is thought from the previous statement that the SSEA is limited to two-electron systems, I note that the applications include many-electron atoms and topics of multiphoton dissociation and association of diatomics, with initial-state wavefunctions being of the closed-shell as well as the open-shell type [1][2][3]. For example, the SSEA calculation of Ψ(t) (including state-specific correlations), allowed the quantitative prediction of the time-resolved creation during photoionization of the profiles of inner-hole autoionizing states of aluminum decaying into a multichannel continuum [6].
The significance of going beyond two-electron systems when attempting to solve the METDSE while taking into account the state-specific self-consistent fields and electron correlations for any type of electronic structures was underlined in [6] as follows: This is clearly a desideratum for experiment as well, if progress is to be characterized by new information on a variety of real systems. For example, only a small part of what one learns about control or about the effects of electron correlation from the experimental and/or theoretical studies of the interaction of the He 1s 2 1 S ground state with a strong electromagnetic pulse, can be transferred

Introduction
The fundamental equation that contains information about the properties and phenomena that are observable when a pulse of radiation interacts with an atomic or molecular state is the METDSE: H is the N-electron field-free Hamiltonian, consisting of one-and two-electron nonrelativistic or relativistic operators. V(ω, t) is the interaction one-electron operator. Apart from time and the photon frequency ω, it contains information about the polarization and intensity of the field, and about the temporal characteristics of the pulse.
The possibility of solving the METDSE for TDMEPs in real systems lies at the heart of many-electron quantum mechanics. For many decades this goal was ignored, not only because of the complexity of the TDMEP, but also because the understanding of properties and phenomena did not require quantitative knowledge of the many-electron Ψ(t) of Equation (1). Instead, the focus was on the possibility of obtaining accurate solutions of the many-electron time-independent Schrödinger equation (METISE) for the calculation of properties of stationary states, including time-independent rates of transition, an activity that remains dominant, especially for N-electron ground states.
Given the countless number of possible many-electron states and spectra and the variety of possible radiation pulses, a generally applicable theory for the ab initio solution of Equation (1) must be able to incorporate in a tractable way the important requirements of advanced theory of electronic structures, of electron correlations, and of the scattering continuum.
It should be recalled that most publications and books related to this subject ignore the serious and challenging issue of properly solving the TDMEP. In other words, no N-electron wave functions are calculated and used. Instead, these publications are exhausted at the level of either pure formalism or of phenomenology. In such categories are many discussions based on two-or three-or four-level models, where the N-electron matrix elements and energies are taken to be parameters.
In other types of investigations, the complexity of the many-electron problem is bypassed by conveniently replacing the true electronic Hamiltonian with a solvable independent electron model (IEM). IEM-type methods of solving Equation (1) are computationally convenient when applied to closed-shell (single determinantal) electronic structures. (When it comes to open-shell initial states, in most cases electron couplings are such that the zero-order description in terms of a single determinant is inadequate.) They may yield gross information for certain observables, analogous to the IEM that is used for the approximate solution of the METISE for closed-shell ground states. However, they cannot handle questions regarding the effects of electron correlations, of double (or multiple) excitations, and of interchannel couplings in the discrete and the continuous spectrum. Hence, unless one has a good understanding from experience or formal analysis of the situations for which the IEM may be reasonably accurate and meaningful, predictions for arbitrary systems based on IEM-type approaches are either practically impossible or unreliable.

The State-Specific Expansion Approach (SSEA)
For one-electron systems, the solution of Equation (1) for various problems has been achieved numerically, either via direct integration on a space-time grid or via a method of expansion over some type of one-electron basis functions. The grid method has also proven practical in certain cases involving closed-shell N-electron states, where, following the proposal and results of Kulander in the late 1980s [18,19], the formalism is simplified by assuming that the field-driven integration involves only one active electron at a time (single active electron (SAE) model). Space-time numerical integration over a grid in all dimensions of the partial differential Equation (1) has also been implemented in the case of the two-electron closed-shell 1 S state. (e.g., see [20], where differences between the SAE and the full two-electron results are attributed to electron correlations).
The full grid approach is practically inapplicable to arbitrary N-electron electronic structures, i.e., beyond He or H 2 , even with the help of supercomputers. This limitation holds not only for the METDSE, but also for the eigenvalue problem of the METISE, as was recognized in the 1960s.
Evidently, the theory for solving the METDSE for arbitrary electronic structures must follow a different path. This fact determined the research program we initiated in the early 1990s. We concluded that a robust and general approach would be one that invokes the fundamental expansion principle in terms of stationary states with complex time-dependent coefficients, and capitalizes on knowledge of, and experience with, calculating the stationary state-specific N-electron wave functions of the discrete and of the continuous spectrum that solve the METISE [1][2][3]21,22]: k labels stationary states of either the discrete or the continuous spectrum. The latter are energy-normalized. q stands for the coordinates of the N electrons.
In our SSEA work on the solution of Equation (1), the electron-field coupling operator in V(ω, t) has been used either in the electric dipole approximation or as the full electric operator of the multipolar Hamiltonian [23] constrained by dipole selection rules ( [24] and its references).
The critical question is how to nonperturbatively solve the METDSE using H(t) and the many-electron H of Equations (1) and (2). For this purpose, we write the formal solution of Equation (1) as follows (here, for simplicity, I use only one-channel for the continuum): where |m > are the stationary states of the discrete spectrum and |ε > are those of the continuous spectrum, including possible resonances. The scattering wave functions are energy-normalized according to < ε|ε > = δ(ε − ε ). The complex coefficients represent the projection of the wavepacket |Ψ(t) > onto the corresponding stationary states. The objective of the SSEA is to represent the |m > and |ε > of Equation (3) by the state-specific wavefunctions that are relevant to the problem of interest, and to correctly compute the time-dependent coefficients. The latter task is accomplished by solving the coupled equations that result from substituting Equation (3) into the METDSE. The calculation is transparent for each trial set of state-specific wavefunctions, and is carried out by adjusting parameters such as the time step or the range of the free-electron energy and the number of points (solutions) in it, so as to achieve good convergence. This endeavor requires calculating energies and bound-bound, bound-free, and free-free N-electron matrix elements [1,2,21,22,24].
One of the advantages of constructing expansion (3) in terms of wavefunctions that are calculated separately for each state, in terms of their own zero-order and correlation function spaces, is that in the process of solving, one can monitor and evaluate with transparency and economy the dependence of the evolution of Ψ(t) on each |m > and |ε > . In other words, unlike other methods that are often used for the approximate solution of Equation (1), the SSEA is not a "one-shot" algorithm that is applied in the same way to different cases. Instead, the methodology involves systematic expansion in terms of state-specific N-electron wavefunctions for the discrete and the continuous spectrum, according to the physics of each problem and the observed convergence of the calculation.
Furthermore, the SSEA is not limited to special electronic structures that are described by a single determinant, nor does it ignore the variety of effects of electron correlations, multiply excited resonances, and interchannel couplings.
Once the expansion coefficients of Equation (3) have been obtained from the solution of Equation (1), it is straightforward to use them, together with the stationary wavefunctions, in order to unambiguously compute the dynamics of interest as a function of time during and at the end of the pulsed interaction. This simple and direct way of quantitatively recognizing the correspondence of the stationary states with the field-induced dynamics is a hallmark of the SSEA. Its utility has been demonstrated in a variety of applications to difficult and computationally demanding TDMEPs, where, among other things, insight and experimentally relevant quantitative information has been gained by using the SSEA coefficients, as, for example, in cases of photoelectron angular distributions, above threshold ionization, the contribution of Rydberg states in multiphoton ionization, and the relative time delay in the emission of photoelectrons from different subshells (see [1][2][3][4][5][6] and Section 6 of this paper). The solution of the METDSE according to the SSEA takes into account the following: The contribution of perturbed or unperturbed Rydberg levels, and of doubly or even multiply and inner-hole excited states. * The contribution of the continuous spectrum of energy-normalized, channel-dependent scattering states, without or with resonance states. * The interchannel coupling.
It is important to add that the structure of the SSEA allows the use of relativistic wavefunctions as well.

Molecules
As regards molecules, the SSEA has been implemented only for the purpose of treating special themes involving vibrational states and dissociation continua in diatomics, where the smooth continuum has been described in terms of either box-normalized basis sets [1,25] or numerical energy-normalized scattering states [26]. It has not yet been applied to problems of time-resolved electron dynamics in molecules, for which descriptive discussions can be found (e.g., in [27]). Given our limited human resources, our prospects of tackling such problems in molecules at the level of rigor that has been demonstrated for atoms via the SSEA are meager. However, the following comment is relevant.
Tracking the hyperfast motions of electrons in molecules while discerning their effects from those caused by the other degrees of freedom is not easily amenable to rigorous computational analysis by accurately solving the METDSE, especially when the molecule is excited by a pulse of relatively short wavelength to a multitude of interacting states and channels. Of course, approximate treatments can, in principle, produce rough answers. As regards possible SSEA-type approaches, where the inclusion of excited states must be done carefully and economically in terms of separately optimized function spaces, the area is essentially open. Here, one fundamental difficulty lies in the proper and reliable calculation of electronic scattering wavefunctions. For diatomics, incorporating such wavefunctions into the SSEA expansion can, in principle, be done as a symmetry-adapted superposition of orbital configurations containing the scattering orbital obtained in the potential of the (N-1)-electron ion. Work and references regarding the calculation of the scattering orbital can be found, e.g., in [28]. However, significant progress is still needed toward the goal of treating polyatomics, with respect to both the electronic excited discrete and resonance states and the accurate representation and use of multichannel continua for different geometries of potential energy surfaces.

The Meaning of the Term "Electron Correlations"
As the title of our paper [14] shows, the concept of electron correlation was at the core of our initial considerations and proposals aiming to connect the prospects of observable effects of electron dynamics to spectroscopy with attosecond resolution. At that time (2001), published statements regarding the possible connection of attosecond pulses to electron dynamics had not yet identified specific concepts and effects whose interrogation would require attosecond resolution. The general argument regarding the relevance of attosecond pulses simply invoked the classical picture of a 1 s Bohr electron of hydrogen going around a circle with a period of about 150 as.
In order to explore the prospects of practical spectroscopic use of attosecond pulses for issues of time-resolved effects of electron dynamics on spectra, we identified as suitable candidates the doubly excited states (DESs) in the continuous spectrum of helium [1][2][3]14,15]. Helium is attractive to both theoreticians and experimentalists for practical reasons. By now, there are a number of publications in the area of attophysics with results on time-resolved spectra involving DESs of helium.
In recent years, in the growing literature on themes of the nascent field of attosecond physics [29], the term electron correlation is often used to describe a feature that is inextricably linked to electron dynamics studied with high temporal resolution. In order to avoid misunderstanding as to the meaning of this term, which is used in our work, I comment as follows: The term is not the same as "electron interactions", whose meaning is related to the obvious fact that there are two-electron operators of the atomic (molecular) Hamiltonian, nonrelativistic or relativistic. In other words, if, say, a calculation is carried out in the absence of the Coulomb operator 1/r 12 , the result represents an artificial situation. In our theoretical and computational work on the electronic structure and properties of atoms and molecules, electron correlation is defined and computed in terms of wave functions, as outlined below.
The term electron correlation was coined by Wigner in his theory of the electron gas for metals in 1934 [30]. It was later picked up and first used in the literature of quantum chemistry in 1952 by Taylor and Parr [31], in their configuration-interaction study of the 1s 2 1 S state of helium. They defined it as the part of the exact solution of Equation (2) that remains after calculating the single-configuration HF wave function and energy. The same definition was later emphasized by Löwdin [32] in his review article on electron correlation, which dealt with solving the eigenvalue of Equation (2) for ground states.
The above definition is translated to the following standard equations of quantum chemistry: where Ψ and E satisfy Equation (2). Equations (4) are valid for those states for which the HF equations are solvable. Furthermore, they are useful for those states for which the HF wave function is overwhelmingly dominant in the zero-order description of the exact eigenfunction. This normally occurs in many ground states of atoms and many closed-shell molecules in equilibrium. Therefore, in all those cases, and especially for ground states, the MEP becomes equivalent to the electron correlation problem defined by Equation (4).
When it comes to complicated cases of electronic structures with open subshells or of excited states, including resonance (autoionizing) states, where the mixing coefficients of a couple of configurations are large, the description of electron correlation in terms of Equation (4) has to be modified. Now, it is more meaningful and much more economical to compute the zero-order wave function as a self-consistent multiconfigurational expansion based on the state-specific "Fermi-sea", and to proceed with calculating the most important remaining electron correlation terms for the property of interest, according to the state-and property-specific (SPS) theory [21,22].
The concept of the "Fermi-sea" was introduced in the early 1970s in order to describe, in zero-order, classes of electronic structures with heavy configurational mixings, and to allow the development of an efficient theory of electron correlation for ground and excited states, and of field-induced properties and processes. Its nature, its background, and its utility are discussed in [21,22].
Accordingly, the form in which the state-specific stationary wave functions, Ψ n , are defined, computed, and used, is [21,22]: Ψ 0 n is, in general, a self-consistently optimized multiconfigurational wavefunction containing the Fermi-sea configurations, and X i n are the remaining symmetry-adapted correlation configurations containing either one, two, or three, etc., virtual excitations. In practice, the virtual (correlation) orbitals are represented by analytic orbitals that are optimized variationally.
The SPS criterion for an economic calculation is that a 0 must be very close to unity. The number and nature of the correlation terms in the sum of Equation (5) then depend on the property of interest and on the level of accuracy required for meaningful comparison with the experiment in [21,22].
Here I should add that the wavefunctions of electronic excited states are often highly mixed due mainly to valence-Rydberg-continuum state interactions. Therefore, their treatment requires a special advanced theory of perturbed spectra and channel mixings, which is applicable (and has been applied) to arbitrary electronic structures [1,2,21,22,33,34].
The main argument of the SPS theory [21,22] is that in order to accurately compute the matrix element of an operator that represents a particular property other than the total energy, what matters most is the "optimal matching" of the necessarily approximate wavefunction with the characteristics of the operator, and not whether the total energy has been obtained to a very good approximation.
Of course, for some operators the degree of convergence of the two calculations (first the total energy and then another property) is normally in harmony, especially for ground states. In general, however, the aforementioned optimal matching is not guaranteed by the criterion of the accuracy of the total energy, since non-negligible contributions to the value of the property may come from components of the wavefunction that do not contribute significantly to the total energy. These are configurations that have small coefficients, whilst their orbital radials are evidently not optimized with respect to the matrix element of interest. Yet, they are responsible for cumulative effects involving constructive or destructive contributions, with quantitative consequences.
The preceding observations suggest that, in general, when it comes to calculating properties and phenomena, the treatment of the omnipresent MEP need not rely exclusively on huge calculations that first obtain a very accurate total energy, especially since these are most often impractical for problems involving field-induced nonstationary states. Instead, in addition to the basics of the energy criterion, it is rewarding to also focus on and develop an understanding of the characteristics of the N-electron wavefunction(s) with respect to the state-specific self-consistent fields and the corresponding state-specific electron correlations [3,21,22]. For this purpose, it is essential for the theory not only to be formally correct and general, applicable to states of all kinds, but also to be implementable in terms of adequately optimized and numerically accurate function spaces [3,21,22].

Other Computational Schemes for Solving the METDSE for Problems of Atomic (Molecular) States Interacting with Ultrashort and/or Strong Pulses
As was already stated in the previous sections, the aim of the present paper is to present elements of the theory and methods of the SSEA, and of quantitative results that the Athens team has produced in connection with problems that fall in the domain of attophysics. Nevertheless, it may be of interest to the reader to cite other methods that have been proposed and used to solve the METDSE. Much of the text that follows is taken from [3].
As new experimental possibilities for the study of phenomena resulting from the interaction of atomic and molecular states with short and/or strong pulses were emerging, the introduction in 1994 of the SSEA [1,2] was linked to the need at that time to develop theory and many-electron methods of computation that are suitable for nonperturbatively solving the METDSE, beyond the cases of one-and two-electron ground states interacting with ultrashort and/or strong pulses. Regardless of the type of physical problem, in pursuing such a goal, the structure of the formalism must allow, in a practical way, the calculation of effects due to the state-specific features of the pertinent electronic structures and electron correlations. This requirement is analogous to, but much more complex and demanding than, that of many-electron problems of quantum chemistry on the energy axis.
Subsequent publications by a number of other groups have also presented methods aiming to provide solutions of the METDSE for many-electron systems. A recent review of a number of such methods was given by Ishikawa and Sato [35]. Its subsections are titled "Direct solution of the TDSE", "Multiconfiguration self-consistent-field (MCSCF) approach", "Time-dependent configuration interaction (TDCI) method", "Time-dependent multiconfiguration Hartree-Fock (TDMCHF) method", "Time-dependent complete active-space self-consistent field (TD-CASSCF) method", "Time-dependent coupled-cluster method", and "R-matrix approach". In a note added in proof, the authors stated that their review did not cover the SSEA.
I comment briefly on a few of these methods, starting with the time-dependent R-matrix (TDRM) theory [36][37][38][39]. The TDRM was first outlined by Burke and Burke in 1997 [36]. The fact that it was then applied only to a toy model, multiphoton ionization of a one-dimensional model, is testimony to the difficulties and high demands that such calculations have when it comes to problems involving real polyelectronic atoms or molecules. In a recent publication [39], the TDRM method was applied to a novel problem that is discussed in Section 6, namely the determination of the relative time delay in the (2s,2p) photoionization of neon.
Given the experience with advanced electron correlation methods for the ground and low-lying discrete stationary states and the convenience of available computer programs of computational quantum chemistry that are applicable to solving the METISE, in recent years many research groups have entered the area of solving the METDSE for certain problems by converting known time-independent methods into time-dependent ones in terms of expansions with time-dependent coefficients. For example, one is the TDMCHF (e.g., [40][41][42]), and another is the TDCI (e.g., [43,44]).
Even when the initial-state electronic structure is favorable for their formal implementation (i.e., single-determinantal closed-shell ground states), methods such as those in [40][41][42][43][44][45] express the content of the SSEA only in an approximate way. However, as I already mentioned in Section 3.2, for the time being approximate treatments may indeed be necessary for the study of molecules. They are constrained by the necessary truncation of the basis function space and by the corresponding uncertainties and inaccuracies in descriptions of the excited states in the discrete and, especially, the continuous spectrum, including Rydberg series (unperturbed or perturbed) and multiply excited resonances. In other words, when it comes to excited states, CI methods using a common basis set produce excited roots that need not constitute a good representation of excited electronic states of the discrete spectrum, or of the multichannel continuum, especially if the function space of virtual excitations is truncated significantly.
For example, consider the "TDCI singles" method described in [44]. It is formulated as an expansion over a set of virtual excitations from a zero-order Hartree-Fock determinant, which includes only one-electron excitations. Yet, even for simple closed-shell initial states (let alone open-shell, heavily mixed wavefunctions), ignoring the presence of doubly excited states, either in the physically significant portion of the spectrum or in the contribution of the character of the excited-state wavefunctions, introduces into the method and into the corresponding results a drastic and uncontrollable uncertainty. Even in the simple case of one-photon transition amplitudes between discrete states, it has been known since the early 1970s that, although the matter-field interaction is represented by a one-electron operator, orbital relaxation (upon excitation) and electron correlations imply that both singly and doubly (at least) excited configurations must be included in the CI wavefunctions; see the discussion on the first-order theory of oscillator strengths (FOTOS) in [21,22]. In fact, this is mandatory when valence-Rydberg-scattering state mixing is important in a portion of the spectrum that may contribute to the time-dependent dynamics.
In general, for strong fields and nonperturbative situations, the degree of reliability of CI-type approximations to the SSEA will be limited, regardless of the availability of hugely increased computer power, unless the contribution of the multichannel continuous spectrum is taken into account consistently and in a way that is mathematically justifiable. Here, one must remember that when a molecule is subjected to strong fields, channels of dissociation may compete with those of ionization.
Finally, I refer to coupled-channel methods such as those published in [46,47]. As in the SSEA, where the coupled equations can reach huge numbers due to the energy-discretization of the multichannel continuum in terms of energy-normalized scattering states (e.g., [6]), coupled-channel methods such as those of [46,47] also have to deal with such issues. In this context, I note that the recent work of Majety, Zielinski, and Scrinzi [47] deals with the four-electron beryllium, for which they stated that a multichannel description is needed. It is worth pointing out that the four-electron, multiconfigurational, multichannel problem was first solved about 25 years ago in terms of the SSEA [1]. It had to do with the interaction of laser pulses with the ground state of Li − , i.e., with a state of a negative ion where near-degeneracy effects, i.e., 2s 2 − 2p 2 , must be included in zero-order, as in the case of Be. The demonstration of the theory involved the continuous spectrum with two open channels, whose thresholds are the Li 1s 2 2s 2 S and 1s 2 2p 2 P o states.

Time-Resolved Many-Electron Physics of Hyperfast Electron Rearrangements
In recent publications by Kaldun and Blättermann et al. [7,8], by Krausz [9], and by Isinger et al. [10], there is discussion and presentation of results of new experiments on time-resolved hyperfast electron dynamics in many-electron systems, made possible through the application of modern techniques of spectroscopy with attosecond resolution. These measurements and the corresponding implications are directly related to, and in harmony with, our proposals and theoretical results that were produced in the '00s, based on time-dependent formalisms that allow the practical use of many-electron wavefunctions [1][2][3][4][5][6][11][12][13][14][15][16][17]. In the following paragraphs, I recall how these compare with the results and discussions in [7][8][9][10][11].

Time-Resolved Coherent Excitation and Decay of Strongly Correlated Autoionizing States
Over many decades, information on properties of resonance (autoionizing) states inside the continuous electronic spectra of real N-electron atoms and molecules has been obtained, to various degrees of accuracy, on the energy axis (e.g., energy positions, partial and total energy widths, photoionization cross-section in the energy region where a resonance is present, perturbation by an external field). This area of research has seen significant experimental and theoretical progress. Regarding theory, the fundamentals of resonance-state theory are settled. However, when it comes to practical implementation for problems involving complex excited electronic structures and multichannel continua, what is needed in addition to the fundamentals is the possibility of reliable calculations. Indeed, this direction contains many possibilities for application to complex situations.
In the mid-1990s, a new direction in the investigation of N-electron autoionizing (resonance) states was initiated. It was based on the many-electron calculation of the N-electron Ψ(t) via the SSEA on the time axis [12,13], a possibility that provides a different type of insight into their properties.
For the treatment of autoionization, the Hamiltonian is field-free and time-independent, H. The initial state is an N-electron wavepacket, Ψ 0 , which is square-integrable at t = 0 and has energy E 0 = < Ψ 0 |H|Ψ 0 > . The main goal in [12,13] was to calculate the time-resolved survival probability against electron emission, |< Ψ 0 |Ψ(t) >| 2 , from which information such as the lifetime, or the energy width corresponding to exponential decay of the decaying state, can also be obtained. The temporal resolution with which |< Ψ 0 |Ψ(t) >| 2 is obtained can be refined down to the attosecond scale.
The prototypical cases that were studied fall into two categories of processes of dissipation into the electronic continuum: (1) One-electron tunneling from a core-excited shape resonance, He − 1s2p 2 4 P, according to . This decay process is due mainly to the state-specific self-consistent field changing from initial to final state, and can be (was) calculated as such.
(2) Two-electron rearrangement. This is the "mechanism" that normally dominates nonrelativistic or relativistic autoionization. For example, an ab initio calculation produced the time-resolved decay curve and the lifetime, 3.5 × 10 −14 s, of the doubly excited state Ca KL3p 6 3d5p 3 F o [13].
The results provided, for the first time, quantitative results for real N-electron systems regarding the time-resolved regimes of both the exponential decay and, most importantly, the nonexponential decay (NED) [12,13]. I note that the states were not chosen haphazardly. Instead, their choice was based on the theoretical results and proposal published in 1977 [48], according to which, for NED of an isolated decaying state to acquire a magnitude of possible physical significance, the energy, E, of the state must be very close to threshold and the ratio R = E/Γ must be close to unity. (Γ is the width of the isolated narrow resonance state.) For a recent application of the results of [48] to radioactivity and carbon dating, see [49].
The insight and computational experience gained from the ab initio calculation of the time-dependent wave functions of autoionizing states, including the effects of electron correlation and of their time-resolved decay [12,13], served as input for the theoretical and computational work we initiated immediately after the important publications by Hentschel et al. [50] and Paul et al. [51] concerning the generation of attosecond pulses and trains of attosecond pulses.
Specifically, in 2002, we published the proposal-which was supported by solutions of the METDSE from first principles-that the then emerging prospects of new electronic spectroscopic techniques that can interrogate electron dynamics on the attosecond time scale could find fertile ground in the study of strongly correlated excited states in the continuum and the ensuing electron rearrangements [14,15]. It was pointed out that, whereas applying femtosecond spectroscopy provides information about nuclear motion, where the "motion" is described classically in terms of positions on potential energy surfaces, in the case of electrons in atoms or molecules, the time-resolution of electronic motion on the attosecond scale can be described in terms of properties of N-electron wave functions and their superpositions.
The understanding and treatment of such dynamics engages the following concepts and computable quantities:

*
The main symmetry-adapted orbital configurations obtained in state-specific self-consistent fields, the corresponding physically relevant transition amplitudes, and the "geometries" of the electron distributions in correlating configurations, and * The dominant state-specific electron correlations as determined by the degree of configurational mixing [3,14].
As we wrote in [15], "The problem was formulated around the following two questions, the first one being more general.
"(1) Is it possible to probe the electronic structure in such a way as to extract from it meaningful information on dynamics that can be resolved within the attosecond regime? If so, what is an observable process that can reflect the concept of 'motion' of electrons?
(2) Considering the first question, given the quantum mechanical formalism of stationary states and of corresponding eigenfunctions of 3N coordinates, how are we to probe and reveal the possible connection of electronic structure to atomic dimensions, to hypershort time scales and to electron correlation?" ( [15], p. 2).
The case for which the first calculations were carried out was the coherent preparation and decay of the strongly correlated low-lying 1 P o DESs of helium. Their theoretical preparation was effected by letting two femtosecond pulses interact either with the ground state, 1s 2 1 S, or with the 1s2s 1 S metastable state, so as to simultaneously excite two 1 P o states, labeled by 2s2p and 2p3d configurations and undergoing hyperfast "electron correlation beats": He 1s 2 1 S or, 1s2s 1 S femtosecond pulse hν 1 → femtosecond pulse hν 2 He 2s2p + 1sεp) 1 P o He 2p3d + 1sεp) 1 P o .
This system was chosen so as to define "bright", "dark", "gray", and "bath" electron configurations with respect to excitation, all playing roles in the overall physics [14].
I note that, according to the SPS theory of resonance states [21], the labels "2s2p" and "2p3d" correspond to two-electron wavefunctions that include the localized electron correlations to a high degree of accuracy.
One of our proposals, which was supported by quantitative results, was that hyperfast electron rearrangements can be recognized in terms of time-dependent geometries of the electron distributions as the dominant orbital configurations labeling the states correlate [3,14,15] (e.g., see [15], Figure 6). Evidently, a similar picture can exist in molecules, where one may expect that, tracked on attosecond scales, the main features of the electron distribution undergo some degree of oscillation from one nucleus to another.
In summarizing the work of [14], we wrote: "The determination of the period of population transfer between the two main configurations is one of the crucial results of the present theory, for it suggests that a laser-generated probe in the attosecond range, which will excite (de-excite) the bright or dark configuration state to a higher (lower) state that can be registered, will allow the recording of the motion of two electrons as they occupy different configurations at different times. We note that the geometry of such configurations can be calculated in terms of appropriate mean values and of angular and radial probability densities." ( [14], p. L276).
The possibility of time-resolving, on hyperfast scales, the consequences of electron correlations and rearrangements in atoms and molecules is now an area of research that attracts interest both theoretically and experimentally. For example, the significance of the possibility of determining results of orbital rearrangements in molecules in future measurements with attosecond resolution was argued by Krausz in a recent review [9]. I quote from Section 6 of [9], which is titled "Future prospects: exploring electronic motions in complex systems": However, there is another observable of crucial importance: the expectation value of the electrons' spatial coordinates. Dynamic changes in the electrons' atomic-scale distribution ... atomic-scale rearrangements of electron charge distributions occurring on the electronic timescale remain inaccessible in complex systems, such as biological molecules and solid-state nanostructures, to this day. Making them perceivable requires attosecond-duration multi-keV photon (X-ray) or electron pulses, permitting diffraction imaging with picometer spatial and attosecond temporal resolution.
From the point of view of theory, fulfilling this prospect will require quantitatively time-resolving the occupation of orbital configurations in different geometries, analogous to the results for intra-atomic electron rearrangements [15].

Photoionization with Attosecond Resolution. Measurements on Neon and on Helium Agree with Theoretical Predictions
I now turn to two new phenomena that were originally revealed during the '00s in studies of photoionization of helium and neon. Thus far, their understanding has been possible not only descriptively, but also quantitatively, thanks to achievements at both the experimental and the theoretical levels.
As regards experiment, the development and implementation of novel pump-probe techniques have allowed the measurement of quantities that describe time-dependent electron dynamics with attosecond resolution. As regards theory, the same quantities were understood and predicted accurately, thanks to the development and implementation of formalisms that allow the calculation and use of the Ψ(t) that is appropriate for the solution of the problem. The solution of the METDSE has been obtained either nonperturbatively via the SSEA, or to first order in perturbation theory. In the latter case, we showed that the results of interest can be obtained from an analytic formula that uses N-electron wavefunctions and matrix elements [16].
Only a description of the essentials is given here. The main purpose is to point out that the experimental findings on helium [7,8] and neon [10,11] are in agreement with our theoretical results of 2007 [16] and 2010 [11] (see also the conclusion of [17]). Complete information can be found in the original publications [7,8,10,11,16], including the supporting online material of [11]. The same themes were also briefly reviewed in [3], which is based on an invited talk at the 2015 Sanibel conference.

Time-Resolved Generation of the Photoionization Cross-Section Profiles of Doubly Excited and Inner-Hole Autoionizing N-Electron States on Hypershort Time Scales
In the mid-'00s, in continuation of the initial work on attosecond electron dynamics [14,15] discussed in Section 5, we completed two projects along the same direction. The corresponding publications of 2007 [6,16] contain theory and quantitative predictions pertaining to hyperfast electron processes of excitation and decay of autoionizing states, initiated by the interaction of femtosecond EUV pulses with the closed-shell He 1s 2 1 S state [16] and the open-shell Al KL3s 2 3p 2 P o state [6].
In the aluminum case, photoionization by the pulsed interaction creates the 2s-hole autoionizing states (1s 2 2s2p 6 3s 2 3p) 3 P o and 1 P o . Because of the open valence shell, the different electron couplings result in different self-consistent fields and electron correlations for the two states. Hence, the dynamics are channel-dependent. For reasons of economy, I will not discuss this case here. The reader is referred to the original paper [6] and the review [2].
In the helium case, the pulsed interaction causes the two-electron excitation above the ionization threshold, He 1s 2 1 S → (2s2p) 1 P o (A).
The label (2s2p) 1 P o represents the zero-order bound configuration of the resonance state inside the 1sεp 1 P o continuous spectrum at about 60.1 eV above the ground state. The energy-dependent stationary description of the resonance is a superposition of the correlated localized component corresponding to (2s2p) 1 P o , with the energy-normalized scattering components 1sεp 1 P o .
The helium project extended the theory and calculations of [14] in order to rigorously treat, from first principles, a novel problem of time-dependent many-electron physics, namely tracking the creation of a resonance state (scheme A) as this is recorded by the time-dependent appearance of the cross-section resonance profile inside the continuous spectrum. The much more demanding work on Al [6] also produced results on this phenomenon.
The work was inspired by the 2005 publication of Wickenhauser et al. [52] who, using a model system that assumes data from (super) Coster-Kronig transitions with lifetimes of ∼ 400 attoseconds, "investigated the feasibility of observing the buildup of a Fano resonance in the time domain by attosecond streaking techniques" ( [52], abstract), and produced results on the phenomenology of such a possibility.
The theory and calculations in [16] were developed and carried out within two frameworks. One involved nonperturbatively solving the METDSE via the SSEA, which is valid for weak as well as strong fields. The other was based on a new time-dependent approach, whereby the METDSE is solved to first order in perturbation theory, while employing the Fano description of resonance states on the real energy axis. In both cases, the calculation is done ab initio, using accurate wavefunctions and matrix elements.
The results were obtained for a sin 2 pulse of duration 450 a.u (11 fs) and central energy of about 60.1 eV. (The lifetime of the He (2s2p) 1 P o resonance is about 17 fs.) They were depicted as snapshots of P(ε, t) ≡|< scattering state(ε)|Ψ SSEA (t) >| 2 , where ε is the energy above the ionization threshold, at various femtosecond intervals. They showed quantitatively how, as the two electrons in the bound orbitals 2s and 2p are allowed to interact with the scattering continuum, the interference of the bound and scattering components evolves until the formation of the resonance asymmetric profile on the energy axis is completed. Upon completion of this profile, the agreement with the experimental values and the earlier theoretical results on the energy axis (energy position, width, and the Fano q parameter for the profile) is excellent [16].
The aforementioned analytic treatment, which is valid for weak fields only, produced a practical analytic formula (corrected for a misprint in [16]) that uses N-electron matrix elements and incorporates the pulse characteristics. It was demonstrated in [16] that when the field is weak, the numerical results from applying this formula are identical to those from the nonperturbative SSEA calculation ( [16], Figure 2). The formula can be used for appropriate predictions in many-electron atoms, provided its parameters are known from computation or experiment. For example, its application shows that Auger states of atoms in the middle of the periodic table, whose width is sufficiently large, about 2 eV or more, are suitable candidates for observing the formation of their stationary profiles with attosecond duration.
The predictions for Al have not yet been tested experimentally. However, those for helium were indeed subjected to a novel pump-probe experiment with a time delay of sub-femtosecond precision [7,8]. The findings and analysis of [7,8] confirmed the theory of [16], although the experiment did not cover the full duration of the generation of the complete resonance profile of the (2s2p) 1 P o state on the energy axis, which was calculated to be about 180 femtoseconds (i.e., about 10 lifetimes) [16].
Given the continuing experimental progress on refined measurements on hypershort time scales, we expect that novel experiments such as those of [7,8] (see also [53]) can be developed further, in order to allow the study of even more intriguing time-resolved processes involving multiply excited or inner-hole-excited states inside the continuous spectrum. Such states are normally unstable (autoionizing). However, due to symmetry restrictions, it is also possible for a few of them to belong to the nonrelativistic discrete spectrum in spite of their very high excitation. For example, as a prototypical case we have treated aspects of the time-dependent dynamics of the preparation of the bound triply excited state of a negative ion, He − 2p 3 4 S o , which is 59.33 eV above the helium ground state, using two pulses with time delay [54]. The wavelengths of the two pulses resonantly connecting the states are in the infrared (10,080 Å) and soft x-ray (323 Å) regions. The ladder of excitation is [54], where 1s2s2p 4 P 0 is metastable and 4 P(ε) > is a shape resonance labeled by the zero-order (1s2p 2 ) 4 P.

Introduction
In 2008, experimentalists at the department of attosecond physics of the Max Planck Institute of Quantum Optics in Garching, Germany, headed by F. Krausz, were in the process of carrying out a spectacular pump-probe experiment of fundamental significance. This involved the implementation of their attosecond streaking technique for the photoionization of neon. Their goal was to measure whether, upon photoemission of the electrons from the 2s and 2p subshells of neon by an attosecond pulse with photon energy centered around 106 eV, the electrons are ejected absolutely simultaneously, as the conventional wisdom was at the time, or with a time difference, i.e., with a relative time delay (TD).
In view of that experimental project, V. Yakovlev of the Garching institute and the Athens team started a collaboration with the purpose of developing a many-electron theory that could determine whether such a TD could indeed occur in neon photoionization, and, if so, what its sign and magnitude ought to be. This possibility would be unique in the study of atomic physics, since there are no other collision processes (e.g., other atoms) that could interfere.
At that time we were in the dark, since this new problem of time-dependent many-electron quantum mechanics had never been considered and understood theoretically, let alone solved quantitatively. The only certainty was that if this were a real phenomenon, the corresponding information ought to be contained in the METDSE describing the synchronized pump (EUV attosecond pulse) plus probe (IR laser pulse) experiment, in terms of the Hamiltonian H tot (t) = H Ne + V int (t). H Ne is the Hamiltonian of the field-free neon and V int (t) = V IR (t) + V EUV (t − ∆t), where ∆t is the controlled time delay between the application of the two pulses in the attosecond streaking experiment [11].
The theoretical treatment and results were published in [11], including its SOM (see also [17] and below), together with the experimental findings of Schultze et al. [11]. Both the experiment and the theory introduced into atomic and molecular physics a new and observable concept in attosecond electron dynamics-namely the relative TD in the photoemission from atoms-and methods in terms of which quantitative results can be obtained.
After the 2010 publication [11], a large number of publications, experimental as well as theoretical, have dealt with aspects of this problem, not only for neon but also for other cases. This body of work was recently reviewed comprehensively by Pazourek, Nagele, and Burgdörfer [55]. Section V of [55], titled "Time-resolved photoionization of many-electron atoms", contains a discussion related to the seminal work on neon. As regards theory, even time-independent treatments have been implemented, reducing the understanding of a phenomenon where the concept of time is inevitably involved (e.g., quantum mechanical time of free N-electron wavepacket formation), to the calculation of photoionization Eisenbud-Wigner-Smith (EWS) time delays in terms of energy derivatives of scattering phases.

Elements of the Theoretical Model and Results of Schultze et al. on Neon
The (2s,2p) photoionization of neon by the attosecond pump pulse centered at about 106 eV with a width of about 14 eV involves the dipole transitions 2s → εp, 2p → (εs,εd) and 1 P o final states in the continuum. This picture excludes two-electron excitations. In the presence of the probe V IR (t), states with many angular momenta in the continuum are mixed. These states are labeled as ψ ,L 2s (ε) = (1s 2 2s2p 6 )ε 1 L and ψ ,L 2p (ε) = (1s 2 2s 2 2p 5 )ε 1 L . If we were to use the SSEA with H tot (t), the corresponding Ψ SSEA (t) would be constructed in the following form: (6) According to the SSEA methodology, the upper limits of the above summations, max and max , as well as the upper limits of the integrals ε max , are increased until full convergence of the time-dependent complex coefficients is achieved. The bound wavefunctions of the ground state, Φ gs ( 1 S), and the ionic cores are computed by including the few single and double orbital electron correlations, which contribute according to the arguments and methods of the SPS theory.
The complex coefficients of the SSEA are probability amplitudes with physical meaning, containing quantum mechanical information about transitions to each energy-normalized scattering state to which they correspond. By asserting that photoionization is a half-collision process, it is reasonable to invoke the concept of time-delay developed in scattering theory and associated with the names Eisenbud, Wigner, and Smith (EWS) [55][56][57]. Accordingly, the time delay for each channel is obtainable from d dε arg[b 1,1 2s (ε, t)] and d dε arg[b 2,1 2p (ε, t)] in the energy region fixed by the EUV excitation energy. For neon, as well as for other systems, a quantitative understanding of the degree of sensitivity of these quantities (energy derivatives of scattering phases) on the wavefunctions and the resulting complex coefficients is a crucial element of theoretical work in this new area of research.
The theoretical treatment in [11] did not apply the SSEA with the Hamiltonian containing both EUV and IR interactions, since this was computationally prohibitive for our computer facilities. Instead, the problem was separated into two parts, which were put together by Yakovlev and collaborators in Garching.
The second part was conceived and carried out in Garching. It involved implementing the one-electron Coulomb-Volkov approximation (CVA) to the estimation of the streaking delay caused by the interaction of the outgoing wavepackets with the IR laser pulse.
The first part had to do with the many-electron solution of the METDSE via the SSEA using only the EUV pump interaction, with the purpose of determining the relative TD for the one-photon processes, Ne 1s 2 2s 2 2p 6 1 S atto sec ond pulse → hν ≈ 106 eV Ne + 1s 2 2s2p 6 2 S + εp or, Ne + 1s 2 2s 2 2p 5 2 P o + εs, εd (C).
Systematic calculations using the parameters of the experimental attosecond pulse determined the effects on the SSEA scattering coefficients and on the corresponding TD, of the state-specific self-consistent fields, of the important state-specific electron correlations (according to the SPS theory [21,22]), of interchannel coupling, and of resonance states corresponding to two-electron excitations from the 2s subshell [3,11,17].
The theory and results of [11] (see below) corroborated the essence of the experimental findings. Specifically, it was demonstrated that upon photoabsorption, assumed to occur instantaneously, the emergence of the wavepackets (of different kinetic energies) representing the electrons initially occupying the neon 2p and 2s HF orbitals takes place with a relative TD. In agreement with the experiment, the 2p electron emerges a few attoseconds after the 2s electron does.
The SSEA result for the spectrally averaged relative TD of the electron wavepackets emitted according to scheme (C), including doubly excited resonances (see below) was 6.38 as. It was impressive that the same number (6.37 as) was obtained independently from the simulation of the streaking measurement in terms of the CVA ( [11], page 6 of the SOM).
On the other hand, the contribution of the IR-laser perturbation on the emitted electrons, deduced according to the CVA [11], was 2.3 as. (The text in [11] reads: "Thus, for the current experimental parameters, the small deviation between the exact electron's motion and the modeled via the CVA gives rise to a two-attosecond discrepancy in the relative delay. Accepting this small discrepancy, many-electron models were applied to investigate the effects of electron correlation" (p. 1661).
Therefore, as regards the simulation of the streaking measurements of [11], the theoretical prediction for the spectrally averaged relative TD for neon was This pump + probe number, or just the SSEA one, the latter representing what was called in [11] the "Wigner-Smith" TD for the process of photoionization, was in quantitative disagreement with the measurement of [11], which gave 21 ± 5 as.
In the next subsections, I focus on three themes related to [11]. The first has to do with the possible role of doubly excited resonances, whose existence at the time of the experiment was unknown. The second has to do with the standard question of the degree of agreement between the theoretical and experimental results. The third has to do with information that only theory can provide, namely, how to understand this phenomenon in neon in terms of the interplay between interelectronic interactions (i.e., wavefunctions) and dynamics (i.e., transition amplitudes).

The Relative Time Delay in the Neon Problem: The Role of Doubly Excited Resonances
The quantitative disagreement between theory and experiment persisted in the course of the investigations on the relative TD of neon in Athens and Garching. Both experimentalists and theoreticians were confident that their results were obtained correctly. In this context, the Athens team considered the possibility that double excitations reaching resonances that were unknown at the time might be responsible for at least some of the observed discrepancy.
Given the energetics and the nature of the electronic structure, we concluded that good candidates ought to be states that are created by double excitations from the 2s subshell [17]. Although we established the existence of a series of resonances with a double hole in the 2s 2 subshell, of which one or two have their positions at the experimentally accessible energy region of the photoelectron, the conclusion as to their influence was negative because of their small width in relation to the spectral averaging over the large width of the attosecond pulse (see the conclusion of [17]).
By inference, in other cases of photoionization, the relative TD might indeed be affected by the presence of a very wide scattering resonance, or of closely lying resonances, at the energy where the photoelectron is emitted. For example, the dependence of time delays on resonance states was examined recently by Sabbar  The quantitative discrepancy between the experimental and theoretical results reported in [11] evidently presented a scientific challenge. On page 465 of [3], I commented as follows: "The enigma stemming from the reported difference between the experimental and the theoretical findings in [11] remains to date, in spite of subsequent theoretical investigations by a number of researchers. It is challenging to both theory and experiment. In this context, it is worth noting that the same theoretical number, [to that of the SSEA], i.e., 6.4 as, was reported recently by Feist et al. [60] who carried out very large R-matrix type calculations, including a series of resonance states near and far from the critical region of 106 eV. Their results confirm the original theoretical predictions which were made in [11,17]." The "enigma" stated above was resolved recently with the publication by Isinger et al. [10] of new measurements with high spectral resolution, which allowed the disentanglement of the direct 2s ionization from the two-electron process of ionization plus excitation to a 3p state [10]. The qualitative conclusions of Isinger et al. [10] as to the sequence of emission of electrons from the 2s and 2p subshells of neon confirm the original experimental and theoretical predictions of Schultze et al. [11].
However, the new experimental number [10] is about 10 as, much closer to the theoretical result first predicted via the theory and the SSEA calculation in [11] than to the experimental value [11].
Here, a parenthetical clause is necessary. The authors of [10] present the discrepancy between the theoretical and experimental TD as a "puzzle" that came to the fore after calculations that followed the publication of [11] produced such a result. (The theoretical work in [10] used the many-body perturbation theory.) In fact, this was already known and reported explicitly, with specific numbers, in the seminal publication [11], which demonstrated, conceptually and quantitatively, the presence of the relative TD in neon, both experimentally and theoretically. (See the paragraph from [3] quoted above.) The calculations by other groups that appeared after the publication of [11], (quoted in [10]), simply confirmed its theoretical finding, which was obtained based on the analysis outlined in Section 6.2.2. It is gratifying to see that, once the phenomenon was understood and presented in [11], later large-scale calculations by Feist et al. [60], Moore et al. [39], and others, produced a number that is the same as, or very close to, the theoretical numbers (6.4 as and 8.7 as) reported in [11] "Contrary to our initial anticipation, the results for this system did not show a strong dependence on details of electron correlation beyond the MCHF level, or on interchannel coupling. For example, whereas the final result due to the pump step was 6.4 as, the result without interchannel coupling, i.e., using just the dipole matrix elements connecting the initial state-specific orbitals, calculated at the MCHF level, to their corresponding final state channel, was 4 as (see page 1661 of [11])." In other words, we found that the cause of the bulk of the relative (2s,2p) TD in neon is already produced at the level of state-specific HF-MCHF wavefunctions for initial and final states, even though the hole-filling electron correlation, 2p 2 − 2sd, normally important for transition processes is present only in the case of photoionization from the 2s subshell. d is a one-electron function of d symmetry representing the sum of the d orbitals of the discrete and the continuous spectrum calculated in the term-specific N-1-electron HF potential.
Nevertheless, the fact that the theoretical results (6.4 as, or 8.7 as with streaking) differed from the experimental value (21 ± 5 as) was later assessed in some publications as being caused by omitted electron correlations of some unknown type.
For example, soon after the publication of [11], Sansone [61], commenting on its results, wrote: "... more refined models including multi-electron correlation, such as the state-specific expansion approach. Even though these models turn out to be in qualitative agreement with the experimental observation (the photoionization from the 2p shell is delayed with respect to that from the 2s shell), they do not offer proper agreement with the experiment as the estimated delay of 5 as is three standard deviations smaller than the measured delay. The discrepancy clearly indicates that our understanding and capability of modelling multi-electron correlation effects is still far from being complete and that subtle effects leading to attosecond time delay cannot yet be satisfactorily described. In this contest, the experimental data represents a benchmark for more sophisticated and advanced models including electron correlations." In the above quotation, the "blame" for the discrepancy between theory and experiment that was reported in [11] was placed on the ostensible lack of "understanding and capability" of handling electron correlation effects.
Along a similar line of reasoning on the significance of electron correlations in determining the relative TD in neon is the statement made in the commentary by Day [62] on the important advances in physics of the past few years, where this TD is attributed to "subtle electron-electron correlations".
On the other hand, the authors of the recent publication [10] have stated that according to their calculations, the effect of electron correlation on the neon TD is rather small compared to the result based on the zero-order orbital approximation, thereby confirming our conclusion from the original theoretical work [11].
Of course, as in many other properties or phenomena, only when the important electron correlation effects are accounted for is the theoretical understanding reliable.
In conclusion, our calculations and analysis of the neon TD problem [11] showed that the phenomenon is due mainly to the different state-specific self-consistent fields, and that electron correlation in the bound states and interchannel coupling in the continuous spectrum act as correction factors. This is why, in the time-independent framework of calculating "Wigner time delays", algorithms such as the random phase approximation and limited many-body perturbation theory have produced similar numbers for the closed-shell neon 1 S state.
However, in general, for N-electron systems whose electronic structures and spectra are quite different than those of neon, the contributions to the energy derivatives of the phases of the complex coefficients of the scattering channel in the photoelectron energy of interest, and therefore to the relative TD, may indeed depend significantly on electron correlations and interchannel couplings, or on the presence of broad resonances.

Conclusions
Thanks to significant developments of the last two to three decades, both theory and experiment can now provide the same or complementary information and insight into attosecond-resolved electron dynamics associated with hyperfast electron rearrangements upon pulsed excitation of atoms and molecules into the continuum.
As regards experiment, the novel techniques that have been developed and used for the study of such phenomena are discussed in the cited works [7,8,10,11,29].
As regards theory and computation, the formal construction, the credible solution, and the quantitative understanding of the related new types of TDMEPs require that the METDSE is treated and solved from a many-electron point of view, for closed-as well as for open-shell atomic or molecular states, rather than in terms of models or one-electron potentials. To this purpose, a fundamental requirement is for theory to account for the details of the electronic structures and the scattering continua of the states participating in the dynamics.
The SSEA, whose main features were presented in Section 3, satisfies the above requirements [1][2][3][4][5][6][14][15][16]21,22]. Indeed, because of its structure, its implementation for each problem of interest can be easily understood, regardless of the wavelength(s) of the applied pulse(s), which can range from the IR to the x-ray region.
In Sections 4-6, I discussed various types of TDMEPs whose solutions, with attosecond resolution, have provided quantitative information and good insight into hyperfast electron dynamics and their relation to state-specific self-consistent fields, electron correlations, interchannel couplings, and resonances. The results were obtained in terms of many-electron wave functions Ψ(t), calculated by solving the METDSE nonperturbatively via the SSEA. In addition, when the pulse is chosen to be short and relatively weak, it is possible to formulate and calculate perturbatively the time-resolved preparation and evolution of states in the continuum in terms of practical analytic formulas containing information about the characteristics of the pulses and the time delay between two pulses [16,54,63].
The cases that were brought to attention as examples were: (1) The ab initio determination of the details of the time-resolved decay of autoionizing states via one-electron tunneling or two-electron rearrangements, thereby identifying quantitatively the regimes of nonexponential decay (NED) for states very close to threshold. I point out that the predictions for NED were made not only for the case of very long times [12,13] but also for the case of very short times after the preparation of the N-electron decaying state [64]. The degree to which the preparation step affects the NED for each particular N-electron isolated decaying state might be a challenging question for attosecond physics. (2) The time-resolved coherent preparation and decay of strongly correlated doubly excited and inner-hole excited autoionizing states. (Prototypical cases of electron dynamics in helium and aluminum.) (3) The time-resolved identification of electron correlation beats and correlation-induced geometrical rearrangements of electrons inside the atom, which occur on the attosecond scale. For molecules, such rearrangements may take also place interatomically. (4) The time-resolved generation of the energy profiles of doubly excited (helium) and inner-hole excited (aluminum) autoionizing states during photoionization by hypershort radiation pulses. The published figures explicitly and quantitatively depicted for the first time the time-dependent effects on these profiles of the interference between the localized and scattering components of these states. The information that has been obtained from the results of our treatment of these prototypical TDMEPs agrees with what has been produced via novel experiments capable of attosecond resolution [7,8,10,11].
In accordance with the original proposal and results [14,15], in all the above examples the hyperfast imaging of electron dynamics reveals the significance of both the state-specific self-consistent orbitals and their dominant electron correlations.

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