Pith. sign in

REVIEW 4 major objections 5 minor 1 cited by

Parametric mapping of the efficiency$\unicode{x2013}$instability relation in plasma-wakefield accelerators

T0 review · 4 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read Simulations show the efficiency–instability relation in plasma-wakefield accelerators is a lower bound, not an exact equality, and reveal the field that reaches it.

desk verdict A solid, useful PIC parameter-mapping that reframes Lebedev's efficiency-instability relation as a lower bound and gives a practical optimal-field design rule, but the 'universal' phrasing overreaches the part of parameter space actually simulated. read the letter →

arxiv 2501.10602 v2 pith:C6HP5HDB submitted 2025-01-17 physics.acc-ph

classification physics.acc-ph
keywords plasma-wakefieldaccelerationbeam-breakupinstabilityefficiency–instabilityrelationoptimalbeamloadingblowoutregimeparticle-in-cellsimulationtransversemitigation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Plasma-wakefield accelerators can only power a collider if the accelerated bunch stays transversely stable while energy is efficiently transferred from the driver. This paper tests a previously proposed formula that ties power-transfer efficiency to the strength of the beam-breakup instability. Using a particle-in-cell parameter scan over all relevant variables, it finds that the formula is a lower bound: the instability is never weaker than predicted, and with the right accelerating field it comes close to that floor. The field that reaches the floor grows with wake radius, and choosing the wrong field can raise the instability by orders of magnitude.

What carries the argument

The central object is $\eta_t = -F_t/F_r$, the ratio of the transverse deflecting force to the ion focusing force at small offset, which measures the strength of the beam-breakup instability. The analytic benchmark is the predicted relation $\eta_{t,\mathrm{pred}} = \eta_p^2/[4(1-\eta_p)]$. The argument maps $\eta_t$ over three independent variables: the normalized wake radius $R_b k_p$, the normalized accelerating field $E_z/E_0$, and the bunch length, with the bunch current profile shaped so that the accelerating field is uniform. A single long simulation per grid point yields many bunch lengths by treating each longitudinal position as the tail of a shorter bunch, and the instability is evaluated at the tail particle.

What would settle it

Run the same parameter grid but deliberately include configurations the paper excludes: trailing bunches whose tails lie beyond half a plasma wavelength behind the driver, at very low $E_z/E_0$, with no energy spread and no ion motion. If any such configuration yields $\eta_t < \eta_p^2/[4(1-\eta_p)]$ at the tail particle, the claimed lower bound is not universal.

Watch

Extended reading notes

Core claim

Under the assumptions of a perfectly beam-loaded uniform accelerating field, no energy spread, and no ion motion, the paper establishes that the normalized transverse deflection force satisfies $\eta_t \ge \eta_p^2/[4(1-\eta_p)]$ rather than the equality proposed previously. The bound is approached only at one accelerating field for each wake radius, fitted as $E_{z,\mathrm{opt}}/E_0 \approx 0.23(1-0.78\eta_p^{1.86})(R_b k_p)^{1.5}$ for efficiencies up to about 60%. Deviations from this field, especially toward higher fields, can increase $\eta_t$ by orders of magnitude. The paper also shows that the amplitude-growth formula from the prior work remains a good description once $\eta_t$ is known, and that uncorrelated energy spread or ion motion can damp the growth.

Load-bearing premise

The lower-bound claim rests on restricting the simulations to perfectly flat, optimally loaded wakes and to bunch tails satisfying the cutoff in Eq. (7), so if an excluded configuration—say a very low accelerating field with a bunch spanning several plasma wavelengths—beats the inequality, the universal lower bound as stated would fail.

Editorial extensions

If this is right

  • Designers can treat $\eta_p^2/[4(1-\eta_p)]$ as the best achievable instability for a given efficiency, not an inevitable cost; any larger instability indicates non-optimal operating parameters.
  • For each wake radius, operating at the fitted optimal field keeps the instability within a small factor of the bound, and deviating to higher fields is considerably worse than deviating to lower fields.
  • The amplitude-growth formula from Ref. [18] remains usable for planning, since simulated oscillation growth tracks it closely once $\eta_t$ is known.
  • Mitigation by uncorrelated energy spread or ion motion works without changing the lower-bound mapping, but the required spread can exceed the sub-percent level allowed at the end of a collider stage.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the lower bound holds across the excluded low-field configurations, then the main practical cost of high efficiency shifts from an unavoidable instability to the precision with which the accelerating field can be set.
  • The fitted scaling $E_{z,\mathrm{opt}}/E_0 \propto (R_b k_p)^{1.5}$ is empirical; a closed-form derivation from the wakefield structure would test whether the exponent is exact.
  • The single-step static wake measurement used to define $\eta_t$ could be checked against full-length, self-consistent simulations, since wake evolution during acceleration could shift the lower-bound map.
  • Combining a large early-stage energy spread with ion motion may provide strong damping while keeping the final energy spread acceptable, because only the final-stage constraint matters for colliders.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper uses the quasi-static PIC code HiPACE++ with SALAME optimal beam loading to scan the transverse beam-breakup instability parameter ηt over a grid of normalized wake radii Rbkp, loaded accelerating fields Ez/E0, and bunch lengths, under the assumptions of no energy spread and no ion motion. From single-step static wake evaluations the authors extract the drive-to-trailing-bunch efficiency ηp and compare ηt to the analytic relation ηt,predicted = ηp^2/[4(1−ηp)] of Lebedev et al. They conclude that this relation is a lower bound, ηt ≥ ηp^2/[4(1−ηp)], that for each Rbkp there is an optimal accelerating field, fitted as Ez,opt/E0 ≈ 0.23(1−0.78ηp^1.86)(Rbkp)^1.5, at which the lower bound is approached, and that deviations from this field can increase ηt by orders of magnitude. Two long simulations at Rbkp = 1.59 and ηp = 0.4 are used to check that the amplitude-growth formula of Lebedev et al. reproduces the simulated oscillation growth, and additional simulations explore mitigation by uncorrelated energy spread and ion motion.

Significance. If the main claim holds, the paper materially sharpens the efficiency–instability relation in plasma-wakefield accelerators, changing a proposed equality into an inequality and providing a practical design rule (Eq. 9) for choosing the accelerating gradient that minimizes transverse instability. The strengths are the broad parametric scan (20 values of Rbkp and 40 values of Ez/E0), the use of SALAME to make beam loading exact, the longitudinal-slicing trick that extracts many bunch lengths from one simulation, and the fact that the central comparison uses the external Lebedev relation as a benchmark, so there is no definitional circularity in the lower-bound claim. The paper also makes falsifiable predictions (Eqs. 8–10). However, the universality of the inequality is qualified by the bunch-length cutoff in Eq. (7), the empirical fits carry no uncertainty quantification, and the amplitude-growth check covers only two moderate operating points; these caveats reduce the strength of the conclusions as currently stated.

major comments (4)
  1. [Section III, Eq. (7); Section VII, Eq. (10)] The paper's central claim, Eq. (10), is asserted as a universal lower bound under the stated assumptions, but it is only tested for configurations satisfying the cutoff ∆z < max(πRb/2, λp/2). For low Ez/E0, particularly at small Rbkp, the SALAME-loaded current profile needed to flatten the wake becomes several plasma wavelengths long; the paper explicitly describes these solutions as 'physically possible, but unusable in practice' and excludes them from the evaluation of ηt. Thus the observation that 'the lowest value never goes below 1 across all parameters' applies only to the surviving subset. If an excluded, perfectly loaded long-bunch configuration yields ηt/ηt,predicted < 1, Eq. (10) is false as stated. The authors should either extend the diagnostics to evaluate ηt over the full bunch length for the low-field corner, or qualify Eq. (10) to apply only to bunches satisfying Eq. (7). This is load-bearing because the advertised 'universal' lower bound is the central claim.
  2. [Section IV, Eqs. (8) and (9)] Equations (8) and (9) are empirical fits to the same simulation data that the paper later uses to support the predictive claims. No error bars, fit residuals, number of fit points, or cross-validation are reported. In particular, the statement that Eq. (9) predicts the optimal Ez/E0 'to an accuracy of 14% or better' measures the deviation on the fitting set, not on independent data. The paper should report the fit uncertainties and validate the formulas on a held-out subset of the simulations (or on additional simulations) before presenting Eq. (9) as a design rule.
  3. [Section V, Fig. 7] The amplitude-growth check of Eq. (3) is performed for only two operating points, both at Rbkp = 1.59 and ηp = 0.4, with ηt/ηt,predicted = 1.6 and 3.1. The extreme high-ratio points (ηt/ηt,predicted up to 10^3–10^5) that motivate the 'orders of magnitude' claim are never propagated over a full acceleration length. The static single-step measurement of ηt is the basis for the entire parameter map, and while the two moderate cases support the use of Eq. (3) in that regime, they do not establish that the exponential growth saturates or follows Eq. (3) at the extreme ratios. The authors should either run at least one long simulation at a strongly off-optimal point or explicitly restrict the 'orders of magnitude' language to the ηt scaling rather than the amplitude growth.
  4. [Sections II–V (overall simulation methodology)] No resolution convergence study, grid-convergence check, or data/code archive is reported for the PIC simulations that underpin the lower-bound claim. The main conclusion rests on a large number of simulations, but the reader cannot assess whether the numerical resolution (grid spacing, number of particles, time step, or box size) affects the measured ηt, particularly near the optimal valley where ηt/ηt,predicted is close to unity and in the high-ratio corner. Please report a convergence study for representative points in the parameter space and make the simulation input files and analysis scripts available or provide a clear data-availability statement.
minor comments (5)
  1. [Abstract and Ref. [18]] The abstract cites the Lebedev relation as 'Phys. Rev. Accel. Beams 21, 059901 (2018)', while the reference list entry [18] gives Phys. Rev. Accel. Beams 20, 121301 (2017); these bibliographic data are inconsistent and should be corrected.
  2. [Fig. 3] The color map indicating Ez/E0 is not accompanied by a visible colorbar in the printed figure; the reader cannot map the colors to numerical values without cross-referencing Fig. 5.
  3. [Footnote 1 (Section IV)] The statement 'We exclude the two highest Rbkp values to improve the fit' is given without justification; please specify the excluded values and quantify how their inclusion degrades the fit of Eq. (9).
  4. [Section VI.A, Fig. 8] The horizontal axis in Fig. 8 is labeled simply 't' without units; using µηt or kp z would be clearer and consistent with Eq. (3).
  5. [Section VI.B, Fig. 9] The text reports that with ion motion the measured ηt/ηt,predicted increases to 4.5 and 4.6 (from 1.6 and 3.1), yet Fig. 9 uses the no-ion-motion values of ηt in the comparison; this choice should be explained, since the ion-motion-modified ηt is not used in the plotted growth rates.

Circularity Check

1 steps flagged · score 4.0 of 10

Lower-bound claim is externally benchmarked and not circular; the main issue is Eq. (9), an in-sample fit presented as a prediction with accuracy quoted on the fitting set.

  1. fitted input called prediction [Section IV (Results), Eqs. (8)-(9) and Fig. 6 footnote 1; echoed in Section VII (Conclusions)]
    "This power-law scaling is confirmed by performing a fit of Ez/E0 versus Rbkp and ηp, giving Ez,opt/E0 ≈ 0.23(1 − 0.78ηp^1.86)(Rbkp)^1.5. ... The optimal values of Ez/E0 found from simulation do not deviate by more than 10% at low efficiency and 14% at high efficiency from the predicted values given by Eq. 9. 1 We exclude the two highest Rbkp values to improve the fit."

    Eq. (9) is a least-squares fit to the values of Ez/E0 at which each simulation in the grid minimizes ηt/ηt,predicted; the 'predicted values' are therefore evaluations of the fitted curve on the same dataset used to determine its coefficients. The stated 10%/14% agreement is an in-sample residual, not an out-of-sample prediction, and the two largest-Rbkp points that would worsen the agreement are explicitly excluded. Eq. (9) also inherits its power-law form from Eq. (8), another fit to the same scan, so the claimed optimal-field formula is an interpolation of the simulation minima rather than a derivation from Eq. (2) alone.

full rationale

The central result, Eq. (10), compares PIC-measured ηt with the externally published Lebedev formula Eq. (2), and no parameter of Eq. (2) is fitted in this paper; the inequality is therefore not true by construction, and the lower-bound claim has independent content within the simulated subset. The principal circularity burden is confined to the empirical optimal-field formula: Eq. (9) is fit to the same simulation minima it is said to 'predict,' its accuracy is quoted on the fitting set after removing the two highest-Rbkp points, and its scaling is inherited from Eq. (8), another fit to the same grid. This is a moderate fitted-input-called-prediction issue, but it does not force the main inequality. A scope caveat, not itself a circular step, is that the lower bound is established only for bunch lengths satisfying the Eq. (7) cutoff; the paper explicitly excludes low-field long-bunch solutions that are 'physically possible, but unusable in practice,' so configurations consistent with the stated assumptions but outside the cutoff are not tested. No load-bearing self-citation chain appears; the self-citations are contextual references rather than the argument's foundation.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central claim relies on the simulation setup, the optimal-loading assumption, the Eq. 7 cutoff, and the empirical fits listed above. No new physical entities are introduced.

free parameters (4)
  • 0.028 prefactor in Eq. 8 = 0.028
    Least-squares fit of eta_t/eta_t,pred to (Rb kp)^3 (Ez/E0)^-2 for Ez below the optimum.
  • Power-law exponents in Eq. 8 = 3 and -2
    Reported scaling of instability ratio with wake radius and accelerating field below the optimum.
  • Equation 9 coefficients = 0.23, 0.78, 1.86, and exponent 1.5
    Fit of the optimal Ez/E0 versus Rb kp and eta_p; the fit excludes the two highest Rb kp values.
  • Bunch-length cutoff in Eq. 7 = max(pi Rb/2, lambda_p/2)
    Hand-chosen limit on z-locations used to compute deflection force; excludes very low-field long-bunch configurations.
assumptions (5)
  • domain assumption Optimal beam loading produces a perfectly uniform accelerating field (SALAME algorithm, Sec. III)
    The entire scan is restricted to flat loaded fields, so the lower-bound claim is conditional on this loading scheme.
  • domain assumption No ion motion and cold plasma in the main parameter scan (Secs. I and III)
    The inequality is stated under this assumption; ion motion later changes the effective eta_t.
  • domain assumption The quasi-static PIC code HiPACE++ accurately models the wake and deflection forces (Sec. III)
    No convergence study or experimental validation is provided.
  • domain assumption Small transverse offset of 0.01/kp keeps the response linear and eta_t offset-independent (Sec. II)
    Underpins the definition and measurement of eta_t.
  • ad hoc to paper The z-location cutoff in Eq. 7 can be used to define the tail without excluding relevant physics
    Very low field and long bunch cases are excluded; the lower-bound universality is not tested there.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Parametric mapping of the efficiency$\unicode{x2013}$instability relation in plasma-wakefield accelerators." pith.science (2026). https://pith.science/paper/C6HP5HDB

@misc{pith2026250110602,
  author       = {Pith},
  title        = {Pith review of: Parametric mapping of the efficiency$\unicodex2013$instability relation in plasma-wakefield accelerators},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/C6HP5HDB}},
  note         = {Machine review of arXiv:2501.10602}
}
abstract

High efficiency is essential for plasma-wakefield accelerators to be a cost-effective alternative in high-power applications, such as a linear collider. However, in a plasma-wakefield accelerator the beam-breakup instability can be seeded by a transverse offset between the driver and trailing bunch. This instability, which rapidly increases the oscillation amplitude of the trailing bunch, grows with higher power-transfer efficiency from the driver to the trailing bunch [V. Lebedev et al., Phys. Rev. Accel. Beams 21, 059901 (2018)]. In this paper, we use particle-in-cell simulations to investigate the efficiency$\unicode{x2013}$instability relation that constrains the driver-to-trailing-bunch power-transfer efficiency in beam-driven plasma accelerators. We test the relation using a grid of simulations across all parameters that affect the beam-breakup instability, assuming a uniform accelerating field (optimal beam loading) and no ion motion. We find that the previously proposed efficiency$\unicode{x2013}$instability relation represents a lower limit on the strength of the instability for a given efficiency. For each normalized wake radius, only a certain accelerating field reaches this lowest value of the transverse instability; deviating from this point can increase the growth rate by several orders of magnitude. Lastly, we highlight how the oscillation-amplitude growth of the trailing bunch can be reduced or damped with an initial uncorrelated energy spread and the presence of ion motion.

Figures

Figures reproduced from arXiv: 2501.10602 by the authors.

Figure 1
Figure 1. FIG. 1. Plasma-electron density (blue) and beam-electron [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Grid schematic illustrating the parameter space (selected points only shown). In one dimension, we scan over the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. The strength of the instability divided by the [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: FIG. 4. The evolution of [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. The strength of the instability compared to the pre [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]
Figure 6
Figure 6. Figure 6: shows the value of Ez/E0 for which ηt/ηt,predicted is minimized at each Rbkp for various power-transfer efficiencies. The expected relation can be found using Eq. 8 by considering where ηt/ηt,predicted ∝ (Rbkp) 3/(Ez/E0) 2 intersects the value 1 (i.e., close to the min…
Figure 7
Figure 7. Figure 7: FIG. 7. Initial charge density (top) for the bunches and plasma in the blowout regime immediately after the plasma wake is [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Oscillation-amplitude growth of the trailing bunch with optimal (left) and less-optimal (right) operating points at [PITH_FULL_IMAGE:figures/full_fig_p008_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. Tail oscillation-amplitude growth comparison between the case with optimal (a) and non-optimal field value (b) at [PITH_FULL_IMAGE:figures/full_fig_p009_9.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. PWFA linear collider improvements -- from previous concepts to HALHF

    physics.acc-ph 2025-09 conditional novelty 1.0 of 10

    A summary of design changes from the 2013 PWFA-LC concept to the HALHF 2.0 baseline, containing no new results.

Reference graph

Works this paper leans on

29 extracted references · 25 canonical work pages · cited by 1 Pith paper

  1. [1]

    V. I. Veksler, Coherent principle of acceleration of charged particles, CERN Symposium on High-Energy Accelerators and Pion Physics (1956), p. 80

  2. [2]

    Tajima and J

    T. Tajima and J. M. Dawson, Laser electron accelerator, Phys. Rev. Lett. 43, 267 (1979)

  3. [3]

    P. Chen, J. M. Dawson, R. W. Huff, and T. Katsouleas, Acceleration of electrons by the interaction of a bunched electron beam with a plasma, Phys. Rev. Lett. 54, 693 (1985)

  4. [4]

    Aicheler, P

    M. Aicheler, P. Burrows, M. Draper, T. Garvey, P. Le- brun, K. Peach, N. Phinney, H. Schmickler, D. Schulte and N. Toge, A Multi-TeV linear collider based on CLIC technology, CLIC Conceptual Design Report (CERN, Geneva, 2013)

  5. [5]

    Behnke, J

    T. Behnke, J. E. Brau, B. Foster, J. Fuster, M. Harrison, J. M. Paterson, M. Peskin, M. Stanitzki, N. Walker and H. Yamamoto, The International Linear Collider Tech- nical Design Report - Volume 1: Executive Summary, (International Linear Collider, 2013)

  6. [6]

    Litos, E

    M. Litos, E. Adli, W. An, C. I. Clarke, C. E. Clayton, S. Corde, J. P. Delahaye et al., High-efficiency acceleration of an electron beam in a plasma wakefield accelerator, Nature 515, 92 (2014)

  7. [7]

    C. A. Lindstrøm, J. M. Garland, S. Schr¨ oder, L. Boulton, G. Boyle, J. Chappell, R. D’Arcy et al., Energy-spread preservation and high efficiency in a plasma-wakefield ac- celerator, Phys. Rev. Lett. 126, 014801 (2021)

  8. [8]

    C. A. Lindstrøm, J. Beinortait˙ e, J. Bj¨ orklund Svensson, L. Boulton, J. Chappell, S. Diederichs, B. Foster et al., Emittance preservation in a plasma-wakefield accelera- tor, Nat. Commun. 15, 6097 (2024)

Show all 29 references
  1. [9]

    Blumenfeld, C

    I. Blumenfeld, C. E. Clayton, F.-J. Decker, M. J. Hogan, C. Huang, R. Ischebeck, R. Iverson et al., Energy dou- bling of 42 GeV electrons in a metre-scale plasma wake- field accelerator, Nature 445, 741 (2007)

  2. [10]

    J. B. Rosenzweig, B. Breizman, T. Katsouleas, and J. J. Su, Acceleration and focusing of electrons in two- dimensional nonlinear plasma wake fields, Phys. Rev. A 44 R6189 (1991)

  3. [11]

    Foster, R

    B. Foster, R. D’Arcy, and C. A. Lindstrøm, A hy- brid, asymmetric, linear Higgs factory based on plasma- wakefield and radio-frequency acceleration, New J. Phys. 25, 093037 (2023)

  4. [12]

    D. H. Whittum, W. M. Sharp, S. S. Yu, M. Lampe, and G. Joyce, Electron-hose instability in the ion-focused regime, Phys. Rev. Lett. 67, 991 (1991)

  5. [13]

    T. J. Mehrling, R. A. Fonseca, A. Martinez de la Ossa, and J. Vieira, Mitigation of the hose instability in plasma-wakefield accelerators, Phys. Rev. Lett. 118, 174801 (2017)

  6. [14]

    Martinez de la Ossa, T

    A. Martinez de la Ossa, T. J. Mehrling, and J. Oster- hoff, Intrinsic stabilization of the drive beam in plasma- wakefield accelerators, Phys. Rev. Lett. 121, 064803 (2018)

  7. [15]

    Raimondi and A

    P. Raimondi and A. Seryi, Novel final focus design for future linear colliders, Phys. Rev. Lett. 86, 3779 (2001)

  8. [16]

    Mehrling, J

    T. Mehrling, J. Grebenyuk, F. S. Tsung, K. Floettmann and J. Osterhoff, Transverse emittance growth in staged laser-wakefield acceleration, Phys. Rev. ST Accel. Beams 15 (2012) 111303

  9. [17]

    W. K. H. Panofsky and M. Bander, Asymptotic Theory of Beam Break-Up in Linear Accelerators, Rev. Sci. In- strum. 39, 206 (1968)

  10. [18]

    Lebedev, A

    V. Lebedev, A. Burov, and S. Nagaitsev, Efficiency ver- sus instability in plasma accelerators, Phys. Rev. Accel. Beams 20, 121301 (2017)

  11. [19]

    Aksoy, D

    A. Aksoy, D. Schulte, and O. Yavas, Beam dynamics sim- ulation for the Compact Linear Collider drive-beam ac- celerator, Phys. Rev. ST Accel. Beams14, 084402 (2011)

  12. [20]

    Efficiency versus instability in plasma 11 accelerators

    S. S. Baturin, On “Efficiency versus instability in plasma 11 accelerators”, arXiv:2010.01046 (2020)

  13. [21]

    Diederichs, C

    S. Diederichs, C. Benedetti, A. Huebl, R. Lehe, A. Myers, A. Sinn, J.-L. Vay, W. Zhang, M. Th´ evenet, HiPACE++: A portable, 3D quasi-static particle-in-cell code, Comput. Phys. Commun. 278, 108421 (2022)

  14. [22]

    W. Lu, C. Huang, M. Zhou, W. B. Mori, and T. Kat- souleas, Nonlinear theory for relativistic plasma wake- fields in the blowout regime, Phys. Rev. Lett. 96, 165002 (2006)

  15. [23]

    Diederichs, C

    S. Diederichs, C. Benedetti, E. Esarey, J. Osterhoff, and C. B. Schroeder, High-quality positron acceleration in beam-driven plasma accelerators, Phys. Rev. Accel. Beams 23, 121301 (2020)

  16. [24]

    Tzoufras, W

    M. Tzoufras, W. Lu, F. S. Tsung, C. Huang, W. B. Mori, T. Katsouleas, J. Vieira, R. A. Fonseca, and L. O. Silva, Beam loading in the nonlinear regime of plasma-based acceleration, Phys. Rev. Lett. 101, 145002 (2008)

  17. [25]

    Barov, J

    N. Barov, J. B. Rosenzweig and M. C. Thompson, Energy loss of a high charge bunched electron beam in plasma, Phys. Rev. ST Accel. Beams 7, 061301 (2004)

  18. [26]

    W. An, M. Zhou, N. Vafaei-Najafabadi, K. A. Marsh, C. E. Clayton, C. Joshi, W. B. Moriet al.Strategies for miti- gating the ionization-induced beam head erosion problem in an electron-beam-driven plasma wakefield accelerator, Phys. Rev. ST Accel. Beams 16, 101301 (2013)

  19. [27]

    C. A. Lindstrøm, Self-correcting longitudinal phase space in a multistage plasma accelerator, arXiv:2104.14460 (2021)

  20. [28]

    J. B. Rosenzweig, A. M. Cook, A. Scott, M. C. Thomp- son, and R. B. Yoder, Effects of ion motion in intense beam-driven plasma wakefield accelerators, Phys. Rev. Lett. 95, 195002 (2005)

  21. [29]

    T. J. Mehrling, C. Benedetti, C. B. Schroeder, E. Esarey, and W. P. Leemans, Suppression of beam hosing in plasma accelerators with ion motion, Phys. Rev. Lett. 121, 264802 (2018)

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.