{"id":"9f29d683-fbb7-491b-bf3e-72b16a90cca2","arxiv_id":"2505.22992","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"If all positrons form positronium before annihilation in Type Ia supernova ejecta, the 511 keV gamma-ray line flux decreases by about 70 percent and energy deposition increases by up to 2 percent near 100 days.","lead":"TARDIS-HE, a new open-source gamma-ray transport module for the TARDIS supernova code, shows that if positronium forms during positron annihilation, the 511 keV line flux in Type Ia supernova models drops by about 70 percent while energy deposition rises by about 2 percent. The study provides a tool and sensitivity estimate for interpreting future gamma-ray observations of nearby supernovae.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The §2.3/Appendix A sampling of the 3γ continuum uses the Ore-Powell number distribution F(x) for fixed-energy packets; the emergent spectrum should follow xF(x), so the simulated continuum is biased soft, likely inflating the 2% energy-deposition enhancement and affecting the 70% line reduction.","rationale":"I identify the 3γ continuum sampling as the most load-bearing concern because the paper's headline numbers are computed with a fixed f_p and 75% branching—assumptions the paper explicitly labels as extreme and explores as parameters—so the physical realism of f_p is not the central claim. However, if the Monte Carlo emission sampling is inconsistent with the equal-energy packet scheme, the quantitative results are wrong even under the stated assumptions. Section 2.2 sets equal packet energies, and Section 2.3/Appendix A samples the 3γ photon energy from F(x). In Lucy's equal-energy packet scheme, the packet frequency should be sampled from the energy distribution xF(x). The resulting bias to lower photon energies directly impacts the claimed 2% deposition enhancement, because the enhancement is driven by photoabsorption of the low-energy continuum. The 70% line reduction could also shift because the continuum level near 511 keV is misestimated. This is a concrete, testable implementation issue rather than a matter of physical modeling choices. The reader's weakest assumption about fixed f_p is valid but explicitly acknowledged in Section 5 and framed as a parameter study; the sampling issue is not acknowledged. Therefore, I recommend retaining the CONDITIONAL verdict, with the condition being a verification or correction of the 3γ sampling. My agreement with the reader is partial: their concern is complementary but not the most load-bearing; the internal sampling inconsistency is more fundamental.","tokens_in":19168,"tokens_out":28779,"duration_ms":290722,"concrete_test":"Set up a TARDIS-HE run with f_p=1 in an optically thin, static, homogeneous shell with the 3γ channel only and no Compton/photoabsorption interactions, and compare the binned emergent energy spectrum to the analytic xF(x) and F(x) curves from Eq. A1. If the simulation matches F(x), resample packet frequencies from the normalized xF(x) distribution, rerun the ddt_N100 day-100 model, and check whether the 70% line reduction and the 2% deposition enhancement change by more than ~1σ Monte Carlo noise.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 2.2 fixes the comoving packet energy to E_total/N_packet (Eq. 7), so TARDIS-HE uses equal-energy indivisible packets in the Lucy (2005) scheme. For such packets, the emission frequency must be sampled from the energy spectrum of the source, not the photon number spectrum. Appendix A samples the 3γ photon energy from the Ore-Powell distribution F(x) in Eq. A1, which is the number distribution of photons per unit x. Because each packet carries equal energy, the simulated energy spectrum becomes proportional to F(x), whereas the true 3γ continuum energy spectrum is proportional to xF(x). Since x ∈ (0,1], this over-represents low-energy photons relative to high-energy photons. The paper's 2% deposition increase is attributed to the low-energy continuum being more readily photoabsorbed (Section 4.1), so this bias likely overestimates that increase. It also modifies the continuum level just below 511 keV, which enters the measured line flux after continuum subtraction, so the 70% reduction could shift. This is an internal sampling inconsistency, independent of the (explicitly parametric) choice of f_p and 75% branching.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents TARDIS-HE, a new time-dependent Monte Carlo transport module for gamma rays built on the TARDIS radiative-transfer framework. The code uses indivisible energy packets, homologous expansion, and three opacity channels (Compton scattering, pair production, photoabsorption), and it reads decay chains and transition data from ENSDF to produce spectra and deposition curves. Positronium formation is included in a parameterized way: a fraction f_p of annihilation packets is assumed to form positronium, and 75% of those are assumed to decay through the three-photon channel. The authors validate the spectra against Summa et al. (2013) and the deposition against several published codes, then apply the code to four SN Ia models. The central quantitative result is that for f_p = 1 and a 75% three-photon branching, the 511 keV line flux is reduced by approximately 70% at around 95-100 days and the energy deposition is increased by up to 2%, compared with models without positronium.","tokens_in":19483,"tokens_out":10763,"duration_ms":134258,"significance":"If the implementation is correct, TARDIS-HE would be a useful open-source addition to the SN Ia gamma-ray modeling toolkit, and the paper gives a transparent parameter study of positronium's effect on line and continuum observables, including a discussion of COSI detectability. The code description is detailed enough that the emission and opacity treatments can be checked, and the comparison against multiple deposition codes is a strength. However, the headline numbers are directly tied to the assumed constant values of f_p and the three-gamma branching, and the three-gamma sampling appears inconsistent with the equal-energy indivisible-packet scheme, so the claimed 2% deposition enhancement is not yet established. The spectral validation is also weakened by a model-mass mismatch. These issues are fixable with a corrected sampling scheme and a same-model comparison, but they must be addressed before the quantitative conclusions can be accepted.","major_comments":[{"comment":"The three-gamma continuum sampling is inconsistent with the equal-energy indivisible-packet scheme. Equation (7) fixes every packet's rest-frame energy to E_tot/N_packet, so in such a scheme the emission frequency must be drawn from the source's energy spectrum. For positronium three-gamma decay, the Ore-Powell function F(x) in Eq. (A1) is the number distribution of photon energies; the physical continuum energy spectrum is proportional to xF(x). Sampling x from F(x) instead produces an emergent energy spectrum proportional to F(x), which is biased toward low photon energies. Because §4.1 attributes the 2% deposition increase to low-energy continuum photons being more readily photoabsorbed, this bias likely inflates that number, and it also shifts the continuum level subtracted under the 511 keV line. In addition, the manuscript does not state how many packets are created per three-gamma decay; since a three-photon decay has total energy 2 m_e c^2 while each equal-energy packet has energy E_tot/N_packet, the energy bookkeeping for this channel must be specified explicitly. The authors should sample from the energy-weighted distribution (or assign variable packet energies) and rerun the positronium cases.","section":"Appendix A; §2.2–2.3, Eq. (7)"},{"comment":"The spectral comparison against Summa et al. (2013) is not a clean validation because the two runs use different 56Ni masses: the text states that the ddt N100 model has 0.6 M_sun of 56Ni in Summa et al., while TARDIS-HE is run with 0.67 M_sun after normalization, and the late-phase flux difference is attributed to this mass difference. This means the comparison does not test whether the code reproduces the published absolute line and continuum fluxes. The authors should run the same 56Ni mass (or explicitly rescale the comparison) and quantify the residual difference.","section":"§3.1, Fig. 3"},{"comment":"The central quantitative claims — approximately 70% reduction in the 511 keV line flux and up to 2% increase in energy deposition — are reported without statistical or systematic uncertainties. With a finite packet count, the line and continuum fluxes carry Monte Carlo Poisson errors, and the specutils Gaussian-fit plus continuum-subtraction procedure adds a systematic component. The authors should report these uncertainties or, at minimum, demonstrate convergence with respect to packet number and show that the line-measurement method does not drive the quoted values.","section":"§4.1, Fig. 5"}],"minor_comments":[{"comment":"The text says positronium decays by two-gamma or three-gamma emission 'in the ratio of 1:4.5', which is inconsistent with the 75% three-gamma branching used throughout the paper; the statistical ratio from the 1:3 singlet-to-triplet spin weights is 1:3, and the lifetime ratio for ortho- to para-positronium is much larger than 4.5. Please check and correct this statement.","section":"§2.3"},{"comment":"The phrase 'We find that full positronium formation can reduce the 511 keV line flux by approximately 70%' should be framed as conditional on the assumed input f_p = 1 and the assumed 75% three-gamma branching, since these quantities are not derived in this work; the current wording could be read as a prediction rather than a parameter-study result.","section":"Abstract and §5"},{"comment":"There are several typos and formatting inconsistencies, including 'radiaoctive' (Introduction), 'desposition' (§2.6 and Fig. 1), 'contin num opacity' (§4.1), and inconsistent rendering of 'tardis-he' (sometimes italic, sometimes roman). A careful copy edit is needed.","section":"Throughout"},{"comment":"The statement that 'the ortho-Ps continuum is independent of the composition of the ejecta' is true only for the intrinsic emission spectrum; after transport through the ejecta the observed continuum is modulated by Compton scattering and photoabsorption, so the sentence should be qualified.","section":"§4.2, last paragraph"}],"recommendation":"major_revision","confidential_remarks":"The Appendix A sampling issue is the key technical concern: for equal-energy packets, drawing emission frequencies from the number distribution rather than the energy distribution biases the three-gamma continuum soft and likely changes both headline numbers. I would ask for a rerun with energy-weighted sampling and for a clean same-mass validation against Summa et al. The paper's scope as an upper-limit parameter study is acceptable if these technical issues are resolved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two quick things to know. The new TARDIS-HE module is a real community asset: time-dependent, three-dimensional, equal-energy-packet gamma-ray transport inside TARDIS, with a detailed method section, validation against Summa et al. spectra and the Blondin et al. deposition code comparison, and a convergence check. Applied to four modern SN Ia models, the positronium sensitivity study gives clean numbers: 511 keV line down by ~70% at f_p=1, deposition up by ~2% at ~100 days. The second thing: those numbers are conditional on inputs (f_p=1, 75% three-photon branching), and the authors say so. The 70% is basically the 75% branching ratio moderated by radiative transfer; the paper is a sensitivity study, not a prediction, and reads that way.\n\nWhat is actually good: the code is described well enough to reproduce, the validation is honest, and the 511/1238 keV line-ratio comparison across delayed detonation, merger, deflagration, and double detonation models is useful for COSI-era planning. Citation choices are appropriate: Ore & Powell, Leising & Clayton, Milne et al., Lucy (2005).\n\nNow the soft spots, in proportion. The one that matters: Appendix A samples the 3γ continuum from the Ore-Powell number distribution F(x), but in the Lucy equal-energy packet scheme the emission frequency should be sampled from the energy distribution xF(x). With fixed packet energy, sampling F(x) biases the simulated continuum soft. Since the 2% deposition increase is attributed to low-energy photoabsorption, that number is likely inflated, and the continuum shape near and below the line in Figure 5 is distorted. The 70% line reduction is largely protected because the bias at 511 keV is small, so the main conclusion survives, but the deposition number and low-energy continuum need to be redone. Also missing: no Monte Carlo uncertainties anywhere, and no code link in the text despite the open-source claim—both easy fixes.\n\nWho this is for: COSI observers, SN Ia modelers, and people building gamma-ray transport codes. It deserves a serious referee. I would insist on correcting the 3γ sampling (or convincing me that packets are weighted differently than Eq. 17 suggests), and on adding packet-count error bars. The qualitative picture—positronium reshapes the 511 keV line and adds a modest deposition bump—will survive those changes.","headline":"TARDIS-HE is a genuinely useful new gamma-ray transport tool and the positronium study deserves a serious referee, but the 3γ continuum sampling is biased soft and the 2% deposition number needs a recheck before anyone quotes it.","tokens_in":19991,"tokens_out":14402,"would_cite":true,"duration_ms":149734,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Positronium formation can cut the 511 keV line flux of Type Ia supernovae by about 70 percent around 100 days after explosion.","keywords":["Type Ia supernovae","positronium","511 keV annihilation line","gamma-ray spectra","Monte Carlo radiative transfer","energy deposition","TARDIS-HE","gamma-ray line ratios"],"falsifier":"Measure the 511 keV line and the continuum shape of a nearby Type Ia supernova around 100 days after explosion with a high-resolution gamma-ray spectrometer: if the line flux matches the no-positronium prediction or no ortho-positronium continuum appears, the assumed full positronium formation and fixed 75 percent three-photon branching are not realized in real ejecta.","tokens_in":19021,"feed_emoji":"💥","tokens_out":6290,"duration_ms":63247,"temperature":0.7,"pith_summary":"Type Ia supernovae are powered by radioactive decay, and their gamma-ray lines are readouts of the explosion; but one microphysical step is usually skipped: the positrons from 56Co decay can first bind with electrons into positronium, and when positronium decays, most of the time it emits three photons instead of two, smearing the sharp 511 keV annihilation line into a continuum. Using a new open-source Monte Carlo transport module coupled to the radiative transfer code TARDIS, the paper varies the positronium fraction from 0 to 100 percent and finds that full positronium formation cuts the 511 keV line flux by about 70 percent and increases deposited energy by up to 2 percent near 100 days. The paper argues that these are measurable, non-dominant but real effects that should be included when interpreting gamma-ray observations, and that missions such as the Compton Spectrometer and Imager could constrain positronium formation in nearby SNe Ia.","feed_headline":"Positronium cuts the 511 keV supernova line by ~70 percent","feed_subtitle":"A new gamma-ray transport code shows the effect also boosts heating and shifts line ratios, so precise models should include it.","key_machinery":"The central machinery is the indivisible energy-packet Monte Carlo transport in TARDIS-HE: gamma-ray packets representing radioactive decay radiation are moved through spherical, homologously expanding ejecta shells and interact through Compton scattering, pair production, and photoabsorption. The positronium channel is handled separately: a chosen positronium fraction $f_p$ diverts that share of annihilation packets into positronium, and 75 percent of those decay through the three-photon ortho-positronium channel, with photon energies sampled from the Ore–Powell distribution. That branching is what moves flux out of the 511 keV line into the continuum below it and shifts energy deposition.","core_discovery":"The central discovery is that including positronium formation in gamma-ray transport changes predicted observables in a specific, quantifiable way. With a positronium fraction of unity and 75 percent of positronium decaying through the three-photon channel, the 511 keV line flux falls by roughly 70 percent at about 95 days in the delayed-detonation model, while the continuum below 511 keV brightens; energy deposition rises by up to 2 percent at late times because the softer three-photon continuum is more likely to be photoabsorbed. The paper also reports that the 511 keV to 1238 keV line ratio decreases in all four explosion models considered when positronium forms, and that low-energy 56Ni lines at 158 and 270 keV can distinguish explosion scenarios, disappearing into the continuum for the violent merger geometry.","pith_inferences":["The 70 percent line reduction is an upper-end scenario attached to the paper's assumed 100 percent positronium formation and fixed 75 percent three-photon branching; a self-consistent model with a time- and density-dependent positronium fraction will likely place the effect somewhere between the paper's 0 and 70 percent endpoints.","The same line-versus-continuum redistribution should apply to other positron-producing transients, such as ejecta powered by 44Ti or 57Ni decay, where the neglected positronium channel would alter gamma-ray diagnostics similarly.","Because the ortho-positronium continuum is independent of ejecta composition while the Compton-scattered continuum is not, separating these two components in an observed spectrum could isolate the positronium contribution without requiring detailed abundance information.","A direct extension would replace the fixed 75 percent three-photon branch with a density- and temperature-dependent ortho-para conversion treatment and test whether the energy-deposition excess grows or shrinks at late times."],"forward_implications":["The 511 keV line is not a clean measure of the 56Co positron yield; interpreting its flux without correcting for positronium redistribution would bias estimates of the radioactive mass.","The 511 keV to 1238 keV line ratio decreases with positronium formation in all four explosion models studied, offering an observational diagnostic for the positronium fraction that is less sensitive to absolute flux calibration.","Low-energy 56Ni lines at 158 and 270 keV survive in delayed-detonation, deflagration, and double-detonation models but are lost in the continuum for the violent merger model, so gamma-ray spectra can discriminate between explosion scenarios.","Energy deposition rises by up to 2 percent near 100 days, which feeds into late-time bolometric light curves and the UV-optical-IR emission powered by radioactive decay.","Future observations of late-phase spectra of a nearby Type Ia supernova could detect the ortho-positronium continuum and constrain the positronium fraction, refining models of radioactive energy transport.","If the claims hold, gamma-ray line-flux measurements used as explosion diagnostics must account for positronium formation even when the effect is not the dominant shaping process."],"supporting_citations":[{"why":"Supplies the three-photon decay energy distribution used to sample the ortho-positronium continuum.","marker":"Ore & Powell 1949"},{"why":"Provides the two-photon versus three-photon branching and positronium lifetimes adopted for the decay treatment.","marker":"Crannell et al. 1976"},{"why":"Source of the 75 percent three-photon decay probability assumed for positronium.","marker":"Leising & Clayton 1987"},{"why":"Defines the gamma-ray interaction processes and the role of positronium in SNe Ia that the module implements.","marker":"Milne et al. 2004"},{"why":"Supplies the indivisible energy-packet Monte Carlo scheme on which TARDIS-HE transport is built.","marker":"Lucy 2005"},{"why":"Provides the comparison gamma-ray spectra and the delayed-detonation diagnostics the code is checked against.","marker":"Summa et al. 2013"},{"why":"Supplies the multi-code energy-deposition comparison that validates TARDIS-HE deposition.","marker":"Blondin et al. 2022"},{"why":"Provides the delayed-detonation ejecta model used for the positronium line-flux and energy-deposition results.","marker":"Seitenzahl et al. 2013"}],"fun_headline_variants":["Positronium cuts 511 keV line by 70% in SNe Ia","Supernova gamma-ray line dims 70% with positronium","Positronium alters supernova spectra and boosts heating","New code shows positronium reduces 511 keV flux by 70%","Type Ia supernovae: positronium reshapes gamma-ray spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes a fixed fraction of positron annihilations goes through positronium and that 75 percent of those positronium atoms decay into three photons, but in a real supernova both numbers depend on local density, temperature, and ionization, so the 70 percent line reduction is tied to those assumed values.","fun_headline_variants_meta":{"raw":{"variants":["Positronium cuts 511 keV line by 70% in SNe Ia","Supernova gamma-ray line dims 70% with positronium","Positronium alters supernova spectra and boosts heating","New code shows positronium reduces 511 keV flux by 70%","Type Ia supernovae: positronium reshapes gamma-ray spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000361,"raw_usage":{"total_tokens":2004,"prompt_tokens":1054,"completion_tokens":950,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":670,"completion_tokens_details":{"reasoning_tokens":853}},"tokens_in":670,"tokens_out":950,"duration_ms":9457,"temperature":1.0,"reasoning_tokens":853,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:55:51.646135+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 511 keV line and the continuum shape of a nearby Type Ia supernova around 100 days after explosion with a high-resolution gamma-ray spectrometer: if the line flux matches the no-positronium prediction or no ortho-positronium continuum appears, the assumed full positronium formation and fixed 75 percent three-photon branching are not realized in real ejecta.","supporting_citations":[{"cited_title":"J., Joyce, G., Ramaty, R., & Werntz, C","cited_arxiv_id":null,"evidence_quote":"Provides the two-photon versus three-photon branching and positronium lifetimes adopted for the decay treatment."},{"cited_title":"D., & Clayton, D","cited_arxiv_id":null,"evidence_quote":"Source of the 75 percent three-photon decay probability assumed for positronium."},{"cited_title":"2013, A&A, 554, A67, doi: 10.1051/0004-6361/201220972","cited_arxiv_id":null,"evidence_quote":"Provides the comparison gamma-ray spectra and the delayed-detonation diagnostics the code is checked against."}],"review_version":1}