A general approach to few-cycle intense laser interactions with complex atoms

A general {\it ab-initio} and non-perturbative method to solve the time-dependent Schr"odinger equation (TDSE) for the interaction of a strong attosecond laser pulse with a general atom, i.e., beyond the models of quasi-one-electron or quasi-two-electron targets, is described. The field-free Hamiltonian and the dipole matrices are generated using a flexible $B$-spline $R$-matrix method. This numerical implementation enables us to construct term-dependent, non-orthogonal sets of one-electron orbitals for the bound and continuum electrons. The solution of the TDSE is propagated in time using the Arnoldi-Lanczos method, which does not require the diagonalization of any large matrices. The method is illustrated by an application to the multi-photon excitation and ionization of Ne atoms. Good agreement with $R$-matrix Floquet calculations for the generalized cross sections for two-photon ionization is achieved.

💡 Research Summary

The paper presents a fully ab‑initio, non‑perturbative computational framework for solving the time‑dependent Schrödinger equation (TDSE) of a many‑electron atom interacting with an intense, few‑cycle attosecond laser pulse. The authors combine two powerful numerical techniques: a flexible B‑spline R‑matrix (BSR) method for constructing the field‑free atomic Hamiltonian and dipole matrix elements, and an Arnoldi‑Lanczos time‑propagation algorithm that avoids explicit diagonalization of large Hamiltonian matrices.

In the BSR part, the atomic space is divided into an inner region, where electron‑electron correlation is treated explicitly, and an outer region, where a single electron moves in the long‑range Coulomb potential. Within each region, term‑dependent one‑electron orbitals are generated from B‑splines and assembled into a non‑orthogonal basis set. This non‑orthogonal, term‑specific basis allows simultaneous representation of bound and continuum states for all electrons, thereby incorporating correlation effects without resorting to model potentials or frozen‑core approximations. The resulting Hamiltonian and dipole matrices retain the full multi‑electron structure and are ready for time‑dependent calculations.

For the time evolution, the authors employ the Arnoldi‑Lanczos algorithm. Starting from an initial wavefunction, the algorithm builds a Krylov subspace of dimension m by repeatedly applying the time‑dependent Hamiltonian. Within this subspace the exponential propagator is approximated by a small Hessenberg matrix, which can be exponentiated efficiently. Because the method never requires the full Hamiltonian to be diagonalized, memory consumption scales linearly with the size of the basis, and the computational cost is dominated by matrix‑vector products, which are highly amenable to parallelisation. The algorithm is stable for the non‑Hermitian Hamiltonians that arise when the laser field is included in the length gauge, and the time step Δt and Krylov dimension m can be tuned to achieve a desired balance between accuracy and speed.

To validate the approach, the authors apply it to neon (Ne) atoms, focusing on multiphoton excitation and ionization processes driven by an intense few‑cycle pulse. They compute generalized two‑photon ionization cross sections and compare them with results obtained from an independent R‑matrix Floquet calculation. The agreement is excellent across a range of photon energies, demonstrating that the BSR‑Arnoldi‑Lanczos scheme captures both the correct magnitude and the subtle energy dependence of the multiphoton ionization probability. The authors also present three‑photon ionization yields and discuss the role of intermediate resonances that are naturally included in the full‑dimensional wavefunction.

Beyond the specific neon example, the paper outlines several broader implications. First, the use of non‑orthogonal term‑dependent orbitals eliminates the need for separate orthogonalisation procedures and enables a compact representation of highly correlated states. Second, the Arnoldi‑Lanczos propagator makes it feasible to simulate attosecond dynamics in systems with dozens of electrons, a regime previously inaccessible due to prohibitive matrix sizes. Third, the method is readily extensible: spin‑orbit coupling, relativistic corrections, and even nuclear motion can be incorporated by augmenting the BSR Hamiltonian, while the time‑propagation algorithm remains unchanged.

The authors acknowledge current limitations, such as the omission of explicit spin‑orbit interaction and the need for careful selection of Krylov subspace parameters to avoid loss of orthogonality over long propagation times. Nevertheless, the presented framework constitutes a significant step toward a universal, high‑accuracy tool for strong‑field and attosecond physics, capable of addressing complex atoms, ions, and eventually molecules under realistic laser conditions.