Evaluation of QED cross sections in strong magnetic fields

This paper presents a formalism and computational techniques for calculating tree-level QED scattering cross sections in strong magnetic fields exceeding the Schwinger limit, implementing the results in an open-source Python package to support future studies of magnetar magnetospheric dynamics.

The Gist

Imagine the universe as a giant, bustling kitchen. Usually, the ingredients (particles like electrons and photons) move around freely, bumping into each other and reacting in predictable ways. This is how physics works in our everyday world.

But in the neighborhood of a Magnetar—a type of neutron star with a magnetic field so strong it could wipe a credit card from halfway across the galaxy—the kitchen rules change completely. The magnetic field is so intense that it acts like a giant, invisible grid or a set of train tracks that forces the particles to move in very specific, rigid lanes.

This paper, written by Olavi Kiuru, is essentially a new cookbook for this extreme kitchen. Here is a breakdown of what the author did, using simple analogies:

1. The Problem: The "Traffic Jam" of Physics

In normal physics, we have a standard rulebook (Quantum Electrodynamics, or QED) that tells us how particles interact. But when you turn on a magnetar's magnetic field, the old rulebook stops working. The field is so strong (over 44 trillion times stronger than a fridge magnet) that it forces electrons and positrons (anti-electrons) to get stuck in specific energy "rungs" on a ladder, called Landau levels.

Think of it like this: In a normal room, you can walk anywhere. In a magnetar's room, you are forced to walk only on specific, invisible stepping stones. If you try to jump off a stone, the physics gets weird and nonlinear. Scientists knew this was happening, but they didn't have a complete list of how particles bounce off each other on these stepping stones. Without this list, they couldn't accurately simulate what happens inside these stars.

2. The Solution: A New "Feynman Map"

The author developed a new mathematical method to calculate exactly what happens when these particles collide.

• The Old Way: Previous studies were like looking at the kitchen through a keyhole. They only looked at one specific collision at a time or assumed the particles were stuck on the very bottom stepping stone (the "Lowest Landau Level").
• The New Way: Kiuru created a comprehensive map that works for all the stepping stones, from the bottom to the top. He used a technique called the Furry Picture, which is like upgrading the camera lens to see the entire kitchen at once, including the invisible magnetic tracks the particles are forced to follow.

3. The Recipe: Calculating the Collisions

The paper details the "recipes" for every possible collision that can happen in this environment without involving a photon traveling between particles (which is a bit like calculating the direct handshake between two people, rather than them passing a note through a third person).

The main collisions calculated include:

• Synchrotron Radiation: An electron sliding down a magnetic track and spitting out a photon (like a car skidding and throwing off sparks).
• Pair Creation: A high-energy photon hitting the magnetic field and splitting into an electron and a positron (like a billiard ball hitting a cushion and suddenly splitting into two balls).
• Compton Scattering: A photon bouncing off an electron (like a tennis ball hitting a racket).

4. The "Secret Sauce": Handling the Resonances

One of the hardest parts of this math is dealing with resonances. Imagine pushing a child on a swing. If you push at exactly the right time, the swing goes super high (resonance). In magnetars, particles can get stuck in these "swing" states where the math blows up to infinity.

The author figured out how to add a tiny bit of "friction" (decay rates) to the math to stop the numbers from exploding. This is like realizing that even a perfect swing has a little bit of air resistance, so it never goes to infinity. This makes the calculations realistic and usable for computer simulations.

5. The Result: An Open-Source Toolkit

The most practical outcome of this paper isn't just the math on the page; it's that the author turned all these complex calculations into a free, open-source Python package.

Think of this as the author taking his new cookbook, writing it into a computer program, and handing it to every astrophysicist in the world. Now, anyone can plug in the conditions of a magnetar and instantly get the correct "collision rates" to simulate how the plasma (the hot, charged gas) behaves around these stars.

Why Does This Matter?

Magnetars are cosmic laboratories. They are the only places in the universe where we can test physics under these extreme conditions. By understanding how particles interact here, we can:

1. Decode the Light: Explain the strange double-peaked light spectra we see coming from these stars.
2. Simulate the Storms: Better model the violent storms and particle cascades that happen in their atmospheres.
3. Test the Laws of Physics: See if our current understanding of the universe holds up under the most extreme pressure imaginable.

In short: This paper provides the missing instruction manual for how matter behaves in the most magnetic places in the universe, turning a chaotic, confusing mess of equations into a clear, usable tool for scientists to understand the most extreme stars in the cosmos.

Technical Summary

1. Problem Statement

Quantum Electrodynamics (QED) becomes nonlinear when magnetic field strengths exceed the critical Schwinger limit ($B_Q \approx 4.41 \times 10^{13}$ G). This regime is realized in the magnetospheres of magnetars (highly magnetized neutron stars), where fields can reach $30 B_Q$.

• The Gap: While QED processes like pair creation and synchrotron radiation are known to be crucial for magnetar plasma dynamics, current simulations lack a complete set of scattering cross sections. Existing studies often focus on specific processes, assume particles are confined to the Lowest Landau Level (LLL), or rely on approximations that limit their applicability across different field strengths ($b \ll 1, b \approx 1, b \gg 1$).
• The Goal: To calculate the cross sections for all tree-level 1-to-2, 2-to-1, and 2-to-2 particle QED scattering processes (excluding those with photon propagators) in a constant, strong background magnetic field, valid for all Landau levels and field strengths.

2. Methodology

The author employs a formalism originally developed for magnetic catalysis but adapted for astrophysical scattering. The approach treats the background magnetic field non-perturbatively while treating the interaction with the dynamical photon field perturbatively (Furry picture).

A. Theoretical Framework

• Background Field: A constant, homogeneous magnetic field $\mathbf{B} = (0,0,B)$ is assumed in the Landau gauge ($A^\mu = (0, 0, Bx, 0)$).
• Wave Functions: The author uses Sokolov-Ternov (ST) wave functions for electrons and positrons. These are eigenfunctions of the magnetic moment operator and preserve spin states under Lorentz boosts along the field, unlike the older Johnson-Lippmann (JL) wave functions.
• Propagator: The fermion propagator is derived using the Schwinger proper-time method but explicitly projected onto Landau levels.
  • Rationale: While the proper-time integral form is elegant, the Landau level projection is essential for calculating cross sections involving resonances. It allows for the explicit inclusion of decay rates (imaginary parts of the self-energy) to regulate infinite peaks caused by on-shell intermediate states.
• Cross Section Calculation:
  • Calculations are performed in a box of size $L$ with periodic boundary conditions, later taking the limit $L \to \infty$ to ensure gauge invariance.
  • The density of states for electrons/positrons is modified due to the lack of a well-defined $x$-momentum in the Landau gauge. The cross section is averaged over the unknown transverse momentum $p_y$.
  • Two reference frames are utilized: the Electron Rest Frame (ERF) and the Center-of-$z$-Momentum (CozM) frame.

B. Computational Techniques

• Integration: The scattering amplitudes involve integrals over transverse coordinates ($x, y$). These are solved using Gaussian integrals and properties of Hermite and Laguerre polynomials.
• Recursive Relations: To handle excited Landau levels for external particles, the author derives a recursive ansatz using the properties of Hermite polynomials and their derivatives with respect to momentum. This allows the generalization of results from the LLL approximation to arbitrary Landau levels.
• Decay Rates: To handle resonances in 2-to-2 scattering, the author calculates the imaginary part of the electron self-energy at 1-loop order. This provides a physically motivated decay rate ($\Gamma$) that shifts the pole of the propagator ($m \to m - i\Gamma/2$), replacing ad-hoc regularization methods.

3. Key Contributions

1. Comprehensive Formalism: The paper provides a unified, systematic framework for calculating QED cross sections in strong magnetic fields that is valid for all three regimes: weak ($b \ll 1$), critical ($b \approx 1$), and strong ($b \gg 1$).
2. Explicit Spin and Polarization Dependence: Unlike many previous works that average over spins, this formalism retains explicit dependence on the spin states of fermions and polarization modes (Ordinary O-mode and Extraordinary X-mode) of photons.
3. Generalization beyond LLL: The derivation explicitly handles transitions between arbitrary Landau levels ($n_i \to n_f$), rather than restricting calculations to the Lowest Landau Level (LLL) approximation, which is often insufficient for high-energy astrophysical scenarios.
4. Resonance Regulation: The paper rigorously implements the decay rate derived from the fermion self-energy to regulate resonances in the cross sections, offering a more robust physical description than previous phenomenological approaches.
5. Open-Source Implementation: All derived cross sections are implemented in a publicly available Python package (SFQED), making these complex calculations accessible for magnetar plasma simulations.

4. Results

The author successfully derived the scattering amplitudes and cross sections for the following processes:

• $O(\alpha_e)$ Processes (3-particle interactions):

  • Synchrotron Radiation (SR): $e^\pm \to e^\pm + \gamma$.
  • One-Photon Pair Creation (PC): $\gamma \to e^+ + e^-$.
  • One-Photon Pair Annihilation (PA): $e^+ + e^- \to \gamma$.
  • Synchrotron Self-Absorption (SSA): $e^\pm + \gamma \to e^\pm$.
  • Note: The SSA cross section contains a $\delta$-function enforcing energy conservation, requiring specific integer conditions on the Landau levels.

• $O(\alpha_e^2)$ Processes (4-particle interactions):

  • Compton Scattering: $e^\pm + \gamma \to e^\pm + \gamma$. The full analytical expressions for the scattering amplitudes (s-channel and u-channel) are provided, including contributions from all Landau levels.
  • Two-Photon Pair Creation: $\gamma + \gamma \to e^+ + e^-$.
  • Two-Photon Pair Annihilation: $e^+ + e^- \to \gamma + \gamma$.
  • Note: Møller and Bhabha scattering (involving photon propagators) are noted as future work.

Key Findings on Scaling:

• Pair Creation/Annihilation: The decay rates and cross sections for one-photon pair creation are amplified by a factor of $b$ (magnetic field strength) in strong fields, while one-photon annihilation is suppressed by $b^{-1}$.
• Compton Scattering: The cross section depends heavily on the Landau level transitions and the polarization modes. The results match known LLL approximations but extend validity to excited states.

5. Significance

• Astrophysical Modeling: This work fills a critical gap in magnetar physics. By providing a complete set of cross sections valid for strong fields, it enables more accurate simulations of magnetospheric plasma dynamics, particle cascades, and emission spectra (e.g., explaining the double-peak structure in magnetar spectra).
• Theoretical Advancement: The paper demonstrates the utility of the "magnetic catalysis" formalism in high-energy astrophysics and clarifies the validity of the Furry picture in strong fields. It also addresses the Ritus-Narozhny conjecture regarding the validity of perturbation theory in strong fields.
• Future Directions: The formalism is designed to be easily generalized to higher-order processes (loops) and 1-to-3 particle decays (e.g., photon splitting), which are vital for understanding low-energy particle cascades in magnetars where standard pair creation is kinematically forbidden.

In summary, this paper provides the rigorous mathematical and computational foundation necessary to move beyond simplified approximations in strong-field QED, directly impacting our understanding of the most extreme magnetic environments in the universe.