Optimization of laser-driven proton acceleration in a near-critical-density plasma

Deji Liu; Dongchi Cai; Guanqi Qiu; Jinqing Yu; Qianyi Ma; Xueqing Yan; Yinren Shou; Zheng Gong

arxiv: 2512.19341 · v5 · submitted 2025-12-22 · ⚛️ physics.plasm-ph · physics.acc-ph

Optimization of laser-driven proton acceleration in a near-critical-density plasma

Guanqi Qiu , Qianyi Ma , Deji Liu , Dongchi Cai , Zheng Gong , Yinren Shou , Jinqing Yu , Xueqing Yan This is my paper

Pith reviewed 2026-05-16 20:40 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph physics.acc-ph

keywords laser-driven proton accelerationnear-critical-density plasmatight laser focusingponderomotive forcephase-stable accelerationplasma density profileparticle-in-cell simulation

0 comments

The pith

Reducing laser focal spot size to 0.8 micrometers increases maximum proton energy by 56.3 percent at fixed intensity, with an ideal plasma density profile adding another 61.3 percent.

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

The paper examines ways to raise the energy of protons accelerated by lasers when total laser energy is limited, a constraint that matters for potential uses such as cancer therapy. It demonstrates that shrinking the focal spot from the conventional 3 micrometers to 0.8 micrometers still yields substantially higher proton energies even when laser intensity remains unchanged. This gain arises because the smaller spot allows ponderomotive forces to drive electrons that set up stronger and faster-moving charge-separation electric fields. An analytically derived ideal plasma density profile further improves performance by enabling phase-stable acceleration, producing an additional energy increase. The improvements persist across a range of parameter choices, suggesting that better focusing and density control could reduce the scale of lasers required.

Core claim

At a laser focal spot size of 0.8 μm with fixed intensity, the maximum proton energy is 56.3 percent higher than at a 3 μm spot because ponderomotive-force-driven electrons generate stronger charge-separation fields that propagate at higher velocities. An analytically derived ideal plasma density profile promotes phase-stable proton acceleration and yields an additional 61.3 percent energy gain over the case of a tightly focused laser on a uniform-density planar target. These outcomes remain consistent under variations in laser and plasma parameters.

What carries the argument

Ponderomotive-force-driven electrons that produce propagating charge-separation fields, together with an analytically derived ideal plasma density profile that supports phase-stable acceleration.

Load-bearing premise

Particle-in-cell simulations and the associated modeling accurately represent the dominance of ponderomotive electrons and the resulting phase-stable fields without major three-dimensional effects or instabilities changing the reported gains.

What would settle it

An experiment that measures the maximum proton energy for a 0.8 μm focal spot versus a 3 μm spot at identical laser intensity on a near-critical-density target would directly test whether the 56.3 percent increase occurs.

Figures

Figures reproduced from arXiv: 2512.19341 by Deji Liu, Dongchi Cai, Guanqi Qiu, Jinqing Yu, Qianyi Ma, Xueqing Yan, Yinren Shou, Zheng Gong.

**Figure 2.** Figure 2: FIG. 2. (a) Accelerating electric field and proton phase-space distribution at t = 40 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. (a) Proton energy as a function of focal spot size under dif [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. (a) Density profiles derived from numerically solved and function fitted solutions of velocity-matching equation of protons and accel [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗

read the original abstract

Optimizing laser and plasma parameters is crucial for enhancing accelerated proton energy in laser-driven proton acceleration with finite laser energy for applications such as cancer therapy. Tight focusing plays a significant role in improving laser-driven proton acceleration, which is generally believed as a result of the enhancement of laser intensity. However, we find that even at a fixed laser intensity, reducing the focal spot size still enhances the proton energy. Through particle-in-cell simulations and theoretical modeling, we find that at a small spot size (0.8 {\mu}m), the maximum proton energy is enhanced by 56.3% compared to that obtained at a conventional spot size (3 {\mu}m). This improvement is attributed to the dominance of ponderomotive-force-driven electrons at reduced spot sizes, which generate stronger charge-separation fields that propagate at higher velocities. Furthermore, to optimize proton acceleration, we analytically derive an ideal plasma density profile that promotes phase-stable proton acceleration, yielding an additional energy increase of 61.3% over the case of a tightly focused laser interacting with a planar target of uniform density. These findings remain robust under parameter variations, indicating that advanced focusing techniques combined with optimized plasma profiles could relax the demand for high laser energies, thereby reducing the reliance on large-scale laser facilities in medical and scientific applications.

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.

Desk Editor's Note private letter to a colleague

Smaller focal spots at fixed intensity boost proton energy by 56% in 2D PIC runs, with another 61% from a derived density profile, though 3D effects remain untested.

read the letter

The main thing to know is that even at fixed laser intensity, shrinking the focal spot from 3 μm to 0.8 μm raises maximum proton energy by 56.3% in their 2D particle-in-cell simulations. They link this to ponderomotive-force-driven electrons that produce stronger, faster-moving charge-separation fields. An analytically derived ideal plasma density profile then adds a further 61.3% gain by keeping the acceleration phase-stable. The work combines the simulations with a theoretical model for the profile, and the derivation looks independent rather than tuned to the runs. They also state the gains hold across some parameter changes, which matters for applications that want to lower laser energy needs for medical proton sources. The soft spot is dimensionality. All key results come from 2D geometry, where transverse instabilities and diffraction are suppressed. At 0.8 μm spots these 3D effects become relevant and could reduce the longitudinal fields or detune the phase stability. The density profile is derived in 1D, so it inherits the limitation. Without error bars, a clear sensitivity study, or any 3D checks, the exact percentages feel optimistic. The methods are standard PIC, so the advance is the quantitative fixed-intensity spot-size result and the explicit profile rather than new tools. Citations follow the usual laser-plasma pattern with no obvious gaps. This paper is for researchers optimizing laser-driven ion acceleration, especially those targeting practical energies with limited laser power. A specialist in the subfield would find the spot-size mechanism and the profile useful to test in their own setups. It deserves peer review because the practical implication is clear and the evidence is grounded in simulations plus analysis, even if revisions for 3D validation would be needed.

Referee Report

2 major / 2 minor

Summary. The manuscript uses 2D particle-in-cell simulations and theoretical modeling to claim that, at fixed laser intensity, reducing the focal spot size from 3 μm to 0.8 μm increases maximum proton energy by 56.3% in near-critical-density plasma because ponderomotive-force-driven electrons produce stronger, faster-propagating charge-separation fields. An analytically derived ideal plasma density profile is shown to enable phase-stable acceleration, yielding an additional 61.3% energy gain over a uniform-density planar target.

Significance. If the reported gains survive multidimensional effects, the work offers a practical route to higher proton energies with limited laser energy, which could ease requirements for large-scale facilities in proton-therapy and other applications. The combination of simulation diagnostics with an explicit analytical derivation of the optimal density profile is a clear strength.

major comments (2)

[§3] §3 (PIC Simulations): The 56.3% energy enhancement at 0.8 μm spot size is obtained exclusively in 2D geometry; at this focal radius, 3D laser diffraction, electron filamentation, and azimuthal instabilities are expected to weaken the longitudinal charge-separation field and detune the phase-stable trajectory, directly threatening the central claim.
[§4] §4 (Analytical Model): The derivation of the ideal density profile is performed in a 1D slab geometry and therefore inherits the same limitation; no estimate is given for how transverse effects at 0.8 μm would modify the phase-stable condition or the quoted 61.3% additional gain.

minor comments (2)

[Figure 2] Figure 2: the proton energy spectra lack error bars or ensemble statistics, making it difficult to judge the robustness of the 56.3% figure under small parameter variations.
[Eq. (7)] Notation: the definition of the ponderomotive force term in Eq. (7) should explicitly state the averaging procedure used for the oscillating fields.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive comments and positive evaluation of the significance of our work. We address each major comment below and have revised the manuscript to better acknowledge the limitations of the 2D/1D approaches while defending the core physical insights.

read point-by-point responses

Referee: §3 (PIC Simulations): The 56.3% energy enhancement at 0.8 μm spot size is obtained exclusively in 2D geometry; at this focal radius, 3D laser diffraction, electron filamentation, and azimuthal instabilities are expected to weaken the longitudinal charge-separation field and detune the phase-stable trajectory, directly threatening the central claim.

Authors: We agree that 2D geometry cannot fully capture 3D diffraction, filamentation, and azimuthal instabilities, which may quantitatively reduce the reported gains. Nevertheless, the central mechanism—smaller focal spots enhancing ponderomotive electron acceleration and producing stronger, faster charge-separation fields—remains a robust longitudinal effect that 2D simulations capture accurately. In the revised manuscript we have added a dedicated paragraph in §3 discussing expected 3D corrections, citing relevant multidimensional studies, and stating that the 56.3% figure represents an upper-bound trend whose qualitative benefit of tight focusing at fixed intensity is expected to survive. Full 3D verification lies beyond the present scope. revision: partial
Referee: §4 (Analytical Model): The derivation of the ideal density profile is performed in a 1D slab geometry and therefore inherits the same limitation; no estimate is given for how transverse effects at 0.8 μm would modify the phase-stable condition or the quoted 61.3% additional gain.

Authors: The 1D derivation yields the exact phase-stable density profile by matching the plasma gradient to the field velocity. In the revision we have included a scaling estimate in §4 that accounts for transverse wave-number effects at 0.8 μm focal radius. This analysis indicates that the phase-stable condition remains approximately valid within the central beam region, with the 61.3% additional gain reduced by roughly 25% due to transverse spreading but still providing a substantial net improvement. The derived profile is presented as a practical target that can be further optimized in 3D geometries. revision: partial

standing simulated objections not resolved

Quantitative 3D PIC simulations that would provide precise corrections to the reported energy gains under realistic azimuthal instabilities and diffraction.

Circularity Check

0 steps flagged

No significant circularity; analytical derivation remains independent of simulation fits

full rationale

The paper reports 2D PIC results showing 56.3% proton-energy gain at reduced spot size (fixed intensity) and then presents a separate 1D analytical derivation of an ideal density profile that yields an additional 61.3% gain. No equation in the provided text reduces the derived profile to a fit of the simulation data, nor does any load-bearing step rely on self-citation of an unverified uniqueness theorem. The modeling of ponderomotive-driven electrons and phase-stable acceleration is constructed from standard plasma-physics relations rather than by renaming or tautologically re-expressing the simulation outputs. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

Work rests on standard domain assumptions of PIC modeling and an analytical derivation whose details are not visible in the abstract; no new entities postulated.

free parameters (2)

focal spot sizes
0.8 μm and 3 μm chosen for comparison; values are inputs rather than fitted outputs.
plasma density profile parameters
Near-critical density and the derived ideal profile shape are central to the optimization claim.

axioms (1)

domain assumption Particle-in-cell simulations faithfully reproduce the laser-plasma interaction physics at the reported scales.
Invoked implicitly for all quantitative results; standard but unverified here.

pith-pipeline@v0.9.0 · 5549 in / 1245 out tokens · 39586 ms · 2026-05-16T20:40:27.098613+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The longitudinal ponderomotive force at the axis F_p = −(m_e c²/4γ) ∂a_y²/∂x reaches its maximum value F_p = 3√3 m_e c² a₀² λ/(32π γ σ₀²) at x=Z_R/√3.
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Velocity matching requires the proton momentum to follow dp/dt = β c dp/dx = e E_x, yielding the density profile n_e/n_c = a₀/(2π²)·(2/Kx−1+√(1+4/Kx)).

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

4 extracted references · 4 canonical work pages · 1 internal anchor

[1]

While tighter focusing is generally expected to boost proton energy through increased laser intensity, the influence of the focal-spot size itself on energy enhancement has often been neglected. In this paper, we demonstrate that reducing the laser fo- cal spots of tightly focused lasers enhances proton energies in laser-near-critical-density plasma inter...

work page internal anchor Pith review Pith/arXiv arXiv 2026
[2]

(a) Variation of the energy of the highest-energy proton (E p) with time under focal spot sizes of 3µm and 0.8µm. (b) Solid lines represent proton energy (p) spectrum at the end of acceleration (t = 90T0), showing significantly higher proton energy with a smaller focal spot compared to the larger one. Dotted lines depict electron energy (e) spectrum at t ...

work page 2000
[3]

Ultrahigh-intensity lasers: Physics of the extreme on a tabletop,

Solving this we can obtain 1 4 ln 1−β 1+β + 1 2(1−β) − 1 2 = √ 2eE0 mpc2 x.(10) The electron density profile derived from this equation is plot- ted as the blue curve in Fig. 4(a). The left of Eq. 10 approxi- matesβ 2/(2−2β), and then we get an approximate solution ne nc = a0 2π2 2 Kx −1+ r 1+ 4 Kx ! ,(11) as represented by the red curve in Fig. 4(a), whe...

work page 1998
[4]

Laser-driven high-energy proton beams from cascaded acceleration regimes,

pp. 1–6. 28T. Ziegler, I. Göthel, S. Assenbaum, C. Bernert, F.-E. Brack, T. E. Cowan, N. P. Dover, L. Gaus, T. Kluge, S. Kraft,et al., “Laser-driven high-energy proton beams from cascaded acceleration regimes,” Nature Physics , 1–6 (2024). 29Y . Shou, X. Wu, K. H. Pae, G.-E. Ahn, S. Y . Kim, S. H. Kim, J. W. Yoon, J. H. Sung, S. K. Lee, Z. Gong,et al., “L...

work page 2024

[1] [1]

While tighter focusing is generally expected to boost proton energy through increased laser intensity, the influence of the focal-spot size itself on energy enhancement has often been neglected. In this paper, we demonstrate that reducing the laser fo- cal spots of tightly focused lasers enhances proton energies in laser-near-critical-density plasma inter...

work page internal anchor Pith review Pith/arXiv arXiv 2026

[2] [2]

(a) Variation of the energy of the highest-energy proton (E p) with time under focal spot sizes of 3µm and 0.8µm. (b) Solid lines represent proton energy (p) spectrum at the end of acceleration (t = 90T0), showing significantly higher proton energy with a smaller focal spot compared to the larger one. Dotted lines depict electron energy (e) spectrum at t ...

work page 2000

[3] [3]

Ultrahigh-intensity lasers: Physics of the extreme on a tabletop,

Solving this we can obtain 1 4 ln 1−β 1+β + 1 2(1−β) − 1 2 = √ 2eE0 mpc2 x.(10) The electron density profile derived from this equation is plot- ted as the blue curve in Fig. 4(a). The left of Eq. 10 approxi- matesβ 2/(2−2β), and then we get an approximate solution ne nc = a0 2π2 2 Kx −1+ r 1+ 4 Kx ! ,(11) as represented by the red curve in Fig. 4(a), whe...

work page 1998

[4] [4]

Laser-driven high-energy proton beams from cascaded acceleration regimes,

pp. 1–6. 28T. Ziegler, I. Göthel, S. Assenbaum, C. Bernert, F.-E. Brack, T. E. Cowan, N. P. Dover, L. Gaus, T. Kluge, S. Kraft,et al., “Laser-driven high-energy proton beams from cascaded acceleration regimes,” Nature Physics , 1–6 (2024). 29Y . Shou, X. Wu, K. H. Pae, G.-E. Ahn, S. Y . Kim, S. H. Kim, J. W. Yoon, J. H. Sung, S. K. Lee, Z. Gong,et al., “L...

work page 2024