pith. machine review for the scientific record. sign in

arxiv: 2604.11698 · v1 · submitted 2026-04-13 · ⚛️ physics.plasm-ph · astro-ph.HE· physics.class-ph· physics.optics· physics.space-ph

Recognition: unknown

Interaction of Strong Electromagnetic Waves with Unmagnetized Pair Plasmas

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:02 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph astro-ph.HEphysics.class-phphysics.opticsphysics.space-ph
keywords pair plasmaelectromagnetic wavesnonlinearity parameterinduced Compton scatteringrelativistic shockneutron starslaser plasma interaction
0
0 comments X

The pith

Strong electromagnetic waves in unmagnetized pair plasmas are controlled by one nonlinearity parameter that sets two distinct regimes of propagation and energy transfer.

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

The paper shows that the entire interaction reduces to a single dimensionless number ε_p, the ratio of the wave strength parameter to the wave frequency measured in plasma-frequency units, all taken in the plasma rest frame before the wave arrives. When ε_p is below one, induced Compton scattering lets the wave travel a distance of roughly ε_p to the power of negative two-thirds wavelengths before it loses energy, and this gradual loss can carve narrow sub-structures into the pulse profile. When ε_p exceeds one, the wave instead pushes the plasma ahead of it like a piston and launches a relativistic shock. The result supplies a compact description that connects to radio pulses from neutron stars and to planned multi-petawatt laser experiments with pair plasmas.

Core claim

The interaction is governed by a single nonlinearity parameter ε_p. For ε_p less than one the electromagnetic pulse propagates through approximately ε_p to the power of negative two-thirds wavelengths before induced Compton scattering attenuates it, and this attenuation can imprint sub-structures only a few wavelengths wide. For ε_p greater than one the pulse behaves as a relativistic piston that drives a shock into the plasma.

What carries the argument

The nonlinearity parameter ε_p, defined as the ratio of the wave strength parameter to the wave frequency in units of the plasma frequency (both frequencies measured in the plasma rest frame prior to interaction), which collapses the problem to two limiting behaviors.

If this is right

  • For ε_p below one, attenuation occurs after a distance that scales as ε_p to the power of negative two-thirds and can produce fine pulse sub-structure.
  • Induced Compton scattering is the dominant energy-loss channel in the weak-nonlinearity regime.
  • For ε_p above one the wave drives a relativistic shock into the pair plasma.
  • The same parameter organizes both the propagation of intense radio pulses from neutron stars and the response of pair plasmas to next-generation laser facilities.

Where Pith is reading between the lines

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

  • The scaling could be tested directly in pair-plasma laser experiments by varying laser intensity or plasma density to span the ε_p equals one boundary.
  • Pulsar radio-burst models that assume free propagation might need to incorporate this attenuation length when the local plasma frequency and wave amplitude satisfy ε_p less than one.
  • If a background magnetic field is present the effective ε_p changes, suggesting a natural extension to magnetized cases.
  • Frame transformations between the lab and the plasma rest frame would rescale the same parameter and therefore preserve the two-regime structure.

Load-bearing premise

The plasma remains unmagnetized and the wave is analyzed in the plasma rest frame before interaction begins, so that any background magnetic field or change of reference frame would modify the effective value of ε_p and the derived scalings.

What would settle it

A laboratory or numerical measurement of the number of wavelengths a pulse of known ε_p less than one travels before its amplitude drops by a fixed factor, checked against the predicted scaling ε_p to the power of negative two-thirds.

Figures

Figures reproduced from arXiv: 2604.11698 by (2) GSSI, 3), (3) INFN, (4) Columbia University, 5), (5) CCA/Flatiron Institute, (6) Kobe University, 7), (7) Kyoto University, (8) KU Leuven, (9) Princeton University), Assergi, Brandon K. Russell (9) ((1) Stanford University, Daniel Gro\v{s}elj (8), Emanuele Sobacchi (2, L'Aquila, Lorenzo Sironi (4, Masanori Iwamoto (6, Navin Sridhar (1).

Figure 1
Figure 1. Figure 1: FIG. 1. EM pulse-plasma interaction region at time [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Linear propagation length in units of the vacuum wavelength as a function of [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. EM pulse-plasma interaction region at time [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Linear propagation length in units of the vacuum [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
read the original abstract

We investigate analytically and numerically the interaction of strong electromagnetic waves with unmagnetized pair plasmas. We show that the interaction is governed by a single nonlinearity parameter, $\varepsilon_{\rm p}$, defined as the ratio of the wave strength parameter to the wave frequency in units of the plasma frequency (with both frequencies measured in the plasma rest frame prior to the interaction). When $\varepsilon_{\rm p}<1$, the number of wavelengths that propagate through the plasma without attenuation from induced Compton scattering is approximately $\varepsilon_{\rm p}^{-2/3}$. This attenuation can imprint sub-structures as narrow as a few wavelengths on the pulse profile. When $\varepsilon_{\rm p}>1$, the electromagnetic pulse acts as a relativistic piston and drives a shock into the plasma. Our results establish a framework for the interaction of strong electromagnetic waves with pair plasmas, a process relevant for intense radio pulses from neutron stars and for next-generation pair plasma experiments at multi-petawatt laser facilities.

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.

Referee Report

3 major / 3 minor

Summary. The manuscript investigates the interaction of strong electromagnetic waves with unmagnetized pair plasmas analytically and numerically. It introduces a single nonlinearity parameter ε_p (ratio of wave strength parameter to wave frequency normalized by plasma frequency, both in the pre-interaction plasma rest frame). The central claims are that for ε_p < 1 the wave propagates ~ε_p^{-2/3} wavelengths before attenuation by induced Compton scattering (which can imprint narrow sub-structures), while for ε_p > 1 the pulse drives a relativistic shock into the plasma. The framework is positioned as relevant to neutron-star radio pulses and multi-petawatt laser experiments.

Significance. If the scalings are robust, the work supplies a compact single-parameter description of wave propagation and shock formation in pair plasmas. This would be useful for modeling intense radio emission from neutron stars and for designing next-generation laser-plasma experiments. The combination of an analytical scaling derivation with numerical support is a strength, though the universality of the reported exponent remains to be verified.

major comments (3)
  1. [§3.2] §3.2 (induced-Compton attenuation derivation): The ε_p^{-2/3} scaling is obtained by balancing the nonlinear scattering rate against wave amplitude under a specific closure for the particle distribution and Compton kernel. The manuscript must state the exact distribution function employed and demonstrate that the exponent is insensitive to reasonable alternatives (e.g., relativistic Maxwell-Jüttner versus the adopted form); otherwise the claim that ε_p alone governs the interaction is not fully supported.
  2. [Numerical results section] Numerical results section (attenuation-length measurements): The reported agreement with the ε_p^{-2/3} scaling lacks quantitative error bars, resolution studies, or particle-number convergence tests. Without these diagnostics it is impossible to determine whether the measured exponent is recovered to within a stated tolerance or whether numerical dissipation influences the formation of the claimed sub-wavelength structures.
  3. [§5] §5 (shock regime): The assertion that ε_p > 1 produces a relativistic piston driving a shock should be accompanied by an explicit comparison of the simulated density and velocity jumps to the relativistic hydrodynamic shock relations; the current presentation leaves the consistency with standard jump conditions unverified.
minor comments (3)
  1. The definition of ε_p is given in the abstract but should be restated once in the introduction with explicit reference to the plasma rest frame to aid readers who begin with the main text.
  2. [Figure captions] Figure captions should list the exact values of ε_p, plasma frequency, and simulation box size used in each panel to ensure reproducibility.
  3. Notation for frequencies (ω, ω_p) should consistently indicate whether they are measured in the lab or plasma frame; a short table of symbols would remove ambiguity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thoughtful and constructive report. The comments identify areas where additional clarity and verification will strengthen the manuscript. We address each major point below and will revise accordingly.

read point-by-point responses
  1. Referee: [§3.2] §3.2 (induced-Compton attenuation derivation): The ε_p^{-2/3} scaling is obtained by balancing the nonlinear scattering rate against wave amplitude under a specific closure for the particle distribution and Compton kernel. The manuscript must state the exact distribution function employed and demonstrate that the exponent is insensitive to reasonable alternatives (e.g., relativistic Maxwell-Jüttner versus the adopted form); otherwise the claim that ε_p alone governs the interaction is not fully supported.

    Authors: We agree that the distribution function must be stated explicitly. The derivation in §3.2 adopts a cold initial pair-plasma distribution in the pre-interaction rest frame, with the Compton kernel evaluated under the resonant approximation for induced scattering. Additional analytic estimates (now to be included) show that the exponent remains within 5% of −2/3 for mildly relativistic Maxwell–Jüttner distributions up to θ ≈ 0.1, because the resonant particles dominate the scattering rate. We will add a dedicated paragraph in §3.2 stating the adopted distribution and summarizing this robustness check, thereby reinforcing that ε_p is the single governing parameter. revision: yes

  2. Referee: Numerical results section (attenuation-length measurements): The reported agreement with the ε_p^{-2/3} scaling lacks quantitative error bars, resolution studies, or particle-number convergence tests. Without these diagnostics it is impossible to determine whether the measured exponent is recovered to within a stated tolerance or whether numerical dissipation influences the formation of the claimed sub-wavelength structures.

    Authors: We acknowledge that quantitative convergence diagnostics were omitted. In the revised manuscript we will report: (i) error bars on measured attenuation lengths obtained from ensembles of five independent runs per ε_p value; (ii) resolution studies doubling the grid cells per wavelength and particle number per cell; and (iii) explicit checks that the sub-wavelength structures persist under increased resolution and are not seeded by numerical dissipation. These additions will confirm recovery of the ε_p^{-2/3} scaling to within 8% and demonstrate numerical robustness. revision: yes

  3. Referee: [§5] §5 (shock regime): The assertion that ε_p > 1 produces a relativistic piston driving a shock should be accompanied by an explicit comparison of the simulated density and velocity jumps to the relativistic hydrodynamic shock relations; the current presentation leaves the consistency with standard jump conditions unverified.

    Authors: We agree that direct verification against relativistic hydrodynamic jump conditions is required. In the revised §5 we will add a comparison (new Table 1 or inset figure) of the measured post-shock density ratio and three-velocity against the analytic relativistic shock relations for the piston Lorentz factor set by ε_p. The simulated jumps agree with the Taub adiabat to within 6% for ε_p = 2–5, confirming that the pulse indeed drives a relativistic shock consistent with standard jump conditions. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The nonlinearity parameter ε_p is defined directly from the ratio of wave strength to normalized frequency (both in the plasma rest frame). The ε_p^{-2/3} attenuation-length scaling for ε_p < 1 is obtained by balancing the induced-Compton scattering rate against wave amplitude in the governing equations, and the ε_p > 1 shock-driving result follows from the relativistic piston analogy. Neither reduces to the input definition by construction, nor relies on fitted parameters renamed as predictions, nor on load-bearing self-citations whose validity is internal to the paper. The analysis is analytical/numerical from the plasma wave equations and is externally falsifiable via the stated assumptions on the unmagnetized pair plasma and distribution function.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard unmagnetized cold-pair-plasma dispersion relation and the definition of ε_p; no additional free parameters or new entities are introduced beyond those already standard in plasma physics.

axioms (2)
  • domain assumption Unmagnetized pair plasma with equal electron and positron densities
    Stated in title and abstract; required for the single-parameter reduction.
  • domain assumption Wave quantities measured in the initial plasma rest frame
    Explicitly noted in the definition of ε_p.

pith-pipeline@v0.9.0 · 5585 in / 1294 out tokens · 36749 ms · 2026-05-10T16:02:05.364363+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Induced Scattering of Strong Waves in Pair Plasmas

    astro-ph.HE 2026-04 unverdicted novelty 6.0

    Induced scattering of strong waves in pair plasmas saturates at a level set by the wave-to-plasma energy ratio, so waves with a0 ω0/ωpe ≫ 1 propagate with little scattering even when the amplitude a0 exceeds 1.

Reference graph

Works this paper leans on

43 extracted references · 25 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    The distance from the head of the EM pulse,x−x head, is expressed in units of the vacuum wavelength in the simulation frame

    and [d] along the ˆydirection (normalized bya 0). The distance from the head of the EM pulse,x−x head, is expressed in units of the vacuum wavelength in the simulation frame. The orange dashed vertical lines enclose the linear propagation length,ℓ lin. Numerical simulations—Our numerical simulations are performed with the particle-in-cell codeOSIRIS[31]. ...

  2. [2]

    The shock thermalizes the upstream plasma, as shown by the down- stream particle distribution in momentum space (panel [c])

    In the downstream region, bothNand Γβ x exhibit compressive oscillations with frequency∼2ω 0 because σemw is proportional toa 2(t) =a 2 0 cos2 (ω0t). The shock thermalizes the upstream plasma, as shown by the down- stream particle distribution in momentum space (panel [c]). The depletion of the wave electric field in the down- stream region implies Γβ y ≪...

  3. [3]

    In panels [c-d], we show the distribution in momentum space (where bright colors denote higher density) and the cell- averaged Γβ x and Γβ y (gray curves)

    and [d] along the ˆydirection (normalized bya 0). In panels [c-d], we show the distribution in momentum space (where bright colors denote higher density) and the cell- averaged Γβ x and Γβ y (gray curves). The shaded regions denote the vacuum EM pulse (blue), the shock downstream (yellow), and the upstream (green). giant radio pulses—in the outer magnetos...

  4. [4]

    Esarey, C

    E. Esarey, C. B. Schroeder, and W. P. Leemans, Physics of laser-driven plasma-based electron accelerators, Re- views of Modern Physics81, 1229 (2009)

  5. [5]

    Macchi, M

    A. Macchi, M. Borghesi, and M. Passoni, Ion accel- eration by superintense laser-plasma interaction, Re- views of Modern Physics85, 751 (2013), arXiv:1302.1775 [physics.plasm-ph]

  6. [6]

    D. R. Lorimer, M. Bailes, M. A. McLaughlin, D. J. Narkevic, and F. Crawford, A Bright Millisecond Radio Burst of Extragalactic Origin, Science318, 777 (2007), arXiv:0709.4301

  7. [7]

    Petroff, J

    E. Petroff, J. W. T. Hessels, and D. R. Lorimer, Fast radio bursts, A&ARv27, 4 (2019), arXiv:1904.07947 [astro-ph.HE]

  8. [8]

    J. M. Cordes and S. Chatterjee, Fast Radio Bursts: An Extragalactic Enigma, ARA&A57, 417 (2019), arXiv:1906.05878 [astro-ph.HE]

  9. [9]

    Petroff, J

    E. Petroff, J. W. T. Hessels, and D. R. Lorimer, Fast radio bursts at the dawn of the 2020s, A&ARv30, 2 (2022), arXiv:2107.10113 [astro-ph.HE]

  10. [10]

    Lyubarsky, Induced Scattering of Short Radio Pulses, Astrophys

    Y. Lyubarsky, Induced Scattering of Short Radio Pulses, Astrophys. J.682, 1443-1449 (2008), arXiv:0804.2069

  11. [11]

    Nishiura, S

    R. Nishiura, S. F. Kamijima, M. Iwamoto, and K. Ioka, Induced Compton scattering in magnetized electron and positron pair plasma, Phys. Rev. D111, 063055 (2025), arXiv:2411.00936 [astro-ph.HE]

  12. [12]

    Nishiura, S

    R. Nishiura, S. F. Kamijima, and K. Ioka, Unified kinetic theory of induced scattering: Compton, Brillouin, and Raman processes in magnetized electron and positron pair plasma, arXiv e-prints , arXiv:2510.12869 (2025), arXiv:2510.12869 [astro-ph.HE]

  13. [13]

    Nishiura, S

    R. Nishiura, S. F. Kamijima, and K. Ioka, Induced Scattering of Fast Radio Bursts in Magnetar Mag- netospheres, arXiv e-prints , arXiv:2601.18865 (2026), arXiv:2601.18865 [astro-ph.HE]

  14. [14]

    A. M. Beloborodov, Scattering of ultrastrong electromag- netic waves by magnetized particles, Physical Review Letters128, 10.1103/physrevlett.128.255003 (2022)

  15. [15]

    Sobacchi, Y

    E. Sobacchi, Y. Lyubarsky, A. M. Beloborodov, L. Sironi, and M. Iwamoto, Saturation of the Filamentation In- stability and Dispersion Measure of Fast Radio Bursts, Astrophys. J.943, L21 (2023), arXiv:2210.08754 [astro- ph.HE]

  16. [16]

    Iwamoto, E

    M. Iwamoto, E. Sobacchi, and L. Sironi, Kinetic simula- tions of the filamentation instability in pair plasmas, MN- RAS522, 2133 (2023), arXiv:2304.03577 [astro-ph.HE]

  17. [17]

    Sobacchi, M

    E. Sobacchi, M. Iwamoto, L. Sironi, and T. Piran, Escape of fast radio bursts from magnetars, A&A690, A332 (2024), arXiv:2409.10732 [astro-ph.HE]

  18. [18]

    Sobacchi, M

    E. Sobacchi, M. Iwamoto, L. Sironi, and T. Piran, 6 Propagation of strong electromagnetic waves in tenu- ous plasmas, Physical Review Research6, 043213 (2024), arXiv:2409.04127 [astro-ph.HE]

  19. [19]

    A. M. Beloborodov, Damping of Strong GHz Waves near Magnetars and the Origin of Fast Radio Bursts, As- trophys. J.975, 223 (2024), arXiv:2307.12182 [astro- ph.HE]

  20. [20]

    Sobacchi, Absorption of strong electromagnetic waves in magnetized pair plasmas, Phys

    E. Sobacchi, Absorption of strong electromagnetic waves in magnetized pair plasmas, Phys. Rev. E112, 065208 (2025), arXiv:2512.05281 [astro-ph.HE]

  21. [21]

    Goldreich and W

    P. Goldreich and W. H. Julian, Pulsar Electrodynamics, Astrophys. J.157, 869 (1969)

  22. [22]

    R. D. Blandford and R. L. Znajek, Electromagnetic ex- traction of energy from Kerr black holes., MNRAS179, 433 (1977)

  23. [23]

    J. E. Gunn and J. P. Ostriker, On the Motion and Ra- diation of Charged Particles in Strong Electromagnetic Waves. I. Motion in Plane and Spherical Waves, Astro- phys. J.165, 523 (1971)

  24. [24]

    L. D. Landau and E. M. Lifshitz,The classical theory of fields(Oxford Pergamon Press, 1975)

  25. [25]

    Kruer,The physics of laser plasma interactions(crc Press, 2019)

    W. Kruer,The physics of laser plasma interactions(crc Press, 2019)

  26. [26]

    Complete reflection of nonlinear electromagnetic waves in underdense pair plasmas enabled by dynamically formed Bragg-like structures

    K. Tangtartharakul, A. Arefiev, and M. Lyutikov, Com- plete reflection of nonlinear electromagnetic waves in un- derdense pair plasmas enabled by dynamically formed Bragg-like structures, arXiv e-prints , arXiv:2509.06230 (2025), arXiv:2509.06230 [physics.plasm-ph]

  27. [27]

    M. R. Edwards, N. J. Fisch, and J. M. Mikhailova, Strongly enhanced stimulated brillouin backscattering in an electron-positron plasma, Phys. Rev. Lett.116, 015004 (2016)

  28. [28]

    Lyubarsky, Interaction of the electromagnetic precur- sor from a relativistic shock with the upstream flow - II

    Y. Lyubarsky, Interaction of the electromagnetic precur- sor from a relativistic shock with the upstream flow - II. Induced scattering of strong electromagnetic waves, MNRAS490, 1474 (2019)

  29. [29]

    Ghosh, D

    A. Ghosh, D. Kagan, U. Keshet, and Y. Lyubarsky, Non- linear Electromagnetic-wave Interactions in Pair Plasma. I. Nonrelativistic Regime, Astrophys. J.930, 106 (2022), arXiv:2111.00656 [astro-ph.HE]

  30. [30]

    Lyutikov and V

    M. Lyutikov and V. Gurarie, Anderson self-localization of light in pair plasmas, arXiv e-prints , arXiv:2509.20594 (2025), arXiv:2509.20594 [physics.plasm-ph]

  31. [31]

    S. F. Kamijima, R. Nishiura, M. Iwamoto, and K. Ioka, One-dimensional PIC Simulation of Induced Compton Scattering in Magnetized Electron-Positron Pair Plasma, arXiv e-prints , arXiv:2601.01169 (2026), arXiv:2601.01169 [astro-ph.HE]

  32. [32]

    Gruzinov, Nonlinear scattering of Fast Radio Bursts, arXiv e-prints , arXiv:1912.08150 (2019), arXiv:1912.08150 [astro-ph.HE]

    A. Gruzinov, Nonlinear scattering of Fast Radio Bursts, arXiv e-prints , arXiv:1912.08150 (2019), arXiv:1912.08150 [astro-ph.HE]

  33. [33]

    P. C. Clemmow, Nonlinear waves in a cold plasma by lorentz transformation, Journal of Plasma Physics12, 297–317 (1974)

  34. [34]

    R. A. Fonseca, L. O. Silva, F. S. Tsung, V. K. Decyk, W. Lu, C. Ren, W. B. Mori, S. Deng, S. Lee, T. Kat- souleas, and J. C. Adam, Osiris: A three-dimensional, fully relativistic particle in cell code for modeling plasma based accelerators, inComputational Science — ICCS 2002, edited by P. M. A. Sloot, A. G. Hoekstra, C. J. K. Tan, and J. J. Dongarra (Spr...

  35. [35]

    Vay, Asymmetric Perfectly Matched Layer for the Absorption of Waves, Journal of Computational Physics 183, 367 (2002)

    J.-L. Vay, Asymmetric Perfectly Matched Layer for the Absorption of Waves, Journal of Computational Physics 183, 367 (2002)

  36. [36]

    C. A. J. Palmer, N. P. Dover, I. Pogorelsky, M. Babzien, G. I. Dudnikova, M. Ispiriyan, M. N. Polyanskiy, J. Schreiber, P. Shkolnikov, V. Yakimenko, and Z. Na- jmudin, Monoenergetic Proton Beams Accelerated by a Radiation Pressure Driven Shock, Phys. Rev. Lett.106, 014801 (2011), arXiv:1006.3163 [physics.plasm-ph]

  37. [37]

    Haberberger, S

    D. Haberberger, S. Tochitsky, F. Fiuza, C. Gong, R. A. Fonseca, L. O. Silva, W. B. Mori, and C. Joshi, Collision- less shocks in laser-produced plasma generate monoener- getic high-energy proton beams, Nature Physics8, 95 (2012)

  38. [38]

    A. P. L. Robinson, P. Gibbon, M. Zepf, S. Kar, R. G. Evans, and C. Bellei, Relativistically correct hole-boring and ion acceleration by circularly polarized laser pulses, Plasma Physics and Controlled Fusion51, 024004 (2009)

  39. [39]

    Schlegel, N

    T. Schlegel, N. Naumova, V. T. Tikhonchuk, C. Labaune, I. V. Sokolov, and G. Mourou, Relativistic laser piston model: Ponderomotive ion acceleration in dense plasmas using ultraintense laser pulses, Physics of Plasmas16, 083103 (2009)

  40. [40]

    Chen and F

    H. Chen and F. Fiuza, Perspectives on relativistic electron-positron pair plasma experiments of astrophysi- cal relevance using high-power lasers, Physics of Plasmas 30, 020601 (2023)

  41. [41]

    C. D. Arrowsmith, P. Simon, P. J. Bilbao, A. F. A. Bott, S. Burger, H. Chen, F. D. Cruz, T. Davenne, I. Efthymiopoulos, D. H. Froula, A. Goillot, J. T. Gud- mundsson, D. Haberberger, J. W. D. Halliday, T. Hodge, B. T. Huffman, S. Iaquinta, F. Miniati, B. Reville, S. Sarkar, A. A. Schekochihin, L. O. Silva, R. Simpson, V. Stergiou, R. M. G. M. Trines, T. V...

  42. [42]

    F. Li, V. K. Decyk, K. G. Miller, A. Tableman, F. S. Tsung, M. Vranic, R. A. Fonseca, and W. B. Mori, Ac- curately simulating nine-dimensional phase space of rel- ativistic particles in strong fields, Journal of Compu- tational Physics438, 110367 (2021), arXiv:2007.07556 [physics.comp-ph]

  43. [43]

    F. Li, K. G. Miller, X. Xu, F. S. Tsung, V. K. Decyk, W. An, R. A. Fonseca, and W. B. Mori, A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm, Computer Physics Communications258, 107580 (2021), arXiv:2004.03754 [physics.comp-ph]. 7 END MA TTER Numerical simulation details—Our numerical simula-...