The theory of measurement projection operators in grid-based time-dependent wave-packet calculations involving electronic continua in atoms and molecules is discussed. A hierarchy of projection operators relevant in their individual restricted configuration spaces is presented. At asymptotically large distances from the scattering or interaction center the projection operators involve plane waves only. To reach this asymptotic regime, however, large propagation times and large boxes may be required. At somewhat smaller distances from the scattering center, the projection operators are expressed in terms of analytical single-center Coulomb scattering waves with incoming wave boundary conditions. If propagation of the wave packet to these asymptotic regimes is impeded, the projection operators involve the exact scattering states which are not readily available in the wave-packet calculation and hence must be supplied by an additional, typically very demanding, calculation. The present approach suggests an exact way of analyzing the timely problem of the one-electron continuum in nonperturbative calculations. A key feature is that the propagated wave packet includes every interaction of the full Hamiltonian. The practicality of the proposed method is illustrated by the nontrivial example of strong-field ionization of the molecular hydrogen ion. Finally, the extension of the presented ideas to single and double ionization of two-electron systems is discussed.

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.

We present experimental results and theoretical analysis of two-photon ionization of He by a superposition of the 7th to the 13th harmonic of a Ti:sapphire laser. Solving the time-dependent Schrödinger equation for He in a coherent polychromatic field, the ${\mathrm{H}\mathrm{e}}^{+}$ yield is calculated. From this yield the number of ${\mathrm{H}\mathrm{e}}^{+}$ ions produced has been estimated and found in reasonable agreement with its measured value. The present results establish the feasibility of a second-order autocorrelation measurement of superposition of harmonics, and thus they represent the precursor towards the direct temporal characterization of attosecond pulse trains.