arxiv: 2605.05543 · v1 · submitted 2026-05-07 · ⚛️ physics.plasm-ph

Recognition: unknown

An Electromagnetic Particle-Particle Method for Relativistic Electron Bunch Dynamics from Early Expansion to Long-Range Transport

Yinjian Zhao , Yibo Liang , Yanan Zhang , Xiaochun Ma , Hui Liu

Authors on Pith no claims yet

Pith reviewed 2026-05-08 04:32 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph

keywords relativistic electron bunchesparticle-particle methodretarded fieldsmagnetospheric transportmesh-free simulationbunch divergenceelectromagnetic interactionsself-field effects

0 comments

The pith

The electromagnetic particle-particle method provides a mesh-free way to accurately simulate relativistic electron bunch dynamics from early self-field expansion to long-range geomagnetic transport.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Particle-mesh methods cannot retain exact pairwise electromagnetic interactions at scales below the mesh cell size. This limitation hinders accurate modeling of dense relativistic electron bunches, where self-interactions cause significant divergence. The work extends the electromagnetic particle-particle model to many particles by incorporating retarded electromagnetic fields, an improved calculation of time delays between particles, and a relativistic integrator for particle motion. A two-stage strategy is used to connect the initial phase dominated by the bunch's own fields with the later phase controlled by the external geomagnetic field. This approach would matter to readers interested in simulating the complete dynamics of such bunches in space environments where both close and distant forces play roles.

Core claim

The central claim is that the electromagnetic particle-particle method, which combines Liénard-Wiechert fields, an improved retarded-time evaluation procedure, and a relativistic particle pusher, together with a two-stage strategy to couple the dense early self-field-dominated evolution to the later long-range geomagnetic-field-controlled transport, provides a practical mesh-free approach for accurately simulating long-range transport of relativistic electron bunches when short-range electromagnetic interaction is important.

What carries the argument

The extended electromagnetic particle-particle (EM-PP) model that computes exact pairwise interactions via retarded-time fields and uses two-stage coupling for different regimes.

If this is right

It retains exact pairwise electromagnetic interactions without mesh-induced approximations.
It accurately models the initial expansion and divergence driven by self-fields in dense bunches.
It enables the transition to transport dominated by the geomagnetic field over long distances.
It offers a feasible alternative when both short-range and long-range effects must be treated together.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If successful, the method could be extended to study how bunch properties affect radiation or other secondary effects not directly simulated here.
One could test the scaling of accuracy with increasing particle numbers in controlled numerical experiments.
This framework might connect to problems in beam physics or space weather modeling where similar scale-bridging is needed.

Load-bearing premise

That the improved retarded-time evaluation procedure and the two-stage coupling can be stably extended to many particles without introducing uncontrolled numerical errors or prohibitive computational cost.

What would settle it

A side-by-side comparison of the method's predicted final bunch size and divergence after long-range transport with observational data from relativistic electron bunches in the magnetosphere would determine if the simulation matches reality.

Figures

Figures reproduced from arXiv: 2605.05543 by Hui Liu, Xiaochun Ma, Yanan Zhang, Yibo Liang, Yinjian Zhao.

**Figure 1.** Figure 1: , the retarded time satisfies |η| ≡ |r − w(tr)| = c(t − tr), (1) where η ≡ r − w(tr) (2) is the vector from the retarded source position to the field point. O x z y 𝒓 w(tr) 𝜼 Retarded position Particle trajectory Present position view at source ↗

**Figure 3.** Figure 3: Particle velocity sampling. 2.6 Geomagnetic Field Model and Injection Geometry The background magnetic field is prescribed and modeled as a static centered dipole field. The origin of the laboratory frame is placed at the center of the Earth, and the dipole moment is aligned with the laboratory z axis, M = M0ez, (77) with M0 = 8.6 × 1022 A m2 . (78) For numerical convenience, the geomagnetic field is evalu… view at source ↗

**Figure 2.** Figure 2: Particle spatial sampling view at source ↗

**Figure 4.** Figure 4: Trajectory projections of the singleparticle benchmark in the prescribed dipole geomagnetic field: x–y (top), x–z (middle), and y–z (bottom). and the relativistic particle pusher independently of the collective self-field effects considered in the subsequent bunch calculations. Having verified the single-particle transport, the simulations are next extended to the full relativistic electron bunch, for w… view at source ↗

**Figure 6.** Figure 6: , where the transverse extent has increased substantially relative to the initial state. To examine the effect of increasing particle number, the same calculation is repeated for Np = 800 and Np = 1600. A representative initial distribution for the Np = 1600 case is shown in view at source ↗

**Figure 8.** Figure 8: Particle distribution of the relativistic view at source ↗

**Figure 9.** Figure 9: Time evolution of the transverse and longitudinal RMS bunch sizes for the cases with Np = 400, 800, and 1600. Here σxy denotes the transverse RMS size and σz denotes the longitudinal RMS size. slower growth stage. This behavior indicates that the initial inter-particle potential energy is pro13 view at source ↗

**Figure 10.** Figure 10: Mean electron energy and its time derivative for the cases with Np = 400, 800, and 1600. The threshold view at source ↗

**Figure 11.** Figure 11: Evolution of the adaptive time step for the cases with Np = 400, 800, and 1600 during the early fully coupled calculation. The time step remains within the prescribed bounds ∆tmin = 10−15 s and ∆tmax = 10−12 s, shown by the dashdotted lines. for Np = 400, 800, and 1600, which suggests that the behavior is physical rather than numerical. At the same time, the mean particle energy rises rapidly, indicating… view at source ↗

**Figure 13.** Figure 13: Mean electron energy and its time derivative for the case with Np = 1600 and E0 = 10 MeV. The upper panel shows the evolution of the mean energy, while the lower panel shows the corresponding energy growth rate. The dashed vertical line marks view at source ↗

**Figure 14.** Figure 14: Mean electron energy and its time derivative for the case with Np = 1600 and E0 = 100 MeV. The upper panel shows the evolution of the mean energy, while the lower panel shows the corresponding energy growth rate. The dashed vertical line marks view at source ↗

**Figure 15.** Figure 15: Long-range evolution of the relativistic electron bunch for the case view at source ↗

**Figure 16.** Figure 16: Long-range evolution of the relativistic electron bunch for the case view at source ↗

**Figure 17.** Figure 17: Long-range evolution of the relativistic electron bunch for the case view at source ↗

**Figure 18.** Figure 18: Comparison of the final particle distributions for the E0 = 10 MeV case with and without inter-particle interaction. The red markers denote the calculation without inter-particle interaction, whereas the blue markers denote the calculation with inter-particle interaction. The three panels show the final distributions in the z–y, z–(x − x0), and (x − x0)–y planes, respectively. For the long-range transp… view at source ↗

read the original abstract

Particle-mesh methods, such as the particle-in-cell (PIC) method, cannot retain exact pairwise interaction at sub-cell scales. For dense nonneutral relativistic electron bunches, this makes it difficult to accurately capture the inter-particle electromagnetic interaction and the associated bunch divergence. In this work, the previously developed electromagnetic particle-particle (EM-PP) model for relativistic two-particle interaction is extended to many-particle electron bunch transport in the Earth's magnetosphere. The method combines the Li\'enard--Wiechert fields, an improved retarded-time evaluation procedure, and a relativistic particle pusher, and adopts a two-stage strategy to couple the dense early self-field-dominated evolution to the later long-range geomagnetic-field-controlled transport. The method provides a practical mesh-free approach for accurately simulating long-range transport of relativistic electron bunches when short-range electromagnetic interaction is important.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper extends a two-particle EM-PP model to relativistic bunches via Liénard-Wiechert fields and a two-stage dense-to-long-range handoff, but supplies no validation data or scaling tests to back its practicality claims.

read the letter

The core contribution is taking an earlier electromagnetic particle-particle approach that worked for two relativistic particles and stretching it to many-particle electron bunches in the magnetosphere. They keep exact pairwise Liénard-Wiechert interactions for the early dense phase, add an improved retarded-time solver, and then switch to geomagnetic-field transport once the bunch spreads out. That two-stage coupling is the part that looks new relative to standard PIC work cited in the abstract. It directly targets the problem that mesh methods lose sub-cell pairwise forces in nonneutral bunches, which can matter for divergence and radiation-belt modeling. The mesh-free premise is a clean way to preserve those interactions without artificial smoothing. The description is straightforward and avoids overclaiming the method as a full replacement for everything. The main weaknesses sit in the lack of evidence. The abstract and available text give no error norms, no convergence plots, no comparison against known two-particle solutions or PIC runs, and no timing or scaling numbers for N greater than a few dozen. Without those, it is impossible to judge whether the handoff between stages preserves trajectories or fields without jumps, or whether the all-pairs retarded evaluations stay feasible before cost blows up. The stress-test point about O(N^2) scaling and possible truncation errors in the transition is reasonable given what is shown; nothing in the paper description contradicts it. If the full manuscript contains solid benchmarks and controlled-error checks across the handoff, those would change the picture. This is aimed at space-plasma and beam-transport researchers who already care about exact short-range EM in dense relativistic flows and are willing to trade mesh simplicity for pairwise accuracy. A reader looking for a practical alternative to PIC in that niche could get value once the numerics are demonstrated. I would send it to peer review so referees can check the implementation details and any tests that exist in the complete manuscript. It is not ready as-is, but the underlying idea is worth a proper look.

Referee Report

3 major / 2 minor

Summary. The paper extends a prior electromagnetic particle-particle (EM-PP) model for relativistic two-particle interactions to many-particle electron bunches. It combines Liénard-Wiechert fields, an improved retarded-time evaluation procedure, and a relativistic particle pusher within a two-stage strategy that couples the early dense self-field-dominated phase to later long-range geomagnetic-field-controlled transport. The central claim is that this mesh-free approach accurately captures short-range pairwise electromagnetic interactions that affect bunch divergence during long-range transport in the Earth's magnetosphere.

Significance. If the accuracy and practicality claims hold, the method could fill a gap left by particle-mesh approaches such as PIC, which lose exact pairwise interactions at sub-cell scales; this would be relevant for modeling nonneutral relativistic bunches in space-plasma and accelerator contexts. The work builds on standard Liénard-Wiechert theory and a relativistic pusher, but currently offers no machine-checked proofs, reproducible code, or falsifiable predictions to strengthen the assessment.

major comments (3)

[Abstract] Abstract: the claim that the method 'provides a practical mesh-free approach for accurately simulating long-range transport' is unsupported; the text supplies no validation data, error metrics, convergence tests, or comparisons against known analytic solutions or other codes.
[Method description (two-stage strategy)] Two-stage coupling: no explicit error bounds, conservation checks, or continuity analysis are given for the handoff between the early all-pair self-field phase and the later geomagnetic phase; this is load-bearing for the claim that divergence accumulated early is preserved without discontinuities.
[Implementation and performance considerations] Scalability: the early-stage all-pair retarded-time evaluations remain O(N²); for realistic bunch sizes the text provides no cost analysis, truncation-error bounds, or controlled approximations, leaving open whether the exact-pairwise premise can be retained without prohibitive cost or uncontrolled numerical error.

minor comments (2)

[Numerical implementation] Clarify the precise algorithmic steps of the 'improved retarded-time evaluation procedure' with pseudocode or an equation reference to allow reproducibility.
[Particle pusher] Add a brief discussion of how the relativistic pusher is integrated with the Liénard-Wiechert fields to ensure consistency with the Lorentz force.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We are grateful to the referee for the thorough review and valuable suggestions. We address each of the major comments in detail below and will make the necessary revisions to strengthen the manuscript.

read point-by-point responses

Referee: [Abstract] Abstract: the claim that the method 'provides a practical mesh-free approach for accurately simulating long-range transport' is unsupported; the text supplies no validation data, error metrics, convergence tests, or comparisons against known analytic solutions or other codes.

Authors: We acknowledge that additional validation is required to support the claim in the abstract. Although the method is grounded in exact Liénard-Wiechert theory and the relativistic pusher, we will revise the manuscript to include explicit validation data, error metrics from comparisons to analytic two-particle solutions, convergence tests with varying particle numbers and timesteps, and limited comparisons to other simulation approaches for benchmark cases. This will better substantiate the accuracy for long-range transport. revision: yes
Referee: [Method description (two-stage strategy)] Two-stage coupling: no explicit error bounds, conservation checks, or continuity analysis are given for the handoff between the early all-pair self-field phase and the later geomagnetic phase; this is load-bearing for the claim that divergence accumulated early is preserved without discontinuities.

Authors: We agree that the handoff between stages requires more rigorous justification. In the revised version, we will add a section detailing the continuity of particle positions, velocities, and fields at the transition time, along with numerical conservation checks (e.g., total energy) and error bounds derived from the retarded-time evaluation accuracy and the assumption that self-fields become subdominant. revision: yes
Referee: [Implementation and performance considerations] Scalability: the early-stage all-pair retarded-time evaluations remain O(N²); for realistic bunch sizes the text provides no cost analysis, truncation-error bounds, or controlled approximations, leaving open whether the exact-pairwise premise can be retained without prohibitive cost or uncontrolled numerical error.

Authors: The all-pairs O(N²) complexity is a known characteristic of the exact EM-PP method and restricts it to moderate particle numbers. For the electron bunch sizes relevant to the magnetospheric transport problem studied here, the approach is computationally feasible. We will add a performance analysis subsection with scaling benchmarks, estimates of truncation errors in the improved retarded-time procedure, and discussion of the regime where the exact pairwise treatment is essential versus where approximations might be applied. revision: yes

Circularity Check

0 steps flagged

No circularity detected; derivation uses standard Liénard-Wiechert fields and relativistic pusher with explicit two-stage coupling.

full rationale

The paper's central method extends a prior two-particle EM-PP model by combining established Liénard-Wiechert potentials, an improved retarded-time solver, and a standard relativistic particle pusher. The two-stage handoff from self-field-dominated early evolution to geomagnetic transport is described as a practical coupling strategy without any fitted parameters renamed as predictions, self-definitional loops, or load-bearing uniqueness theorems imported solely via self-citation. All load-bearing components rest on externally verifiable electromagnetic theory and numerical integration techniques that do not reduce to the paper's own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The method relies on established relativistic electrodynamics without introducing new free parameters or invented entities in the abstract description.

axioms (2)

standard math Liénard-Wiechert fields give the exact electromagnetic fields of arbitrarily moving point charges
Invoked as the core field calculation for the particle-particle interactions.
domain assumption A relativistic particle pusher can accurately integrate trajectories under time-retarded electromagnetic fields
Assumed for the many-particle extension and two-stage transport phase.

pith-pipeline@v0.9.0 · 5455 in / 1238 out tokens · 58984 ms · 2026-05-08T04:32:14.883587+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

39 extracted references · 5 canonical work pages

[1]

Frontiers in Astronomy and Space Sciences , volume=

A mission concept to determine the magnetospheric causes of aurora , author=. Frontiers in Astronomy and Space Sciences , volume=. 2020 , publisher=

2020
[2]

Frontiers in Astronomy and Space Sciences , volume=

Relativistic particle beams as a resource to solve outstanding problems in space physics , author=. Frontiers in Astronomy and Space Sciences , volume=. 2019 , publisher=

2019
[3]

Journal of Geophysical Research: Space Physics , volume=

Modification of the loss cone for energetic particles in the Earth's inner magnetosphere , author=. Journal of Geophysical Research: Space Physics , volume=. 2022 , publisher=

2022
[4]

1980 , publisher=

The application of artificial electron beams to magnetospheric research , journal=. 1980 , publisher=

1980
[5]

Physics of Plasmas , volume=

Three-dimensional particle-particle simulations: Dependence of relaxation time on plasma parameter , author=. Physics of Plasmas , volume=. 2018 , publisher=

2018
[6]

Physical Review E , volume=

Electron-proton relaxation in hot-dense plasmas with a screened quantum statistical potential , author=. Physical Review E , volume=. 2024 , publisher=

2024
[7]

Journal of Propulsion and Power , volume =

Yinjian Zhao and Joseph Wang and Hideyuki Usui , title =. Journal of Propulsion and Power , volume =. 2018 , doi =

2018
[8]

New Jersey , volume=

Introduction to Electrodynamics Prentice Hall Upper Saddle River , author=. New Jersey , volume=
[9]

Journal of Experimental and Theoretical Physics , volume=

Relativistic modeling of ultra-short electron pulse propagation , author=. Journal of Experimental and Theoretical Physics , volume=. 2019 , publisher=

2019
[10]

Nature communications , volume=

Sub-phonon-period compression of electron pulses for atomic diffraction , author=. Nature communications , volume=. 2015 , publisher=

2015
[11]

Journal of Aerosol Science , volume=

Electrospray plume evolution: Influence of drag , author=. Journal of Aerosol Science , volume=. 2023 , publisher=

2023
[12]

The 36th International Electric Propulsion Conference, University of Vienna, Austria, IEPC-2019 , volume=

A particle-particle simulation model for droplet acceleration in colloid thrusters , author=. The 36th International Electric Propulsion Conference, University of Vienna, Austria, IEPC-2019 , volume=

2019
[13]

Physics of Plasmas , volume=

A binary collision Monte Carlo model for electron-ion temperature relaxation , author=. Physics of Plasmas , volume=. 2018 , publisher=

2018
[14]

Physics Letters A , volume=

Investigation of effective impact parameters in electron--ion temperature relaxation via particle--particle coulombic molecular dynamics , author=. Physics Letters A , volume=. 2017 , publisher=

2017
[15]

Molecular-Dynamics Simulations of Electron-Ion Temperature Relaxation in a Classical Coulomb Plasma , author =. Phys. Rev. Lett. , volume =. 2008 , month =. doi:10.1103/PhysRevLett.101.135001 , url =

work page doi:10.1103/physrevlett.101.135001 2008
[16]

IEEE Transactions on Plasma Science , year=

Comparison Between Electrostatic PP and PIC Simulations of Electron Bunch Expansion , author=. IEEE Transactions on Plasma Science , year=
[17]

R. W. Hockney and J. W. Eastwood , title =
[18]

Birdsall and A

Charles K. Birdsall and A. Bruce Langdon , title =
[19]

Plasma Sources Science and Technology , volume=

Vlasov simulations of electric propulsion beam , author=. Plasma Sources Science and Technology , volume=. 2024 , publisher=

2024
[20]

Physics of Plasmas , volume=

Grid-based Vlasov simulation of collisionless plasma expansion , author=. Physics of Plasmas , volume=. 2021 , publisher=

2021
[21]

and Porazik, Peter and Greklek-Mckeon, Michael and Amin, Kailas and Shaw, David and Kaganovich, Igor D

Powis, Andrew T. and Porazik, Peter and Greklek-Mckeon, Michael and Amin, Kailas and Shaw, David and Kaganovich, Igor D. and Johnson, Jay and Sanchez, Ennio , TITLE=. Frontiers in Astronomy and Space Sciences , VOLUME=. 2019 , DOI=

2019
[22]

John David Jackson , title =
[23]

2025 , doi =

An electromagnetic particle-particle model on solving relativistic binary collision , journal =. 2025 , doi =

2025
[24]

A Particle-Particle Simulation Model for Droplet Acceleration in Colloid Thrusters , year=

Yinjian Zhao and Joseph Wang , booktitle=. A Particle-Particle Simulation Model for Droplet Acceleration in Colloid Thrusters , year=
[25]

2023 , issn =

Electrospray plume evolution: Influence of drag , journal =. 2023 , issn =. doi:https://doi.org/10.1016/j.jaerosci.2022.106079 , author =

work page doi:10.1016/j.jaerosci.2022.106079 2023
[26]

Journal of Propulsion and Power , volume =

Zhao, Yinjian and Wang, Joseph and Usui, Hideyuki , title =. Journal of Propulsion and Power , volume =. 2018 , doi =

2018
[27]

2017 , issn =

Investigation of effective impact parameters in electron–ion temperature relaxation via Particle–Particle Coulombic molecular dynamics , journal =. 2017 , issn =. doi:https://doi.org/10.1016/j.physleta.2017.05.066 , author =

work page doi:10.1016/j.physleta.2017.05.066 2017
[28]

Physics of Plasmas , volume =

Zhao, Yinjian , title =. Physics of Plasmas , volume =. 2018 , month =

2018
[29]

and Gremillet, L

Pérez, F. and Gremillet, L. and Decoster, A. and Drouin, M. and Lefebvre, E. , title =. Physics of Plasmas , volume =. 2012 , month =

2012
[30]

Lavell, M. J. and Kish, A. J. and Sexton, A. T. and Masti, R. L. and Mohammad, I. and Kim, M. J. and Srinivasan, A. and Jarvis, K. and Scullin, W. and Shaw, J. G. and Sefkow, A. B. , title =. Physics of Plasmas , volume =. 2024 , month =

2024
[31]

and Lehe, R

Zhao, Y. and Lehe, R. and Myers, A. and Thévenet, M. and Huebl, A. and Schroeder, C. B. and Vay, J.-L. , title =. Physics of Plasmas , volume =. 2022 , month =

2022
[32]

Derouillat, A

Smilei : A collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation , journal =. 2018 , issn =. doi:https://doi.org/10.1016/j.cpc.2017.09.024 , author =

work page doi:10.1016/j.cpc.2017.09.024 2018
[33]

and Andersen, Hans C

Swope, William C. and Andersen, Hans C. and Berens, Peter H. and Wilson, Kent R. , title =. The Journal of Chemical Physics , volume =. 1982 , month =

1982
[34]

Pushing the Frontier in the Design of Laser-Based Electron Accelerators with Groundbreaking Mesh-Refined Particle-In-Cell Simulations on Exascale-Class Supercomputers , year=

Fedeli, Luca and Huebl, Axel and Boillod-Cerneux, France and Clark, Thomas and Gott, Kevin and Hillairet, Conrad and Jaure, Stephan and Leblanc, Adrien and Lehe, Rémi and Myers, Andrew and Piechurski, Christelle and Sato, Mitsuhisa and Zaim, Neïl and Zhang, Weiqun and Vay, Jean-Luc and Vincenti, Henri , booktitle=. Pushing the Frontier in the Design of La...
[35]

Griffiths , title =

David J. Griffiths , title =
[36]

D. J. Griffiths ,. 1999 ,

1999
[37]

Higuera, A. V. and Cary, J. R. , title =. Physics of Plasmas , volume =. 2017 , month =

2017
[38]

2018, ApJS, 235, 21, doi: 10.3847/1538-4365/aab114

Ripperda, B. and Bacchini, F. and Teunissen, J. and Xia, C. and Porth, O. and Sironi, L. and Lapenta, G. and Keppens, R. , title =. doi:10.3847/1538-4365/aab114 , year =

work page doi:10.3847/1538-4365/aab114
[39]

Energy loss due to binary collisions in a relativistic plasma , author =. Phys. Rev. A , volume =. 1979 , month =

1979