🔬 optics

Meshfree versus grid-based Schrödinger solvers for modeling the interactions between free-electron wave packets and light

This paper introduces a mesh-free numerical framework that solves the time-dependent Schrödinger equation for free-electron wave packets interacting with light, offering an efficient and scalable alternative to grid-based methods that is up to 800 times faster while eliminating the need for spatial meshing and avoiding eikonal approximations.

Original authors: Mitja Funk, Sebastian Merk, Marlis Hochbruck, Caroline Lasser, Nahid Talebi

Published 2026-07-03

📖 3 min read☕ Coffee break read

CC BY 4.0

Original authors: Mitja Funk, Sebastian Merk, Marlis Hochbruck, Caroline Lasser, Nahid Talebi

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

Imagine you are trying to predict how a swarm of tiny, invisible bees (electrons) will dance when they fly through a complex, shifting wind tunnel (electromagnetic light fields). Scientists have long needed a way to simulate this dance to build better microscopes and quantum devices.

For a long time, the standard way to do this was like using a giant, rigid grid (like a chessboard) laid over the entire room. To see how the bees move, the computer has to check every single square on the board, even the empty ones where no bees are flying. This is accurate, but it's incredibly slow and computationally expensive, especially when the bees are moving fast or the wind is changing rapidly.

The New Approach: A Mesh-Free "Swarm" Method

This paper introduces a new, much faster way to solve the same problem. Instead of using a rigid grid, the authors use a "mesh-free" method.

Think of it this way: Instead of checking every square on a chessboard, the computer tracks a swarm of individual, glowing balloons (Gaussian wave packets). Each balloon represents a small part of the electron's "cloud."

How it works: The computer calculates the path of each balloon individually as it flies through the wind. Because the math of quantum mechanics is "linear" (like mixing paints), you can calculate the path of 100 balloons separately and then just add them all together to get the final picture of the whole swarm.
The Result: This method doesn't need to check empty space. It only focuses on where the "bees" actually are.

The "Thawed" Balloon

The authors use a specific type of balloon called a "Thawed Gaussian Wave Packet." Imagine a balloon that can stretch, shrink, and wiggle as it flies, adapting to the wind currents around it. This allows the simulation to capture the complex, wavy nature of the electron without getting bogged down in unnecessary calculations.

The Race: Old Grid vs. New Swarm

The researchers tested this new method against the old, trusted grid method using two different scenarios:

The "Near-Field" Dance: An electron flying past a tiny gold rod (like a nanorod) that is vibrating with light. This is like a bee flying close to a vibrating speaker.
The "Free-Space" Dance: An electron flying through two crossing laser beams in empty space. This is like a bee flying through a complex pattern of intersecting wind streams.

The Findings

Accuracy: The new "swarm" method produced results that were virtually identical to the old "grid" method. The patterns of energy and momentum the electrons gained were the same.
Speed: This is the big winner. The new method was up to 800 times faster than the grid method.
- Analogy: If the old grid method took 14 minutes to simulate a specific interaction, the new method did it in just 1 second.
Efficiency: The new method is so efficient that it can handle complex, changing environments without needing massive computing power.

Why This Matters

Currently, many scientists use shortcuts (approximations) to simulate these interactions because the full calculation is too slow. This new tool removes the need for those shortcuts. It allows scientists to simulate the full, complex dance of electrons and light with high precision and incredible speed.

In short, the authors have built a high-speed, flexible calculator that lets us watch the quantum dance of electrons and light in real-time, opening the door to designing better electron microscopes and controlling electron beams with unprecedented precision.

Technical Summary: Meshfree versus Grid-based Schrödinger Solvers for Modeling Free-Electron Wave Packet Interactions with Light

Problem Statement
The interaction between free-electron wave packets and electromagnetic fields is central to ultrafast electron microscopy and quantum nanophotonics, enabling the engineering of electron wave packets for applications such as energy combs and momentum-state superpositions. While existing theoretical descriptions often rely on eikonal or no-recoil approximations, these fail to capture the full dynamical evolution of electron wave packets in regimes involving electron recoil, strong-field effects, or dynamics beyond the rotating-wave approximation.

Numerical approaches based on the minimal-coupling Hamiltonian are required for a complete treatment. However, standard grid-based methods coupling the time-dependent Schrödinger equation with Maxwell's equations face significant computational bottlenecks. This arises from the substantial mismatch between the spatial and temporal discretization scales required for the Schrödinger equation versus Maxwell's equations, making the Schrödinger solver the limiting factor in hybrid simulations.

Methodology
The authors present a mesh-free numerical framework that directly solves the time-dependent single-particle Schrödinger equation for arbitrary electromagnetic potentials. The method is based on the Thawed Gaussian Wave Packet (TGWP) approximation.

Mathematical Formulation: The wave function $\psi(t, r)$ is approximated as a sum of $N$ Gaussian wave packets:
$\psi(t, r) \approx \sum_{j=1}^N c_j e^{iS_j(t)/\hbar} g_j(t, r)$
where $c_j$ are time-independent complex coefficients, and $g_j(t, r)$ are time-dependent Gaussians defined by their centers in position ( $q_j$ ) and momentum ( $p_j$ ), width matrices ( $C_j$ ), and normalization/phase factors ( $\gamma_j$ ).
Initialization: The initial Gaussian wave packet is represented as an over-complete set of Gaussians. The expansion coefficients are determined using quasi-Monte Carlo quadrature with Sobol' points weighted by a Gaussian function.
Evolution: Due to the linearity of the Schrödinger equation, each Gaussian is evolved semi-classically and independently. The parameters ( $q_j, p_j, Q_j, P_j, S_j$ ) follow classical equations of motion derived from the local harmonic approximation of the potential $V(t, r)$ around the packet's center. The width matrices are parametrized using the Hagedorn parametrization ( $C_j = P_j Q_j^{-1}$ ). The equations of motion are integrated using a second-order Störmer-Verlet integrator.
Benchmarking: The TGWP method is benchmarked against a grid-based reference solver using the split-step Fourier method with Strang splitting. Both methods are applied to two distinct physical scenarios:
1. Quasi-static near-field: Interaction with the optical near-field of a gold nanorod (PINEM configuration).
2. Stimulated Compton scattering: Interaction with two superposed, inclined laser pulses in free space.

Key Results
The study compares the TGWP approximation against the grid-based solver across various parameter sets, including different electric field strengths, electron kinetic energies, and laser wavelengths.

Accuracy:
- Near-field (PINEM): The TGWP method reproduces the transient deformation of the wave packet and the final momentum space distribution with high fidelity. While the method produces artificial oscillations in regions of negligible wave function intensity (contributing to $L^2$ error), the spectral envelope, sideband positions, and intensity redistribution match the grid-based solution closely. Peak height errors remained below 1.2%, and width errors below 2.3%.
- Stimulated Compton Scattering: The agreement was excellent, with the TGWP method accurately capturing the coherent redistribution of the wave function in momentum space and the oscillatory interference structures. The regularized error remained below 1% throughout propagation.
- Robustness: The method remained accurate even when parameters deviated from ideal phase-matching conditions, correctly reproducing physical trends in off-resonant regimes.
Performance:
- The mesh-free approach demonstrated significant computational speedups. In the near-field case, a speedup of approximately 100x was observed on comparable hardware.
- In the stimulated Compton scattering case, the speedup reached approximately 800x (TGWP taking ~1 second vs. ~14 minutes for the split-step method per parameter set).
- The method scales efficiently by eliminating the need for spatial meshing, which is the primary bottleneck in grid-based coupling with Maxwell solvers.

Significance and Claims
The paper establishes the TGWP method as a versatile, efficient, and scalable computational tool for quantum electron optics. Its primary contributions are:

Beyond Approximations: It enables the simulation of free-electron-light interactions beyond the eikonal and no-recoil approximations, capturing full minimal-coupling dynamics.
Efficiency: By removing the spatial mesh requirement, the method overcomes the discretization scale mismatch that plagues hybrid Schrödinger-Maxwell simulations, offering speedups of up to two orders of magnitude.
Validation: The method is rigorously validated against a grid-based benchmark for both localized near-field interactions and free-space scattering, showing excellent quantitative agreement.

The authors conclude that this framework provides a solid foundation for modeling complex time-dependent electromagnetic environments. They suggest that, combined with optimization algorithms and machine learning, this approach could facilitate the inverse design of electromagnetic environments for the on-demand shaping and coherent control of electron wave packets. The work also positions TGWP as a promising candidate for resolving full Maxwell–Schrödinger dynamics, particularly given its prior successful application to magnetic (minimal coupling) Hamiltonians.

Technical Summary: Meshfree versus Grid-based Schrödinger Solvers for Modeling Free-Electron Wave Packet Interactions with Light

More like this