Multielectron Tunneling Ionization

• Multielectron tunneling ionization is defined as the removal of two or more electrons via strong-field quantum tunneling, exhibiting both sequential and correlated pathways.
• It involves shake-up effects, under-the-barrier correlations, and multiorbital contributions that extend beyond the single-active-electron model.
• Advanced computational and experimental methods, such as TD-MCSCF, TDDFT, and COLTRIMS, validate predictions and elucidate complex electron dynamics.

Multielectron tunneling ionization refers to the removal of two or more electrons from an atom or molecule by strong-field-induced quantum tunneling, with multielectron correlations, multi-orbital character, or collective dynamics playing a fundamental—sometimes essential—role. Unlike the single-active-electron (SAE) scenario, where only the highest-occupied orbital participates and electrons are ejected independently, multielectron tunneling encompasses sequential and nonsequential ionization pathways, shake-up and correlation effects, and complex angular and momentum distributions in the ionization products. Advances in theory, computation, and experiment have led to a comprehensive understanding of this phenomenon, spanning single atoms, small molecules, and complex polyatomic systems.

1. Theoretical Foundations: Sequential and Correlated Models

The original theoretical description of strong-field ionization employs the SAE tunneling picture, formalized by the Ammosov–Delone–Krainov (ADK) rate: $$W_{\rm ADK}(F)=C_{n^* l^*}^2\left(\frac{2F}{\pi\kappa^3}\right)^{1/2}\left(\frac{2\kappa^3}{F}\right)^{2n^*-|m|-1}\exp\left(-\frac{2\kappa^3}{3F}\right)$$ where $F$ is the instantaneous field, $\kappa=\sqrt{2I_p}$, and $C_{n^* l^*}$ is an angular-momentum constant. In the sequential SAE model, electrons tunnel one by one, their ionization rates exponentially dominated near the pulse maximum.

For quantitative corrections in very strong fields, the “augmented quantum tunneling” (AQT) rate incorporates empirical saturation corrections: $W_{\rm AQT} = W_{\rm ADK} \exp(-\alpha\kappa^3/F)$.

The fully correlated classical ensemble approach instead propagates all electrons under Newtonian dynamics with full electron-electron and electron-nucleus couplings plus the laser field. Here, recollision, nonsequential emission, and explicit correlation effects are naturally included. Emission times and asymptotic electron momenta are extracted from the classical trajectories. These two approaches bracket the primary theoretical frameworks for interpreting experiments on multielectron tunneling in strong fields (Wang et al., 2012).

2. Multielectron Dynamics and Correlation Mechanisms

Beyond sequential tunneling, multielectron effects arise from shake-up, under-the-barrier correlation, multiorbital contributions, and collective or correlated electron escape.

Shake-up scenarios, as realized in the METI model, predict that as one electron tunnels, the rapid change in the ionic potential projects the residual electron(s) onto ionic excited states with probability $$P_\text{exc}^j = |\langle\phi_\text{core}^j|\phi_\text{core}^0\rangle|^2$$. This mechanism, directly observed in experiments with circularly-polarized pulses on argon, enhances ionization into excited ionic states beyond single-orbital ADK expectations (Bryan et al., 2010).

Correlation-assisted escape and under-the-barrier interactions can be captured analytically using the R-matrix formalism, where higher-order tunneling amplitudes account for ionization into “deeper” channels (those with higher and wider tunneling barriers) via in-barrier electron–electron (or electron–ion) correlation. Exact evaluation of spatial and momentum integrals beyond the standard saddle-point approximation captures the transfer of angular momentum and the resulting nontrivial structures in the photoelectron angular distributions (Pisanty, 2013).

Nonsequential and collective tunneling regimes emerge under ultrashort pulses: if pulse durations are shorter than the saturation time for single-electron emission, correlated two-electron tunneling ('collective tunneling') is possible. Electron–electron repulsion suppresses, but does not prohibit, such events; the breakdown of the sequential model and the appearance of a diagonal ridge in two-electron coordinate or momentum densities are direct signatures (Tyurin et al., 29 Apr 2025).

3. Advanced Many-Electron Theory and Computational Methods

Contemporary approaches systematically include multielectron and correlation effects by working at the level of the full multielectron wavefunction or reduced density matrices.

Time-dependent multiconfiguration self-consistent-field (TD-MCSCF) and optimized coupled-cluster (TD-OCC) methods propagate all-electron wavefunctions expanded in time-dependent determinant (or cluster) spaces. The time-dependent variational principle (TDVP) ensures gauge invariance and Ehrenfest’s theorem. Electron correlation, orbital relaxation, and multichannel ionization are captured on equal footing. Gauge-invariant observables such as the $k$-electron ionization probabilities and time-resolved photoelectron spectra are obtained from real-space grids with absorbing boundaries (Sato et al., 2022).

Time-dependent density functional theory (TDDFT), solved in real time with real-space grids, resolves both total and orbital-resolved ionization yields, enabling identification of ionization from inner Kohn–Sham orbitals as a direct probe of multielectron effects. Photoelectron angular distributions (PADs) computed by integrating the radial flux through a spherical surface elucidate the spatial signatures of multiorbital participation (Yahel et al., 2017).

Integral-representation many-electron weak-field asymptotic theory (IR-ME-WFAT) permits robust evaluation of ionization rates in molecules by explicitly including key many-electron matrix elements in the tunneling amplitude through Gaussian-basis-friendly integrals. This approach, applicable to arbitrary geometries and quantum-chemistry methods, matches experiment and unambiguously resolves important orientation controversies (e.g., CO field orientation dependence) (Wahyutama et al., 2024).

4. Experimental Probes and Observables

Modern experiments deploy coincidence detection, COLTRIMS “reaction microscopes,” and intensity-selective scanning to reconstruct the full angular and momentum correlation of the ejected electrons and residual ions.

Elliptical or circular polarization is used to suppress recollision, enabling unambiguous identification of tunneling times and emission direction (“attoclock” protocols) (Wang et al., 2012). The mapping of release times to momentum components leverages the rotating laser field and drift momentum relation $\mathbf{p}(\infty) = -\mathbf{A}(t_i)$.

Orientation-resolved molecular frame PADs and KER spectra allow assignment of events to specific orbitals and sequential or nonsequential channels. As shown in CO, covariance analysis of ion sum-momentum distributions and their orientation-resolved asymmetries provides clear evidence for multiorbital sequential tunneling and the minor role of the linear Stark shift (Wu et al., 2012).

Ionization yield deconvolution and normalization protocols, especially “conserved probability of ionization” (CPI), enforce probability conservation and enable absolute comparison of experimental yields with theory, circumventing focal-geometry artifacts (Bryan et al., 2010). Experimental deviations from sequential-ADK predictions validate the necessity of including multielectron shake-up and correlated effects.

5. Multielectron Effects in Molecular Systems

Molecules display several unique avenues for multielectron tunneling not present in atoms:

• Enhanced ionization of diatomics at critical internuclear separation $R_c$ due to field-induced inner-barrier lowering, with unambiguous demonstration that the up-field center is the active tunneling site (COLTRIMS with ellipticity-probe) (Wu et al., 2013).
• Multiorbital contributions: Multiple occupied molecular orbitals participate in strong-field ionization depending on orientation, as shown by orientation-dependent PADs and switching of “dominant” orbital yield at high intensity or certain angles in polyatomic and substituted aromatics (Yahel et al., 2017).
• Induced dipole and Stark effects: Multielectron polarization of the parent cation modifies both the effective tunneling barrier and the electron’s trajectory post-tunneling. Quantitative 3D TDSE and semiclassical analyses in CO show MEP-induced symmetry restoration and fringe wash-out in PMDs (Abu-samha et al., 2024).
• Ultrashort pulse regime: Collective tunneling becomes more prominent as the relative timing of pulse envelopes and ionization saturation times aligns, as observed in double-ionization suppression and correlated coordinate densities in Br$^-$ (Tyurin et al., 29 Apr 2025).

6. Nonadiabatic, Coherence, and Application-Driven Perspectives

Recent advances highlight the interplay of multielectron tunneling with ionic coherence, subcycle nonadiabatic dynamics, and applications to coherent light-matter control:

• Nonadiabatic corrections to the ionization rate, including saddle-point and subcycle effects, improve quantitative accuracy in capturing the timing and coherence of tunnel-ionization bursts at the attosecond scale (Yuen, 11 Jan 2026).
• Density-matrix SFI frameworks formalize the equivalence of reduced quantum state propagation and multielectron wavefunction approaches, with direct implications for modeling strong-field-induced ionic coherence and gain in air-lasing and mid-IR pumping scenarios (e.g., $\mathrm{CO}_2^+$ $C$–$B$ population inversion) (Yuen, 11 Jan 2026).
• Control paradigms: By engineering tunnel-ionization coherence (via field parameters, pulse duration, alignment), one can aim to steer attosecond electron emission and photochemical outcome dynamically.

7. Outlook and Open Questions

Multielectron tunneling ionization is now understood as a multi-faceted process involving a spectrum of mechanisms:

• Sequential SAE-like tunneling is often valid for long pulses in atoms and simple molecules.
• Deviations due to shake-up, correlation, or multiorbital involvement are pronounced whenever the field reaches over-barrier intensities, the pulse duration drops below saturation time, or recollision is suppressed.
• Current theories, including TD-MCSCF, TD-OCC, TDDFT, R-matrix, and IR-ME-WFAT, offer systematically improvable descriptions and are validated across atomic and molecular systems.
• Open issues remain regarding the detailed classical mechanism of “release” in regimes where ADK/AQT underestimates observed emission times or yields, the full quantification of under-the-barrier correlation pathways, and the universal extension to large, strongly correlated, or multireference systems.

Further investigation into experimental regimes combining attosecond temporal resolution, complex molecular targets, and state-resolved detection will continue to refine the physical picture, driving both fundamental and application-driven advances in ultrafast strong-field science (Wang et al., 2012, Yuen, 11 Jan 2026, Wahyutama et al., 2024).