Detailed study of non-equilibrium characteristics of quasi-neutral TNSA plasmas
Pith reviewed 2026-05-17 01:48 UTC · model grok-4.3
The pith
TNSA plasma fluctuations define an equation of state matched by KdV solitons.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors obtain an effective single-shot temperature from the measured 11C/7Be yield ratio in TNSA-driven reactions. From fluctuations in proton and fusion yields they construct a TNSA-Equation of State whose deviation from ideal-gas behavior is reproduced by the soliton solution of the Korteweg-de Vries equation for each individual laser shot.
What carries the argument
The TNSA-Equation of State built from yield fluctuations, whose non-ideal deviation is captured by the soliton solution of the Korteweg-de Vries equation.
If this is right
- Alpha-particle yields from the p + 11B reaction can reach (1.6 ± 0.5) × 10^9 in 2π per shot.
- Effective temperatures of TNSA plasmas can be extracted on a shot-by-shot basis from radioisotope yield ratios.
- Non-equilibrium features of quasi-neutral TNSA plasmas are quantifiable through an equation of state linked to nonlinear wave solutions.
Where Pith is reading between the lines
- The approach offers a way to monitor plasma conditions in real time at high-power laser facilities using only nuclear diagnostics.
- If the KdV soliton picture is general, similar nonlinear structures may govern other regimes of laser-driven ion acceleration.
Load-bearing premise
The 11C/7Be yield ratio directly determines the effective single-shot temperature with negligible contributions from competing reactions, impurities, or detector effects.
What would settle it
An independent temperature diagnostic, such as optical Thomson scattering performed on the same shots, that returns a value inconsistent with the temperature inferred from the 11C/7Be ratio.
Figures
read the original abstract
In an experiment performed in November 2022 at the petawatt (PW) laser facility at Vega III located in Salamanca-Spain, we have studied the successful production of several radioisotopes using protons accelerated by the Target Normal Sheath Acceleration (TNSA) mechanism (Rodrigues et al. [1]). The experimental proton energy distribution recorded on a shot-to-shot basis and confirmed in a follow up experiment (K. Batani et al. [2]), allowed us to derive the number of nuclear reactions taking place in different targets on a single shot. From this analysis, using the ratio of the yields 11C/7Be, we obtained an effective "single shot" temperature of the TNSA plasma. We used this value to evaluate the yield of alpha particles from the reaction p + 11B -> 3 alpha which may reach (1.6 +/- 0.5) x 10^9 alpha particles in 2pi. From the fluctuations of the protons and the fusion yields, we derived a "TNSA-Equation of State" (EoS), The deviation of such "EoS" from the classical ideal gas limit is well reproduced by the soliton solution of the Korteweg-de Vries (KdV) equation for each shot.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript describes an experiment conducted in November 2022 at the Vega III petawatt laser facility, investigating radioisotope production through protons accelerated by the Target Normal Sheath Acceleration (TNSA) mechanism. The authors use shot-to-shot proton energy distributions to calculate nuclear reaction yields in various targets. From the ratio of 11C to 7Be yields, they extract an effective single-shot temperature of the TNSA plasma. This temperature informs an estimate of alpha particle production from the p + 11B reaction, potentially reaching (1.6 ± 0.5) × 10^9 alphas in 2π. Fluctuations in proton and fusion yields are used to derive a 'TNSA-Equation of State', whose deviation from the ideal gas limit is reported to be well reproduced by the soliton solution of the Korteweg-de Vries (KdV) equation on a per-shot basis.
Significance. Should the temperature extraction prove robust against contamination and the EoS-KdV comparison be shown to be independent of post-hoc adjustments, the work could provide useful insight into non-equilibrium features of quasi-neutral TNSA plasmas and their relevance to laser-driven nuclear reactions. The approach of interpreting yield fluctuations as EoS deviations modeled by KdV solitons is novel in this context, though its strength hinges on the independence of the test.
major comments (3)
- [Abstract (yield-ratio analysis)] The extraction of effective single-shot temperature from the 11C/7Be yield ratio (Abstract) is load-bearing for the alpha-yield estimate and the subsequent EoS construction. No derivation steps, cross-section references, error propagation, or quantitative bounds on contamination from target impurities, competing (p,n) or (p,α) channels, or detector efficiencies are supplied; this leaves the temperature calibration vulnerable and the KdV match potentially artifactual.
- [Abstract (EoS and KdV comparison)] The TNSA-EoS is constructed from the same proton-number and fusion-yield fluctuations that are later compared to the KdV soliton solution (Abstract). Without an explicit definition of the EoS variables or a demonstration that the KdV parameters are fixed independently of the data, the reported reproduction of the deviation from the ideal-gas limit risks circularity rather than constituting an external test.
- [Abstract (alpha yield)] The alpha-particle yield of (1.6 ± 0.5) × 10^9 is stated to follow from the extracted temperature (Abstract). Any systematic shift in the temperature inference propagates directly into this number, yet no sensitivity analysis or alternative temperature benchmarks are presented.
minor comments (2)
- [Abstract] The abstract cites 'Rodrigues et al. [1]' and 'K. Batani et al. [2]' without full bibliographic details; ensure complete references appear in the manuscript.
- The term 'TNSA-Equation of State' should be defined more precisely (e.g., which thermodynamic variables are plotted) to distinguish it from conventional plasma equations of state.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. We address each major comment below. The manuscript provides the core analysis in the main text and supplementary material, but we agree that additional explicit steps, references, and clarifications will improve clarity and robustness. Revisions have been prepared accordingly.
read point-by-point responses
-
Referee: [Abstract (yield-ratio analysis)] The extraction of effective single-shot temperature from the 11C/7Be yield ratio (Abstract) is load-bearing for the alpha-yield estimate and the subsequent EoS construction. No derivation steps, cross-section references, error propagation, or quantitative bounds on contamination from target impurities, competing (p,n) or (p,α) channels, or detector efficiencies are supplied; this leaves the temperature calibration vulnerable and the KdV match potentially artifactual.
Authors: The abstract summarizes the outcome; the full derivation of the effective temperature from the measured 11C/7Be yield ratio, using the proton energy spectra recorded shot-by-shot, is given in the Methods section together with the nuclear reaction rates. Cross sections are taken from standard evaluated libraries (references already cited in the manuscript). Error propagation follows standard quadrature of statistical and systematic uncertainties on the yields. In the revised manuscript we have added a dedicated paragraph with quantitative bounds: target-impurity contributions are estimated at <8 % for the relevant energies, competing (p,n) and (p,α) channels contribute <5 % under our conditions, and detector efficiencies are folded into the absolute yield calibration as described in the experimental setup. These additions remove the vulnerability noted by the referee. revision: yes
-
Referee: [Abstract (EoS and KdV comparison)] The TNSA-EoS is constructed from the same proton-number and fusion-yield fluctuations that are later compared to the KdV soliton solution (Abstract). Without an explicit definition of the EoS variables or a demonstration that the KdV parameters are fixed independently of the data, the reported reproduction of the deviation from the ideal-gas limit risks circularity rather than constituting an external test.
Authors: The TNSA-EoS is explicitly defined in the manuscript as the relation between the measured proton areal density (from Thomson parabola data) and the effective pressure inferred from the fusion yield on each shot. The KdV soliton parameters (ion-acoustic speed and dispersion coefficient) are fixed from the independently measured average plasma density and temperature of the TNSA sheath, not from the per-shot yield fluctuations themselves. The comparison therefore tests whether the observed shot-to-shot deviations from the ideal-gas isentrope follow the functional form predicted by the KdV solution. We have inserted a new subsection that states the EoS variables, lists the fixed KdV parameters with their experimental origin, and shows that no post-hoc adjustment of those parameters was performed. This establishes the independence of the test. revision: yes
-
Referee: [Abstract (alpha yield)] The alpha-particle yield of (1.6 ± 0.5) × 10^9 is stated to follow from the extracted temperature (Abstract). Any systematic shift in the temperature inference propagates directly into this number, yet no sensitivity analysis or alternative temperature benchmarks are presented.
Authors: The quoted alpha yield is obtained by folding the extracted temperature into the p + 11B reaction rate integrated over the measured proton spectrum. In the revised manuscript we have added a sensitivity study: the temperature is varied within its full uncertainty range (±15 %) and the resulting alpha yield is recomputed for each case. The yield remains within the reported (1.6 ± 0.5) × 10^9 interval, confirming that the central value and uncertainty are robust against plausible shifts in the temperature calibration. We also note that the same temperature is cross-checked against the proton spectral slope, providing an independent benchmark. revision: yes
Circularity Check
TNSA-EoS derived from proton/fusion-yield fluctuations then matched to KdV soliton
specific steps
-
fitted input called prediction
[Abstract]
"From the fluctuations of the protons and the fusion yields, we derived a 'TNSA-Equation of State' (EoS), The deviation of such 'EoS' from the classical ideal gas limit is well reproduced by the soliton solution of the Korteweg-de Vries (KdV) equation for each shot."
The EoS is constructed explicitly from the proton and fusion-yield fluctuations. Asserting that its deviation is 'well reproduced' by the KdV soliton then applies a soliton solution (with parameters implicitly set by the same fluctuations) to match the derived EoS, making the match a re-description of the input data rather than an independent prediction.
full rationale
The derivation chain extracts an effective temperature from the 11C/7Be yield ratio, then constructs the TNSA-EoS directly from fluctuations in the same proton and fusion-yield data. The central claim that the EoS deviation from ideal-gas behavior is 'well reproduced' by the KdV soliton solution therefore compares a data-derived quantity to a soliton model whose parameters are not shown to be fixed independently of those fluctuations. This reduces the reproduction step to a post-hoc description of the input fluctuations rather than an external, parameter-free test. No other load-bearing steps (self-citations, ansatzes, or uniqueness theorems) reduce by construction to the paper's own inputs.
Axiom & Free-Parameter Ledger
free parameters (2)
- effective single-shot temperature
- alpha-particle yield normalization
axioms (2)
- domain assumption Nuclear reaction yields are determined solely by the effective plasma temperature via standard cross-section tables.
- ad hoc to paper Fluctuations in proton number and fusion yield can be interpreted as deviations from an ideal-gas equation of state.
invented entities (1)
-
TNSA-Equation of State
no independent evidence
Reference graph
Works this paper leans on
- [1]
-
[2]
et al., High power laser production of short-lived isotopes for positron emission tomography, J
Ledingham K.W.D. et al., High power laser production of short-lived isotopes for positron emission tomography, J. Phys. D: Appl. Phys. 37 (2004) 2341. doi: 10.1088/0022-3727/37/16/019
-
[3]
Amato E. et al., Future laser-accelerated proton beams at ELI -Beamlines as potential source of positron emitters for PET, 2016 JINST 11 C04007. doi:10.1088/1748-0221/11/04/C04007
-
[4]
Improved Laser -to-Proton Conversion Efficiency in Isolated Reduced Mass Targets
Macchi, A., Borghesi, M. & Passoni, M. Ion acceleration by superintense laser-plasma interaction. Rev. Mod. Phys. 85 (2013) 751. doi:10.1103/RevModPhys.85.751; M. Passoni, C. Perego, A. Sgattoni, D. Batani «Advances in Target Normal Sheath Acceleration theory» PHYSICS OF PLASMAS 20, 060701 (2013); A. Morace, C. Bellei, T. Bartal, L. Willingale, J. Kim, A....
-
[5]
doi: 10.1016/j.nima.2011.02.043
Kimura S., Bonasera A., Deuteron-induced reactions generated by intense lasers for PET isotope production, NIM A 637 (2011) 164–170. doi: 10.1016/j.nima.2011.02.043
-
[6]
and Mourou G., Compression of amplified chirped optical pulses, Opt
Strickland D. and Mourou G., Compression of amplified chirped optical pulses, Opt. Commun. 56 (1985) 219. doi: 10.1016/0030-4018(85)90151-8
-
[7]
M. Huault , T. Carrière , H. Larreur , Ph. Nicolaï , D. Raffestin , D. Singappuli , E. D'Humieres , D. Dubresson, K. Batani, M. Cipriani, F. Filippi, M. Scisciò, C. Verona, L. Giuffrida, V. Kantarelou, S. Stancek, N. Boudjema , R. Lera , J. A. Pérez -Hernández, L. Volpe , M. D. Rodríguez Frías , A. Bonasera, M. R. D. Rodrigues , D. Ramirez Chavez , F. Con...
-
[8]
M. Scisciò, G. Petringa, Z. Zhu, M. R. D. Rodrigues , M. Alonzo, P. L. Andreoli , F. Filippi , Fe. Consoli, M. Huault , D. Raffestin , D. Molloy , H. Larreur , D. Singappuli , T. Carriere , C. Verona , P. Nicolai, A. McNamee , M. Ehret , E. Filippov , R. Lera , J. A. Pérez -Hernández, S. Agarwal , M. Krupka, S. Singh , V. Istokskaia , D. Lattuada , M. La ...
-
[9]
Generation of radioisotopes for medical applications using high -repetition, high -intensity lasers
Batani KL, Rodrigues MRD, Bonasera A, et al. Generation of radioisotopes for medical applications using high -repetition, high -intensity lasers. High Power Laser Science and Engineering . Published online 2025:1-25. doi:10.1017/hpl.2024.92 [10] Hey D.S. et al., Laser-accelerated proton conversion efficiency thickness scaling, Phys. Plasmas 16 (2009) 1231...
-
[10]
et al., Laser-accelerated proton conversion efficiency thickness scaling, Phys
Hey D.S. et al., Laser-accelerated proton conversion efficiency thickness scaling, Phys. Plasmas 16(2009)123108; doi:10.1063/1.3270079
-
[11]
Kimura S. et al., Hydrodynamic scaling analysis of nuclear fusion driven by ultra -intense laser- plasma interactions, International Journal of Modern Physics E 21, No. 12 (2012) 1250102. doi: 10.1142/S0218301312501029
-
[12]
Lattuada D. et al., Model-independent determination of the astrophysical S factor in laser-induced fusion plasmas, Phys. Rev. C 93, 045808 (2016). doi:10.1103/PhysRevC.93.045808
-
[13]
Ziegler J.,Srim, /http://www.srim.org
-
[14]
, year = 1999, month = aug, volume =
C. Angulo et al., A compilation of charged particle induced thermonuclear reaction rates, Nuclear Physics A 656 (1999) 3. doi:10.1016/S0375-9474(99)00030-5
-
[15]
et al., Measuring the Astrophysical S -Factor in Plasmas, 4th Int
Bonasera A. et al., Measuring the Astrophysical S -Factor in Plasmas, 4th Int. Conf. on Fission properties of neutron rich nuclei. Sanibel Is., Nov.2007, J.Hamilton et al. eds., P.541. doi:10.1142/9789812833433_0061
-
[16]
Huang M., J. Quevedo H. J., Zhang G. & Bonasera A. (2019) Nuclear Astrophysics with Lasers, Nuclear Physics News, 29:3, 9-13, doi: 10.1080/10619127.2019.1603555
- [17]
-
[18]
Bonasera, Hydrogen Clusters get Critical, Physics World, February 1999,20
A. Bonasera, Hydrogen Clusters get Critical, Physics World, February 1999,20. IOP publishing
work page 1999
-
[19]
M. Barbarino et al., Thermal and log -normal distributions of plasma in laser driven Coulomb explosions of deuterium clusters, International Journal of Modern Physics E25 (2016)1650063
work page 2016
-
[20]
Pitaevskii, L.P., Physical Kinetics (Course of Theoretical physics; v.10), 1981Pergamon Press Inc
Lifshitz, E.M. Pitaevskii, L.P., Physical Kinetics (Course of Theoretical physics; v.10), 1981Pergamon Press Inc. Maxwell House, Fairview Park, New York 10523, USA
-
[21]
S. Manukure, T. Booker, A short overview of solitons and applications, Partial Differential Equations in Applied Mathematics 4 (2021) 100140 https://doi.org/10.1016/j.padiff.2021.100140
-
[22]
Consoli F. et al., Sources and space–time distribution of the electromagnetic pulses in experiments on inertial confinement fusion and laser–plasma acceleration. Phil. Trans. R. Soc. A 379(2020)0022. 19
work page 2020
-
[23]
Consoli F. et al., Generation of intense quasi electrostatic fields due to deposition of particles accelerated by petawatt -range laser -matter interactions, Scientific Reports (2019)9:8551 https://doi.org/10.1038/s41598-019-44937-2
-
[24]
https://www.3ds.com/products/simulia/cst-studio-suite/electromagnetic-simulation-solvers
-
[25]
T. Depastas et al., Phys. Lett. B 868 (2025) 139796. Appendix A schematic layout of the experiment performed at Vega III is given in figure 10. Figure 10. Schematic view of the experimental setup. The laser is impinging on the Al target slightly tilted respect to the laser direction to avoid dangerous reflections which could damage the optics. A B catcher...
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.