pith. machine review for the scientific record. sign in

arxiv: 2603.28473 · v1 · submitted 2026-03-30 · ⚛️ physics.plasm-ph

Recognition: no theorem link

Simulation Design for Velocity-Controlled Spatio-Temporal Drivers in Laser Wakefield Acceleration

Authors on Pith no claims yet

Pith reviewed 2026-05-14 01:36 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords spatio-temporal laser driversvelocity-controlled pulseslaser wakefield accelerationparticle-in-cell simulationssubluminal regimewakefield excitationcontinuous wall injectionspectral construction
0
0 comments X

The pith

Velocity-controlled spatio-temporal laser drivers decouple intensity peak velocity from group velocity for tailored wakefield excitation.

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

The paper develops a simulation workflow for modeling subluminal velocity-controlled spatio-temporal laser pulses in particle-in-cell codes for laser wakefield acceleration. It constructs these drivers as superpositions of exact vacuum solutions in a Maxwell-consistent spectral form that can be discretized in k-space for initialization. The work shows that the spatio-temporal geometry directly couples the longitudinal extent of the high-intensity region to the transverse scale, yielding scaling guidelines for near-resonant excitation. It also identifies geometric constraints that inflate computational cost over long distances and demonstrates that continuous wall injection reproduces the intended vacuum propagation while shrinking the transverse domain.

Core claim

Velocity-controlled ST laser drivers, constructed via a Maxwell-consistent spectral superposition of exact vacuum solutions and discretized for numerical use, excite wakes in the subluminal regime where the ST geometry links longitudinal high-intensity length to transverse beam size; scaling rules for near-resonant drive follow directly, and continuous wall injection preserves the designed propagation characteristics while cutting transverse domain size.

What carries the argument

Maxwell-consistent spectral construction of spatio-temporal pulses as superposition of exact vacuum solutions, discretized in k-space for PIC initialization and combined with continuous wall injection.

Load-bearing premise

The spectral construction stays accurate after k-space discretization and insertion into plasma without creating enough numerical dispersion or nonlinear shifts to change the intended velocity control.

What would settle it

A side-by-side comparison of the simulated intensity-peak velocity over many Rayleigh lengths against the analytically prescribed subluminal value, checking for deviation beyond discretization error.

Figures

Figures reproduced from arXiv: 2603.28473 by Chiara Badiali, Jorge Vieira, Rafael Almeida, Thales Silva.

Figure 1
Figure 1. Figure 1: Spectral design of ST pulses with controllable focal velocities. The desired spectral components (red curve) are defined by the intersection of the vacuum light cone (green surface) with the plane imposed by the spatio-temporal coupling (purple plane). (a) In the superluminal regime (vf = 1.5c), the intersection yields an open hyperbolic trajectory, characteristic of X-waves with infinite spectral extent. … view at source ↗
Figure 2
Figure 2. Figure 2: Two examples of spatio-temporal pulses obtained with the discrete implementation in OSIRIS. In (a), we have an example of the shape of a superluminal pulse with vf = 1.5c. In (b), we have an example of the shape of a subluminal pulse with vf = 0.5c. The blue-red colour-scale represents the electric field, and the black line corresponds to the pulse’s full longitudinal envelope f∥(z) in configuration space.… view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of the injected pulse electric field for three choices of the longitudinal and transverse repetition lengths used in the spectral reconstruction. (a) k0Trep,z = 700 and k0Trep,(x,y) = 1000, for which the reconstructed pulse is smooth and free of aliasing artefacts. (b) k0Trep,z = 200 and k0Trep,(x,y) = 100, where the longitudinal repetition length is too small, producing spurious replicas along … view at source ↗
Figure 4
Figure 4. Figure 4: Wakefield excitation for different choices of vf and w0. Panels (a) and (b) show the case vf = 0.8c: for (a) k0w0 = 10 the driver is too long and overlaps with most of the first plasma period, whereas for (b) k0w0 = 7.1 the driver length is better matched and a clearer wake structure is obtained. The same trend is more pronounced in panels (c) and (d) for vf = 0.6c, comparing (c) k0w0 = 10 with the tighter… view at source ↗
Figure 5
Figure 5. Figure 5: Waterfall plots of the accelerating electric-field amplitude for different focal velocities: (a) vf = 0.8c, (b) vf = 0.99c, and (c) vf = 1.4c. 7 [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Wakefield excitation for a fixed driver with vf = 0.96c and k0w0 = 17, for different values of a0. The upper panels show the electron density response and the on-axis laser electric-field envelope, while the lower panels show the on-axis accelerating field. (a) a0 = 0.5, (b) a0 = 1, (c) a0 = 2, and (d) a0 = 3. for this elongation by reducing the spot size w0 according to the scaling in equation (18). This … view at source ↗
Figure 7
Figure 7. Figure 7: Schematic of a subluminal ST pulse showing the full laser envelope (black) and the localised high-intensity region (grey line; blue–red shading). Because the envelope propagates at c while the focus advances at vf < c, the intensity peak slips backwards within the envelope between (a) t1 and (b) t2 > t1. Once the peak reaches the trailing edge of the envelope, the driver can no longer sustain an efficient … view at source ↗
Figure 8
Figure 8. Figure 8: Reduced-domain simulations of a subluminal ST driver with and without wall injection. (a) schematic of the injected field and boundary locations. (b,d) early-time snapshots in the co-moving coordinate ξ without and with wall injection, respectively. (c,e) late-time snapshots showing wake degradation without injection and sustained excitation with injection. (f) peak intensity in vacuum versus time, confirm… view at source ↗
read the original abstract

Velocity-controlled spatio-temporal (ST) laser drivers offer a route to tailoring laser-plasma interactions by allowing the velocity of the intensity peak to be controlled independently of the envelope group velocity. In this work, we present a simulation-design workflow for PIC modelling of subluminal velocity-controlled ST pulses in OSIRIS based on a Maxwell-consistent spectral construction expressed as a superposition of exact vacuum solutions, and we describe its discrete k-space representation for numerical initialisation. We then examine wakefield excitation with velocity-controlled drivers, showing how the ST geometry couples the effective longitudinal extent of the high-intensity region to the transverse scale and deriving scaling guidelines for near-resonant excitation in the subluminal regime. Finally, we discuss the geometric constraints that make long-distance simulations costly, including focus-envelope slippage and strong transverse expansion, and we show that continuous wall injection can reproduce the intended vacuum propagation while substantially reducing the transverse domain size. Together, these results provide practical guidelines for accurate and computationally efficient PIC simulations of velocity-controlled ST drivers in wakefield-relevant regimes.

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

2 major / 2 minor

Summary. The paper outlines a PIC simulation workflow in OSIRIS for subluminal velocity-controlled spatio-temporal laser pulses in laser wakefield acceleration. It constructs the drivers via superposition of exact vacuum Maxwell solutions, provides a discrete k-space initialization, derives scaling guidelines that couple the longitudinal high-intensity extent to the transverse scale for near-resonant wake excitation, and demonstrates continuous wall injection to reproduce vacuum propagation while shrinking the transverse domain. The central claim is that these elements together enable accurate, computationally efficient modeling of velocity-controlled drivers.

Significance. If the velocity control survives discretization and plasma coupling, the work supplies concrete, practical guidelines for a class of tailored drivers that could improve control over wakefield phase velocity and injection. The Maxwell-consistent spectral construction and the wall-injection technique are genuine strengths that reduce ad-hoc parameter tuning and domain size, respectively.

major comments (2)
  1. [initialization and wakefield-excitation sections] The manuscript asserts that the discrete k-space representation preserves the intended intensity-peak velocity after insertion into plasma, yet provides no quantitative vacuum-to-plasma comparison (e.g., measured peak-velocity drift versus propagation distance) that would bound numerical dispersion from finite k-space truncation or grid staggering. This comparison is load-bearing for the scaling guidelines in the subluminal regime.
  2. [scaling guidelines paragraph] The derived scaling guidelines for near-resonant excitation assume the ST geometry maintains its vacuum velocity throughout the interaction length; without reported error bars on velocity preservation or a test case showing that dispersion-induced drift remains << resonant detuning, the guidelines risk being applied outside their validity range.
minor comments (2)
  1. [spectral construction] Notation for the spectral components and the definition of the effective longitudinal extent should be made explicit in a single equation block to avoid ambiguity when readers implement the initialization.
  2. [wall-injection discussion] The abstract states that continuous wall injection 'reproduces the intended vacuum propagation' but does not specify the metric used (e.g., peak-position error or envelope fidelity); a short quantitative table would strengthen the claim.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed review. We have addressed both major comments by adding quantitative comparisons and validation tests in the revised manuscript. These additions directly support the velocity preservation claims and the applicability of the scaling guidelines.

read point-by-point responses
  1. Referee: [initialization and wakefield-excitation sections] The manuscript asserts that the discrete k-space representation preserves the intended intensity-peak velocity after insertion into plasma, yet provides no quantitative vacuum-to-plasma comparison (e.g., measured peak-velocity drift versus propagation distance) that would bound numerical dispersion from finite k-space truncation or grid staggering. This comparison is load-bearing for the scaling guidelines in the subluminal regime.

    Authors: We agree that a direct quantitative vacuum-to-plasma comparison strengthens the central claim. The spectral construction is exact in vacuum by design, but to bound discretization effects we have added new simulations and a dedicated figure comparing the intensity-peak velocity evolution in vacuum versus plasma over propagation distances matching the wakefield interaction length. The measured drift remains below 0.1 % and is negligible relative to the resonant detuning, thereby confirming that the discrete k-space initialization preserves the intended velocity to the accuracy required by the scaling guidelines. revision: yes

  2. Referee: [scaling guidelines paragraph] The derived scaling guidelines for near-resonant excitation assume the ST geometry maintains its vacuum velocity throughout the interaction length; without reported error bars on velocity preservation or a test case showing that dispersion-induced drift remains << resonant detuning, the guidelines risk being applied outside their validity range.

    Authors: We thank the referee for this important clarification. In the revised manuscript we have inserted error bars on all reported velocity measurements and added an explicit test-case subsection that quantifies dispersion-induced drift versus resonant detuning across the subluminal parameter space. The results show that the drift remains at least an order of magnitude smaller than the detuning for the geometries considered, thereby establishing the validity range of the scaling guidelines and removing the risk of misapplication. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained from vacuum solutions

full rationale

The paper constructs velocity-controlled ST drivers via superposition of exact vacuum Maxwell solutions expressed in a spectral representation, then discretizes this for PIC initialization in OSIRIS. Scaling guidelines for near-resonant wakefield excitation follow directly from the geometric coupling between longitudinal extent and transverse scale in this independent vacuum construction. Continuous wall injection is shown to reproduce the same vacuum propagation properties while shrinking the domain. No step reduces a claimed prediction or uniqueness result to a fitted parameter, self-citation chain, or ansatz imported from the authors' prior work; the central claims remain externally falsifiable against vacuum benchmarks and plasma simulations without circular reduction to inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the validity of the Maxwell-consistent spectral construction for vacuum solutions and the assumption that geometric coupling between longitudinal and transverse scales can be derived without additional fitted parameters beyond standard plasma scaling.

axioms (1)
  • domain assumption Maxwell-consistent spectral construction expressed as superposition of exact vacuum solutions can be discretized in k-space for numerical initialization
    Invoked as the basis for the PIC initial condition in the workflow description.

pith-pipeline@v0.9.0 · 5483 in / 1302 out tokens · 60536 ms · 2026-05-14T01:36:37.294798+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    High energy laser plasma accelerators,

    T. Tajima, “High energy laser plasma accelerators,”Laser and Particle Beams, vol. 3, no. 4, pp. 351–413, 1985

  2. [2]

    Physics of laser-driven plasma-based electron accelerators,

    E. Esarey, C. B. Schroeder, and W. P. Leemans, “Physics of laser-driven plasma-based electron accelerators,”Rev. Mod. Phys., vol. 81, pp. 1229–1285, Aug 2009

  3. [3]

    Efficient laser wakefield accelerator in pump depletion dominated bubble regime,

    V. Horn´ y, P. G. Bleotu, D. Ursescu, V. Malka, and P. Tomassini, “Efficient laser wakefield accelerator in pump depletion dominated bubble regime,”Phys. Rev. E, vol. 110, p. 035202, Sep 2024

  4. [4]

    Acceleration of electrons in the plasma wakefield of a proton bunch,

    E. Adli, A. Ahuja, O. Apsimon, R. Apsimon, A.-M. Bachmann, D. Barrientos, F. Batsch, J. Bauche, V. Berglyd Olsen, M. Bernardini,et al., “Acceleration of electrons in the plasma wakefield of a proton bunch,”Nature, vol. 561, no. 7723, pp. 363–367, 2018

  5. [5]

    Review of laser-driven ion sources and their applications,

    H. Daido, M. Nishiuchi, and A. S. Pirozhkov, “Review of laser-driven ion sources and their applications,”Reports on progress in physics, vol. 75, no. 5, p. 056401, 2012

  6. [6]

    Femtosecond x rays from laser-plasma accelerators,

    S. Corde, K. Ta Phuoc, G. Lambert, R. Fitour, V. Malka, A. Rousse, A. Beck, and E. Lefebvre, “Femtosecond x rays from laser-plasma accelerators,”Rev. Mod. Phys., vol. 85, pp. 1–48, Jan 2013

  7. [7]

    Prospects for free-electron lasers powered by plasma-wakefield-accelerated beams,

    M. Galletti, R. Assmann, M. Couprie, M. Ferrario, L. Giannessi, A. Irman, R. Pompili, and W. Wang, “Prospects for free-electron lasers powered by plasma-wakefield-accelerated beams,” Nature Photonics, vol. 18, no. 8, pp. 780–791, 2024

  8. [8]

    High-resolutionγ-ray radiography produced by a laser-plasma driven electron source,

    Y. Glinec, J. Faure, L. L. Dain, S. Darbon, T. Hosokai, J. J. Santos, E. Lefebvre, J. P. Rousseau, F. Burgy, B. Mercier, and V. Malka, “High-resolutionγ-ray radiography produced by a laser-plasma driven electron source,”Phys. Rev. Lett., vol. 94, p. 025003, Jan 2005

  9. [9]

    Wakefield acceleration,

    T. Tajima, X. Yan, and T. Ebisuzaki, “Wakefield acceleration,”Reviews of Modern Plasma Physics, vol. 4, no. 1, p. 7, 2020

  10. [10]

    Circumventing the dephasing and depletion limits of laser-wakefield acceleration,

    A. Debus, R. Pausch, A. Huebl, K. Steiniger, R. Widera, T. E. Cowan, U. Schramm, and M. Bussmann, “Circumventing the dephasing and depletion limits of laser-wakefield acceleration,”Phys. Rev. X, vol. 9, p. 031044, Sep 2019

  11. [11]

    Controlling the velocity of ultrashort light pulses in vacuum through spatio-temporal couplings,

    A. Sainte-Marie, O. Gobert, and F. Quere, “Controlling the velocity of ultrashort light pulses in vacuum through spatio-temporal couplings,”Optica, vol. 4, no. 10, pp. 1298–1304, 2017

  12. [12]

    Spatiotemporal control of laser intensity,

    D. H. Froula, D. Turnbull, A. S. Davies, T. J. Kessler, D. Haberberger, J. P. Palastro, S.-W. Bahk, I. A. Begishev, R. Boni, S. Bucht,et al., “Spatiotemporal control of laser intensity,” Nature photonics, vol. 12, no. 5, pp. 262–265, 2018

  13. [13]

    Flying focus: Spatial and temporal control of intensity for laser-based applications,

    D. Froula, J. Palastro, D. Turnbull, A. Davies, L. Nguyen, A. Howard, D. Ramsey, P. Franke, S.-W. Bahk, I. Begishev,et al., “Flying focus: Spatial and temporal control of intensity for laser-based applications,”Physics of Plasmas, vol. 26, no. 3, 2019

  14. [14]

    Space-time wave packets,

    M. Yessenov, L. A. Hall, K. L. Schepler, and A. F. Abouraddy, “Space-time wave packets,” Advances in Optics and Photonics, vol. 14, no. 3, pp. 455–570, 2022. 12 C. Badialiet al

  15. [15]

    Photon acceleration in a flying focus,

    A. Howard, D. Turnbull, A. Davies, P. Franke, D. Froula, and J. Palastro, “Photon acceleration in a flying focus,”Physical review letters, vol. 123, no. 12, p. 124801, 2019

  16. [16]

    Optical shock-enhanced self-photon acceleration,

    P. Franke, D. Ramsey, T. T. Simpson, D. Turnbull, D. Froula, and J. Palastro, “Optical shock-enhanced self-photon acceleration,”Physical Review A, vol. 104, no. 4, p. 043520, 2021

  17. [17]

    Phase-locked laser-wakefield electron acceleration,

    C. Caizergues, S. Smartsev, V. Malka, and C. Thaury, “Phase-locked laser-wakefield electron acceleration,”Nature Photonics, vol. 14, pp. 475–479, 2020

  18. [18]

    Ionization waves of arbitrary velocity driven by a flying focus,

    J. Palastro, D. Turnbull, S.-W. Bahk, R. Follett, J. Shaw, D. Haberberger, J. Bromage, and D. Froula, “Ionization waves of arbitrary velocity driven by a flying focus,”Physical Review A, vol. 97, no. 3, p. 033835, 2018

  19. [19]

    Ionization waves of arbitrary velocity,

    D. Turnbull, P. Franke, J. Katz, J. Palastro, I. Begishev, R. Boni, J. Bromage, A. Milder, J. Shaw, and D. Froula, “Ionization waves of arbitrary velocity,”Physical Review Letters, vol. 120, no. 22, p. 225001, 2018

  20. [20]

    Dephasingless laser wakefield acceleration,

    J. P. Palastro, J. L. Shaw, P. Franke, D. Ramsey, T. T. Simpson, and D. H. Froula, “Dephasingless laser wakefield acceleration,”Phys. Rev. Lett., vol. 124, p. 134802, Mar 2020

  21. [21]

    Programmable-trajectory ultrafast flying focus pulses,

    M. Ambat, J. Shaw, J. Pigeon, K. Miller, T. Simpson, D. Froula, and J. P. Palastro, “Programmable-trajectory ultrafast flying focus pulses,”Optics Express, vol. 31, no. 19, pp. 31354–31368, 2023

  22. [22]

    Dephasingless laser wakefield acceleration in the bubble regime,

    K. G. Miller, J. R. Pierce, M. V. Ambat, J. L. Shaw, K. Weichman, W. B. Mori, D. H. Froula, and J. P. Palastro, “Dephasingless laser wakefield acceleration in the bubble regime,” Scientific Reports, vol. 13, no. 1, p. 21306, 2023

  23. [23]

    Flying focus laser wake field acceleration by donut shape pulse,

    A. Ghasemi, S. Mirzanejhad, and T. Mohsenpour, “Flying focus laser wake field acceleration by donut shape pulse,”Applied Physics B, vol. 130, no. 6, p. 105, 2024

  24. [24]

    Laser wakefield acceleration of ions with a transverse flying focus,

    Z. Gong, S. Cao, J. P. Palastro, and M. R. Edwards, “Laser wakefield acceleration of ions with a transverse flying focus,”Phys. Rev. Lett., vol. 133, p. 265002, 2024

  25. [25]

    Laser wakefield acceleration driven by a discrete flying focus,

    J. R. Pierce, K. G. Miller, F. Li, J. P. Palastro, and W. B. Mori, “Laser wakefield acceleration driven by a discrete flying focus,”Physical Review Accelerators and Beams, vol. 28, no. 10, p. 101303, 2025

  26. [26]

    Electron acceleration in a wakefield generated by structured light,

    A. Liberman, S. Smartsev, A. Golovanov, S. Tata, I. A. Andriyash, S. Benracassa, E. Y. Levine, E. Kroupp, R. Lahaye, C. Thaury,et al., “Electron acceleration in a wakefield generated by structured light,” inLaser Acceleration of Electrons, Protons, and Ions VIII, vol. 13534, pp. 19–25, SPIE, 2025

  27. [27]

    Plasma accelerator for decaying particle,

    C. Badiali, R. Almeida, B. Malaca, R. Fonseca, T. Silva, and J. Vieira, “Plasma accelerator for decaying particle,”arXiv preprint arXiv:2506.23668, 2025

  28. [28]

    OSIRIS: A Three-Dimensional, Fully Relativistic Particle in Cell Code for Modeling Plasma Based Accelerators,

    R. A. Fonsecaet al., “OSIRIS: A Three-Dimensional, Fully Relativistic Particle in Cell Code for Modeling Plasma Based Accelerators,” inComputational Science — ICCS 2002(P. M. A. Sloot, A. G. Hoekstra, C. J. K. Tan, and J. J. Dongarra, eds.), (Berlin, Heidelberg), p. 342, Springer Berlin Heidelberg, 2002

  29. [29]

    Arbitrarily non-paraxial electromagnetic wave-packets in particle-in-cell codes,

    R. R. de Almeida, “Arbitrarily non-paraxial electromagnetic wave-packets in particle-in-cell codes,” master’s thesis, Instituto Superior T´ ecnico, Universidade de Lisboa, Lisbon, Portugal,

  30. [30]

    Maxwell-consistent non-paraxial laser injection enabling prescribed focal velocity (flying focus) in PIC

  31. [31]

    Electromagnetic gaussian beams beyond the paraxial approximation,

    C. Sheppard and S. Saghafi, “Electromagnetic gaussian beams beyond the paraxial approximation,”Journal of the Optical Society of America A, vol. 16, no. 6, pp. 1381–1386, 1999

  32. [32]

    Gaussian beam and pulsed-beam dynamics: complex-source and complex-spectrum formulations within and beyond paraxial asymptotics,

    E. Heyman and L. B. Felsen, “Gaussian beam and pulsed-beam dynamics: complex-source and complex-spectrum formulations within and beyond paraxial asymptotics,”Journal of the Optical Society of America A, vol. 18, no. 7, pp. 1588–1611, 2001

  33. [33]

    Exact solutions for the electromagnetic fields of a flying focus,

    D. Ramsey, A. Di Piazza, M. Formanek, P. Franke, D. H. Froula, B. Malaca, W. B. Mori, J. R. Pierce, T. T. Simpson, J. Vieira, M. Vranic, K. Weichman, and J. P. Palastro, “Exact solutions for the electromagnetic fields of a flying focus,”Phys. Rev. A, vol. 107, p. 013513, Jan 2023. 13 C. Badialiet al

  34. [34]

    Cas proceedings of the 2019 course on high gradient wakefield accelerators, sesimbra, portugal, 11–22 march 2019,

    V. Malka, “Cas proceedings of the 2019 course on high gradient wakefield accelerators, sesimbra, portugal, 11–22 march 2019,” tech. rep., 2020

  35. [35]

    Direct observation of a wakefield generated with structured light,

    A. Liberman, A. Golovanov, S. Smartsev, S. Tata, I. A. Andriyash, S. Benracassa, E. Y. Levine, Y. Wan, E. Kroupp, and V. Malka, “Direct observation of a wakefield generated with structured light,”Nature Communications, vol. 16, p. 10957, dec 2025

  36. [36]

    B. E. A. Saleh and M. C. Teich,Beam Optics, ch. 3, pp. 80–107. John Wiley & Sons, Ltd, 1991

  37. [37]

    Particle-in-cell modelling of laser–plasma interaction using fourier decomposition,

    A. F. Lifschitz, X. Davoine, E. Lefebvre, J. Faure, C. Rechatin, and V. Malka, “Particle-in-cell modelling of laser–plasma interaction using fourier decomposition,”Journal of Computational Physics, vol. 228, no. 5, pp. 1803–1814, 2009

  38. [38]

    Implementation of a hybrid particle code with a pic description in r–z and a gridless description inϕinto osiris,

    A. Davidson, A. Tableman, W. An, F. S. Tsung, W. Lu, J. Vieira, R. A. Fonseca, L. O. Silva, and W. B. Mori, “Implementation of a hybrid particle code with a pic description in r–z and a gridless description inϕinto osiris,”Journal of Computational Physics, vol. 281, pp. 1063–1077, 2015. 14