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

This paper presents fully kinetic PIConGPU simulations demonstrating that laser-driven compression of magnetized cryogenic hydrogen targets produces a dominant non-thermal ion acceleration mechanism via charge-separation double layers, which remains robust under laboratory-scale magnetic fields but is significantly suppressed and altered by kilotesla-scale fields that magnetize hot electrons and extend compression times.

The Gist

The Big Picture: Squeezing a Tiny Ice Cube with Light

Imagine you have a tiny cylinder of frozen hydrogen (like a microscopic ice cube) and you want to crush it in the center to create extreme pressure. To do this, scientists are using powerful lasers. This paper is a computer simulation that acts like a "digital wind tunnel," predicting exactly what happens when these lasers hit the ice.

The researchers are testing two different types of laser pulses:

1. The "Snap" (30 femtoseconds): A super-fast, sharp burst of energy, like a hammer hitting a nail.
2. The "Push" (150 femtoseconds): A longer, sustained shove, like a hand slowly but firmly pressing down on a spring.

They also tested what happens if they add a giant magnetic field to the mix, acting like an invisible cage around the ice.

The Main Discovery: Two Types of Particles

When the lasers hit the hydrogen, they don't just heat it up; they create a strange "traffic jam" of particles. The simulation revealed that the hydrogen splits into two distinct groups, much like a crowd of people reacting to a sudden event:

1. The "Sprinters" (Fast Ions): A small group of particles gets kicked hard and zooms inward at incredible speeds (millions of electron-volts).
2. The "Walkers" (Bulk Ions): The rest of the particles move inward much more slowly, like a crowd shuffling forward.

The "Magic Mirror" Analogy:
The paper explains that the "Sprinters" aren't being pushed by the laser directly. Instead, the laser creates a moving wall of electric charge (a "charge-separation front") that acts like a moving mirror.

• When the laser hits the ice, it pushes electrons away, leaving a gap.
• This gap creates a massive electric field (about 3 trillion volts per meter!).
• As this electric "mirror" moves inward, it bounces the positive hydrogen ions off it.
• Just like a tennis ball bouncing off a racket moving toward you, the ions gain speed. The paper found a simple rule: if the mirror moves at speed $v$, the ball bounces back at speed $2v$.

The Difference Between the "Snap" and the "Push"

The type of laser pulse changes how these "Sprinters" behave:

• The "Snap" (30 fs): Because the laser is so short, the electric mirror moves at a constant speed for a split second. This creates a neat, uniform group of Sprinters, all hitting the center with the exact same speed. It's like a perfectly timed volley of arrows.
• The "Push" (150 fs): Because the laser lasts longer, the electric mirror keeps accelerating as it moves. This means the Sprinters are launched at different speeds over time. Some are slow, some are fast. It's like a stream of water where the speed varies, creating a "sweep" of energies rather than a single sharp group.

The Magnetic Field Experiment: The Invisible Cage

The researchers then turned on a magnetic field to see if it would trap the particles and squeeze the ice harder. They tested fields ranging from what we can build in a lab (20 Tesla) to extreme, theoretical fields (10,000 Tesla).

• The Lab-Scale Field (20 T): This is like a gentle breeze. The particles are moving so fast and are so energetic that they simply ignore the magnetic field. They zoom right through it. The simulation showed zero change in the results.
• The Extreme Field (1,000–10,000 T): This is like a steel cage. At these levels, the magnetic field is strong enough to trap the fast-moving electrons.
  • The Result: When the electrons are trapped, they can't run away to form that "moving mirror" anymore. Without the mirror, the "Sprinters" (the fast ions) disappear. The laser loses its ability to kick the ions inward.
  • The Twist: Even though the "Sprinters" are gone, the magnetic field actually helps the "Walkers" (the bulk ions) stay compressed for twice as long. It's as if the magnetic cage holds the pressure in longer, allowing the slow-moving crowd to squeeze the center more effectively before they bounce back out.

A Surprising Side Effect: The Balloon Effect

You might think a magnetic cage would squeeze everything tighter. However, the simulation showed something counter-intuitive: the outer edge of the hydrogen target actually expanded more when the magnetic field was strong.

The Analogy: Imagine a balloon. If you squeeze the middle, the ends might bulge out. The magnetic field traps the hot electrons, but it also changes how they push against the outer layers of the target. Instead of collapsing neatly, the outer "skin" of the target puffs out further into space.

The "Geometric Trick"

The paper notes a clever way to test this in the real world. The 10,000 Tesla fields used in the simulation are impossible to build for a tiny 15-micron target. However, the physics depends on the ratio of the particle's path to the size of the target.

The authors argue that if you used a much larger target (like a jet of hydrogen 1,000 times bigger), you wouldn't need 10,000 Tesla. You could use a modest 10 Tesla field (which is easy to build) and get the exact same magnetic trapping effect. It's like how a small toy car and a real car can both turn the same way if you adjust the steering wheel speed relative to their size.

Summary

• Lasers create a moving electric wall that bounces ions inward.
• Short lasers create a uniform group of fast ions; long lasers create a mixed group.
• Weak magnets do nothing.
• Super-strong magnets stop the fast ions but help the slow ions stay compressed longer.
• Strong magnets also make the outer edge of the target puff out, rather than shrink.
• Big targets can experience these "super-magnet" effects using normal, lab-sized magnets.

Technical Summary

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

Problem Statement
The paper addresses the fundamental mechanisms governing laser-plasma coupling, hot-electron transport, and ion acceleration in the high-energy-density regime, specifically focusing on interactions with solid-density cryogenic hydrogen targets. Standard radiation-hydrodynamic (rad-hydro) models, which treat plasma as a quasi-neutral fluid, fail to capture intrinsically non-quasi-neutral electrostatic structures and non-thermal particle channels that arise on sub-picosecond timescales. The authors aim to establish a predictive numerical baseline for upcoming experimental campaigns at the Helmholtz-Zentrum Dresden-Rossendorf (HZDR), contrasting the operational DRACO laser ($\tau = 30$ fs) with the upcoming PEnELOPE facility ($\tau = 150$ fs). Furthermore, the study investigates the influence of external axial magnetic fields on these dynamics, seeking to understand thresholds for magnetization and their impact on compression and ion acceleration.

Methodology
The authors employed fully kinetic, two-dimensional, three-velocity-component (2D3V) Particle-in-Cell (PIC) simulations using the open-source code PIConGPU.

• Target: A solid-density cryogenic hydrogen cylinder ($R = 7.5$ µm, $n_0 = 5.24 \times 10^{28}$ m$^{-3}$).
• Laser Configuration: A symmetric three-beam direct-drive geometry with linear polarization ($\lambda = 800$ nm). Two pulse durations were modeled: $\tau = 30$ fs (impulsive) and $\tau = 150$ fs (sustained), at normalized vector potentials of $a_0 = 12.7$ and $a_0 = 22.0$.
• Physics Modules: The simulations utilized the FLYonPIC module for self-consistent non-LTE atomic kinetics and field/collisional ionization. The fully kinetic approach explicitly resolved Maxwell's equations and the Vlasov-Maxwell system, avoiding fluid heat-flux closures.
• Magnetic Field Scan: An external static axial magnetic field ($B_z$) was applied, scanning from 0 T to 10 kT. While 20 T represents a lab-achievable seed field, the kT-scale fields were used as a proxy for magnetized physics on larger targets via a geometric equivalence argument (holding the ratio of electron Larmor radius to target radius constant).
• Diagnostics: Key metrics included spatiotemporal electric fields ($E_r$), ion energy spectra, inward-momentum ratios to define compression time ($t_{comp}$), and local compression ratios at $r=1$ µm.

Key Contributions and Results

1. Kinematic Bifurcation and the $2v_{hb}$ Reflection Law:
   The simulations reveal a robust bifurcation of the ion population into a fast beam (1–5 MeV) and a slower bulk flow (1–100 keV). The authors demonstrate that the fast-ion population originates from direct reflection off a moving, non-quasi-neutral electrostatic double layer (the charge-separation front).

   • A simple reflection scaling, $E_k \approx 2v_{hb}$ (where $v_{hb}$ is the instantaneous velocity of the charge-separation front), quantitatively predicts the fast-ion energy.
   • For the 30 fs driver, the front moves at a nearly constant velocity, producing a monoenergetic fast-ion band.
   • For the 150 fs driver, the front accelerates continuously during the pulse rise, sweeping the fast-ion energy from hundreds of keV to several MeV.
   • This mechanism is identified as a genuinely kinetic feature, lying outside the closure assumptions of rad-hydro models which cannot resolve the coherent double-layer structure.

2. Magnetic Field Thresholds and Effects:

   • 20 T Regime: Fields achievable in current laboratories (20 T) have negligible impact on macroscopic observables because the Larmor radii of MeV hot electrons and protons remain orders of magnitude larger than the target radius.
   • kT-Scale Regime: Fields in the kilotesla range progressively magnetize the MeV hot-electron population ($r_L < R$). This leads to:
     • Quenching of Ion Acceleration: The laser-driven charge-separation mechanism is suppressed as hot electrons become magnetized and locked to field lines, preventing the formation of the strong inward-pointing electrostatic double layer. Consequently, the MeV fast-ion band disappears from the spectra.
     • Extended Compression: The net-inward compression time ($t_{comp}$) for the 30 fs driver more than doubles under strong magnetization (e.g., from ~0.52 ps to ~1.1 ps at $a_0=12.7$). The sustained inward push allows the bulk ion wave to reach the core, creating a second, larger density peak that is absent in the unmagnetized case.
     • Envelope Expansion: Contrary to naive confinement expectations, the outer target envelope expands further with increasing $B_z$. This is attributed to magnetically pinned hot electrons sustaining a late-time outward ambipolar field that expands the ion envelope even as the core compresses.

3. Geometric Equivalence:
   The authors argue that the kT-scale physics observed in the 15 µm target is kinematically equivalent to modest field strengths (~10 T) applied to much larger cryogenic hydrogen jets (e.g., 7.5 mm diameter), based on the conservation of the ratio $r_L/R$. This suggests the reported physics is accessible to larger-scale experiments, though exact plasma similarity (scaling density as $1/R^2$) is not maintained.

Significance and Claims
The paper claims two central findings:

1. Mechanistic Insight: The front-surface ion acceleration in this regime is driven by a non-quasi-neutral electrostatic double layer whose instantaneous velocity sets the fast-ion energies via a $2v_{hb}$ reflection law. This mechanism is distinct from Target Normal Sheath Acceleration (TNSA) and requires fully kinetic treatment to be captured.
2. Magnetization Impact: An applied axial magnetic field, once it brings the MeV hot-electron Larmor radius below the target radius, coherently quenches this acceleration channel. While this suppresses the fast-ion beam, it extends the net-inward compression phase and alters the target envelope dynamics.

The study emphasizes that these findings provide a predictive baseline for DRACO and PEnELOPE campaigns. The authors note limitations inherent to the 2D3V geometry (suppression of longitudinal return currents) and the lack of exact plasma similarity in their geometric mapping, cautioning that absolute field values and compression ratios should be interpreted as indicative rather than directly scalable to 3D experiments without further investigation. Future work is proposed to include 3D simulations, self-generated field studies via tilted incidence, and direct comparison with experimental observables.