{"id":"b4810060-bf0b-47b5-9f30-ac2153bba3c7","arxiv_id":"2501.07134","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A one-zone radioactive heating model with macroscopic nickel-oxygen mixing reproduces the [O III] luminosity of SN 2018ibb and yields a scaling relation for line strength among PISNe.","lead":"This paper proposes that the strong [O III] emission lines seen in the superluminous supernova 2018ibb can be powered by radioactive cobalt decay inside oxygen-rich ejecta, rather than by collision with surrounding gas. A simple one-zone model with two adjustable parameters reproduces the observed line luminosity and ratios, and a scaling relation explains why such lines vary among pair-instability supernovae.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model assumes oxygen absorbs a mass-weighted share of the bolometric luminosity, but in a clumpy two-component ejecta the gamma-ray energy deposited in oxygen is set by covering factor and optical depth, not by MO/M.","rationale":"The reader's weakest assumption correctly identifies macroscopic mixing and the density contrast χ as free, assumed ingredients. The stress-test sharpens this into a specific mechanism: the deposition formula Ld(O) = (MO/M)L bypasses the actual gamma-ray transport in a clumpy two-component medium. The condition δ < 0.5 addresses inter-clump transparency but not the probability that gamma rays encounter oxygen clumps, which is controlled by the filling factor and clump optical depth. Because f = (MO/M)/χ, the same χ that raises the local ionization rate also reduces the geometric covering fraction, so the net effect on oxygen heating is not captured by the current one-zone treatment. This is load-bearing because the absolute [O III] luminosity is directly proportional to Ld(O); if the true deposition is lower by a factor of χ, the model's central match fails. The proposed test is concrete and would settle whether the concern lands. The paper is otherwise a plausible and clearly written demonstration of a possible mechanism, and the scaling relation in Section 4 is a useful independent contribution. Since the reader verdict was already CONDITIONAL, this concern does not change the verdict; it reinforces the need for a transport-based deposition calculation before the radioactive mechanism is confirmed.","tokens_in":6739,"tokens_out":8359,"duration_ms":93040,"concrete_test":"Run a simple Monte Carlo or analytic two-component gamma-ray deposition calculation: place 30 M⊙ of 56Ni in clumps and 15 M⊙ of oxygen in clumps with filling factor f = (MO/M)/χ for χ = 3 and inter-clump velocity separation δ in [0.1, 0.5]; compute the absorbed power in oxygen from the 56Co gamma-ray lines at t ≈ 381 days. If the absorbed power falls below (MO/M)L by more than 30%, the model's energy normalization is not self-consistent and the reproduced [O III] luminosity is an artifact of the assumed mass-proportional energy partition.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Section 2 the deposited power in oxygen is set to Ld(O) = (MO/M)L without a radiative-transfer calculation. The mixing condition Δτγ = 2.2δ < 1 only ensures that the material between nickel and oxygen fragments is transparent; it does not ensure that a 56Co gamma ray actually intersects an oxygen clump and is absorbed there. With oxygen filling factor f = (MO/M)/χ, the covering fraction of oxygen as seen from nickel clumps is roughly f for optically thick clumps, so the absorbed fraction can scale as 1/χ relative to the mass fraction. For the fiducial case MO = 15 M⊙, M = 60 M⊙, and χ = 3, this would reduce Ld(O) by a factor of about 3, from 3×10^42 to ~10^42 erg/s, breaking the [O III] luminosity match. The density contrast χ is used to boost the local deposition rate and ionization, but the same χ lowers the oxygen filling factor and therefore the gamma-ray covering factor; the paper does not compute this coupling. Thus the central luminosity agreement is not an independent test of the mixing hypothesis; it is an assumed energy partition. The acknowledged one-zone limitations for line profiles are secondary; the normalization problem is more fundamental.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"SN 2018ibb is a candidate pair-instability supernova whose nebular spectrum shows strong [O III] 5007, 4959, and 4363 Å emission of debated origin. The paper proposes a one-zone radioactive-deposition model: the ejecta are a homogeneous-on-average sphere of 60 Msun and 1.2e52 erg, with 30 Msun of 56Ni macroscopically mixed into oxygen, and the deposited gamma-ray power in oxygen is assumed to be Ld(O) = (MO/M)L. Oxygen ionization is solved in steady state with free electron temperature Te and density contrast chi. The model reproduces the observed [O III] luminosity and line ratios for MO = 15 Msun with Te about 9000 K and chi about 3, and similar solutions for MO = 10 and 5 Msun. Section 4 derives a scaling y3 proportional to Mni^{3/2} M^{-7/2} to explain the range of [O III]/[O I] ratios among PISNe. The paper explicitly acknowledges that the one-zone model cannot reproduce the [O I] 5577 Å nondetection or the broader [O I] profile.","tokens_in":6972,"tokens_out":8366,"duration_ms":87441,"significance":"If the central claim is correct, the strong [O III] emission of SN 2018ibb would be powered by 56Co decay rather than by circumstellar interaction, and the simple scaling relation would provide a testable explanation for why only a few PISNe show strong [O III]. The model is transparent, the parameter set is small, and the scaling relation is a falsifiable prediction. The manuscript also honestly lists its limitations. However, the quantitative support for the mechanism rests on an assumed energy partition and on Te and chi being fitted to the observations, so the evidence presented is weaker than the concluding claim of confirmation.","major_comments":[{"comment":"The relation Ld(O) = (MO/M)L assumes that oxygen absorbs a mass-weighted share of the total bolometric luminosity, but this is not derived from the clumpy geometry introduced in the same section. With the oxygen filling factor f = (MO/M)chi^{-1} and chi = 3, the covering fraction of optically thick oxygen clumps as seen from 56Co clumps is about f, so the absorbed fraction can be a factor of about chi lower than MO/M. The paper neither computes the coupling between chi and this covering factor nor supplies the clump-size and optical-depth argument needed to justify the mass-weighted partition. Because the [O III] luminosity normalization is the central test of the radioactive mechanism, this missing step is load-bearing.","section":"Section 2"},{"comment":"The agreement with the observed [O III] luminosity and the 4363/5007 ratio is achieved by choosing Te and chi; for each assumed MO the parameters are adjusted to pass through the observed point, so the agreement is a two-parameter fit rather than an independent confirmation. The comparison eta_h approximately eta_h,num in Table 1 is a useful internal consistency check, but it uses the same Ld(O) and the same one-zone ionization balance, so it does not by itself establish the radioactive origin of the emission.","section":"Section 3, Figure 1, Table 1"},{"comment":"The acknowledged failures on [O I] 5577 Å and the [O I] line profile are not merely cosmetic. The model predicts a strong 5577 Å line at +276 d that is not observed, and the observed [O I] doublet is broader than [O III], implying different emission zones. The paper defers these issues to a two-zone model; until that model is specified, the reproduced [O III]/[O I] ratio in Figure 1b cannot be regarded as a physical prediction for a single oxygen component, and the concluding sentence in Section 5 is stronger than the demonstrated result.","section":"Sections 3 and 4"}],"minor_comments":[{"comment":"The wavelength of the [O III] auroral line is given as 4363 Å in the Introduction but as 4359 Å in Figure 1 and Section 3; one consistent value should be used.","section":"Introduction and Section 3"},{"comment":"There are numerous typographical errors, e.g., 'deposion', 'fsctor', 'ioniztion', 'distibution', and 'radiactive'; the manuscript needs a careful proofreading pass.","section":"Throughout"},{"comment":"The filling factor f = (MO/M)chi^{-1} is introduced but then plays no explicit role in the deposition or line-luminosity equations; stating where f enters would make the cost of a high density contrast visible.","section":"Section 2"},{"comment":"The extension of w = 2.3I to O II and O III is stated in one sentence; a brief justification or a reference for the ion case would help the reader assess the ionization balance.","section":"Section 2"},{"comment":"The scaling Mni proportional to E^3 is taken from helium-core PISN models, whereas SN 2018ibb is modeled as an oxygen-core explosion; the paper should state whether the same proportionality is expected for oxygen-core ejecta.","section":"Section 4"},{"comment":"The caption says 'Colored lines show the dependence on the chi parameter for the fixed temperature' and identifies red as 7000 K and green as 10000 K, but the black lines showing the temperature dependence at fixed chi should also be described explicitly in the caption, as they are in the text.","section":"Figure 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal, and the author is admirably explicit about the model's limitations. The main risk is the deposition partition assumption; if the author can add a simple covering-factor estimate or a two-zone demonstration, the conclusion would be considerably strengthened. I see no concerns about citation practice or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short paper, single author, decent physics. The new thing is applying the well-known radioactive heating mechanism to SN 2018ibb with macroscopic mixing, and deriving a scaling relation for [O III] strength among PISNe. The one-zone model fits the observed [O III] luminosity and the 4363/5007 ratio with Te≈9000 K and χ≈3 for 15 M⊙ oxygen, and similar values for 10 and 5 M⊙. The ionization calculation is standard, and the comparison of ηh with the Spencer-Fano number is a nice cross-check. The scaling relation y3 ∝ Mni^{3/2} M^{-7/2} is derived cleanly and is independent of the fitting; it explains naturally why a low-mass, high-56Ni PISN like 2018ibb would show strong [O III] while a more massive He-core PISN would not. That is the most useful part of the paper.\n\nThe soft spots are real but the paper is honest about them. The central reproduction of the [O III] luminosity is a fit to Te and χ; those are the two free parameters, so the agreement in Figure 1 is not a prediction. The model fails on the [O I] 5577 line (predicts it too strong at +276 d) and on the broader [O I] profile; the author says a two-zone model with a cool outer oxygen region is needed, which is a reasonable way out but not demonstrated here. The stress-test note about the energy partition Ld(O)=(MO/M)L also has merit: in a clumpy geometry with χ>1, the gamma-ray absorption in oxygen is not guaranteed to be mass-weighted, and the paper does not compute the covering factor. I would not call it fatal—the clump optical depths are order unity and χ is adjustable—but it adds another layer of freedom to the 'agreement'.\n\nThe scaling relation and the explicit call for macroscopic mixing are worth having in the literature. The paper is clearly written, the equations check out, and the author flags the limitations instead of hiding them. A serious referee can engage with whether the two-zone model would actually produce the observed line profiles. I would send it to review. For a reading group, it is a quick, discussable paper on a topical object.","headline":"A short, honest one-zone model makes a plausible case that 56Co decay, given macroscopic mixing, can power the [O III] emission in SN 2018ibb, but the central luminosity match is a fit to Te and chi, and the two-zone discrepancies are acknowledged; the scaling relation is the cleanest result.","tokens_in":7536,"tokens_out":6733,"would_cite":true,"duration_ms":68418,"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":"With macroscopic mixing of nickel and oxygen, $^{56}$Co decay can power the strong [O III] lines of SN 2018ibb, so circumstellar interaction is not required.","keywords":["SN 2018ibb","pair-instability supernovae","[O III] emission","radioactive energy deposition","56Co decay","macroscopic mixing","nebular phase","superluminous supernovae"],"falsifier":"Monitor the [O III] 5007 luminosity at late times, say from day 500 to 900, and compare its decline with the $^{56}$Co decay input: the radioactive model predicts the line should follow that radioactive input, while a circumstellar-interaction origin would give a different light curve. In parallel, a three-dimensional explosion model without imposed macroscopic mixing and with $\\delta > 0.5$ would yield a deposited oxygen power far below $10^{41}$ erg s$^{-1}$, directly falsifying the mixing assumption.","tokens_in":6477,"feed_emoji":"💥","tokens_out":9937,"duration_ms":91610,"temperature":0.7,"pith_summary":"The paper asks whether the strong [O III] 5007, 4959, and 4363 Å emission of the pair-instability supernova SN 2018ibb, at roughly $10^{41}$ erg s$^{-1}$, can be produced by radioactive heating alone. Its answer is yes, provided the ejecta contain macroscopic mixing between $^{56}$Ni and oxygen: gamma-rays from $^{56}$Co decay deposit energy directly into oxygen matter, ionizing and exciting O III. A one-zone model with realistic ejecta mass, energy, and $^{56}$Ni mass reproduces the observed line luminosity and temperature-sensitive flux ratios for an oxygen mass of 10–15 $M_\\odot$, a temperature near 9000 K, and a density contrast $\\chi \\approx 2$–3. The paper also derives a scaling relation $y_3 \\propto M_{\\mathrm{Ni}}^{3/2} M^{-7/2}$ that explains why the [O III]/[O I] ratio should vary widely among pair-instability supernovae.","feed_headline":"Cobalt decay can power SN 2018ibb's [O III] lines","feed_subtitle":"A simple one-zone model reproduces the oxygen-line luminosity and explains why some pair-instability supernovae shine in [O III].","key_machinery":"The central object is a one-zone, free-expansion model of the ejecta with mass $M = 60\\,M_\\odot$, energy $E = 1.2 \\times 10^{52}$ erg, and $M_{\\mathrm{Ni}} = 30\\,M_\\odot$, in which oxygen occupies a filling factor $f = (M_O/M)\\chi^{-1}$ with density contrast $\\chi$. The key condition is $\\Delta\\tau_\\gamma = \\delta k_\\gamma \\rho v_0 t \\approx 2.2\\delta < 1$, meaning gamma-rays from $^{56}$Co cross oxygen without strong absorption, so the deposited power in oxygen is $L_d(O) = (M_O/M)L$. Ionization balance is solved in steady state, and line emissivities come from three-level balance equations with radiative and collisional transitions. The argument is carried by the scaling $y_3 \\propto M_{\\mathrm{Ni}}^{3/2} M^{-7/2}$, derived from the O II–O III ionization balance combined with the relation $M_{\\mathrm{Ni}} \\propto E^3$; this scaling translates ejecta mass and nickel mass into the O III fraction and therefore into the [O III]/[O I] line ratio.","core_discovery":"The central claim, stated in Section 5, is that observed [O III] emission lines can originate from deposition of the radioactive energy of $^{56}$Co decay into oxygen matter, given macroscopic mixing between $^{56}$Co and oxygen. In the model, oxygen receives a fraction $M_O/M$ of the total radioactive luminosity, and gamma-ray absorption is kept weak by the mixing condition $\\Delta\\tau_\\gamma = 2.2\\delta < 1$ (velocity separation $\\delta < 0.5$); a density contrast $\\chi$ then boosts the volumetric deposition. With $M_O = 15\\,M_\\odot$, $T_e \\approx 9000$ K, and $\\chi \\approx 3$, the computed [O III] doublet luminosity, the 4363/5007 ratio, and the [O III]/[O I] ratio agree with observations, and the thermalization fraction matches a detailed numerical deposition calculation. The favored high oxygen mass, together with the scaling $y_3 \\propto M_{\\mathrm{Ni}}^{3/2} M^{-7/2}$, leads the author to conclude that the radioactive mechanism, not circumstellar interaction, powers the observed [O III] emission.","pith_inferences":["A testable extension the paper does not perform: track the [O III] 5007 luminosity after day 400; the radioactive model predicts it should decline on the $^{56}$Co decay timescale set by the bolometric input, whereas a circumstellar-interaction origin would give a different decline.","The same mixing condition and scaling could be used as a classification tool: a measured [O III]/[O I] ratio combined with an ejecta-mass estimate would constrain the nickel mass and the degree of mixing in other superluminous supernovae.","If three-dimensional PISN explosion simulations resolve macroscopic mixing, they can check the $\\delta < 0.5$ requirement directly; the model implies that mixing must be present at that level or the [O III] line would be far fainter.","The two-zone extension suggested by the paper implies that late-time spectra with better spatial or velocity resolution should show [O III] concentrated at lower velocities than [O I]."],"forward_implications":["If the model is right, the [O III] luminosity of SN 2018ibb is powered by $^{56}$Co decay in oxygen, so no massive circumstellar shell or pulsational pair-instability interaction is needed to explain the lines.","The inferred oxygen mass of 10–15 $M_\\odot$ and temperature near 9000 K place the [O III]-emitting zone in a hot inner region distinct from the cooler outer zone that emits [O I].","The scaling $y_3 \\propto M_{\\mathrm{Ni}}^{3/2} M^{-7/2}$ predicts that pair-instability supernovae with high nickel mass and low ejecta mass show strong [O III]/[O I] ratios, while those with helium shells and lower nickel mass show ratios tens of times smaller.","The agreement between the modeled thermalization fraction and detailed numerical deposition calculations supports the radioactive deposition picture for oxygen masses near 15 $M_\\odot$.","The discrepancy with earlier one-dimensional PISN spectral models is attributed to the lack of macroscopic $^{56}$Ni–oxygen mixing in those models, identifying mixing as a necessary ingredient for radioactivity-powered [O III] emission."],"supporting_citations":[{"why":"Supplies the observed [O III] luminosity, line ratios, bolometric luminosity, and ejecta parameters for SN 2018ibb that the model must reproduce.","marker":"Schulze et al. 2024"},{"why":"Provides the mod60 reference model (60 $M_\\odot$, $1.2 \\times 10^{52}$ erg, 30 $M_\\odot$ of $^{56}$Ni) and the argument that SN 2018ibb exploded from an oxygen core.","marker":"Chugai 2024"},{"why":"Gives the gamma-ray absorption coefficient $k_\\gamma = 0.03$ cm$^2$ g$^{-1}$ used in the mixing condition $\\Delta\\tau_\\gamma = 2.2\\delta$.","marker":"Sutherland & Wheeler 1984"},{"why":"Provides the numerical thermalization fraction $\\eta_{h,\\mathrm{num}}$ used to check the computed deposition efficiency.","marker":"Kozma & Fransson 1992"},{"why":"Supplies the relation $M_{\\mathrm{Ni}} \\propto E^3$ used to derive the scaling $y_3 \\propto M_{\\mathrm{Ni}}^{3/2} M^{-7/2}$.","marker":"Kasen et al. 2011"},{"why":"Provides the average work per ionization $w = 2.3I$ used in the ionization balance.","marker":"Ahlen 1980"},{"why":"Supplies the radiative recombination coefficients used in the steady-state oxygen ionization equations.","marker":"Tarter 1971"}],"fun_headline_variants":["Radioactive cobalt explains SN 2018ibb's oxygen lines","Cobalt decay powers [O III] in pair-instability supernova","Simple model ties cobalt decay to SN 2018ibb's [O III]","SN 2018ibb's O III lines from 56Co decay, model says","Pair-instability SN oxygen lines traced to cobalt decay"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model stands on the assumption that $^{56}$Ni and oxygen are macroscopically mixed so that the gamma-ray optical depth between them satisfies $\\Delta\\tau_\\gamma \\approx 2.2\\delta < 1$ (velocity separation $\\delta < 0.5$), and that a density contrast $\\chi \\approx 2$\\u20133 increases deposition; if real ejecta are not mixed this way, the power deposited in oxygen is too small to produce the observed [O III] luminosity.","fun_headline_variants_meta":{"raw":{"variants":["Radioactive cobalt explains SN 2018ibb's oxygen lines","Cobalt decay powers [O III] in pair-instability supernova","Simple model ties cobalt decay to SN 2018ibb's [O III]","SN 2018ibb's O III lines from 56Co decay, model says","Pair-instability SN oxygen lines traced to cobalt decay"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000999,"raw_usage":{"total_tokens":4201,"prompt_tokens":893,"completion_tokens":3308,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":509,"completion_tokens_details":{"reasoning_tokens":3213}},"tokens_in":509,"tokens_out":3308,"duration_ms":23649,"temperature":1.0,"reasoning_tokens":3213,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:48:48.700209+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Monitor the [O III] 5007 luminosity at late times, say from day 500 to 900, and compare its decline with the $^{56}$Co decay input: the radioactive model predicts the line should follow that radioactive input, while a circumstellar-interaction origin would give a different light curve. In parallel, a three-dimensional explosion model without imposed macroscopic mixing and with $\\delta > 0.5$ would yield a deposited oxygen power far below $10^{41}$ erg s$^{-1}$, directly falsifying the mixing assumption.","supporting_citations":[{"cited_title":"The power deposited in the oxy gen matter by gamma-quanta of 56Co decay is Ld(O) = ( MO/M)L","cited_arxiv_id":null,"evidence_quote":"Provides the mod60 reference model (60 $M_\\odot$, $1.2 \\times 10^{52}$ erg, 30 $M_\\odot$ of $^{56}$Ni) and the argument that SN 2018ibb exploded from an oxygen core."}],"review_version":1}