{"id":"ee25ab9d-b3bb-4e75-9674-cc24fb3d05a0","arxiv_id":"2411.11831","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"High-energy photons from a gamma-ray burst jet can transmute protons in the surrounding stellar envelope into neutrons, potentially enabling neutron-capture nucleosynthesis in the cocoon.","lead":"This paper proposes that the intense gamma-ray flash of a long gamma-ray burst can convert protons in the surrounding stellar material into free neutrons through photon-induced reactions, creating a new environment for making heavy elements. The idea matters because it offers a possible way for collapsar gamma-ray bursts to produce r-process elements without starting with neutron-rich matter.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The one-zone reaction network violates photon number conservation: at the fiducial jet-head density and thickness, Eq. (10) demands roughly 10^4 neutrons per incident 300 MeV photon, so only a thin surface layer of baryons is processed.","rationale":"The reader identified the photon flux extrapolation in Sec. 3.2 as the weakest assumption. I agree that the flux is uncertain, but the most load-bearing problem is internal: the one-zone network treats the entire jet-head volume as uniformly bathed in unattenuated radiation, even though the quoted baryon density and length scale make the region optically thick by many orders of magnitude. This is not a question of an astrophysical input being uncertain; it is an inconsistency between Eq. (10) and the stated parameters, independent of the absolute flux value. The photon-supply argument shows the integrated neutron production cannot exceed the incident photon number, so the homogeneous network in Eqs. (20)-(23) overproduces neutrons by a factor ~n_b sigma_pgamma L ~ 10^4. The paper deserves credit for using the ANL-Osaka cross sections, for being transparent about simplifications, and for phrasing the proposal as an estimate. However, the nucleosynthesis calculations in Sec. 6 inherit these overestimated neutron fractions. The reader's CONDITIONAL verdict remains appropriate: the idea is not falsified, but it is conditional on demonstrating that photon attenuation does not reduce the neutron-rich ejecta to a negligible column. My disagreement with the reader is about which assumption is most load-bearing, not about the overall verdict.","tokens_in":29,"tokens_out":37554,"duration_ms":1064017,"concrete_test":"Compute the total photopion neutron yield per unit area: Y_n = integral Phi(E) [sigma_pgamma(E) / (sigma_pgamma(E) + sigma_att(E))] dE over the 100-500 MeV band, using Eq. (1) with the Sec. 4.7 parameters, and compare with the baryon column n_b L of the jet head. If Y_n/(n_b L) is much less than unity, the cocoon electron fraction after mixing remains near 0.5 and the low-Y_e initial conditions (0.0035, 0.034, 0.334) used in Sec. 6 are not realized. This check uses only the paper's own fiducial parameters and standard photon attenuation cross sections.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that stellar material becomes inundated with neutrons rests on the one-zone network of Sec. 4.7, which uses n_b = 10^24 cm^-3 and a jet-head length scale L ~ 5 x 10^7 cm (Sec. 5.4). At the Delta-resonance peak (E_gamma ~ 300 MeV), sigma_pgamma ~ 3 x 10^-28 cm^2, so the photopion optical depth across the jet head is tau_pgamma = n_b sigma_pgamma L ~ 1.5 x 10^4. The rate r_n = n_p integral Phi sigma dE in Eq. (10), integrated over a volume A L, therefore produces ~tau_pgamma neutrons per photon entering the zone. This is impossible without photon regeneration: the homogeneous network implicitly lets each 300 MeV photon convert ~10^4 protons. In reality, the photon flux attenuates over a length lambda ~ 1/(n_b sigma_tot) ~ 10^2-10^3 cm (Compton and pair production dominate the opacity at these energies), so only a thin skin of the jet head is irradiated. The resulting neutron-rich mass is far smaller than the ~0.1-1 M_sun cone mass used in Sec. 5.5 to set the low electron fractions of Sec. 6. Even if the Sec. 3.2 photon flux is accepted at face value, the r-process scenarios are not the natural outcome of the stated jet-head conditions.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the intense photon flux at the head of a long-GRB jet can drive hadronic photoproduction, notably gamma + p -> pi+ + n, converting initially proton-only stellar material into neutron-rich material in situ. The authors estimate the jet-head photon spectrum by rescaling observed GRB spectra to the source, compute photohadronic rates using ANL-Osaka cross sections, model the proton/neutron/pion populations with a one-zone reaction network, estimate escape and injection timescales, and then run PRISM nucleosynthesis for three assumed density/temperature/electron-fraction trajectories. They conclude that collapsar long GRBs could be a new site for r-process and i-process nucleosynthesis and discuss several observational signatures.","tokens_in":29409,"tokens_out":8796,"duration_ms":91206,"significance":"If the central mechanism operates as claimed, this paper would add collapsar long GRBs to the list of candidate r-process sites and connect jet physics directly to heavy-element nucleosynthesis. The use of the documented ANL-Osaka partial-wave amplitudes for the photoproduction cross sections and the use of a standard reaction network (PRISM) are strengths, as is the concrete enumeration of observable tests such as a neutron precursor, kilonova signatures, gamma-ray lines, and associated neutrinos. The main reservations concern not the cross-section physics but whether the assumed photon field can actually process the claimed baryon mass, and whether the nucleosynthesis calculations are connected to the photoproduction network by a forward model.","major_comments":[{"comment":"The homogeneous reaction network implicitly neglects photon depletion, and this is load-bearing for the central claim. With the fiducial values n_b = 10^24 cm^-3 and L ~ 5 x 10^7 cm, the optical depth across the jet head for a ~300 MeV photon is tau ~ n_b sigma L ~ 1.5 x 10^4 for sigma ~ 3 x 10^-28 cm^2. Integrating Eq. (10) over the whole zone therefore asks each incident photon to convert ~10^4 protons, which violates photon number conservation. In reality the photon flux attenuates on a scale lambda ~ 1/(n_b sigma_tot) ~ 10^2-10^3 cm, with Compton and pair-production opacity comparable or larger, so only a thin skin of the jet-head material is processed. The cone-mass estimate M ~ 0.1-1 M_sun in Secs. 5.5-5.6 and the very low initial electron fractions used in Sec. 6 therefore do not follow from the stated jet-head conditions. Please provide an attenuation-corrected estimate of the processed baryon mass and neutron yield, or explicitly reframe the claim as applying to a thin surface layer rather than to the full jet-head cone.","section":"Sec. 4.7, Eq. (10); Sec. 5.4"},{"comment":"The estimate Phi_gamma ~ 10^33 ph/cm^2/s/keV at 10^6 keV is an order-of-magnitude extrapolation from Figure 2 of Ahlgren et al. (2019), scaled by (d_l/r_jh)^2 under the assumptions that the sub-photospheric non-thermal tail originates at r_jh ~ 10^10 cm, is isotropic, and is not absorbed before reaching the baryons. Because the proposed mechanism fails for flux below Phi_gamma ~ 10^25 ph/cm^2/s/keV at 10^6 keV (Sec. 4.7), the nominal margin of several orders of magnitude does not remove the need for a sensitivity analysis; the true flux at the interaction region could be much lower if the high-energy tail is produced at larger radii, is beamed, or is attenuated by pair production and scattering. Please state the plausible uncertainty range on Phi_gamma and show explicitly at which flux values the neutron-rich outcome is lost.","section":"Sec. 3.2"},{"comment":"The nucleosynthesis simulations are not yet connected to the photoproduction network. The initial electron fractions Ye = 0.334, 0.034, and 0.0035 are imposed from the extreme-conversion formulas in Sec. 5.5, which assume that all baryons in the jet head are converted to neutrons and then mix with envelope protons. Once the photon-attenuation issue of Sec. 4.7 is taken into account, the fraction of converted baryons is far smaller, so the actually attainable Ye will be higher and closer to the quasi-equilibrium range Ye ~ 0.52-0.66 quoted in Sec. 4.8. The statement that the integration of simulations (a)-(c) 'closely aligns with the full range of solar r-process residuals' should be presented as an illustration of possible outcomes for chosen parameters, not as a prediction of the proposed mechanism, unless a forward model linking the jet-head conditions to Ye and the density trajectory is supplied.","section":"Sec. 6 and Sec. 5.5"}],"minor_comments":[{"comment":"There are several typographical errors: 'photodisentigration' in Sec. 1, 'copius' in Sec. 3.2, and 'Relativitiy' in Sec. 8; these should be corrected.","section":"Throughout"},{"comment":"The caption reads 'Cross section ( -barns)', which appears to have a missing mu symbol; please confirm that the units are microbarns.","section":"Fig. 3 caption"},{"comment":"The adopted neutron escape time tau_esc = 5 x 10^-4 s and the injection time tau_inj = 10^-4 s are not clearly derived from the stated geometry: Sec. 4.6 estimates tau_esc ~ L/v with L = r tan(theta), while Sec. 5.4 gives L ~ 5 x 10^7 cm, for which L/c ~ 0.17 s. Please define a single length scale and velocity and provide corresponding timescales, or present these values explicitly as free parameters in a sensitivity scan.","section":"Secs. 4.6-4.7"},{"comment":"The claim that the conclusions are independent of starting composition is not demonstrated, because the network in Eqs. (20)-(23) tracks only protons, neutrons, and charged pions and does not include photodisintegration of 4He or 12C. Please either add the relevant photodisintegration channels or soften the statement.","section":"Sec. 4.7"},{"comment":"Equation (48) treats tau1, tau2, and xi as adjustable parameters; for a work proposing a new mechanism, a brief justification of the adopted values, or a footnote indicating which physical regimes they represent, would improve clarity.","section":"Sec. 5.6"}],"recommendation":"major_revision","confidential_remarks":"This is a provocative and potentially interesting paper, but the photon-attenuation problem in the one-zone network is severe and affects the central quantitative claims and the nucleosynthesis scenarios. I would not reject outright, because the cross-section treatment is sound and the underlying idea is testable, but the manuscript needs a substantial reworking of the processed-mass estimate and a correspondingly softened or better-justified connection to the r-process conclusion before it can be published."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe new thing here is real: photohadronic neutron production as an r-process route in long GRBs, distinct from the usual cosmic-ray/neutrino framing of pion photoproduction. The cross sections come from ANL-Osaka, the nuclear network work is standard, and the authors are transparent about the simplifications. I believe them that the mechanism is worth exploring.\n\nThe soft spots are not minor. The jet-head photon flux is an order-of-magnitude extrapolation from a published figure, scaled back with (d_l/r_jh)^2 under assumptions about isotropy and subphotospheric origin. That alone would justify a conditional verdict. But the bigger problem is the reaction network in Sec 4.7. At the fiducial n_b=1e24 cm^-3, L~5e7 cm, and sigma~3e-28 cm^2 near the Delta resonance, the optical depth is ~1.5e4. A one-zone network that just multiplies the rate by volume lets each photon convert 10^4 protons. Photon number is not conserved. In reality the beam attenuates on a length scale ~10^2-10^3 cm, so only a thin layer is processed. The resulting neutron-rich mass is set by the total photon number above threshold, not by the cone mass of 0.1-1 M_sun used in Sec 5.5. The nucleosynthesis scenarios (a)-(c) are hand-picked to cover a range of Ye and density histories; they show what the network can do, not what the jet head actually produces.\n\nSo the central claim of the paper—that stellar material becomes inundated with neutrons—is not supported by the present calculation. The idea may survive a proper transport treatment: the photon flux is large enough that even a thin processed layer could yield interesting amounts of neutron-rich ejecta, but that needs to be shown with an actual calculation that tracks photon attenuation and the advancing jet head.\n\nWho is this for? People working on GRB nucleosynthesis and photohadronic processes. It deserves a serious referee, because the idea is new and the flaw is fixable. I would send it to review, but with a strong recommendation for major revision: address photon number conservation, propagate the flux uncertainties, and derive the threshold flux from a proper attenuation calculation. As it stands, I would not cite the quantitative yields.","headline":"A genuinely new r-process mechanism for long GRBs, but the one-zone network violates photon number conservation and the neutron yields are likely overestimated by orders of magnitude.","tokens_in":29999,"tokens_out":4678,"would_cite":false,"duration_ms":46402,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Long gamma-ray bursts can manufacture neutrons from photons at the jet head, converting proton-only stellar gas into fuel for heavy-element nucleosynthesis.","keywords":["gamma-ray bursts","hadronic photoproduction","neutron production","r-process nucleosynthesis","collapsars","cocoon","photomeson reactions","electron fraction"],"falsifier":"Compute the actual photon flux at the jet-head interaction radius with a self-consistent sub-photospheric dissipation model; if the flux at $10^{6}$ keV falls below ~$10^{25}$ ph/$cm^{2}$/s/keV, neutron production is slower than the ~$10^{-4}$ s escape timescale and the claimed mechanism fails.","tokens_in":28864,"feed_emoji":"💥","tokens_out":9417,"duration_ms":90008,"temperature":0.7,"pith_summary":"The paper proposes a new way for heavy elements to be made in long gamma-ray bursts: instead of relying on pre-existing neutrons, the burst's own photon flood creates them. At the jet head, where the relativistic jet plows into stellar material, the estimated flux of roughly $10^{33}$ photons per square centimeter per second per keV at 1 MeV makes the reaction gamma + p -> pi+ + n proceed almost instantly, converting even pure-proton stellar gas into neutron-rich matter. Because neutrons are uncharged, they escape the magnetically confined jet head into the surrounding cocoon, where neutron-capture nucleosynthesis can proceed. The paper shows that the resulting cocoon conditions can produce full r-process patterns, including actinides, as well as intermediate neutron-capture patterns, and that the combined outputs can match solar r-process residuals. If this holds, collapsar GRBs become a new class of r-process site that does not require a neutron-star merger.","feed_headline":"Gamma-ray burst jets can turn pure protons into a neutron flood","feed_subtitle":"If correct, collapsar long GRBs become a new source of r-process elements, alongside neutron-star mergers.","key_machinery":"The load-bearing mechanism is the photo-pion transmutation reaction gamma + p -> pi+ + n, balanced against gamma + n -> pi- + p, with cross sections taken from the paper's coupled-channel photon-nucleon reaction model, a unified fit to world photoproduction data. The photon flux is modeled as a piecewise power law of the observed GRB spectral form, normalized to A = $10^{37}$ ph/$cm^{2}$/s/keV at the pivot energy and scaled to the jet head. A single-zone reaction network tracks protons, neutrons, escaped neutrons, and charged pions, competing the neutron-production rate against neutron decay, pion decay, and the ~$10^{-4}$ s neutron escape timescale while protons and pions remain magnetically confined. This interplay sets the electron fraction of the cocoon, the parameter that ultimately controls whether the outflow produces an r-process, an i-process-like pattern, or proton-rich nucleosynthesis.","core_discovery":"This paper's central claim is that hadronic photoproduction can manufacture neutrons in situ in collapsar gamma-ray bursts. Using a coupled-channel photon-nucleon model to compute cross sections and a single-zone reaction network at the jet head, it finds that with the estimated photon flux Phi_gamma ~ $10^{33}$ ph/$cm^{2}$/s/keV at ~$10^{6}$ keV, neutrons and pions are produced nearly instantaneously from a starting composition of pure protons. The charged products stay confined by the jet's magnetic field, while neutrons escape on ~$10^{-4}$ s timescales into the cocoon; the resulting electron fractions can be extremely low, down to Ye ~ 0.0035 in one simulated case. Nucleosynthesis calculations then yield weak and robust r-process patterns, fission recycling with actinide production, and an i-process-like cold pattern, whose summed abundance curve reproduces the full range of solar r-process residuals. The conclusion is that collapsar long GRBs are a viable new site for neutron-capture nucleosynthesis.","pith_inferences":["If this is right, the standard division of r-process labor between short GRBs and compact-object mergers is incomplete; galactic chemical evolution models would need a long-GRB channel whose yield is sensitive to jet Lorentz factor and envelope density.","Because the same photons that create neutrons also create pions that decay into neutrinos, a stacking search for high-energy neutrinos from long GRBs could independently confirm the jet-head photon flux before any nucleosynthesis signature is observed.","The bottleneck is the sub-photospheric non-thermal tail; the mechanism could be tested by radiative-transfer simulations of jet-head dissipation that compute the actual 1 MeV photon flux in the envelope frame, or by measuring the prompt spectrum's high-energy extension across many bursts."],"forward_implications":["Collapsar long GRBs become candidate r-process sites that do not require pre-existing neutrons or a neutron-star merger; even a helium or carbon envelope works because photodissociation first frees nucleons.","The cocoon can reach electron fractions far below those of disk outflows, down to Ye ~ 0.03 for Gamma ~ 10 and even lower in the paper's longer-lived scenario, which is the regime that produces actinides and fission recycling.","The combined abundance patterns from the paper's three density and temperature evolutions can reproduce the solar r-process residuals across the first, second, third, and lead peaks, giving a specific elemental fingerprint to search for in GRB-associated ejecta.","Observable signatures follow from the same reactions: red kilonova light from lanthanide and actinide production, the 2.6 MeV 208Tl decay line, a possible neutron precursor, and high-energy neutrinos from pion decay."],"supporting_citations":[{"why":"Supplies the observed sub-photospheric non-thermal spectrum and the 10^6 keV flux that is scaled back to estimate the jet-head photon flux.","marker":"Ahlgren et al. (2019)"},{"why":"Provides the finalized coupled-channel reaction amplitudes used to compute the photoproduction cross sections.","marker":"Kamano et al. (2013)"},{"why":"Gives the methodology and partial-wave amplitudes from which the paper computes the cross sections for the hadronic processes.","marker":"Kamano et al. (2019)"},{"why":"Establishes the direct, resonant, and multi-pion photomeson reaction set that the neutron-production channels belong to.","marker":"Mücke et al. (1999)"},{"why":"Provides the observed GRB spectral form used for the piecewise power-law photon flux normalization.","marker":"Band et al. (1993)"},{"why":"Supplies the standard GRB jet phases, radii, and baryon-loading constraints that set the jet head environment.","marker":"Piran (2004)"},{"why":"Provides envelope density profiles and cocoon simulations used to set baryon densities and the inverse-Compton contribution.","marker":"De Colle et al. (2022)"},{"why":"Supplies the reaction network used to simulate neutron-rich nucleosynthesis from the model's electron fractions.","marker":"Sprouse et al. (2021)"},{"why":"Defines the collapsar scenario that places the jet inside a massive stellar envelope.","marker":"Woosley (1993)"}],"fun_headline_variants":["Neutron flood from pure protons in gamma-ray burst jets","Collapsar jets produce neutrons for r-process nucleosynthesis","Photon-rich jets churn out neutrons from protons","Pure proton jets yield neutron seas for heavy elements"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands on the assumed photon flux at the jet head, estimated by scaling an observed GRB spectrum from cosmological distance back to r ~ $10^{10}$ cm; if that flux is even a few orders of magnitude lower than ~$10^{33}$ ph/$cm^{2}$/s/keV at 1 MeV, neutron production cannot keep pace with neutron escape.","fun_headline_variants_meta":{"raw":{"variants":["Neutron flood from pure protons in gamma-ray burst jets","Collapsar jets produce neutrons for r-process nucleosynthesis","Photon-rich jets churn out neutrons from protons","Pure proton jets yield neutron seas for heavy elements"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000194,"raw_usage":{"total_tokens":1330,"prompt_tokens":900,"completion_tokens":430,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":516,"completion_tokens_details":{"reasoning_tokens":364}},"tokens_in":516,"tokens_out":430,"duration_ms":4672,"temperature":1.0,"reasoning_tokens":364,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:05:25.878337+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the actual photon flux at the jet-head interaction radius with a self-consistent sub-photospheric dissipation model; if the flux at $10^{6}$ keV falls below ~$10^{25}$ ph/$cm^{2}$/s/keV, neutron production is slower than the ~$10^{-4}$ s escape timescale and the claimed mechanism fails.","supporting_citations":[],"review_version":1}