Theoretical Analysis and PIC Simulations of Electromagnetic Wakefields Excited by Relativistic Beams in Magnetized Plasmas

This study develops a fully causal three-dimensional Green function formalism and validates it with EPOCH particle-in-cell simulations to demonstrate how external axial magnetic fields and plasma density parameters enhance and hybridize electromagnetic wakefields excited by relativistic electron beams in magnetized plasmas.

The Gist

The Big Picture: Surfing on a Magnetic Wave

Imagine you are trying to push a heavy boat (a beam of electrons) through a calm lake (a plasma). Usually, the boat creates a wake behind it, like the waves a speedboat leaves. In the world of particle physics, scientists want to use these "wakes" to push other particles forward, giving them a massive energy boost. This is called a plasma wakefield accelerator.

This paper asks a specific question: What happens if we put the lake inside a giant, invisible magnetic tunnel?

The authors, Ali Asghar Molavi Choobini and Mehran Shahmansouri, built a mathematical model and ran computer simulations to see how adding a magnetic field changes the shape, strength, and behavior of these waves.

The Two Tools They Used

To solve this puzzle, the team used two different methods, like checking a map with both a compass and a GPS:

1. The Mathematical Map (Green's Function): They developed a new, complex set of equations. Think of this as a perfect, theoretical recipe that predicts exactly how the water (plasma) should ripple when a boat (electron beam) passes through, especially when a magnetic field is pulling the water sideways.
2. The Computer Movie (PIC Simulations): They used a powerful computer code called EPOCH to create a 3D movie of the event. They simulated millions of tiny particles interacting to see if the "movie" matched their "recipe."

What They Discovered

Here are the main findings, explained through analogies:

1. The Magnetic Field Acts Like a "Stiffener"
In a normal lake, the water ripples at a certain speed. But when they added the magnetic field, it was like turning the water into a stiff gel.

• The Result: The waves started vibrating much faster (higher frequency).
• The Analogy: Imagine plucking a loose guitar string versus a tight one. The tight string (the magnetized plasma) vibrates faster and with more energy. The magnetic field made the "restoring force" (the force that tries to snap the water back to calm) much stronger.

2. The Waves Get a "Sidekick" (Hybrid Motion)
Normally, the wake moves mostly forward and backward. But with the magnetic field, the water starts swirling sideways too.

• The Result: The forward motion and the sideways motion became linked. You can't have one without the other anymore.
• The Analogy: Think of a dancer. Without the magnetic field, they just march forward. With the field, they are forced to march forward and spin in circles at the same time. The paper calls this a "hybrid" mode.

3. The "Focusing" Effect Gets Stronger
One of the goals of these accelerators is to keep the beam of particles from spreading out (like a flashlight beam that gets too wide).

• The Result: The magnetic field created much stronger "focusing" forces. It acted like a pair of invisible hands squeezing the beam back together.
• The Analogy: Without the magnet, the wake is like a gentle breeze. With the magnet, the wake acts like a vacuum cleaner hose, pulling particles tightly toward the center.

4. The Shape of the Boat Matters
They tested different shapes for the "boat" (the electron beam).

• Sharp vs. Smooth: If the boat had a sharp, sudden edge (like a square block), it created wild, choppy waves with lots of ringing. If the boat was smooth and rounded (like a teardrop), the waves were smoother and calmer.
• The Finding: Sharper edges on the beam create stronger, more energetic waves, but they also create more "noise" (oscillations) behind the beam.

5. Speed and Density

• Speed: If the boat was moving slowly, the waves were messy and weak. But once the boat reached "ultra-relativistic" speeds (close to the speed of light), the waves settled into a perfect, universal pattern. It didn't matter how much faster they went after that; the wave pattern stayed the same.
• Density: If the water was thicker (higher plasma density), the initial wave was huge and powerful, but it died out (damped) very quickly. If the water was thinner, the wave lasted longer but was weaker.

The Bottom Line

The paper proves that by adding an external magnetic field, scientists can fundamentally change how plasma wakes behave.

• They can make the waves stronger and faster.
• They can create a tighter focus for the particles.
• They can mix the forward and sideways movements into a single, powerful hybrid wave.

The authors confirmed that their mathematical "recipe" perfectly matched their computer "movie." This means they now have a reliable tool to design future accelerators that use magnetic fields to get better results, provided they can control the density of the plasma and the shape of the electron beam.

Note: The paper focuses entirely on the physics of how these waves are created and shaped. It does not discuss using these results for medical treatments, specific future machines, or clinical applications; it is purely about understanding the mechanics of the wakefields themselves.

Technical Summary

1. Problem Statement

Plasma Wakefield Acceleration (PWFA) is a promising technology for achieving high-energy particle acceleration in compact devices. However, conventional models often rely on simplified one-dimensional approximations or assume unmagnetized plasmas. In reality, external magnetic fields can significantly alter plasma dynamics, yet a unified, fully causal, and three-dimensional theoretical framework describing the coupled longitudinal and radial wakefields in magnetized cold plasmas has been lacking. The specific challenges addressed include:

• Understanding how an external axial magnetic field ($B_0$) couples longitudinal charge-separation dynamics with cyclotron-driven transverse electron motion.
• Determining how magnetization modifies effective restoring forces, wake amplitudes, and the radial focusing-defocusing structure.
• Validating theoretical predictions against realistic, high-resolution numerical simulations across a wide range of parameters (density, magnetic field strength, beam profile, and Lorentz factor).

2. Methodology

The study employs a dual approach combining rigorous analytical theory with high-fidelity numerical simulations.

A. Theoretical Framework

• Governing Equations: The authors start from linearized Maxwell-fluid equations in the presence of a magnetized plasma dielectric tensor.
• Green's Function Formalism: A fully causal, three-dimensional Green-function solution is derived. This approach avoids simplified 1D models and explicitly accounts for the full 3D plasma response.
• Mathematical Techniques:
  • Fourier Transform: Applied to the co-moving coordinate ($\xi = v_b t - z$) to transform wave equations into spectral space.
  • Hankel Transform: Used for radial variables to handle the cylindrical symmetry of the beam.
  • Dielectric Tensor: The magnetized plasma response is captured via the dielectric tensor $\boldsymbol{\epsilon}$, which includes the cyclotron frequency ($\omega_c$) and plasma frequency ($\omega_p$).
• Beam Models: Analytical solutions are derived for two specific charge distributions:
  1. A finite-radius, ultra-short "top-hat" beam slice.
  2. An idealized point-like charge (Green's function source).

B. Numerical Simulations

• Code: Extensive 3D Particle-in-Cell (PIC) simulations were performed using the EPOCH code.
• Setup:
  • Domain: Fully 3D Cartesian grid ($600 \times 320 \times 320$ cells).
  • Plasma: Cold, uniform, quasi-neutral electron-ion medium with an imposed static axial magnetic field ($B_0$).
  • Driver: Relativistic electron beams with varying Lorentz factors ($\gamma$), transverse radii ($r_d$), and current profiles (Gaussian, top-hat, exponential, etc.).
  • Parameters: Plasma densities ranging from $10^{17}$ to $10^{20} \text{ m}^{-3}$ and magnetic fields up to 20 T.
• Validation: Simulation results (wakefield components $E_z, E_r, B_\theta$) were directly compared against the analytical Green-function solutions.

3. Key Contributions

• Unified Hybrid Eigenmode: The study identifies and characterizes a hybrid electrostatic-electromagnetic eigenmode where longitudinal charge separation is strongly coupled to cyclotron-induced transverse motion. This mode does not exist in unmagnetized limits.
• Causal 3D Green's Function: Development of a rigorous, causal 3D Green-function formalism that captures the complete electromagnetic response, including transverse plasma currents induced by the magnetic field.
• Quantitative Validation: The paper provides the first direct numerical confirmation of several theoretical predictions, including the density-dependent amplification of wake amplitudes and the emergence of high-frequency radial oscillations under strong magnetization.
• Parameter Sensitivity Maps: Systematic analysis of how wakefields respond to variations in plasma density, magnetic field strength, beam Lorentz factor, and longitudinal current profiles.

4. Key Results

Effect of External Magnetic Field ($B_0$)

• Hybridization: Magnetization couples longitudinal and radial wake components, which behave as independent modes in unmagnetized plasmas.
• Frequency Shift: The oscillation frequency shifts from $\omega_p$ to the hybrid frequency $\sqrt{\omega_p^2 + \omega_c^2}$.
• Transverse Current: In unmagnetized plasmas, an axisymmetric beam induces no transverse current. With $B_0 > 0$, a significant transverse plasma current is generated, leading to strong radial focusing forces.
• Amplitude Enhancement: Increasing $B_0$ enhances the initial wake amplitude and creates pronounced focusing-defocusing cycles. The radial wakefield, zero in the unmagnetized case, grows significantly with $B_0$.
• Optimal Regime: Maximum sensitivity to magnetization occurs at moderate fields ($B_0 \approx 1\text{--}2$ T), where cyclotron-plasma coupling is maximized.

Effect of Plasma Density ($n_e$)

• Amplitude Scaling: Increasing plasma density substantially amplifies the initial wake amplitude due to stronger electrostatic restoring forces.
• Damping: Higher densities lead to accelerated damping of higher-order oscillations, resulting in a shorter wake coherence length but a stronger initial accelerating gradient.

Effect of Beam Parameters

• Lorentz Factor ($\gamma$): As $\gamma$ increases, the wakefield converges to a universal ultrarelativistic structure. Longitudinal fields require very high $\gamma$ ($>10^3$) to stabilize, while radial fields converge slightly earlier.
• Beam Radius ($r_d$): Narrow beams produce highly localized, sharply peaked longitudinal fields with weak radial components. Wide beams generate broader, slower-decaying longitudinal fields with strong, extended radial focusing structures.
• Current Profile: Sharper longitudinal profiles (e.g., top-hat, exponential) with abrupt edges excite stronger high-k plasma modes, leading to higher peak fields and more total wake energy compared to smooth profiles (e.g., Gaussian).

Point-Like vs. Finite Beams

• Point-like drivers excite fundamental Green-function responses, revealing the intrinsic spatial structure of the plasma eigenmode. Finite beams are shown to be convolutions of this point-source response with the beam density profile.

5. Significance and Implications

• Design of Next-Gen Accelerators: The framework provides a robust predictive tool for designing magnetized plasma-based accelerators. It demonstrates that external magnetic fields can be used as a tunable control parameter to simultaneously enhance accelerating gradients, focusing strengths, and beam stability.
• Beam Quality and Stability: The ability to strengthen focusing forces via magnetization offers a pathway to mitigate beam breakup instabilities and preserve emittance in high-energy accelerators.
• Theoretical Benchmark: The work establishes a new conceptual benchmark for understanding wakefield excitation, bridging the gap between simplified 1D theories and complex 3D reality.
• Practical Optimization: By identifying optimized operating regimes (e.g., specific $B_0$ and $n_e$ combinations), the study guides the development of compact, high-gradient accelerators for applications in high-energy physics, free-electron lasers, and medical technologies.

In conclusion, this paper successfully validates that external magnetic fields fundamentally reshape plasma wakefield dynamics, creating a hybrid mode that offers superior control over acceleration and focusing compared to conventional unmagnetized schemes.