🔬 atomic physics

Temperature and mean axial momentum vs. laser intensity of electrons released from O $_2$ by an 800 nm ultrashort pulsed laser

This paper presents a semi-empirical model that modifies theoretical kinetic energy spectra with two adjustable parameters to describe the thermalized temperature and mean axial momentum of electrons released from O $_2$ by an 800 nm ultrashort pulsed laser as a function of peak intensity, thereby providing essential initial conditions for electrodynamic fluid simulations.

Original authors: Edward L. Ruden

Published 2026-07-16

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Edward L. Ruden

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 a world where light doesn't just shine; it punches. When a laser beam is squeezed tight enough and fired fast enough, it can rip electrons right out of the atoms they are glued to. This happens in a flash, faster than a blink, creating a chaotic soup of free-floating particles called plasma. Scientists are fascinated by this because these plasma "filaments" act like invisible wires in the sky, conducting electricity and even beaming out radio waves and terahertz radiation. To understand how these invisible wires work, we need to know two things about the electrons inside them: how hot they are (their temperature) and how fast they are moving in the direction the laser is pointing (their momentum). If we get these numbers wrong, our predictions about how the plasma behaves will be as accurate as guessing the weather by looking at a cloud from a mile away.

This paper is a detective story about those electrons, specifically the ones ripped from oxygen molecules (O2) by a powerful, 800-nanometer laser pulse. The author, E. L. Ruden, is trying to build a better map for computer simulations. The goal is to predict the "initial conditions" of the electron soup right after the laser passes but before the electrons have a chance to calm down and mix together. The paper takes a known theoretical model, which acts like a rough sketch of how electrons should behave, and tweaks it with two "adjustable knobs" to match real-world experimental data. The main finding is that the electrons don't just follow the simple rules of the laser's push; they often crash back into their parent atoms (a process called rescattering) or recombine with them. This crash-and-burn behavior creates a "plateau" of high-energy electrons that the simple theory missed, and it also suppresses the number of very low-energy electrons. By accounting for these crashes, the author creates a new, semi-empirical model (called SFA2) that accurately predicts the electron temperature and their forward momentum across a wide range of laser intensities.

The Laser, The Atom, and The Crash

Think of an oxygen molecule as a tiny, two-story house with a very stubborn tenant: an electron. When a super-strong laser pulse hits this house, the electric field of the light acts like a giant, invisible hand trying to pull the tenant out. In the simplest version of this story (a theory called SFA0), the electron is pulled out, and then the laser just pushes it away, like a surfer riding a wave. The theory predicts exactly how fast the electron should go and how much energy it should have.

However, the real world is messier. Sometimes, after the electron is pulled out, the laser's hand swings back, and the electron gets thrown back toward the house. It might crash into the parent atom (the house) and bounce off, or it might get sucked back in and recombine. The paper suggests that these "rescattering" events are the missing piece of the puzzle.

The author found that the simple theory (SFA0) made two big mistakes when compared to actual data:

It predicted a huge crowd of very slow, low-energy electrons that simply don't exist in the real measurements.
It missed a group of very fast, high-energy electrons that appear in the data.

To fix this, the author introduced two "adjustable parameters," which are like fine-tuning knobs on a radio.

Knob 1: The Ceiling
The simple theory predicted a "surge" of low-energy electrons. The author realized this surge is likely because those electrons crash back into their parent atoms and get stuck (recombine) or bounce off in a way that removes them from the low-energy count. To fix the model, the author placed a "ceiling" on the number of low-energy electrons allowed in the simulation. Any electron predicted to be below a certain energy (1.60 eV) is capped, effectively removing the fake surge. This new version is called SFA1.

Knob 2: The Energy Multiplier
Even with the ceiling, the model still didn't get the temperature right for the high-energy electrons. The data showed that electrons were hotter (more energetic) than even the "ceiling" model predicted. This is because the electrons that do crash into the parent atom and bounce off (rescatter) get a massive energy boost, like a pinball hitting a bumper. To account for this, the author added a second knob: an energy multiplier (a factor of 2.067). This stretches the energy scale of the model to match the "hotter" reality seen in experiments. This final, tweaked version is called SFA2.

The Results: A Hotter, Faster Crowd

When the author ran the numbers with these two adjustments, the model finally matched the experimental data. The paper shows that for laser intensities where the Keldysh parameter (a number that tells us how "strong" the field is) is between 0.82 and 1.30, the new model predicts the electron temperature very accurately.

The paper also calculated the "mean axial momentum," which is the average speed of the electrons moving in the same direction the laser is traveling. This is crucial because this forward motion is what creates the electrical current in the plasma filament. The author found that the electrons moving forward are significantly faster than the simple theory predicted. This is because the electrons that get a "rescatter" boost don't just go faster; they also get a kick in the forward direction.

The paper explicitly notes that this model works best for linearly polarized light (where the electric field wiggles in one flat plane). While the author tried to apply the same logic to circularly polarized light (where the field spins like a propeller), there wasn't enough data to be sure, so the results for that case are left as a "preliminary" guess for future study.

How Sure Are We?

The author is quite confident in the temperature predictions for the specific range of laser intensities mentioned (0.82 ≤ γ0 ≤ 1.30), because this range is based on interpolating between actual data points. However, the paper is more cautious about the momentum calculations, especially for lower laser intensities (γ0 < 1). The author admits that if the laser is so strong that it strips away half the oxygen molecules in the path, the simple math might get a bit wobbly.

Furthermore, the paper suggests that the "surge" of low-energy electrons seen in the simple theory is likely due to electrons recombining with their parents, but it stops short of saying this is the only reason. It suggests that the "rescatter" mechanism is the key to the high-energy "plateau," but acknowledges that other complex effects, like quantum interference between the two atoms in the oxygen molecule, are not included in this model.

In short, this paper doesn't claim to have solved the entire mystery of laser-plasma interaction. Instead, it offers a much better "semi-empirical" map—a blend of theory and real-world tuning—that helps scientists predict how hot and how fast electrons will be when a laser zaps oxygen. This map is essential for building accurate computer simulations of the glowing filaments that form in the air, which could one day help us understand everything from lightning to new types of wireless communication.

Technical Summary: Temperature and Mean Axial Momentum vs. Laser Intensity of Electrons Released from O₂ by an 800 nm Ultrashort Pulsed Laser

Problem Statement
This paper addresses the need for accurate initial conditions—specifically thermalized electron temperature ( $T$ ) and mean post-optical axial momentum ( $\langle p_{fz} \rangle$ )—for electrodynamic fluid simulations of laser filaments in air. While previous work has determined electron density and internal laser intensity based on filament conductivity measurements, a robust model linking these parameters to the peak laser intensity ( $I_0$ ) of a focused 800 nm ultrashort pulsed laser (USPL) interacting with oxygen ( $O_2$ ) was lacking. The challenge lies in reconciling theoretical strong-field approximation (SFA) models with empirical electron kinetic energy spectra, which exhibit features (such as low-energy surges and high-energy plateaus) not predicted by standard non-adiabatic SFA theories due to complex electron-ion interactions like rescattering and recombination.

Methodology
The authors present a semi-empirical model built upon a reformulated non-adiabatic Strong Field Approximation (SFA0), originally formulated by Li et al. and Luo et al. The model assumes a time-dependent but steady-state laser pulse envelope and treats the ionization process quasiclassically, considering only the most probable tunneling path for a given instantaneous ionization time $t_0$ .

To bridge the gap between SFA0 theory and experimental data, the authors introduce two phenomenological modifications to the theoretical electron kinetic energy spectrum $S_0(U)$ :

SFA1 (Ceiling Model): To address a theoretical surge in electron population at low kinetic energies ( $U < U_c$ ) that is absent in experimental data, a "ceiling" is imposed on the spectrum at a critical energy $U_c = 1.60$ eV. This accounts for the suppression of low-energy electrons due to recombination with the parent ion.
SFA2 (Energy Rescaling): To match the higher thermalized temperatures observed in 4 $\pi$ steradian averaged data, an energy multiplier $\zeta = 2.067$ is applied to the SFA1 spectrum. This accounts for the energy boost electrons receive from elastic rescattering off the parent ion within an optical cycle.

The mean axial momentum $\langle p_{fz} \rangle$ is estimated using classical kinematics of rescattering. The model calculates the time of electron return to the parent ion ( $t_1$ ) and the subsequent momentum components after elastic scattering at an angle $\phi$ . The total momentum is derived by integrating over the pulse profile and solid angles, utilizing empirical temperature data ( $T_x$ , $T_+$ , $T_*$ ) recorded in different directions to weight the angular distribution.

Key Contributions

Refined Spectral Models: The paper extends the standard SFA0 by introducing SFA1 and SFA2, which successfully fit empirical kinetic energy spectra for $O_2$ exposed to 800 nm linearly polarized light across a range of Keldysh parameters ( $\gamma_0$ ).
Thermalized Temperature ( $T$ ) vs. Intensity: The authors provide a semi-empirical relationship for the thermalized electron temperature $T$ as a function of peak intensity $I_0$ (or $\gamma_0$ ). The model distinguishes between temperatures derived from x-direction emission ( $T_x$ ), 2 $\pi$ plane averages ( $T_+$ ), and 4 $\pi$ steradian averages ( $T_*$ ).
Mean Axial Momentum ( $\langle p_{fz} \rangle$ ) Estimation: A method is developed to estimate $\langle p_{fz} \rangle$ by combining the SFA2 thermalized temperature with classical rescattering kinematics. This includes a correction factor $\alpha(\gamma_0)$ to account for the spatiotemporal averaging of the laser pulse intensity, distinguishing between the properties of a pulse peaking at $I_0$ and the intrinsic properties at a specific intensity.
Polarization Analysis: The model compares linear ( $\epsilon=0$ ) and circular ( $\epsilon=1$ ) polarization. While SFA0 predicts higher temperatures for circular polarization, the empirical data (limited to linear polarization) suggests that rescattering effects in the linear case lead to comparable or higher effective temperatures, negating the anticipated advantage of circular polarization for microwave emission.

Results

Spectral Fits: The SFA2 model ( $S_2(U)$ ) provides a good fit to experimental spectra averaged over 4 $\pi$ steradians for $0.82 \leq \gamma_0 \leq 1.30$ . The model successfully reproduces the "plateau" in the high-energy tail observed in experiments, attributed to rescattering.
Temperature Trends: The thermalized temperature $T_2$ (from SFA2) is found to be a useful approximation down to $\gamma_0 = 0.47$ , though accuracy decreases for $\gamma_0 > 1.40$ where the multiphoton ionization regime dominates and SFA assumptions break down.
Momentum Calculations: The calculated mean axial momentum $\langle p_{fz} \rangle$ is shown to be significantly influenced by rescattering. For $\gamma_0 = 1$ , the rescattering contribution ( $\langle p_{fz} \rangle_b$ ) is approximately 0.80 $\times 10^{-27}$ kg-m/s, while the non-rescattering contribution ( $\langle p_{fz} \rangle_a$ ) is 1.66 $\times 10^{-27}$ kg-m/s. The paper notes that for $\gamma_0 < 1$ , significant ionization depletion of neutral $O_2$ may introduce errors in the momentum estimates derived from $T_x$ data.
Correction Factors: The analysis reveals that the intrinsic momentum at peak intensity is higher than the pulse-averaged momentum, with a correction factor $\alpha(\gamma_0) \approx 0.76$ for $\gamma_0 = 1$ . This factor decreases further for $\gamma_0 > 1$ .

Significance and Claims
The primary significance of this work, as claimed by the authors, is the provision of initial conditions ( $T$ and $\langle p_{fz} \rangle$ ) for electrodynamic fluid simulations of laser filaments in air. These conditions are derived from a synthesis of published experimental data and a modified SFA theory that accounts for electron rescattering and recombination.

The authors modestly frame their contribution as a tool for data analysis and simulation initialization rather than a fundamental breakthrough in strong-field theory. They acknowledge limitations, including:

The model is most accurate for $0.82 \leq \gamma_0 \leq 1.30$ .
The treatment of $\epsilon=1$ (circular polarization) is preliminary due to a lack of spectral data.
Effects such as Coulomb focusing, quantum interference between the two atoms in $O_2$ , and atomic resonances at high $\gamma$ are not included.
Significant errors may exist for $\gamma_0 < 1$ due to neutral depletion effects not fully captured in the current momentum estimation.

The paper serves as a complementary component to a broader study (referenced as [1] and [2]) aimed at determining the physical properties of air filaments formed by USPLs, specifically linking conductivity measurements to electron density and temperature.

The Laser, The Atom, and The Crash

The Results: A Hotter, Faster Crowd

How Sure Are We?

More like this