# Theory of ultrafast autoionization dynamics of Fano resonances

###### Abstract

We study atomic autoionization processes in the time domain. With the emerging attosecond extreme vacuum ultraviolet and soft x-ray pulses, we first address how to characterize the time evolution of the decay of a discrete state into a degenerate continuum. A short pump beam generates a number of resonance states in a series and the nearby background continuum, and the resultant wave packet evolves with time until the full decay of the bound states. Taking the resonance series embedded in the continuum in beryllium atom as an example, the time evolution of the autoionizing wave packet in energy domain and in coordinate space is calculated and analyzed, where Fano profiles build up in the photoelectron energy during the process. A proposed pump-probe scheme assumes that the probe beam ionizes the inner electron in the wave packet. The lifetimes of the resonances and the photoelectron energy distribution can be obtained from the ionization yield versus the time delay of the probe.

###### pacs:

32.80.Zb,32.80.Fb## I Introduction

Scattering processes – photoionization, electron collision, etc. – have been important tools for studying physical systems since the early days of the development of atomic theories. Autoionization, first observed by Beutler beutler in the photoabsorption of rare gas atoms, is a phenomenon where the energy spectra is characterized by asymmetric resonance peaks, which were explained by Fano in his seminal paper in 1961 fano . In this paper, asymmetric resonance in photoabsorption occurs in the energy region where a discrete state is embedded in the continuum. Photoabsorption populates both the discrete state and the continuum states. The interaction of the discrete state with the continuum results in the ”autoionization” of the former, where the discrete state interferes with the directly populated continuum states to form asymmetric Beutler-Fano profile. Today Fano resonances have been used to describe a wealth of physical phenomena in atomic, molecular and optical systems, condensed matter systemswick14 ; wick15 , nanostructuresmirosh as well as many other fields. Essential to Fano’s theory is that each resonance can be characterized by three parameters, the ”unshifted” resonance position , the shape parameter and the resonance width . The resonance width characterizes the decay time or the lifetime of the discrete state. Since the typical resonance width of an atomic or molecular system is in the order of 10 eV, which corresponds to the lifetime of as short as a few femtoseconds, almost all measurements on these resonances have been carried out in the energy domain, using high-resolution spectroscopy. Today a large set of experimental and theoretical data of these resonances have been collected, and they provide important structural information of the systems.

In the past decade, with the advent of laser technology, single attosecond pulses (SAP) are becoming availablekrausz . Such pulses opened up a new research platform for studying electronic dynamics, such as the autoionization dynamics. Indeed, the first SAP experiment in atomic system was to determine the lifetime of the Auger decay of an -shell hole in Kr created by extreme vacuum ultraviolet (XUV) attosecond pulse by Drescher et al drescher . For Auger decay, the resonance profile is mostly symmetric, and thus only the decay width or the lifetime of the inner-shell hole was determined. In Ref. drescher , the lifetime is determined by creating the inner-shell hole in the presence of an IR laser. By changing the time delay between the XUV pulse with the IR, photoelectron spectra can be used to ”deduce” the lifetime of the inner-shell hole – based on the approximate ”streaking” theorykitzler . Such streaking theory has been generalized to Fano resonances by Zhao and Lin zhao as well as by others wickenhauser ; smirnova . However, the photoelectron spectra from such XUV+IR setup are very complicated and only the lifetime of the Fano resonance can be extracted.

Due to the insufficient fluence of the SAP’s so far, existing measurements using attosecond pulses for Auger decay or autoionization have all been carried out using the XUV+IR pump-probe setup where the two pulses actually overlap in timedrescher ; gilbertson ; wang . Rigorously speaking, the IR is not probing the resonances created by the attosecond XUV pulse in these experiments since the photoelectron is created in the IR field. With improved technology in attosecond pulse generation and the recent advances in X-ray free-electron lasers (XFEL), we envision that intense attosecond pulses will become available in the near future, such that a true time-dependent investigation of the autoionization process will be possible. In this article, we analyze the time evolution a Fano resonance, or the autoionization dynamics of the system in the field-free condition before the discrete state fully decays into the continuum. For example, if the discrete state has a lifetime of 20 fs while the pump pulse is only 1.5 fs, we ask what is the time evolution of the wave packet associated with the decay process after the pump pulse, and how the energy profile eventually evolves into what is characterized by the Fano resonance shape.

We comment that the time evolution of a Fano resonance has been studied in several calculations by solving the time-dependent Schrödinger Equation of helium or a model atom mercouris ; lindroth ; wickenhauser ; cavalieri , with or without the IR probe. These studies did not address how to retrieve the time evolution of the resonance profiles. Similarly, autoionization of double Rydberg states in the time scale of picoseconds has been studiedjones . In the meanwhile, the interaction of a doubly excited state with a whole Rydberg series also falls into the similar categorygallagher1 ; gallagher2 . Understanding the time evolution of the two-electron wave packets and the issue of how to probe their evolution lie at the heart of understanding and controlling two-electron dynamics.

In this article, we use two pulses that are separate in time. The probe is always placed strictly after the pump so that the field-free evolution between the pulses can be defined. The physical process is understood roughly as the decaying of the bound state to the continuum due to the interaction between them. In Sec. II, starting with Fano theory which is for stationary states, we construct the time-dependent wave function of the system which describes the time evolution of the electron wave packet after the pump pulse. The wave function is recasted in the basis set of the bound state and the continuum states, where it contains the time information for the bound and the continuum components. The wave packet for the continuum part is shown to evolve to the typical Fano resonance profile at long time. We thus define the continuum wave packet at short time as the time-dependent Fano profile of the autoionization. Since the pump pulse is short, it inevitably spans a broad-band in the energy spectra. Assuming that the background energy spectrum generated by the pump is described by a Gaussian distribution, we generalize the theory to include all the resonances within the bandwidth. In Sec. III we apply our theory to the resonances in beryllium to investigate wave packet dynamics after the pump pulse. In Sec. IV we demonstrate how to probe the time-dependent Fano resonances, using the example in beryllium given in Sec. III. A short summary and future outlook is given in Sec. V.

## Ii Theory

### ii.1 Fano’s theory of resonance

In a multi-electron atomic system, the high-lying discrete states are embedded in the continuous spectrum if these discrete energies are higher than the binding energy. Both the bound and continuum states are described in terms of configurations. According to Ref. fano , with one bound state and its nearby continuum, the system is governed by the total Hamiltonian given by

(1) | ||||

(2) | ||||

(3) |

where and are bound and continuum configurations respectively. The off-diagonal terms represent the mixing strength between the bound and the continuum. The diagonal terms and are bound and continuum energies respectively.

The eigenstates of the atom are obtained by diagonalizing the Hamiltonian in Eq. 1-3 and are given by Fano as

(4) |

where the coefficients are

(5) | ||||

(6) | ||||

(7) |

Here we assume that is a real constant in the narrow energy region for each Fano resonance, provided that the bound and the continuum states are in real representation.

Originally Fano’s theory was applied to photoabsorption processes. The transition amplitude is where is the ground state, is the excited state of energy , and is the photon transition operator. The states and are eigenstates of the atomic Hamiltonian , so the coefficients do not change with time. In Ref. fano , the continuum wave functions are taken to be real standing waves. The absorption cross section is proportional to , which contains the features of Fano profiles. Until recently, all photoabsorption measurements are performed using light sources of hundreds of picoseconds or longer, where they can be considered as monochromatic. Thus, of a single is promoted in each measurement. To trace out a Fano resonance, monochromatic lights are tuned across the resonance while the cross section for each energy point is recorded.

Although Fano’s theory was applied to photoabsorption originally, it can be trivially generalized to study photoelectron angular distributions. While the continuum functions in Eq. 4 are taken to be real standing waves, is an eigenstate of the Hamiltonian. To obtain the scattering amplitude for an electron with momentum , one only needs to project the standing wave onto the momentum eigenstate of the photoelectrons, in energy-normalized form given by

(8) |

Here are the radial wave function of the electron taken as real standing waves as in Fano’s theory, are the total scattering phase shifts, and the ’s are spherical harmonic functions.

Following the notation in Ref. fano , the width of the resonance is defined by . Assuming that is energy-independent near the resonance energy , the -parameter is defined by

(9) |

Without the flat-background and the constant assumptions, will be replaced by the ”modified bound state” defined in Ref. fano which contains a term based on uneven and near . We assume this modification is negligible in general cases and consider the -parameter as defined in Eq. 9. Fano’s theory predicts the resonance profile in a simple form, where the associated parameters , , and can be extracted by fitting the form to the measured profile.

With the advent of available light pulses of a few femtoseconds or shorter, one can use such a short pulse to ionize the atom. After the pump pulse is over, the time-dependent wave function of the atom is given by

(10) |

where in a single pulse, the integral in Eq. 10 covers a bandwidth more than the width of a resonance. Such a resonance can be measured, long after the pump, in a high-resolution spectrometer. This amounts to projecting out the time-dependent wave function over a long time. Thus the whole Fano resonance is generated in a single pulse and its profile is given by .

We now consider the situation where the pump pulse is shorter than the lifetime of a Fano resonance. One would like to ”see” the dynamics of the autoionization. However, such information is not contained in since are eigenstates of the atomic Hamiltonian and the coefficients are constant when the pump field is absent. Thus in Eq. 10, the time information is hidden and not recognized, even though we ”know” that autoionization is happening!

### ii.2 Time-dependent wave function in configuration basis

It is convenient to define the excited wave function as the total wave function excluding the ground state part since the ground state does not participate in autoionization. Instead of using eigenstates as expressed in Eq. 10, we alternatively expand the excited wave function in configuration basis as

(11) |

Since and are not eigenstates of the atomic Hamiltonian, the coefficients and are time-dependent. One expects that at long time (compared to the lifetime), the coefficient will go to zero and all the information of the total wave function will be contained in . However, the time evolution of these coefficients cannot be measured in laboratory – they change in the time scale of femtoseconds or less, while laboratory electronics can only measure changes of the order of picoseconds or longer. Here we define as the ”time-dependent Fano resonance profile”, or in short, ”resonance profile”. Below we discuss how these coefficients behave. The measurement issues will be addressed later.

Assuming that the initial values of and are known, with the aid of Eq. 4-7, the time-dependent coefficients are solved and given by

(12) | ||||

(13) |

where the functions and are defined by

(14) | ||||

(15) |

The superscript denotes the initial values, i.e.,

(16) | ||||

(17) |

The solutions in Eq. 12-15 are calculated exactly for a given set of parameters , , , and , where these parameters can be either obtained by running a separate program such as Time-dependent Schrödinger Equation (TDSE) code or time-dependent perturbation calculation, based on the prior knowledge of the pump, or extracted from the experimental values in the time-integrated spectrum. In the latter case, with defined by Eq. 9 and rewritten as

(18) |

the initial coefficient can be extracted since the value of and the background magnitude of can be obtained from the Fano profile after long time (the phase of can be deduced if photoelectron angular distributions are measured). Equations 12-15 are the backbone of our model.

Note that the initial coefficients and are in general complex values while is real. Within the narrow range of a resonance width, the energy in is close to the bound state energy such that the phases of and are assumed to be the same in the neighborhood of the resonance, and the corresponding , representing the ratio between and , is real. When spans a wide energy range covering many resonances, its phase varies smoothly with , but in each resonance region, the statement above is still true, and each is real.

The decay lifetime of the bound state is defined as . Rigorously speaking, the dynamics of the autoionization is not only about the decay of the bound state, but also its coupling to the continuum (See Eq. 11). The wave packet can propagate from the continuum to the bound or the other way around, or between continuum states of different energies through the bound state. Thus, the lifetime , determined by the interaction matrix , characterizes the time scale of the whole event.

We can generalize Eq. 12-15 to include multiple resonances by applying Fano’s theory for many bound states fano and assuming that higher-order perturbation by is negligible. The excited wave function is

(19) |

with the coefficients given by

(20) | ||||

(21) |

where the label is the index for the th bound state. The -function and the -function, now subscripted by , are calculated independently for each resonance with the associated parameters as given by Eq. 14-15.

To simulate the effect of the finite range of while keeping the calculation simple, the profiles can be further adjusted by multiplying each ”resonance term” in Eq. 21 (in the bracket) by a Gaussian window function to limit the range of each resonance. For an isolated resonance, the effect of the window function is minimum or unnecessary because the window function mainly removes ”the wings” but not the center of the resonance, where the central part is the significant and informative part. On the other hand, when dealing with many resonances, without the window functions, the wings may be overly extended and unnaturally perturb the nearby resonances. In the latter case, multiplying the window functions is an effective and efficient solution.

#### ii.2.1 Flat initial continuum (IC) distributions

To study the qualitative behavior of the wave function, we first assume that the initial profile is flat, i.e. . This initial condition, dubbed IC, is valid near the resonance positions if the pump bandwidth is broad enough relative to the resonance widths so that the populated profile is nearly flat around . For conceptual presentation, only the single-resonance formulae Eq. 12-15 are referred here, while the idea is the same for the many-resonance case. Defining the scaled time by and scaled energy by , the bound and the continuum coefficients under IC are

(22) | ||||

(23) |

where in the scaled energy, q-parameter becomes

(24) |

and the factor disappears because is measured relative to the resonance energy. While the coefficients given by Eq. 22-23 represent the evolution for , their limits at are different from the initial coefficients and that we have defined at . This discontinuity means an immediate jump of the wave function right after the evolution starts. This unphysical behavior is due to the integral of -function in Eq. 14 and -function in Eq. 15 over infinite energy range at very small times, where the contributions from much higher or lower than are misleadingly counted in, while in practice, only the energy range comparable to should be considered. Nevertheless, the large-time behavior is reliable since both - and -functions are compressed more and more into the near-resonance region with increasing time.

In the limit in Eq. 22-23, the bound state decays to zero, but the continuum recovers the general form of Fano profile,

(25) |

We emphasize that Eq. 25 is a mathematical consequence. It cannot be obtained from the measurement since continuum configurations are not eigenstates of the total atomic Hamiltonian. In the special case where and , the system starts with only bound state, and Lorentz profile will be recovered as , i.e.,

(26) |

Fig. 1 shows the spectra at different times for , , , and . Each sub-figure demonstrates how its spectrum morphs into Fano profile. A feature found in the profile is the oscillation in the wings during the evolution, which shrinks its energy period when time increases, and disappears as . This oscillation, due to the phase between the two exponential terms in Eq. 23, is not seen in the conventional time-integrated spectra. Although analytically simple, the IC simplification is conceptually rich and offer an easy visualization for the autoionization process.

#### ii.2.2 Gaussian initial continuum distributions

The condition IC is valid only in the vicinity of a resonance and not applicable to the whole energy range. When a short pump pulse is used to initiate the system, one should cover its whole bandwidth and the state vectors must be properly normalized. Here we consider the normalization where the total probability for the bound component plus the continuum component is set to 1.0, where the ground state is not considered since it is not involved in the autoionization. This normalization is automatically fulfilled by the exact solution once the initial coefficients are normalized. Additionally, we assume that the pump pulse has a Gaussian envelope, resulting in a Gaussian profile in energy. The initial coefficients are given by

(27) | ||||

(28) |

where

(29) |

This initial condition is dubbed IC. The initial continuum state has a Gaussian energy distribution centered at with standard deviation (). The initial coefficients are determined by the pump width and the resonance parameters , (or equivalently, ), and . For , IC is reduced back to IC near the th resonance. The comparison between the profiles under the IC and IC is plotted in Fig. 2 for , where the bandwidth (FWHM) for IC is 15. As seen in Fig. 2, with the modification from IC to IC, the profile changes from a flat background to a Gaussian shape, and the and profiles change similarly. However, the featuring structures are preserved. The modification should be slight if the bandwidth is much wider than , which is true for a sufficiently short pump.

#### ii.2.3 Growth and decay of a bound state

Considering an isolated resonance, in IC, the total bound state is an exponential decay function of time as described by Eq. 22. However, under the normalization condition of IC, which is more realistic, the bound state decay is not necessarily monotonic. For the IC solution, the bound state coefficient in Eq. 12 and its two terms are shown in Fig. 3 in their absolute squares. The two panels are for two different initial bandwidths of the continuum centered at the resonance. The first term, which comes from the initial bound state, is always an exponential decay. The second term, which comes from the initial continuum, goes through a rising period up to about in Fig. 3(a) and in Fig. 3(b) before it decays. This rising period will shorten if the initial bandwidth increases or the pump duration shortens. In the limit of infinite bandwidth, both the first and the second terms will be monotonic decay functions, and the solution in IC will be recovered. Note that even in a situation where the length of the rising period is comparable to the lifetime, such as what Fig. 3(a) shows, the total bound state population will eventually be dominated by the exponential decay. This growth part of the bound state is usually much shorter than the decay lifetime, so to ”watch” this detail structure, a higher temporal resolution is required.

In the special case for , the pump pulse does not directly populate the bound state of the system. The -parameter defined by Eq. 18 is 0. However, this definition is based on the assumption that and the continuum background are constant near the resonance, and it excludes the modification term for the bound state. In an actual situation where no original bound state was populated, the theory needs to be modified to include higher order terms. Nonetheless, even if we neglect the initial population of the original and the modified bound state, the time-dependent probability as given by Eq. 12 will still have contribution from the initial continuum population . In other words, we will still have the ”decay” of the bound state even if it was not populated by the pump pulse. This somewhat awkward result is best understood in the following way. First recall that the initial bound state population depends on the transition operator connecting it to the ground state, while the ”configuration interaction” which connects the bound state to the continuum is always present. This ”configuration interaction” feeds the bound state from the continuum, but the continuum states also draw the bound state component and take it to large distances. If these continuum states were substituted by highly excited bound states, the outer electron would return to recollide with the inner electron(s) and repopulate the bound state. This would correspond to situations studied in Ref. gallagher1 ; gallagher2 ; jones . To initiate the oscillation, the bound state does not have to be excited by the pump pulse. The ever present ”configuration interaction” is the one which is responsible for the oscillation. Such ”interaction” is the property of the Hilbert space but not of the specific excitation mechanism.

## Iii Application to resonances in Be

In atomic systems, autoionization is understood only in the multiple-electron picture where electron-electron correlation is considered. Our theory is applied here to calculate the autoionizing wave packet of the resonance and of the ( to 9) resonance series embedded in the continuum in beryllium atom, generated by the photoionization from the ground state. For a short pump, we assume that the initial continuum distribution can be approximated by a Gaussian distribution. Only the two outer electrons in beryllium are assumed active, and the core electrons are frozen. The bound orbitals are Slater type orbitals where their analytical forms are determined by fitting the numerical Hartree-Fock calculations. The continuum basis functions are calculated with the model potential given in the literature buendia for standing waves joachain . The resonance parameters are taken from the experiment of Wehlitz et al wehlitz , where is 2.789 eV above the ionization threshold, , and eV, or fs. The pump beam has the carrier frequency right at the resonance and the duration of 2 fs, or mean bandwidth of 0.912 eV.

For the resonance case, the excited wave function, Eq. 11, is given by

(30) |

where the -functions are constructed by corresponding bound orbitals and continuum waves, and symmetrized between the two electrons. The one-electron density is defined by

(31) |

where the integral sums over the angular dependence as well as the radial part of one of the electrons. In the large- region, represents the electron density of the autoionizing electron. In this section, energy, time, and distance are in electron Volts (eV), femtoseconds (fs), and Bohr radius (a.u.) respectively unless otherwise specified. The bandwidth in energy and the duration in time are defined in FWHM.

First we show how the wave packet moves in the coordinate space. Fig. 4(a) displays the electron density at , 2 fs, and 4 fs, where the time scale is about the lifetime of the resonance fs. Before the propagation starts (), a dominant spatial distribution covers 0-60 a.u. This distribution is mainly contributed by the continuum, where the percentage of initial bound state is as low as 7%. The spatial oscillation with period of about 10-15 a.u., or wave number of about 0.4-0.6 a.u., matches the energy distribution centered at the resonance energy of 2.789 eV. As time goes on, the electron wave packet moves outward toward large . When time is long compared with the autoionization lifetime, the motion of the electron wave packet can be analyzed more easily. In Fig. 4(b), the long-time behavior of the wave packet are plotted for fs, 80 fs, 120 fs, and 160 fs. The electron is always confined in some spatial range, but it gradually spreads wider when time increases. By recognizing and tracing the dominant two peaks in Fig. 4(b) over time, we find that the corresponding velocities are almost constant. It is similar to a free electron wave packet traveling in space, where its average velocity keeps the same, but its spatial distribution spreads. In this long time regime, the electron moves in a very predictable way to the large distance toward the detector.

Next, we discuss the evolution of the resonance profile , i.e., the time-dependent Fano profile, for the same wave packet in Eq. 30. In Fig. 5(a), the photoelectron energy profile is shown at , 1.5 fs, 3 fs, and 4.5 fs, i.e., near the lifetime of the resonance fs. The profile starts evolving with a Gaussian shape at determined by the pump pulse, and changes quickly within its lifetime. At , just a little longer than the lifetime, the double-peaked shape characteristic of autoionization can be seen, but the shape is not quite ”settled” into the Fano profile. The long-time behavior is shown in Fig. 5(b). At fs to 20 fs, the resonance profile is already quite close to the final Fano profile that can be determined in the standard energy domain measurements. It is also physically sensible that the electron velocities observed in Fig. 4(b) can be mapped onto, although not exactly, the energies of the two highest peaks in Fig. 5(b).

We have discussed an isolated resonance. Practically, a short pump pulse is likely to populate more resonances at once. Here we consider the case including the resonances up to in Be. Even though the state (about 5.277 eV above the ground state nist ) is lower than the threshold and will not show up in the continuum profile, it can in general interact with the continuum beyond the threshold and is included in the present calculation. The parameters for these resonances are generated by setting , eV, and , where is the quantum defect and is the effective quantum number. The first few widths, starting from , are 0.514, 0.181, 0.0833, and 0.0451, in eV’s, and the decay lifetimes are 1.28, 3.64, 7.90, and 14.6 in femtoseconds, respectively. The resonance series stretches between the threshold and the threshold, a range of about 4 eV. Our attention is directed at the time scale and energy scale of this series. The initial continuum distribution is set to be centered at 2.3 eV above the ionization threshold, with bandwidth of 1.825 eV, or pump duration of 1 fs. The bandwidth, covering the resonances of interest, is limited to the energy range above the threshold to avoid unnecessary complications. Each resonance term in Eq. 21 is multiplied by a Gaussian window whose width is 2 to narrow down the interaction to some finite range in the calculation. The resultant resonance profiles are displayed in Fig. 6, along with the experimental photoabsorption cross section values in Ref. wehlitz , plotted by the light gray curve and normalized to fit in the vertical range of our calculation. The resonance series is seen in the experiment from low to high energy, starting with at about 1.5 eV. For the calculation, the figure shows how the profiles evolve in time and are finalized, where the lower (wider) resonances build up quickly and the higher (narrower) resonances slowly. At short times, the continuum must be viewed as a whole distribution which gives no information on individual resonances; only at large times, e.g. fs, the profiles assigned to individual resonances can be seen. Nonetheless, our calculation and the experiment in Ref. wehlitz should not be compared with each other directly. In the present work, a short pump pulse covers a wide energy range; the system propagates coherently with time. On the other hand, the experiment in Ref. wehlitz was done with many independent incoherent measurements, where each measurement was made for a single energy point using a long pulse of at least of hundreds of picoseconds, and gives back a single quantity for the cross section.

We caution the readers that the time-dependent resonance profile shown in Fig. 6 is from theoretical calculations only, i.e., a plot of vs for different time delays after the pump pulse. As emphasized earlier, cannot be measured directly in the laboratory since any physical measurement will take at least many picoseconds which cannot distinguish time evolution that happens on the femtosecond scale. Secondly, the wave functions associated with the coefficients in Eq. 19 are not eigenstates of the atom such that cannot be directly interpreted as measurable probabilities.

## Iv Method of probing the time-dependent autoionization processes

After we have predicted how the wave packet in the autoionization process evolves with time, an experimental scheme where such predictions can be confirmed is essential. Using SAP-pump and IR-probe, the lifetime of Auger decay of an inner-shell hole drescher or of a Fano resonance in helium gilbertson have been ”deduced”. The streaking of an autoionizing electron in the IR laser field has been calculated by Wickenhauser et al wickenhauser . The streaked electron spectra are too complicated to allow the retrieval of the predicted intermediate profiles. Taking beryllium as example, here we propose in the following a way to measure the time dependence of Fano profiles.

The probe scheme is sketched in Fig. 7. After the pump, we define the starting time for the autoionization as such that the pump field is finished for . A probe pulse is shined onto the system with its peak at time delay . The probe has the duration of 1.5 fs, or bandwidth of 1.216 eV, and carrier energy of 40 eV. This probe pulse ionizes the electron in the continuum (binding energy of is 18.21 eV). After the probe pulse is over, the excited wave packet of Eq. 30 changes to

(32) |

We have assumed that the bound states are not changed by the probe pulse, and only a small part of the continuum states are ionized (thus the amplitude now given by is slightly different from ) to generate new states. After the probe pulse, both and continue to change with time as the bound states decay. The newly created states are eigenstates of the atomic (field-free) Hamiltonian and are stationary states, and they can thus be probed by a laboratory detector. Using 40 eV probe pulses, the energy will be centered around 21.8 eV, well-separated from the electron energy of 2.5 eV due to autoionization. Clearly the probability should be independent of time and is proportional to . To determine one can measure the energies of the two electrons in coincidence or by measuring the low-energy electrons in coincidence with Be ions. Either way lower energy group duplicates the original continuum profile at the moment of the probe. By changing the delay time of the probe a copy of can be determined. Due to the finite duration of the probe pulse, the measured copy of will be somewhat smoothed out.

If the probe is of perturbative strength, we can calculate the predicted two-electron spectra quantitatively. In the presence of the external probe field , the coefficient for is given by the first-order time-dependent perturbation theory by

(33) | ||||

(34) | ||||

(35) |

where is the dipole interaction potential for the probe pulse. The interaction is turned on only for . Note that is now measured with respect to the center of the probe, and the total excited wave function at has already gone through the autoionization process for a duration of . In Eq. 33-34, the time-propagator in operates on and gives in return. In the next step in Eq. 35, since the bound states do not participate in the probe process, only the continuum part is written down. The dipole matrix element is assumed to be a constant. The polarization of the probe field is in the -direction, i.e. . The retrieved profile as a function of is defined as . Fig. 8 shows very good agreement between the retrieved profile at different ’s and the original continuum profile at the associated time .

While the calculation given above is in first-order perturbation, the scheme can be applied to intense probe as well as long as the probe field dominantly removes the electron and has little effect on other components in the system. The electron is not affected by the probe, and the continuum profile of at will still be transferred to plainly (with additional peaks possibly from two-photon ionization). The advantage of a stronger probe is that it increases the signal rate. This may be an important factor to realize this pump-probe scheme since the population of double excited states by the pump can be already weak. In the present example, the photon energy is 11.62 eV for the pump and 40 eV for the probe, presuming that the short UV sources are provided for both the pump and the probe.

## V Conclusions

In this paper we study the temporal behavior of the atomic autoionization process. Following the original paper of Fano, we treat autoionization as the decay of a discrete state into a degenerate continuum. Fano’s theory characterizes such a resonance with a shape parameter and a width , in addition to the resonance energy position. For many decades, Fano resonances are measured and characterized in energy domain experiments. With the advent of XUV and soft X-rays with pulse durations of attoseconds to a few femtoseconds, we address the question whether Fano resonances can be probed directly in the time domain within its typical lifetime. In particular, can we talk about the time-dependent Fano profiles before the discrete state is fully decayed and can such profiles be measured experimentally?

To describe time-dependent Fano profiles, it is inconvenient to express the total time-dependent wave function in terms of the eigenstates of the system. Instead, the configuration basis states, including the bound and the continuum, used in Fano’s original theory are much more useful (or more physical) for analyzing the autoionization process, and we are led to define the time-dependent Fano profiles by referring to the wave packets expanded in the continuum configuration basis states. To probe the temporal Fano profiles with good time resolution, we suggest that the decay process be perturbed by a short soft-X-ray pulse which interrupts the autoionization at the time of the probe. The probe in the meantime creates a stationary wave packet that can be measured in the laboratory. We further consider many Fano resonances together if they fall within the bandwidth of the pump pulse.

The pump-probe scheme proposed here, requiring high-energy short pulses for both the pump and the probe, is not yet available today. On the other hand, attosecond technology is undergoing rapid growth, the proposed scheme of measuring autoionization dynamics may become doable within a few years. At the same time, theory for attosecond electron dynamics is still a barren field today. General approaches for understanding electron dynamics for many-electron systems are badly needed at this time. Not all the probe pulses have similar effectiveness in retrieving the time-dependent wave packets. Without the general theoretical framework, it would be difficult to retrieve electron dynamics from the signals measured by the probe pulse.

###### Acknowledgements.

This work was supported in part by Chemical Sciences, Geosciences and Biosciences Division, Office of Basic Energy Sciences, Office of Science, U.S. Department of Energy.## References

- (1) H. Beutler, Z. Physik 93, 177 (1935).
- (2) U. Fano, Phys. Rev. 124, 1866 (1961).
- (3) J. Faist et al, Nature 390, 589 (1997).
- (4) V. Madhaven et al, Science 280, 567 (1998).
- (5) A. E. Miroshnichenko, S. Flach, and Y. S. Kivshar, Rev. Mod. Phys. 82, 2257 (2010).
- (6) F. Krausz and M. Yu. Ivanov, Rev. Mod. Phys. 81, 163 (2009).
- (7) M. Drescher et al, Nature 419, 803 (2002).
- (8) M. Kitzler, N. Milosevic, A. Scrinzi, F. Krausz, and T. Brabec, Phys. Rev. Lett. 88, 173904 (2002).
- (9) Z. X. Zhao and C. D. Lin, Phys. Rev. A 71, 060702 (2005).
- (10) M. Wickenhauser, J. Burgdörfer, F. Krausz, and M. Drescher, Phys. Rev. Lett. 94, 023002 (2005).
- (11) O. Smirnova, V. S. Yakovlev, and A. Scrinzi, Phys. Rev. Lett. 91, 253001 (2003).
- (12) S. Gilbertson et al, International Quantum Electronics Conference (IQEC), Baltimore, Maryland May 31, 2009 (Optical Society of America, Washington, D.C., 2009).
- (13) H. Wang et al, Phys. Rev. Lett. 105, 143002 (2010).
- (14) Th. Mercouris, Y. Komninos, and C. A. Nicolaides, Phys. Rev. A 75, 013407 (2007).
- (15) L. Argenti and E. Lindroth, Phys. Rev. Lett. 105, 053002 (2010).
- (16) S. Cavalieri and R. Eramo, Phys. Rev. A 58, R4263 (1998).?
- (17) S. N. Pisharody and R. R. Jones, Science 303, 813 (2004).
- (18) D. W. Schumacher, B. J. Lyons, and T. F. Gallagher, Phys. Rev. Lett. 78, 4359 (1997).
- (19) B. J. Lyons, D. W. Schumacher, D. I. Duncan, R. R. Jones, and T. F. Gallagher, Phys. Rev. A 57, 3712 (1998).
- (20) E. Buendía, F. J. Gálvez, P. Maldonado, and A. Sarsa, J. Phys. B: At. Mol. Opt. Phys. 39, 3575 (2006).
- (21) C. J. Joachain, Quantum collision theory (Elsevier Science Publishing, 1984).
- (22) R. Wehlitz, D. Lukić, and J. B. Bluett, Phys. Rev. A 68, 052708 (2003).
- (23) NIST Atomic Spectra Database, http://www.nist.gov/physlab/data/asd.cfm.