QED cross sections in strong magnetic fields

This paper presents a novel formalism that resums magnetic field interactions and incorporates finite decay widths of Landau levels to systematically calculate tree-level QED scattering cross sections relevant to magnetar plasma dynamics, with results made available in an open-source Python package.

The Gist

Imagine the universe as a giant, chaotic dance floor. Usually, the dancers (particles like electrons and photons) move around freely, bumping into each other in predictable ways. Physicists have a very precise rulebook for how they dance, called Quantum Electrodynamics (QED). It's like a perfect choreography that has been tested and proven to work incredibly well in our everyday, "empty" space.

But what happens if you turn the music up to deafening levels and flood the dance floor with an invisible, crushing force? That's what happens around Magnetars.

The Setting: The Magnetar Dance Floor

Magnetars are a special type of dead star (neutron star) with magnetic fields so strong they break the rules. Their magnetic fields are millions of times stronger than anything we can create on Earth. In fact, they are so strong that they push physics into a "nonlinear" zone.

Think of it like this: In normal space, if you push a swing, it moves a little. In a Magnetar's magnetic field, if you push a swing, it might suddenly turn into a trampoline, a rollercoaster, or a portal to another dimension. The standard rulebook (QED) stops working because the background "music" (the magnetic field) is too loud and too heavy.

The Problem: Broken Rules and Broken Math

For a long time, scientists trying to simulate what happens inside these stars had to use a "cheat code." They assumed that all the electrons were stuck on the very bottom step of a ladder (called the Lowest Landau Level).

Imagine a staircase where the bottom step is the only one that exists.

• The Old Way: Scientists assumed every electron was stuck on step 1. If an electron tried to jump to step 2, they just ignored it.
• The Reality: In these super-strong fields, electrons do jump to higher steps. They get excited, they vibrate, and they interact with the magnetic field in complex ways. By ignoring the upper steps, the old simulations were like trying to predict a hurricane by only looking at the ground floor. They were getting the dynamics wrong, sometimes by huge margins.

The Solution: A New, Complete Rulebook

The authors of this paper (Olavi Kiuru and his team) decided to throw away the cheat code. They built a brand-new, comprehensive framework to calculate exactly how particles dance in these super-strong magnetic fields.

Here is how they did it, using some simple analogies:

1. The "All-Seeing" Propagator
In physics, a "propagator" is like a map showing how a particle gets from Point A to Point B. In normal space, the map is a straight line. In a Magnetar, the magnetic field is so strong that the particle doesn't just go straight; it spirals, bounces, and interacts with the field constantly.
The authors created a "super-map" that accounts for every possible interaction with the magnetic field at once. Instead of drawing one line, they drew a tangled web of every possible path the particle could take, summing them all up. This ensures they don't miss any of the "dance moves."

2. The "Fuzzy" Steps (Decay Widths)
In the old math, if an electron hit a specific energy level, the math would break and give an answer of "infinity." It was like a song that got stuck on a single note and screamed forever.
The authors realized that in reality, these energy levels aren't sharp, perfect lines; they are "fuzzy" or "blurred" because the particles are constantly decaying and changing. They added a "fuzz factor" (called decay width) to their calculations. This smoothed out the infinite spikes, turning them into realistic, manageable peaks. It's like taking a laser beam and spreading it out so it illuminates a room instead of burning a hole in the wall.

3. The Spin and the Polarization
Particles have a property called "spin" (like a tiny top spinning) and "polarization" (the direction of their vibration). In these strong fields, the direction they spin matters immensely.
The authors didn't just calculate the average; they calculated exactly what happens if an electron is spinning "up" versus "down." They found that the magnetic field acts like a bouncer at a club, letting only certain spins in and blocking others, which drastically changes how particles collide.

The Result: A New Toolkit for the Universe

The team didn't just write equations; they turned their entire discovery into a free, open-source software package (a Python library).

Think of this as giving astrophysicists a new, high-definition video game engine.

• Before: They were playing a game with low-resolution graphics and broken physics.
• Now: They have a tool that simulates the "Magnetar Dance Floor" with perfect accuracy, accounting for every step of the ladder, every spin of the dancer, and every pulse of the magnetic field.

Why Does This Matter?

Magnetars are some of the most energetic objects in the universe. They blast out X-rays and gamma rays that we can detect from Earth. To understand why they blast out this energy, we need to know exactly how the particles inside them are colliding and creating new particles.

If we use the old, broken math, we might think a Magnetar is quiet when it's actually screaming, or that it's creating a storm of particles when it's actually calm. This new paper provides the correct "physics engine" to finally understand the violent, beautiful, and chaotic lives of these cosmic giants.

In short: They took a broken, simplified model of the universe's most extreme environments and replaced it with a complete, accurate, and free-to-use simulation tool, allowing us to finally see the true dance of matter in the strongest magnetic fields in existence.

Technical Summary

1. Problem Statement

The magnetospheres of magnetars (highly magnetized neutron stars) possess magnetic fields ($B$) exceeding the Schwinger limit ($B_Q \approx 4.41 \times 10^{13}$ G). In this regime, Quantum Electrodynamics (QED) becomes nonlinear, and standard vacuum QED approximations fail.

• Current Limitations: Existing studies of magnetar plasma dynamics often rely on:
  • Lowest-Landau-Level (LLL) Approximation: Restricting both external and virtual fermions to the ground state ($n=0$). This fails at high energies where excited Landau levels ($n > 0$) dominate scattering.
  • Vacuum Cross Sections: Using vacuum QED formulas which ignore the qualitative modifications of the energy spectrum and interaction vertices caused by the background field.
  • Incomplete Scope: Previous works typically focus on single processes (e.g., Compton scattering) or make assumptions that limit applicability (e.g., using Johnson-Lippmann spinors that obscure spin dependence).
• Consequence: These approximations lead to interaction rates that may be orders of magnitude too small, resulting in qualitatively incorrect simulations of magnetar plasma dynamics and particle cascades.

2. Methodology

The authors apply a formalism originally developed for studying magnetic catalysis in hot quark-gluon plasma to strong-field QED (SFQED). The core of their approach involves a consistent treatment of the background field to all orders.

• Furry Picture Formalism: The gauge field is split into a non-interacting background part ($A_\mu$) and an interaction part. The fermion propagator is derived from the non-interacting Lagrangian, capturing all interactions with the background field.
• Resummed Fermion Propagator: Instead of perturbative expansions, they use a mixed representation propagator (Fourier transformed in $t$ and $z$, but not $x, y$) that includes a sum over all Landau levels ($n$) of the virtual fermion.
• Inclusion of Decay Widths: To regulate unphysical divergences (resonances) where virtual fermions go on-shell, the authors calculate the imaginary part of the fermion self-energy ($\Sigma$). This finite decay width ($\Gamma$) is incorporated into the dressed propagator, turning infinite resonance peaks into finite-width resonances.
• Spin and Polarization:
  • Fermions: They utilize Sokolov-Ternov (ST) wave functions, which are preferred in astrophysics for their Lorentz transformation properties and independent decay of spin states.
  • Photons: Photons are treated as having two polarization modes (Ordinary/O-mode and Extraordinary/X-mode). The cross sections are calculated for specific polarization states.
• Scope of Processes: The paper provides a systematic analysis of all tree-level QED scattering processes relevant to magnetar plasmas (1-to-2, 2-to-1, and 2-to-2) with no internal photon lines. This includes Synchrotron radiation, Pair creation/annihilation, Compton scattering, and electron-electron/positron scattering.

3. Key Contributions

• Systematic Framework: The first systematic derivation of tree-level QED cross sections in strong magnetic fields without restricting fermions to the LLL.
• Beyond LLL: The method explicitly includes contributions from excited Landau levels in both external and virtual fermions, correcting the high-energy behavior where LLL approximations fail.
• Spin-Dependent Results: Unlike previous works using spin-summed Johnson-Lippmann spinors, this work provides spin-dependent cross sections, revealing that specific spin configurations (e.g., spin-down electrons and spin-up positrons) dominate certain scattering channels.
• Open-Source Implementation: All derived analytical expressions and their numerical implementations are provided in an open-source Python package, allowing for direct integration into Monte Carlo and plasma simulation codes.

4. Key Results

• Compton Scattering:
  • The authors' results match existing literature (e.g., Mushtukov et al.) at low energies but diverge significantly at higher photon energies ($\omega \gtrsim 2.5 m_e$).
  • The LLL approximation is shown to be valid only for $\omega \lesssim 2.5 m_e$. Beyond this, neglecting excited Landau levels in the propagator leads to unphysically transparent photon targets.
  • The inclusion of decay widths resolves infinite resonance peaks, with peak heights inversely proportional to the spin-averaged decay widths.
• Two-Photon Pair Creation:
  • The cross section exhibits a strong dependence on the spin states of the produced fermions.
  • For ground-state collisions, only specific spin combinations (spin-down electron, spin-up positron) are produced.
  • Comparison with Kozlenkov and Mitrofanov shows agreement in the LLL regime but significant deviations when excited Landau levels of external fermions become relevant.
• One-Photon Pair Annihilation:
  • The cross section is non-monotonic with respect to the longitudinal momentum ($p_z$) and the Landau levels of the incoming pair.
  • While the LLL dominates at low momenta ($p_z \lesssim 3m_e$), higher Landau levels dominate at larger momenta, challenging the assumption that excited states can be ignored in relativistic regimes.
• Resonance Structure: The paper clarifies that resonance peaks in cross sections correspond to the opening of new scattering channels involving specific Landau levels. The finite width of these peaks is determined by the fermion self-energy.

5. Significance

• Astrophysical Modeling: The results provide the necessary microphysical input for simulating magnetar magnetospheres. Accurate cross sections are crucial for modeling pair cascades (exponential production of $e^\pm$ pairs), which drive the plasma dynamics and hard X-ray emission of magnetars.
• Correction of Previous Models: By demonstrating that LLL approximations fail for high-energy particles ($p_\perp \gtrsim m_e \sqrt{b}$), the paper corrects interaction rates used in recent simulations, potentially altering predictions for magnetar emission spectra and plasma stability.
• Generalizability: The formalism is robust and can be extended to non-zero temperatures, chemical potentials, and different background field configurations (e.g., linearly polarized laser fields in high-intensity laser physics).
• Resource Availability: The release of the Python package democratizes access to complex SFQED calculations, enabling the broader community to move beyond vacuum approximations in astrophysical simulations.

In summary, this work bridges the gap between theoretical strong-field QED and practical astrophysical simulation by providing a consistent, spin-resolved, and Landau-level-complete framework for calculating scattering cross sections in the extreme environments of magnetars.