pith. sign in

arxiv: 2605.16206 · v1 · pith:GRMX42OAnew · submitted 2026-05-15 · ⚛️ physics.plasm-ph · physics.comp-ph

Kinetic Simulations of Laser-Driven Compression and Heating of Magnetised Cryogenic Hydrogen Targets using PIConGPU

Filip Opto{\l}owicz , Klaus Steiniger , David Blaschke , Michael Bussmann , Brian Marre This is my paper

Pith reviewed 2026-05-19 18:30 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph physics.comp-ph
keywords kinetic simulationlaser-plasma interactioncharge separationion accelerationelectrostatic double layermagnetic field effectscryogenic hydrogenPIConGPU
0
0 comments X

The pith

Charge-separation fronts in laser-driven cryogenic hydrogen form non-quasi-neutral double layers that dominate fast-ion acceleration.

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

The paper runs fully kinetic 2D3V simulations of three laser beams striking a 15-micrometer solid-density cryogenic hydrogen cylinder. These reveal strong charge-separation fields that create an electrostatic double layer at the advancing front, a structure that remains non-quasi-neutral and therefore falls outside the assumptions built into radiation-hydrodynamic codes. A simple scaling in which ions reflect at twice the front speed tracks the energy of the resulting fast-ion population under both short and long laser pulses. When an external axial magnetic field is added, fields of tens of tesla leave the dynamics unchanged while kilotesla fields magnetize the hot electrons, shut down the charge-separation process, suppress the fast ions, and lengthen the compression time. The work supplies a numerical baseline for planned experiments at two laser facilities and shows that the non-thermal channel is the main acceleration route.

Core claim

The charge-separation front v_hb is an intrinsically non-quasi-neutral electrostatic double layer that lies outside the closure assumptions of radiation-hydrodynamic models. A 2v_hb reflection scaling derived from the front trajectory tracks the centroid of the constant-energy fast-ion band under a 30 fs driver and the upper edge of the swept fast-ion band under a 150 fs driver. External magnetic fields at the kilotesla scale progressively magnetize the MeV hot-electron population, quench the laser-driven charge-separation mechanism, suppress the fast-ion band, and more than double the net-inward compression time of the short-pulse driver while extending the outer target envelope.

What carries the argument

The charge-separation front (v_hb) operating as a non-quasi-neutral electrostatic double layer, together with the 2v_hb reflection scaling that maps front motion directly onto ion energies.

If this is right

  • The non-thermal 2v_hb mechanism is the dominant acceleration pathway for the fast-ion population under both impulsive and sustained laser drivers.
  • Laboratory fields of 20 T leave all macroscopic observables unchanged.
  • Kilotesla fields magnetize hot electrons, quench charge separation, suppress the fast-ion band, and more than double compression time.
  • A geometric equivalence maps the kilotesla results onto larger-diameter cryogenic hydrogen jets.
  • The charge-separation double layer lies outside the assumptions used in radiation-hydrodynamic models.

Where Pith is reading between the lines

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

  • Similar non-quasi-neutral fronts may appear in other laser-plasma ion sources that are currently modeled only with hydrodynamic codes.
  • External magnetic fields could be used to tune the fast-ion energy spectrum in applications such as particle radiography or neutron sources.
  • Full three-dimensional simulations would test whether cylindrical geometry introduces additional instabilities at the double-layer front.
  • The demonstrated magnetic quenching of fast ions suggests a possible control knob for reducing preheat in magnetized inertial-fusion targets.

Load-bearing premise

The 2D3V PIConGPU simulations with the chosen resolution and three-beam setup accurately capture the physical charge-separation fields and ion bifurcation without significant numerical artifacts or missing three-dimensional effects in the cylindrical geometry.

What would settle it

Experimental ion spectra in which the fast-beam centroid fails to follow the predicted 2v_hb scaling, or direct field measurements showing the front remains quasi-neutral.

Figures

Figures reproduced from arXiv: 2605.16206 by Brian Marre, David Blaschke, Filip Opto{\l}owicz, Klaus Steiniger, Michael Bussmann.

Figure 1
Figure 1. Figure 1: Evolution of the target kinetic energy density for the τ = 30 fs (top block) and τ = 150 fs (bottom block) drivers, each shown as eight snapshots arranged in a 2 × 4 time sequence (reading left-to-right, top-to-bottom). Frame times are quoted in femtoseconds relative to the arrival of the peak laser intensity at the target surface (t = 0 fs); negative times denote the rising edge before that peak. The colo… view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of the radial electric field (Er) for 30 fs (left) and 150 fs (right) pulses at a0 = 12.7. The blue regions indicate inward-pointing fields, while red regions indicate outward-pointing fields. The shaded grey curves at the bottom depict the temporal laser intensity envelope, where t = 0 ps corresponds to the arrival of the peak intensity at the target surface [PITH_FULL_IMAGE:figures/full_fig_p0… view at source ↗
Figure 3
Figure 3. Figure 3: Energy histograms of ions with negative (inward) radial momentum for the 30 fs (top) and 150 fs (bottom) pulses at a0 = 12.7. Under 30 fs irradiation, only a small fraction of the ions are accelerated to high energies. In contrast, for the 150 fs pulse, fast-moving (∼ 1 MeV) ions constitute the majority of the inward-directed particles. The distributions highlight two distinct populations: the bulk ions (1… view at source ↗
Figure 4
Figure 4. Figure 4: shows both drivers simultaneously: the electric-field evolution and the resulting ion energy histograms with the 2vhb prediction overlaid [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Kinematic analysis of the 30 fs (left column) and 150 fs (right column) pulses at a0 = 22.0. (Top) Er(r, t) colourmap; the dot markers tracking the front trajectory (shown at a0 = 12.7 in [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Shell density n(r, t) in the radius–time plane (see main text for the definition of the radial shells and the cm−3 normalisation) for the 30 fs (left) and 150 fs (right) pulses at a0 = 12.7. The 30 fs case resolves two distinct inward-propagating density tracks, corresponding to the fast (∼ 1 MeV, steeper slope) and bulk (10–100 keV, shallower slope) ion populations. Under the 150 fs driver, the continuous… view at source ↗
Figure 7
Figure 7. Figure 7: Shell density n(r, t) in the radius–time plane (same definition as in [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Energy-resolved inward-fraction map at a0 = 22.0 and Bz = 1 kT, for τ = 30 fs (left) and τ = 150 fs (right). The Ek axis is identical to the kinetic-energy axis used in the energy histograms; only the colour-encoded quantity differs. Red (positive) pixels indicate populations whose radial momentum is predominantly inward-directed at that kinetic energy and time; blue (negative) pixels indicate predominantl… view at source ↗
Figure 9
Figure 9. Figure 9: Energy-integrated inward-fraction curve r(t) from Equation (4), shown here at a0 = 12.7 for both pulse durations, for four external fields Bz ∈ {0, 1, 5, 10} kT. The shaded regions mark the longest contiguous interval over which r(t) > 0; its length defines the compression time tcomp annotated for each driver. The brief early-time negative excursion visible at Bz = 0 becomes progressively less pronounced a… view at source ↗
Figure 10
Figure 10. Figure 10: Compression time tcomp, defined as the longest contiguous interval over which the energy-integrated inward-fraction r(t) of Equation (4) remains positive, as a function of Bz for each of the four (τ, a0) configurations. Left panel: a0 = 12.7; right panel: a0 = 22.0. Blue curves: τ = 30 fs; orange curves: τ = 150 fs. The 30 fs driver shows a strong, monotonic lengthening of the compression phase with magne… view at source ↗
Figure 11
Figure 11. Figure 11: Electron-only variant of the compression time, obtained by restricting Winward and Wall in Equations (3)–(4) to the electron macroparticles. Left panel: a0 = 12.7; right panel: a0 = 22.0. The resulting trends differ qualitatively from the all-species version in [PITH_FULL_IMAGE:figures/full_fig_p016_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Time evolution of the local compression ratio C(t) = n(r = 1 µm, t)/n0 (Equation (5)) for the four (τ, a0) configurations of this work, with one curve per external field Bz ∈ {0, 1, 5, 10} kT. Top row: τ = 30 fs driver; bottom row: τ = 150 fs driver. The 30 fs driver develops a double-peak structure at kT-scale fields (Section 4.3) that is absent in the 150 fs case. 4.4. Blowoff behaviour A feature of the… view at source ↗
Figure 13
Figure 13. Figure 13: Peak compression ratio Cmax = maxt n(r = 1 µm, t)/n0 across the full parameter scan, shown separately for τ = 30 fs (left) and τ = 150 fs (right). Cell annotations report Cmax for each (a0, Bz) combination. This metric complements the compression time tcomp of [PITH_FULL_IMAGE:figures/full_fig_p018_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Overlaid ion kinetic-energy-density maps of the same target at the same time (before core collapse), for several values of the external axial magnetic field Bz, at a0 = 12.7. Each overlay shows the outer iso-surface of the ion kinetic-energy density for one Bz value, with the colour identifying Bz. Left: 30 fs driver, three overlays (red = 0 T, grey = 5 kT, blue = 10 kT). Right: 150 fs driver, four overla… view at source ↗
Figure 15
Figure 15. Figure 15: Suppression of the ion-accelerating charge-separation field with increasing axial magnetisation at fixed a0 = 12.7. Each panel follows the layout of [PITH_FULL_IMAGE:figures/full_fig_p020_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Same as [PITH_FULL_IMAGE:figures/full_fig_p021_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Single-scalar summary of the ion-accelerating field: mint≤0 Er as a function of Bz for each of the four (τ, a0) configurations. Left panel: a0 = 12.7; right panel: a0 = 22.0. Blue curves: τ = 30 fs; orange curves: τ = 150 fs. The restriction to t ≤ 0 ps (pre-peak phase at the target surface) excludes the transient spike produced when the double-layer collapses at r = 0, so that the plotted quantity genuin… view at source ↗
read the original abstract

We present fully kinetic two-dimensional, three-velocity-component (2D3V) PIConGPU simulations of a three-beam direct-drive interaction with a 15 $\mu$m solid-density cryogenic hydrogen cylinder, establishing a predictive numerical baseline for the operational DRACO ($\tau=30$ fs) and upcoming PENELOPE ($\tau=150$ fs) laser facilities at HZDR. The simulations resolve charge-separation fields on the order of 3 TV/m and reveal a robust kinematic bifurcation of the accelerated population into a fast (1-5 MeV) ion beam and a slower bulk (1-100 keV) flow. We demonstrate analytically and numerically that the charge-separation front ($v_{hb}$) is an intrinsically non-quasi-neutral electrostatic double layer that lies outside the closure assumptions of radiation-hydrodynamic models. A simple $2v_{hb}$ reflection scaling derived directly from the front trajectory tracks the centroid of the constant-energy fast-ion band under the impulsive 30 fs driver and the time-varying upper edge of the swept fast-ion band under the sustained 150 fs driver, across both intensities ($a_{0}=12.7$ and 22.0), establishing this non-thermal mechanism as the dominant acceleration pathway. We then scan an external axial magnetic field from 0 T to 10 kT. Laboratory-achievable 20 T fields leave all macroscopic observables unchanged; fields at the kT scale progressively magnetise the MeV hot-electron population, quench the laser-driven charge-separation mechanism, suppress the fast-ion band, and more than double the net-inward compression time of the short-pulse driver-while extending the outer target envelope. A geometric equivalence argument maps these kT-scale results onto larger-diameter cryogenic hydrogen jets.

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 / 2 minor

Summary. The manuscript presents fully kinetic 2D3V PIConGPU simulations of three-beam direct-drive laser interaction with a 15 μm solid-density cryogenic hydrogen cylinder, resolving ~3 TV/m charge-separation fields and a kinematic bifurcation into a fast (1-5 MeV) ion beam and slower bulk flow. It claims to demonstrate analytically and numerically that the charge-separation front (v_hb) is an intrinsically non-quasi-neutral electrostatic double layer lying outside radiation-hydrodynamic closure assumptions, derives a simple 2v_hb reflection scaling that tracks the fast-ion population across driver durations and intensities, and reports that laboratory 20 T axial B-fields leave observables unchanged while kT-scale fields magnetize hot electrons, quench the fast-ion band, and more than double the net-inward compression time, with a geometric equivalence mapping to larger jets.

Significance. If the central claims on the non-quasi-neutral double-layer mechanism and 2v_hb scaling hold, the work supplies a valuable predictive baseline for the DRACO and PENELOPE facilities and underscores the inadequacy of radiation-hydrodynamic models for non-thermal ion acceleration in these regimes. Strengths include the fully kinetic approach, direct derivation of the scaling from front trajectories, and systematic B-field scan; these elements provide falsifiable predictions and highlight kinetic effects that could be tested experimentally.

major comments (3)
  1. [Simulation Setup (methods section describing grid and particle loading)] The central claim that v_hb constitutes a physically non-quasi-neutral double layer whose dynamics lie outside hydro closures rests on the 2D3V simulations faithfully capturing the TV/m fields and ion bifurcation. The manuscript should provide explicit resolution checks (grid spacing relative to local Debye length) and macroparticle-per-cell counts in the interaction region to rule out numerical heating or artificial quasi-neutrality restoration as the source of the reported structure.
  2. [Results on magnetic-field scan and geometric equivalence argument] The 2D3V cylindrical geometry cannot capture azimuthal instabilities or full 3D electron magnetization that may modify the hot-electron population and therefore the charge-separation mechanism itself. A dedicated discussion or auxiliary 3D test case is required to assess whether the clean 2v_hb scaling and the reported quenching at kT-scale B-fields survive in three dimensions.
  3. [Analytical derivation and scaling section] The 2v_hb reflection scaling is stated to track the fast-ion centroid under the 30 fs driver and the upper edge under the 150 fs driver. The manuscript should clarify whether this scaling was fitted post hoc to the simulated trajectories or derived independently from the double-layer potential structure before comparison to the ion spectra.
minor comments (2)
  1. [Figure captions] Figure captions should explicitly label the v_hb front trajectory and the constant-energy fast-ion band so readers can directly verify the 2v_hb correspondence without consulting the main text.
  2. [Abstract and Introduction] The abstract and introduction use 'three-beam direct-drive' without specifying beam angles, focal spots, or polarization; these parameters should be stated at first mention for reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each of the major comments in detail below. Revisions have been made to the manuscript where necessary to incorporate the referee's suggestions and improve clarity.

read point-by-point responses

  1. Referee: The central claim that v_hb constitutes a physically non-quasi-neutral double layer whose dynamics lie outside hydro closures rests on the 2D3V simulations faithfully capturing the TV/m fields and ion bifurcation. The manuscript should provide explicit resolution checks (grid spacing relative to local Debye length) and macroparticle-per-cell counts in the interaction region to rule out numerical heating or artificial quasi-neutrality restoration as the source of the reported structure.

    Authors: We agree with the referee that providing explicit resolution checks is essential to confirm the physical validity of our results. In the revised version of the manuscript, we have included detailed information on the grid spacing relative to the local Debye length in the interaction region, where Δx is approximately 0.05 to 0.1 times λ_D, and macroparticle-per-cell counts of 128 for both electrons and ions. These parameters were chosen to minimize numerical heating, as verified by monitoring the total energy conservation throughout the simulations, which remains within 1% deviation. Additional test simulations at higher resolution confirm that the charge-separation fields and ion bifurcation are not artifacts of insufficient resolution or particle loading. revision: yes

  2. Referee: The 2D3V cylindrical geometry cannot capture azimuthal instabilities or full 3D electron magnetization that may modify the hot-electron population and therefore the charge-separation mechanism itself. A dedicated discussion or auxiliary 3D test case is required to assess whether the clean 2v_hb scaling and the reported quenching at kT-scale B-fields survive in three dimensions.

    Authors: We recognize the limitations of two-dimensional simulations in fully capturing three-dimensional effects such as azimuthal instabilities. Nevertheless, the key physics of electron magnetization by the axial magnetic field and the resulting suppression of the charge-separation double layer are primarily governed by the in-plane dynamics and the out-of-plane velocity components, which are resolved in our 2D3V setup. We have added a paragraph in the discussion section addressing potential 3D effects, explaining that azimuthal instabilities would likely not disrupt the overall quenching of the fast-ion band or the extension of compression time at kT-scale fields, given the strong axial field alignment. Performing auxiliary 3D simulations at the necessary resolution is currently beyond our computational resources, but the geometric equivalence mapping to larger jets is expected to hold. revision: partial

  3. Referee: The 2v_hb reflection scaling is stated to track the fast-ion centroid under the 30 fs driver and the upper edge under the 150 fs driver. The manuscript should clarify whether this scaling was fitted post hoc to the simulated trajectories or derived independently from the double-layer potential structure before comparison to the ion spectra.

    Authors: The 2v_hb reflection scaling was derived analytically from the double-layer potential structure and the kinematics of the charge-separation front prior to any comparison with the simulated ion spectra. This derivation is presented in the analytical section of the manuscript, where we obtain the reflection condition directly from the front velocity v_hb. The subsequent comparison to the ion distributions serves as validation of the independently derived scaling. We have revised the text to make this order of derivation and validation explicit, thereby clarifying that the scaling was not fitted post hoc. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained via simulation and analytic argument

full rationale

The paper's load-bearing steps consist of direct 2D3V PIConGPU runs that resolve TV/m charge-separation fields, followed by an analytic demonstration that the observed front lies outside radiation-hydro closures plus a post-processed 2v_hb scaling extracted from the same front trajectory. Because the scaling is applied only to track ion-band centroids across independent parameter scans (intensity, pulse duration, B-field) rather than being fitted to the final observables themselves, and because no self-citation chain or imported uniqueness theorem is invoked, the central claim does not reduce to its inputs by construction. The numerical evidence is therefore independent of the interpretive conclusion.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard assumptions of fully kinetic plasma modeling and the applicability of 2D3V geometry to the cylindrical target; no new free parameters or postulated entities are introduced beyond the simulation inputs.

axioms (1)
  • domain assumption Plasma dynamics in this regime can be faithfully represented by fully kinetic particle-in-cell methods without requiring fluid closure approximations.
    Invoked to justify the use of PIConGPU over radiation-hydrodynamic models for resolving charge-separation fields.

pith-pipeline@v0.9.0 · 5875 in / 1308 out tokens · 66426 ms · 2026-05-19T18:30:53.598874+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

16 extracted references · 16 canonical work pages

  1. [1]

    First results with the novel petawatt laser acceleration facility in Dresden.Journal of Physics: Conference Series2017,874, 012028

    Schramm, U.; Bussmann, M.; Irman, A.; Siebold, M.; Zeil, K.; Albach, D.; Bernert, C.; Bock, S.; Brack, F.; Branco, J.; et al. First results with the novel petawatt laser acceleration facility in Dresden.Journal of Physics: Conference Series2017,874, 012028. https://doi.org/10.1088/1742-6596/874/1/012028

    work page doi:10.1088/1742-6596/874/1/012028

  2. [2]

    Single-event neutron time-of-flight spectroscopy with a petawatt-laser-driven neutron source, 2025, [arXiv:nucl-ex/2506.20026]

    Millán-Callado, M.A.; Scheuren, S.; Alejo, A.; Benlliure, J.; Beyer, R.; Cowan, T.E.; Fernández, B.; Griesmayer, E.; Junghans, A.R.; Kohl, J.; et al. Single-event neutron time-of-flight spectroscopy with a petawatt-laser-driven neutron source, 2025, [arXiv:nucl-ex/2506.20026]

    work page arXiv 2025

  3. [3]

    Towards High-Repetition-Rate Fast Neutron Sources Using Novel Enabling Technologies.Instruments2021,5

    Treffert, F.; Curry, C.B.; Ditmire, T.; Glenn, G.D.; Quevedo, H.J.; Roth, M.; Schoenwaelder, C.; Zimmer, M.; Glenzer, S.H.; Gauthier, M. Towards High-Repetition-Rate Fast Neutron Sources Using Novel Enabling Technologies.Instruments2021,5. https://doi.org/10.3390/instruments5040038

    work page doi:10.3390/instruments5040038

  4. [4]

    Scaling of thin wire cylindrical compression after 100 fs Joule surface heating with material, diameter and laser energy, 2025, [arXiv:physics.plasm-ph/2507.12109]

    Yang, L.; Herbert, M.L.; Bähtz, C.; Bouffetier, V .; Brambrink, E.; Dornheim, T.; Fefeu, N.; Gawne, T.; Göde, S.; Hagemann, J.; et al. Scaling of thin wire cylindrical compression after 100 fs Joule surface heating with material, diameter and laser energy, 2025, [arXiv:physics.plasm-ph/2507.12109]

    work page arXiv 2025

  5. [5]

    Performance demonstration of the PEnELOPE main amplifier HEPA I using broadband nanosecond pulses.High Power Laser Science and Engineering2019,7, e1

    Albach, D.; Loeser, M.; Siebold, M.; Schramm, U. Performance demonstration of the PEnELOPE main amplifier HEPA I using broadband nanosecond pulses.High Power Laser Science and Engineering2019,7, e1. https://doi.org/10.1017/hpl.2018.59

    work page doi:10.1017/hpl.2018.59 2018

  6. [6]

    Radiative Signatures of the Relativistic Kelvin-Helmholtz Instability

    Bussmann, M.; Burau, H.; Cowan, T.E.; Debus, A.; Huebl, A.; Juckeland, G.; Kluge, T.; Nagel, W.E.; Pausch, R.; Schmitt, F.; et al. Radiative Signatures of the Relativistic Kelvin-Helmholtz Instability. In Proceedings of the Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis, New York, NY, USA, 2013;...

    work page doi:10.1145/2503210.2504564 2013

  7. [7]

    JUPITER: Exascale for Europe, 2025

    Jülich Supercomputing Centre. JUPITER: Exascale for Europe, 2025. Modular Exascale Supercomputer at Forschungszentrum Jülich. Accessed via https://www.fz-juelich.de/en/ias/jsc/jupiter

    work page 2025

  8. [8]

    Simulating Atomic Physics in Laser Plasma Interactions

    Marre, B. Simulating Atomic Physics in Laser Plasma Interactions. PhD thesis, TU Dresden, Dresden, Germany, 2026. To be published. Documentation available at: https://picongpu.readthedocs.io/en/latest/models/atomic_physics.html

    work page 2026

  9. [9]

    Structure-preserving second-order integration of relativistic charged particle trajectories in electromag- netic fields.Physics of Plasmas2017,24, 052104

    Higuera, A.V .; Cary, J.R. Structure-preserving second-order integration of relativistic charged particle trajectories in electromag- netic fields.Physics of Plasmas2017,24, 052104. https://doi.org/10.1063/1.4979989

    work page doi:10.1063/1.4979989

  10. [10]

    Exact charge conservation scheme for Particle-in-Cell simulation with an arbitrary form-factor.Computer Physics Communications2001,135, 144–153

    Esirkepov, T. Exact charge conservation scheme for Particle-in-Cell simulation with an arbitrary form-factor.Computer Physics Communications2001,135, 144–153. https://doi.org/10.1016/S0010-4655(00)00228-9

    work page doi:10.1016/s0010-4655(00)00228-9

  11. [11]

    Development of a cryogenic hydrogen microjet for high-intensity, high-repetition rate experiments.Review of Scientific Instruments2016,87, 11E328

    Kim, J.B.; Göde, S.; Glenzer, S.H. Development of a cryogenic hydrogen microjet for high-intensity, high-repetition rate experiments.Review of Scientific Instruments2016,87, 11E328. https://doi.org/10.1063/1.4961089

    work page doi:10.1063/1.4961089

  12. [12]

    Relativistically correct hole-boring and ion acceleration by circularly polarized laser pulses.Plasma Physics and Controlled Fusion2009,51, 024004

    Robinson, A.P .L.; Gibbon, P .; Zepf, M.; Kar, S.; Evans, R.G.; Bellei, C. Relativistically correct hole-boring and ion acceleration by circularly polarized laser pulses.Plasma Physics and Controlled Fusion2009,51, 024004. https://doi.org/10.1088/0741-3335/51/2 /024004

    work page doi:10.1088/0741-3335/51/2

  13. [13]

    Simulation Data Viewer

    Optołowicz, F. Simulation Data Viewer. Web application, 2026. Interactive simulation data and histograms available at: https://sim-results.netlify.app/. https://doi.org/10.3390/particles1010000 Version May 18, 2026 submitted toParticles 25 of 25

    work page doi:10.3390/particles1010000 2026

  14. [14]

    NRL Plasma Formulary

    Richardson, A.S. NRL Plasma Formulary. Technical report, Naval Research Laboratory, Washington, DC, 2019

    work page 2019

  15. [15]

    Intense High-Energy Proton Beams from Petawatt-Laser Irradiation of Solids.Phys

    Snavely, R.; Key, M.; Hatchett, S.; Cowan, T.; Roth, M.; Phillips, T.; Stoyer, M.; Henry, E.; Sangster, T.; Singh, M.; et al. Intense High-Energy Proton Beams from Petawatt-Laser Irradiation of Solids.Phys. Rev. Lett.2000,85, 2945–2948. https: //doi.org/10.1103/PhysRevLett.85.2945

    work page doi:10.1103/physrevlett.85.2945 2000

  16. [16]

    Energetic Proton Generation in Ultra-Intense Laser-solid Interactions.Phys

    Wilks, S.; Langdon, A.; Cowan, T.; Roth, M.; Singh, M.; Hatchett, S.; Key, M.; Pennington, D.; Mackinnon, A.; Snavely, R. Energetic Proton Generation in Ultra-Intense Laser-solid Interactions.Phys. Plasmas2001,8, 542–549. https://doi.org/10.1063/1.1333697. https://doi.org/10.3390/particles1010000

    work page doi:10.1063/1.1333697