{"id":"18cadc92-c8a3-4da9-a2b5-71abe2dde7f4","arxiv_id":"1908.05151","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A simulation study predicts that 0.1 µm Mylar foils with a preplasma generated by the laser pedestal yield about 25 MeV protons at 40 TW and about 600 MeV at 4 PW.","lead":"The authors combine hydrodynamic and particle-in-cell simulations to model how the nanosecond laser pedestal reshapes thin plastic foils before the main pulse arrives. They predict optimized proton energies of about 25 MeV for a 40 TW laser and about 600 MeV for a 4 PW laser.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 600 MeV headline rests on one 2D PIC maximum-energy point whose 3D translation is unquantified; the paper itself disclaims slope analysis, so the predicted energy should be treated as an upper bound until a 3D check is run.","rationale":"The reader and I identified the same load-bearing concern: the 2D PIC volumetric overestimation compromises the central quoted energies. I examined the strongest claim and found that it is not independently supported; the paper acknowledges the dimensional limitation but does not quantify its effect on the maximum energy, and no code/data release or convergence study is provided. I considered whether the scaling law from Esirkepov (epsilon ~ sqrt(P)) could independently support the 600 MeV claim; that scaling is cited in the introduction (epsilon_p ~ 173*sqrt(P) MeV) and gives ~550 MeV at 10 PW, which is in the same ballpark, so the 600 MeV figure is not implausible as a trend. However, the paper's own optimization argument changes the target from a flat foil to a preplasma-modified near-critical-density foil, so the simple scaling does not apply directly to the simulated case. I also considered the 25 MeV at 40 TW claim; the published experimental comparison (Ref. 57) is for a nanostructured foil at ~5 MeV at similar intensity, so 25 MeV is a factor of 5 above that, but the paper argues for a 7x enhancement from the preplasma, so the factor is internally consistent. The weakness is that this 7x enhancement is itself derived from the 2D maximum energies. Thus the single most load-bearing concern is the unquantified 2D-to-3D maximum-energy error. My concrete test is a targeted 3D simulation of the single best case; that is the only way to settle whether the 600 MeV number is real or an artifact. I agree with the reader's weakest_assumption, so agreement_with_reader is 'agree'. I maintain the CONDITIONAL verdict, because the qualitative trends and methodology are sound enough to publish once the dimensional uncertainty is either quantified or appropriately caveated.","tokens_in":23202,"tokens_out":2077,"duration_ms":17238,"concrete_test":"Run one 3D EPOCH simulation of the best 2D case (0.1 um Mylar, ASE 1.6e10 Wcm-2, main pulse 1.85e22 Wcm-2) with the same imported density profile on a reduced, axisymmetric-equivalent grid (e.g., 30x30 um transverse, 10 nm cells, scaled macroparticles), and compare the maximum proton energy at 400 fs. If the 3D maximum is within 20% of the 2D value, the headline stands; if it drops by 50% or more, the 2D maximum was inflated and the 600 MeV figure should be reported as an order-of-magnitude upper bound only.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim (abstract) is that optimized Mylar foils give ~25 MeV at ~40 TW and ~600 MeV at ~4 PW. The 25 MeV value is an interpolation from the 1.85e20 W/cm2 runs; the 600 MeV value comes from Fig. V.1 extrapolating the 1.85e22 W/cm2 run. Both rely entirely on the 2D EPOCH maximum proton energy. Sec. IV.D states explicitly that the 2D focal-spot representation 'overestimates protons from the highest intensity regions' and that 'no detailed slope analysis is made for the spectra since only a 3D simulation would give a meaningful result.' That admission covers the spectral shape and particle number, but the quoted maximum proton energy is a spectral-extremum statistic: in 2D the infinite transverse line focus over-weights the high-intensity region and sustains longer acceleration, so the maximum can be inflated even when the rescaled spectrum (Eq. IV.1) corrects only the particle number. The 600 MeV prediction is also obtained at 1.85e22 Wcm-2, an intensity regime where the 0.1 um foil is near relativistic transparency and where 2D and 3D transmission/absorption differ qualitatively; the paper's own Fig. IV.16 shows the optimization is a near-critical-density balancing act, which is precisely the regime most sensitive to dimensionality. No convergence study is presented (the paper reports one cell size, 10 nm, and one macroparticle count), and the beam is only 400 fs long; the spectra in Fig. IV.13/IV.15 are plotted without an axis scale, so the maximum point cannot be audited. Without a 3D simulation or a published convergence test, the 600 MeV headline is an unvalidated upper bound, not a prediction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies laser-driven proton acceleration from nanometre-thick Mylar foils under realistic conditions that include a nanosecond ASE pedestal. The authors combine 3D hydrodynamic simulations (3DLINE code) to obtain the preplasma density profiles with 2D particle-in-cell simulations (EPOCH) of the main pulse interaction. They scan foil thicknesses (0.1, 0.3, 0.9 µm), ASE pedestal intensities (0, 1.6e9, 1.6e10 W/cm²), and main-pulse intensities (1.85e20, 1.85e21, 1.85e22 W/cm²), reporting laser energy absorption, reflected/transmitted/trapped fields, electron spectra, and proton spectra. The central quantitative claims are the abstract's predictions of ~25 MeV protons at ~40 TW and ~600 MeV at ~4 PW. The paper identifies the optimal target parameters for each intensity level and discusses the competing mechanisms (TNSA, MVA, relativistic transparency).","tokens_in":23495,"tokens_out":1847,"duration_ms":19439,"significance":"If the central predictions hold, the paper would be a valuable design study for near-term multi-PW laser facilities such as ELI-Beamlines, connecting realistic preplasma modeling to proton-energy optimization in a parameter regime relevant to applications like hadron therapy. The authors are appropriately explicit about many limitations: they state that the 2D focal-spot representation overestimates protons from the highest-intensity regions, that particle numbers are overestimated, that no slope analysis is performed because only a 3D simulation would be meaningful, and that the FEOS equation of state cannot describe the overheated liquid state at 1 ns pedestal timescales. The paper also ships a detailed numerical setup (cell size, macroparticle numbers, solver choice, box size), which is a strength for reproducibility. However, the headline 25 MeV and 600 MeV values rest on single 2D PIC maximum-energy points, and the 3D translation of those values is unquantified; this limits the strength of the quantitative conclusions.","major_comments":[{"comment":"The headline predictions of ~25 MeV at 40 TW and ~600 MeV at 4 PW are taken from the maximum proton energy in 2D PIC runs, but Sec. IV.D explicitly states that the 2D representation of the focal spot as an infinite line focus 'overestimates protons from the highest intensity regions' and that no slope analysis is made because 'only a 3D simulation would give a meaningful result.' The maximum energy is a spectral-extremum statistic, not a particle-number statistic, so the caveat in Eq. (IV.1), which rescales only particle numbers, does not address how the 2D geometry modifies the maximum energy itself. The 600 MeV point is a single extrapolated point at 1.85e22 W/cm² in Fig. V.1, a regime where the 0.1 µm foil is near relativistic transparency and where 2D and 3D transmission/absorption differ qualitatively. The paper should either provide a 3D simulation of at least one case or explicitly reframe the abstract's numbers as 2D upper bounds with a quantitative estimate of the expected reduction in 3D.","section":"Abstract and Sec. IV.D / Fig. V.1"},{"comment":"No convergence study is presented. The paper reports a single cell size (10 nm) and a single macroparticle number (9×10^25 macroparticles per species, ~27 per cell initially). Since the central claim is a quantitative maximum proton energy in a near-critical-density regime, the sensitivity of the maximum energy to cell size and macroparticle number should be demonstrated, at least for the 0.1 µm foil at 1.85e22 W/cm² that produces the 600 MeV point. Without this, the reader cannot judge whether the headline value is numerically converged or is a resolution artifact.","section":"Sec. III and Sec. IV.D"},{"comment":"The proton spectra in Figs. IV.13 and IV.15 are plotted without axis scale or a y-axis label. The text refers to 'maximum proton energy' and 'a valley' and 'an exponentially decaying distribution,' but the reader cannot quantitatively compare the spectra or extract the maximum energies from the figures. Since the abstract's numbers are maximum energies, the paper should provide the numerical values of the maxima (e.g., in the text or in a table) and preferably plot the spectra with labeled axes and a clear indication of the maximum-energy point.","section":"Sec. IV.D, Figs. IV.13 and IV.15"},{"comment":"The hydrodynamic model's equation of state is stated to be unable to describe the overheated liquid phase at 1 ns ASE pedestal timescales (Sec. II.B). This is a potentially load-bearing limitation for the preplasma density profiles that are imported into the PIC simulations. The paper argues that the target mass dynamics at the rear side are still quantitatively correct, citing Ref. 31, but no direct validation or benchmark is provided for the specific Mylar foils and ASE intensities used here. If the preplasma density distributions are wrong by factors of order unity, the subsequent PIC-based claims about optimal thickness and maximum proton energy could shift. The authors should quantify the expected uncertainty in the imported density profiles or provide a sensitivity test.","section":"Sec. II.B and Sec. II.C"}],"minor_comments":[{"comment":"The paper contains several typographical errors and inconsistent notations, e.g., 'hve' in the caption of Fig. II.2, 'its'' instead of 'its', 'Fig' without a period in several places, and '1.85×1021 Wcm−2 and 1.85×1022 Wcm−2' missing a power in the Sec. V summary sentence. A careful proofreading pass is needed.","section":"General"},{"comment":"The derivation of the peak intensity formula (Eq. II.4) from Eq. (II.2) is not shown in detail; the limit m→0 involves a cancellation of m in the denominator, but the intermediate steps are omitted. Adding one line of derivation would improve clarity.","section":"Sec. II.A, Eq. (II.4)"},{"comment":"The factor 1.27×10^5 in the particle-number rescaling is presented as a fixed number, but it uses the assumed 5 µm effective emission length and the focal-spot radius R. The sensitivity of this factor to the choice of L and R should be stated, as the paper itself notes the correction does not produce an equivalent 3D spectrum.","section":"Sec. IV.D, Eq. (IV.1)"},{"comment":"Fig. IV.8 plots 'maximum magnetic field versus maximum electron energy' and shows a linear fit, but the physical basis for expecting a linear scaling is not explained in the text. A brief justification would help the reader assess the trend.","section":"Sec. IV.C and Fig. IV.8"},{"comment":"The reference list contains a duplicate: Refs. 48 and 49 are identical (Tian et al., Phys. Rev. Lett. 109, 115002 (2012)). This should be corrected in revision.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid parametric study but its central predictions are presented as quantitative expectation values while resting on 2D PIC maxima in a near-critical-density regime. The authors are unusually honest about the 2D limitations, but the abstract's 25 MeV and 600 MeV claims are stronger than the evidence supports. I would encourage the editor to request a 3D benchmark case or a clear upper-bound framing. The paper fits the journal's scope and the underlying methodology is sound; the requested changes are within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The genuinely new content is the parameter scan: three foil thicknesses, three main-pulse intensities, and three ASE pedestal levels, all fed from a hydrodynamic preplasma model into EPOCH PIC runs. The specific predictions — roughly 25 MeV at 40 TW and 600 MeV at 4 PW — are not in the earlier literature, and the optimum preplasma scale lengths (~110 nm at 1.85e20 W/cm2, ~45 nm at 1.85e21 W/cm2) are exactly the kind of quantitative output experimental planners can use. The hydro-to-PIC method is not new; the same group did it in Refs. 14–16. But the systematic thickness-versus-contrast study, with the discussion of transparency, MVA, and TNSA as competing mechanisms, is a real step beyond those papers.\n\nCredit where due: the authors are unusually candid about their own limitations. They state that the 2D focal-spot representation overestimates protons from the highest-intensity regions, that spectral slope analysis needs 3D to be meaningful, that their equation of state cannot describe the overheated liquid state on nanosecond pedestal timescales, and that the assumed contrasts are \"extremely optimistic.\" That level of explicit self-limitation is rare and should count in their favor.\n\nThe soft spots are real, though. The 600 MeV headline comes from a single 2D maximum-energy point at 1.85e22 W/cm2. The stress-test note is right: the rescaling in Eq. IV.1 corrects particle number, not the spectral maximum. In 2D, the infinitely long line focus over-weights the highest-intensity region and can sustain acceleration longer, and this is precisely the near-critical-density regime where 2D and 3D behavior differ most. There is no convergence study — one cell size, one macroparticle count — and the proton spectra in Figs. IV.13 and IV.15 have no axis scales, so the maximum point cannot be independently audited. The 25 MeV number is less exposed, since it comes from interpolation around the 1.85e20 W/cm2 cases, but it still inherits the 2D geometry.\n\nI would not call the central argument broken. The qualitative trends — an optimum preplasma scale length, harmful transparency for very thin foils at high intensity, thicker foils being less sensitive to contrast — are plausible and consistent with existing near-critical-density acceleration physics. What is not supported is the abstract's phrasing that these energies \"are expected.\" A single 2D run with acknowledged overestimation supports an upper-bound estimate, not a firm prediction. A 3D spot check or a convergence test would change my assessment substantially.\n\nWho should read this: groups planning foil experiments at ELI-Beamlines or similar facilities, and anyone wanting a compact map of how ASE preplasma shifts the optimum target thickness. It is not a benchmark or a validated scaling law. I would send it to a knowledgeable referee, let them push on the 2D-to-3D translation, and insist the abstract be softened and the axis scales be supplied. I would not cite the 600 MeV number as a prediction in my own work, but I would cite the scan itself once it is cleaned up.","headline":"Useful parameter scan with honest limitations, but the 600 MeV headline is a single 2D maximum-energy point and should be treated as an untested upper bound until a 3D check is done.","tokens_in":24171,"tokens_out":2280,"would_cite":false,"duration_ms":25193,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.38.Kd","52.65.Rr","52.50.Jm"],"model":"deepseek-v4-flash","headline":"The paper predicts that thin Mylar foils with realistic pedestal-formed preplasma deliver about 25 MeV protons from a 40 TW pulse and 600 MeV from 4 PW, and identifies the foil thickness and scale-length that maximize proton energy.","keywords":["laser-driven proton acceleration","thin foil targets","preplasma formation","amplified spontaneous emission","particle-in-cell simulation","hydrodynamic simulation","target normal sheath acceleration","relativistic transparency"],"falsifier":"Run the planned experiment the paper is designed to support: a 30 fs, 810 nm, P-polarized pulse of about 1.31 J focused to a 5 $\\mu$m spot on a 0.1 $\\mu$m Mylar foil with an amplified-spontaneous-emission pedestal near $1.6\\times10^{10}$ W cm$^{-2}$, and measure the cutoff energy of the rear-surface proton spectrum; a maximum below roughly 10 MeV, against the predicted $\\sim$25 MeV, would indicate that the 2D geometry inflates the energy. A cheaper check is a three-dimensional particle-in-cell simulation of the same nominal cases, which the paper itself identifies as the meaningful test of spectral shape.","tokens_in":22945,"feed_emoji":"⚛️","tokens_out":24771,"duration_ms":204120,"temperature":0.7,"pith_summary":"This paper argues that the nanosecond-long amplified-spontaneous-emission pedestal of a petawatt laser pulse cannot be ignored when predicting proton acceleration from ultra-thin foils: it pre-expands and curves the foil into a preplasma, and only simulations that start from that realistic density profile give trustworthy proton energies. The authors combine three-dimensional hydrodynamic simulations of the pedestal's effect with two-dimensional particle-in-cell simulations of the main pulse, and on that basis predict $\\sim$25 MeV protons from a $\\sim$40 TW laser and $\\sim$600 MeV protons from a $\\sim$4 PW laser using 0.1 $\\mu$m-thick Mylar foils. The optimization result is concrete: thinner foils give the highest proton energies in nearly all simulated cases, and an optimum preplasma scale-length near 110 nm boosts the maximum proton energy roughly sevenfold at $1.85\\times10^{20}$ W cm$^{-2}$ compared with the same foil without preplasma. If correct, the predictions offer a realistic, experimentally grounded route toward hundred-MeV-class laser-driven proton sources, the energy range needed for applications such as proton therapy.","feed_headline":"600 MeV protons predicted from 4 PW laser and thin plastic foil","feed_subtitle":"Modeling the nanosecond pedestal shows 0.1-micron foils can push laser protons from 25 toward 600 MeV.","key_machinery":"The load-bearing apparatus is a two-stage simulation chain. First, the 3DLINE hydrodynamic code models the 1 ns amplified-spontaneous-emission pedestal acting on a flat Mylar foil, using one-group radiation diffusion with tabulated opacities, a hybrid ray-tracing plus one-dimensional Helmholtz laser-deposition model, and a Thomas-Fermi-based equation of state, producing the curved, expanded density distributions that the main pulse actually encounters. Those distributions are imported as initial conditions into the EPOCH particle-in-cell code, which follows the 30 fs main pulse at peak intensities from $1.85\\times10^{20}$ to $1.85\\times10^{22}$ W cm$^{-2}$ and records the proton spectra. The paper's central diagnostic is the preplasma scale-length, defined as the distance between the relativistically corrected critical-density surfaces for the full and half peak electric field, plotted against the effective foil thickness (the thickness at the relativistically corrected critical density); this two-parameter view is what reveals the $\\sim$110 nm optimum and the inversion of the preplasma benefit at high intensity.","core_discovery":"The paper's central claim is that realistic preplasma must be part of any predictive simulation of laser-foil proton acceleration, and that when it is included, thin plastic (Mylar) foils deliver substantially higher maximum proton energies than an ideal flat foil would suggest: $\\sim$25 MeV for a $\\sim$40 TW pulse and $\\sim$600 MeV for a $\\sim$4 PW pulse. The acceleration itself follows standard physics: hot electrons leave the target, the resulting sheath field pulls protons from the rear surface (target-normal sheath acceleration, TNSA), and at the highest intensities a near-critical-density target switches to magnetic-vortex acceleration. The quantitative content is new: among the 0.1 $\\mu$m, 0.3 $\\mu$m and 0.9 $\\mu$m foils studied, the 0.1 $\\mu$m foil maximizes proton energy in nearly every case, and a preplasma scale-length of about 110 nm (measured between the relativistically corrected critical-density surfaces at full and half field amplitude) raises the maximum proton energy roughly sevenfold at $1.85\\times10^{20}$ W cm$^{-2}$ compared with a pristine foil. The authors also find that this preplasma benefit erodes as intensity rises: at $1.85\\times10^{22}$ W cm$^{-2}$ an ultra-thin target can become relativistically transparent, and the best proton energies then come from targets with little or no preplasma.","pith_inferences":["If the $\\sim$600 MeV prediction holds at $\\sim$4 PW, it would exceed the simple $\\varepsilon_p \\approx 173\\sqrt{P[\\mathrm{PW}]}$ MeV scaling quoted in the paper's introduction, suggesting the preplasma-modified thin-foil regime is more efficient per watt than flat-foil scaling implies.","The 40 TW case can be tested on existing laser systems, so the hydro-then-PIC pipeline could be validated against measured spectrum shape and divergence, not just the cutoff energy, before a multi-PW campaign is committed.","The paper's own volumetric-correction argument implies that in a real three-dimensional experiment the number of protons per unit energy will be roughly $10^5$ times smaller than the raw 2D spectra suggest, which matters for flux-hungry applications such as isotope production even when the cutoff energy is as predicted.","The same pipeline could be applied to other target materials or structured foils; the finding that the thinnest foil wins suggests probing tens-of-nanometre foils, with the caveat that such foils are more easily destroyed by the pedestal."],"forward_implications":["A $\\sim$40 TW, $\\sim$1.3 J laser focused on a 0.1 $\\mu$m Mylar foil is predicted to produce protons up to $\\sim$25 MeV, roughly seven times the energy the same foil would give without preplasma.","At $\\sim$4 PW ($\\approx1.85\\times10^{22}$ W cm$^{-2}$), the same optimization predicts proton energies near 600 MeV, close to the GeV range of the title and well above the near-100 MeV protons already produced at near-PW facilities.","Foil thickness becomes a dominant design lever: 0.1 $\\mu$m foils outperform 0.3 $\\mu$m and 0.9 $\\mu$m foils in almost every simulated case, so target manufacture should aim at the thinnest foils that survive the pedestal intact.","The optimum preplasma scale-length shrinks as intensity grows (about 110 nm at $10^{20}$ W cm$^{-2}$, about 45 nm at $10^{21}$ W cm$^{-2}$), so lower-intensity facilities should embrace a moderate pedestal while multi-PW facilities should push for higher contrast.","At the highest intensity a too-thin target whose peak density falls below the relativistically corrected critical density becomes transparent and acceleration is quenched; keeping the density near critical instead engages magnetic-vortex acceleration."],"supporting_citations":[{"why":"Supplies the EPOCH particle-in-cell code used for every main-pulse interaction and proton-spectrum result in the paper.","marker":"[20]"},{"why":"Supplies the 3DLINE hydrodynamic code that computes the ASE-pedestal preplasma density profiles later imported into the PIC simulations.","marker":"[23]"},{"why":"Prior work establishing how the ASE pedestal and preplasma modify petawatt-class laser thin-foil interactions, the effect this paper quantifies across a parameter grid.","marker":"[16]"},{"why":"Provides the hybrid ray-tracing and Helmholtz laser-deposition model used inside the hydrodynamic simulations to compute laser absorption.","marker":"[26]"},{"why":"Provides the equation-of-state package that closes the hydrodynamic model for the Mylar target.","marker":"[29]"},{"why":"Supplies the opacity data used in the one-group radiation-diffusion treatment of the preplasma.","marker":"[24]"},{"why":"Gives an experimental proton-energy measurement under similar conditions against which the low-intensity simulation spectra are compared.","marker":"[57]"},{"why":"Documents near-100 MeV protons at near-PW facilities, the benchmark that motivates the GeV-range optimization.","marker":"[2]"}],"fun_headline_variants":["Preplasma boosts thin-foil laser protons to 600 MeV","0.1-micron plastic foil maximizes laser proton energy","Realistic preplasma modeling predicts 600 MeV protons","Thin plastic foils with preplasma yield 600 MeV protons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative predictions rest on two-dimensional particle-in-cell simulations in which the focal spot is effectively an infinitely long line in the third dimension; as the paper itself notes, this geometry overestimates protons from the highest-intensity regions, so the headline $\\sim$25 MeV and $\\sim$600 MeV values could be lower in a real three-dimensional spot.","fun_headline_variants_meta":{"raw":{"variants":["Preplasma boosts thin-foil laser protons to 600 MeV","0.1-micron plastic foil maximizes laser proton energy","Realistic preplasma modeling predicts 600 MeV protons","Thin plastic foils with preplasma yield 600 MeV protons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000337,"raw_usage":{"total_tokens":1915,"prompt_tokens":1046,"completion_tokens":869,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":662,"completion_tokens_details":{"reasoning_tokens":796}},"tokens_in":662,"tokens_out":869,"duration_ms":8246,"temperature":1.0,"reasoning_tokens":796,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:21:05.406441+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the planned experiment the paper is designed to support: a 30 fs, 810 nm, P-polarized pulse of about 1.31 J focused to a 5 $\\mu$m spot on a 0.1 $\\mu$m Mylar foil with an amplified-spontaneous-emission pedestal near $1.6\\times10^{10}$ W cm$^{-2}$, and measure the cutoff energy of the rear-surface proton spectrum; a maximum below roughly 10 MeV, against the predicted $\\sim$25 MeV, would indicate that the 2D geometry inflates the energy. A cheaper check is a three-dimensional particle-in-cell simulation of the same nominal cases, which the paper itself identifies as the meaningful test of spectral shape.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the EPOCH particle-in-cell code used for every main-pulse interaction and proton-spectrum result in the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the 3DLINE hydrodynamic code that computes the ASE-pedestal preplasma density profiles later imported into the PIC simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior work establishing how the ASE pedestal and preplasma modify petawatt-class laser thin-foil interactions, the effect this paper quantifies across a parameter grid."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the hybrid ray-tracing and Helmholtz laser-deposition model used inside the hydrodynamic simulations to compute laser absorption."},{"cited_title":"Faik , author A","cited_arxiv_id":null,"evidence_quote":"Provides the equation-of-state package that closes the hydrodynamic model for the Mylar target."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the opacity data used in the one-group radiation-diffusion treatment of the preplasma."},{"cited_title":"Margarone , author O","cited_arxiv_id":null,"evidence_quote":"Gives an experimental proton-energy measurement under similar conditions against which the low-intensity simulation spectra are compared."},{"cited_title":"Higginson , author R","cited_arxiv_id":null,"evidence_quote":"Documents near-100 MeV protons at near-PW facilities, the benchmark that motivates the GeV-range optimization."}],"review_version":1}