🔬 atomic physics

Qprop-ND: A Nondipole Extension of the Qprop Package for Strong-Field Laser Physics

The paper introduces Qprop-ND, an extension of the Qprop 3.0 package that incorporates nondipole effects related to laser wave propagation to enable accurate simulations of strong-field atomic phenomena where magnetic field effects or laser wavelength size become significant.

Original authors: Haoyuan He, Zhao-Han Zhang, Karen Z. Hatsagortsyan, Christoph H. KeiteL

Published 2026-09-29

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Haoyuan He, Zhao-Han Zhang, Karen Z. Hatsagortsyan, Christoph H. KeiteL

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the microscopic world of atoms, electrons are usually thought of as tiny planets orbiting a nucleus, held in place by an invisible electric tether. When scientists shine a powerful laser on an atom, they can knock these electrons loose, turning the atom into an ion and sending the electron flying off. For decades, physicists have studied this process using a simplified rule of thumb: they assumed the laser light was so uniform across the tiny atom that the electron only felt the electric push and pull, ignoring the magnetic side of the light wave. This assumption worked perfectly for most experiments, allowing researchers to map out how atoms behave under intense light. However, as laser technology has advanced, creating beams that are both incredibly intense and extremely short in wavelength, this simplification has begun to break down. When the light wave is so short that it is comparable in size to the atom itself, or when the electron is moving fast enough to feel the magnetic tug of the light, the old rules no longer tell the whole story. The electron's path shifts in unexpected ways, and the patterns it leaves behind change, hinting at a more complex interaction between matter and light than previously understood.

To explore these subtle shifts, a team of researchers at the Max Planck Institute for Nuclear Physics in Germany has developed a new version of a powerful computer program called Qprop. This software acts as a virtual laboratory, allowing scientists to simulate how an electron behaves when hit by a laser pulse without having to build a physical experiment that might be too expensive or difficult to control. The original version of this program was excellent at handling the standard electric effects of light, but it ignored the magnetic influences that become important in extreme conditions. The new update, named Qprop-ND, adds these missing magnetic pieces to the simulation. By doing so, it allows researchers to see what happens when the laser wave's magnetic field pushes the electron, a force that was previously left out of the calculations. The team tested this new tool by simulating the ionization of hydrogen and helium atoms under various conditions, including the use of ultra-short, high-energy X-ray pulses that are becoming available at modern research facilities.

The results of these simulations reveal a distinct change in how electrons are ejected when the magnetic effects are included. In the standard view, where only the electric field is considered, the electron tends to fly off in two main directions, creating a pattern with two lobes. However, when the new magnetic corrections are turned on, a third lobe appears in the pattern, pointing directly opposite to the direction the laser beam is traveling. This three-part structure is a clear signature of the nondipole effect, proving that the magnetic field of the light is actively reshaping the electron's journey. The researchers found that this effect becomes more pronounced as the laser gets stronger and the wavelength gets shorter. In some cases, the magnetic push is strong enough to suppress the electron's ability to crash back into the atom, a process that usually creates bright flashes of light known as high-order harmonic generation. This suggests that in the most extreme laser environments, the magnetic field acts as a brake, preventing certain interactions that were once thought to be inevitable.

The team also investigated how the shape of the laser light affects these outcomes, comparing straight-line polarization with circular polarization. They discovered that the electron's behavior depends heavily on how the electric field rotates. When the laser spins in a circle, the electron is pushed away from the atom's core more quickly, feeling less of the core's gravitational-like pull. This allows the magnetic drift to dominate the electron's motion at lower laser intensities compared to when the laser vibrates in a straight line. Consequently, the patterns formed by the electrons look different: circular light produces a more compact, symmetrical burst, while straight light creates a wider, more stretched-out shape. These differences are not just minor details; they represent a fundamental shift in how energy is shared between the electron and the ion. The simulations show that the magnetic field can shift the electron's momentum by a specific amount, leaving a clear mark on the final distribution of the ejected particles.

To ensure their new program was accurate, the researchers compared their results against other established calculations and known physical limits. They found that their simulations matched well with previous data for hydrogen-like atoms, even in the most intense fields where the laser strength approaches the limits of what current theory can handle. In cases where the laser was extremely strong, the new program predicted slightly higher chances of ionization than older models, a difference that the authors attribute to the fact that their simulation includes the first layer of magnetic corrections while other models might rely on fully relativistic equations that are much harder to solve. This agreement gives confidence that the new tool is reliable for exploring the "nondipole regime," a territory where the size of the light wave and the speed of the electron make the magnetic force impossible to ignore.

The work presented in this paper does more than just refine a computer code; it provides a necessary bridge for future experiments. As facilities like the European XFEL and the Linac Coherent Light Source begin to deliver even more powerful X-ray beams, scientists will need tools that can accurately predict what happens when these beams hit matter. The Qprop-ND package offers a way to do this, filling a gap in our understanding of how light and matter interact at the smallest scales. By capturing the subtle influence of the magnetic field, the researchers have shown that the electron's path is not just a simple response to an electric push, but a complex dance influenced by the full nature of the light wave. This deeper understanding is essential for interpreting the data that will come from the next generation of laser experiments, helping to unlock new phenomena in the realm of strong-field physics.

Technical Summary of Qprop-ND: A Nondipole Extension of the Qprop Package for Strong-Field Laser Physics

Problem Statement
The interaction of atoms with superintense laser fields and/or hard X-rays often enters a regime where the standard dipole approximation fails. This nondipole regime arises when the laser wavelength becomes comparable to the atomic size or when the electron velocity approaches the speed of light, rendering the spatial dependence of the laser field and the associated magnetic field non-negligible. While the fully relativistic Dirac equation provides a rigorous framework, it presents significant numerical challenges. Conversely, the nonrelativistic time-dependent Schrödinger equation (TDSE) is insufficient without corrections for laser wave propagation effects. Existing tools like the Qprop package (version 3.0) are highly efficient for solving the 3D TDSE under the dipole approximation but lack the capability to simulate strong-field phenomena where magnetic-field-induced Lorentz forces and retardation effects are critical. There is a need for an efficient computational framework capable of accurately predicting photoelectron momentum distributions (PMDs) in these nondipole regimes without resorting to the full complexity of the Dirac equation.

Methodology
The authors present Qprop-ND, an extension of the Qprop 3.0 package that incorporates first-order nondipole corrections ( $\sim 1/c$ ) into the Hamiltonian. The core methodology involves:

Hamiltonian Formulation: The vector potential $\mathbf{A}(z,t)$ is expanded to first order in $z/c$ . This introduces new terms into the minimal-coupling Hamiltonian, specifically coupling the electric field $\mathbf{E}(t)$ with the spatial coordinate $z$ and the vector potential $\mathbf{A}(t)$ . These terms account for the magnetic field interaction and retardation effects.
Wave Function Expansion: Similar to the dipole version, the total wave function is expanded in spherical harmonics $Y_l^m(\Omega)$ . This reduces the 3D TDSE to a set of coupled one-dimensional radial equations for the partial-wave components $\Phi_{lm}(r,t)$ .
Coupling Mechanisms: The nondipole terms introduce couplings between different angular momentum channels $(l, m)$ . Specifically, the Hamiltonian includes terms that couple $(l, m)$ to $(l, m\pm1)$ , $(l\pm2, m\pm1)$ , and $(l\pm1, m)$ . The authors derive explicit matrix representations for these couplings, decomposing the Hamiltonian into dipole, mixed nondipole, and angular nondipole components.
Numerical Propagation: The time evolution is performed using the Crank-Nicolson scheme to ensure numerical stability and norm conservation. To handle the non-diagonal structure introduced by the nondipole terms (which couple different $l$ and $m$ subspaces), the authors employ a splitting method for the time-evolution operator. This involves factorizing the evolution into a sequence of exponentials of sub-blocks, which are then approximated using the Crank-Nicolson method. Special mathematical treatments (involving unitary transformations and block-tridiagonal matrix decompositions) are applied to efficiently invert the resulting matrices.
Spectrum Calculation: Photoelectron spectra are calculated using the $t$ -SURFF (time-dependent surface flux) and $i$ -SURFV (improved surface flux with velocity) methods, which extract momentum information from the probability flux through a surface, avoiding the need to propagate the wave function over vast distances.

Key Contributions

Algorithmic Development: The paper provides the mathematical derivation and numerical implementation of first-order nondipole corrections within the spherical harmonics expansion framework of Qprop.
Code Release: The authors release the Qprop-ND software package (written in C++), which extends the capabilities of the widely used Qprop tool to the nondipole regime.
Validation: The implementation is rigorously benchmarked against established relativistic calculations and nonperturbative TDSE results from the literature.

Results
The performance of Qprop-ND was demonstrated through several illustrative examples:

Hydrogen-like Helium in Ultrastrong XUV Fields: The code was tested against relativistic calculations by Telnov and Chu for the ionization of He $^+$ $^{+}$ in linearly polarized fields.
- Low-Energy Structure: Qprop-ND successfully reproduced the characteristic three-lobe low-energy structure in the PMD, where a central lobe appears in the laser counterpropagation direction, distinct from the two-lobe dipole pattern.
- Quantitative Agreement: For field strengths corresponding to $a_0 \approx 0.1$ and $0.17$, the ionization probabilities matched the benchmark relativistic results within a few percent. At higher fields ( $a_0 \approx 0.25$ ), a systematic overestimation of ionization probability was observed compared to the fully relativistic Dirac equation results. The authors attribute this to the absence of higher-order relativistic effects (kinematics, spin dynamics) in the nonrelativistic TDSE framework, confirming that Qprop-ND captures the dominant nondipole physics but not full relativistic dynamics at these intensities.
X-ray Regime (Atomic Hydrogen): Comparisons were made with nonperturbative TDSE calculations by Dondera and Bachau for photon energies up to 3 keV.
- Angular Distributions: The code reproduced the slight asymmetry in angular distributions predicted by nondipole theory, while confirming that the energy spectra remain largely unaffected by leading-order nondipole corrections in this regime.
- Agreement: The results showed excellent agreement with the benchmark, validating the code's accuracy in the high-frequency, weakly relativistic domain.
Polarization Dependence: The code was used to investigate zero-energy structures (ZES) under linear (LP) and circular (CP) polarization.
- Distinct Dynamics: The study revealed that under equal-intensity conditions, CP fields produce a dominant central lobe in the ZES, whereas LP fields retain stronger side lobes.
- Regime Transition: The analysis identified two regimes: a "Coulomb-influenced" regime and a "laser-dominated" regime. The transition to the laser-dominated regime (where nondipole drift overcomes Coulomb attraction) occurs at a lower normalized vector potential $a_0$ for CP fields compared to LP fields.

Significance
The paper claims that Qprop-ND provides a robust and efficient computational framework for investigating strong-field ionization dynamics in the nondipole regime. By extending the well-established Qprop package, it allows researchers to simulate phenomena such as above-threshold ionization, high-order harmonic generation, and ionization stabilization where magnetic field effects are non-negligible. The code serves as a reliable tool for interpreting experimental data from modern facilities (e.g., XFELs, ELI) where laser intensities and frequencies push interactions beyond the dipole approximation. The authors emphasize that while Qprop-ND is not a fully relativistic solver, it accurately captures the leading-order nondipole physics essential for understanding the interplay between laser propagation effects and atomic dynamics in the weakly relativistic regime.

More like this