{"id":"eeeeb44d-290a-4799-88a7-f5b92f551d49","arxiv_id":"2505.13825","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"TransFit solves the time-dependent radiative diffusion equation with flexible density and heating profiles and fits the light curves of SN 1993J and SN 2011kl.","lead":"This paper presents TransFit, a numerical model that quickly computes supernova light curves by tracking heat diffusion through expanding debris. It aims to combine the speed of simple analytic fits with the accuracy of full radiation simulations for the upcoming flood of survey data.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dimensionless source term in Eq. (13) omits the ε_heat0 amplitude factor, so the central PDE as written cannot produce the reported MNi dependence unless an unstated normalization is imposed.","rationale":"The reader's weakest assumption is the gray-opacity and diffusion-approximation validity, which is a legitimate but secondary concern about realism and parameter degeneracy. The more load-bearing issue is internal to the derivation: the dimensionless source term in Eq. (13) lacks the heating-amplitude factor that must arise from the non-dimensionalization. If the equation as written is used, the model cannot respond to MNi or engine luminosity, contradicting the paper's own parameter study and the fits to SN 1993J and SN 2011kl. If the factor is implicitly absorbed by choosing u0 = ρ0 tdiff ε_heat0, then ETh,in and the heating amplitude are no longer independent, yet the paper treats them as independent inputs. Either way, the published formulation is not self-consistent, and the central claim that TransFit solves the generalized energy-conservation equation is not supported by the equations as presented. This is a correctness risk that needs to be resolved by the authors, but it does not by itself demand rejection because it may be a typographical omission in the non-dimensionalization rather than a fundamental flaw. The reader's verdict of CONDITIONAL remains appropriate, now with an additional specific condition: correct or clarify the source-term normalization and demonstrate that the code's parameter dependencies match the physical equations. The missing Sedona validation and absence of code/uncertainties also remain conditions for full acceptance, as the reader noted.","tokens_in":36963,"tokens_out":9508,"duration_ms":83016,"concrete_test":"Re-derive Eq. (11) from Eq. (1) using Eqs. (4)-(10) without assuming u0 = ρ0 tdiff ε_heat0, and verify whether the factor ε_heat0 ρ0 tdiff/u0 appears in S. Then run the published algorithm for two models differing only in MNi (e.g., 0.1 and 0.5 M_sun) with all other parameters fixed; if the code literally implements Eq. (13), the predicted bolometric light curves must be identical, revealing the inconsistency. If the code supplies the amplitude factor elsewhere, the paper must state that normalization explicitly and justify treating ETh,in and MNi as independent in Tables 1-2 and Figures 3-7.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that TransFit solves the generalized energy-conservation equation with time-dependent diffusion and heating. The dimensionless reduction from Eq. (1) to Eq. (11) requires multiplying both sides by ρ0 tdiff/u0. Doing this yields a source term S = (ρ0 tdiff ε_heat0/u0) η_ej(x) ξ_heat(x) f_heat(y) (Rmax/R0). The paper's Eq. (13) defines S without the factor ε_heat0 ρ0 tdiff/u0. If Eq. (13) is taken literally, the heating amplitude ε_heat0 (set by MNi or engine luminosity) never enters the PDE; varying MNi would leave the dimensionless solution e(x,y) unchanged and hence Lbol unchanged, contradicting Figure 3. The factor can be hidden only by adopting the normalization u0 = ρ0 tdiff ε_heat0, which is stated only in the Arnett-comparison section (Section 4.1), not in the general derivation. That normalization would tie ETh,in to MNi through Eq. (15), whereas the paper treats ETh,in and MNi as independent inputs (e.g., Table 1 for SN 1993J and Figure 7). Thus either the published equation is not the correct reduction of the stated physics, or an unstated parameter coupling exists. This directly undercuts the claim that TransFit solves the stated energy-conservation problem, independent of the gray-opacity and gamma-ray-leakage concerns raised by the reader.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces TransFit, a one-dimensional time-dependent radiation-diffusion model for supernova and transient light curves. The model solves a dimensionless energy-conservation partial differential equation for the internal-energy density of homologously expanding ejecta, with radioactive or central-engine heating, using a Crank-Nicolson finite-difference scheme. The authors explore how the light-curve morphology depends on ejecta mass, kinetic energy, opacity, initial thermal energy, progenitor radius, density profile, nickel mixing, and magnetar parameters, compare the model with Arnett and one-zone semi-analytic models, and apply it to SN 1993J and SN 2011kl. The abstract and introduction also claim validation against Monte Carlo radiative transfer codes such as Sedona.","tokens_in":37303,"tokens_out":10087,"duration_ms":94510,"significance":"If the derivation and validation were sound, TransFit would fill a practical niche: a fast finite-difference solver that relaxes the self-similar temperature-profile assumption of Arnett-like models while remaining inexpensive enough for survey-scale fitting. The numerical scheme in Appendix B is a standard Crank-Nicolson discretization and is likely stable, and the parameter study reproduces several expected qualitative trends, including the effects of nickel mixing, gray opacity, ejecta mass, and magnetar spin-down. The comparison with Arnett and one-zone models usefully illustrates when Arnett's law fails. However, as written the central dimensionless equation is internally inconsistent, and the promised Sedona validation is absent, so the paper's main claims currently rest on unverified or incorrect derivations.","major_comments":[{"comment":"The dimensionless source term as printed omits the heating amplitude. Substituting Eqs. (7)-(10) into Eq. (1) yields a source term S = (rho0 tdiff epsilon_heat0 / u0) eta_ej(x) xi_heat(x) f_heat(y) (Rmax/R0), whereas Eq. (13) defines S without the prefactor rho0 tdiff epsilon_heat0 / u0. Taken literally, the partial differential equation no longer depends on epsilon_heat0 (and hence on MNi or engine luminosity), so varying MNi would leave e(x,y) and therefore Lbol unchanged, contradicting Figure 3. The normalization u0 = rho0 tdiff epsilon_heat0 is introduced only in §4.1 near Eq. (50), not in the general derivation, and adopting it makes ETh,in dependent on epsilon_heat0 through Eq. (15), which is inconsistent with Section 3.4 and Table 1, where ETh,in and MNi are treated as independent inputs. The authors must either include the explicit prefactor in Eq. (13) or state the normalization and its consequences for the free parameters; the current text is internally inconsistent.","section":"§2, Eq. (13)"},{"comment":"The abstract and Section 1 state that validation comparisons with sophisticated Monte Carlo radiative transfer simulations such as Sedona demonstrate excellent agreement across a variety of transient scenarios, but no Sedona comparison appears anywhere in the body, the appendices, or the figures. The actual comparisons are limited to the Arnett and one-zone analytic models and to two observed objects. Since this is a central validation claim for the framework, the authors need either to add genuine Sedona comparison runs or to substantially soften the wording so that the claim matches the content of the paper.","section":"Abstract and §1"},{"comment":"The paper describes these applications as fits, but no inference procedure is presented. Tables 1 and 2 give point values for Mej, MNi, EK, R0, ETh,in, Pi, and Bd, yet there is no likelihood function, no prior distributions, no search or MCMC algorithm, no convergence check, and no reported uncertainties or goodness-of-fit statistics. The text says that 'identical priors' are used without defining them. Moreover, because the model parameters are adjusted to reproduce the same light curves that are then shown as agreement, the comparison is at best an illustrative model match and cannot serve as validation. The authors should either describe an actual fitting methodology with uncertainties or explicitly present these as manually tuned illustrative models.","section":"§5, Tables 1 and 2"}],"minor_comments":[{"comment":"The integrand appears as 4πr^2 a u dr, but u already denotes the internal-energy density; the factor a is likely a typographical artifact and should be removed.","section":"Eq. (15)"},{"comment":"The optical depth is written as tau0 = kappa0 rho0 R0 I_tau, with kappa0, while the opacity elsewhere is denoted kappa; the notation should be made consistent.","section":"Eq. (22)"},{"comment":"The legend does not display the xheat values used in the mixing-radius comparison, which makes the qualitative trend difficult to evaluate directly from the figure.","section":"Figure 2"},{"comment":"It is unclear whether central-engine heating enters through the source term S(x,y) or through the inner boundary condition f_ib(y) of Eq. (27); the text says the source term is 'incorporated into this inner boundary condition,' but Section 3.1 also motivates a radial heating profile for centrally concentrated heating. The authors should clarify how the two descriptions are reconciled to avoid double counting.","section":"§2.2 and §3.1"},{"comment":"The introduction promises flexibility for 'compositionally dependent opacity variations,' but Section 3.3 implements a single, time-independent gray opacity; this should either be clarified as a future extension or the introductory claim should be adjusted.","section":"§1 and §3.3"},{"comment":"There are several small typographical errors, including 'the the' in the captions of Figures 4 and 5 and 'bolometric ligtcurve' in Section 5.2; these should be corrected in a revision.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The missing prefactor in Eq. (13) is the most serious issue; it suggests the numerical code may use a normalization different from the one stated in the paper, and the discrepancy must be resolved before the model can be assessed. The absent Sedona comparison is also a central unmet promise. If both are fixed, the paper could become publishable, but in its current form the central claims are not verifiable from the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere's my read on TransFit. The paper is a credible attempt to give transient astronomers a fast numerical light-curve solver that improves on Arnett-like models by solving the time-dependent diffusion equation with flexible density and heating profiles. The parameter scans in Figs. 2-8 are useful and reproduce known trends; the comparison against Arnett and one-zone models is instructive, especially the demonstration that Arnett's law breaks down for centrally concentrated heating. The fits to SN 1993J and SN 2011kl show the intended applications, though they are qualitative.\n\nThe soft spots are real. The stress-test is correct: Eq. (13) is missing the (rho0 tdiff eps_heat0/u0) amplitude factor. Working through the non-dimensionalization of Eq. (1), the source term must carry that prefactor. As printed, Eq. (13) makes the solution independent of eps_heat0, so varying MNi would do nothing, contradicting Figure 3. The only way it works is to impose u0 = rho0 tdiff eps_heat0, which is stated only later in the Arnett section, and which ties ETh,in to MNi. The paper treats ETh,in and MNi as independent (Table 1). So either the equation is wrong or there's an unstated parameter coupling. This is load-bearing and needs to be fixed or clarified.\n\nSecond, the abstract and intro claim Sedona validation, but I cannot find a Sedona comparison anywhere in the body. That claim should be substantiated or removed.\n\nThird, the fits have no uncertainties, and the gray opacity and gamma-ray leakage assumptions are untested. That's not fatal for a forward-modeling paper, but it limits the inferred-parameter claims.\n\nNo code is released, so the source-term issue can't be checked. That's an easy fix to require.\n\nThe core physics and numerics are standard, and the idea is useful. The paper deserves peer review, but the authors need to correct the normalization, show or remove the Sedona claim, release the code, and add residual/uncertainty analysis. If they do, this will be a practical tool for LSST/ZTF fitting.\n\nRecommendation: send to peer review with major revision.","headline":"A useful fast light-curve solver with a correct core idea, but the printed PDE is missing the heating-amplitude factor and the reported Sedona validation is absent; fixable, not rejectable.","tokens_in":37891,"tokens_out":38155,"would_cite":false,"duration_ms":266530,"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":"TransFit claims that solving a time-dependent dimensionless diffusion equation computes supernova light curves with Monte Carlo-level accuracy at semi-analytic speed, including early shock-cooling phases.","keywords":["supernovae","light curves","radiative diffusion","Arnett model","shock cooling","magnetar spin-down","time-domain surveys","nickel-56 heating"],"falsifier":"Fit a supernova that has an independent ejecta-mass constraint, for example from nebular spectroscopy, with TransFit using opacities of 0.1 and 0.4 cm$^2$ g$^{-1}$; if the inferred ejecta mass and kinetic energy shift by more than the observational uncertainties, the single-opacity assumption is the practical bottleneck. A second check is to compare the predicted rise time against a frequency-dependent time-dependent radiative-transfer calculation for an event with strong line blanketing.","tokens_in":36743,"feed_emoji":"💥","tokens_out":7374,"duration_ms":70859,"temperature":0.7,"pith_summary":"TransFit claims that solving a generalized, time-dependent energy-conservation equation—rather than assuming a temperature profile of fixed shape as Arnett-like models do—yields supernova light curves nearly as accurate as Monte Carlo radiative transfer but fast enough for survey-scale fitting. The framework follows the diffusion of radiation through homologously expanding ejecta, allows the heating source to be spatially distributed (radioactive nickel or a central engine), and evolves the temperature profile in time. If correct, it lets observers fit the full bolometric light curve, including early shock-cooling and the diffusion-delayed rise, for thousands of transients instead of only the peak and tail. The paper demonstrates this with fits to SN 1993J and SN 2011kl.","feed_headline":"Supernova light-curve fits now run in under a second","feed_subtitle":"TransFit models early shock cooling and radioactive peaks in one framework, ready for survey-scale fitting.","key_machinery":"The central object is the dimensionless energy-diffusion equation $\\partial e/\\partial y = x^{-2}\\partial_x[D(x,y)\\partial_x e] + S(x,y)$, with $D(x,y)=x^2\\eta_{ej}(x)^{-1}(R_{\\max}/R_0)$ as the diffusion coefficient and $S(x,y)=\\eta_{ej}(x)\\xi_{heat}(x)f_{heat}(y)(R_{\\max}/R_0)$ as the heating source; here $x=r/R_{\\max}$ is the comoving radius and $y=t/t_{diff}$ is time in units of the characteristic diffusion time. This equation carries the argument because it replaces the fixed-shape single mode of the Arnett model with a full numerical solution: a Crank-Nicolson finite-difference scheme on a tridiagonal grid whose inner boundary condition can represent either a reflecting core or a central-engine luminosity and whose outer boundary uses the Eddington surface condition. The bolometric luminosity is then read directly from the diffusive flux at the outer grid cell.","core_discovery":"The paper's central claim is that a dimensionless diffusion equation for the internal energy density, $\\partial e/\\partial y = x^{-2}\\partial_x[D(x,y)\\partial_x e] + S(x,y)$, solved on a comoving grid, captures the light-curve signatures that semi-analytic models miss. Unlike Arnett-type solutions based on separation of variables and a single diffusion timescale, TransFit lets a radiative diffusion wave propagate through an evolving energy profile. That evolution is what produces the early 'dark phase' before diffusion breakout, the shock-cooling dip, and the sensitivity to the nickel mixing radius; it also makes the transition from shock cooling to radioactive or magnetar power self-consistent. The authors report agreement with Monte Carlo radiative transfer simulations across a range of transient scenarios, and reproduce the double-peaked SN 1993J light curve and the magnetar-powered SN 2011kl light curve.","pith_inferences":["If the single gray-opacity assumption is the weakest link, re-fitting a sample with slightly different opacity values would show how much of the measured scatter in ejecta mass and nickel mass is an opacity artifact rather than astrophysics.","The same normalized diffusion equation with a different heating function could be applied to kilonovae and fast blue optical transients, giving a uniform cross-class fitting tool.","The explicitly evolving temperature profile could be mapped to filter-specific magnitudes more consistently than the paper's bolometric-plus-effective-temperature treatment, which would strengthen multi-band fits.","The duration of the pre-breakout dark phase could serve as an early alert diagnostic, since it is sensitive to ejecta mass and opacity within days of discovery."],"forward_implications":["Early-time behavior—the shock-cooling dip, the diffusion-delayed dark phase, and the rise to the nickel-powered peak—can be modeled in one self-consistent framework rather than patched together.","Arnett's law, which equates peak luminosity with the instantaneous heating rate, breaks down for centrally concentrated heating; TransFit implies that nickel masses and engine parameters inferred from Arnett fits can be systematically biased.","Because a single model evaluation costs a fraction of a second, Bayesian parameter estimation becomes practical for the thousands of transients expected from wide-field time-domain surveys.","The spatial distribution of heating (nickel mixing radius or central engine) becomes a fit parameter rather than a fixed input, and the model handles radioactive and magnetar power sources with the same machinery."],"supporting_citations":[{"why":"Defines the semi-analytic diffusion model and characteristic timescale that TransFit is compared against and improves on.","marker":"Arnett 1980"},{"why":"Supplies the Monte Carlo radiative-transfer benchmark used for validating TransFit light curves.","marker":"Kasen et al. 2006"},{"why":"Shows Arnett's law breaks down for centrally concentrated heating, an effect TransFit reproduces and quantifies.","marker":"Khatami & Kasen 2019"},{"why":"Provides the eigenfunction solution of the Arnett model with proper boundary conditions, which TransFit contrasts with an evolving temperature profile.","marker":"Pinto & Eastman 2000"},{"why":"Provides the SN 1993J bolometric light-curve data used in the validation fit.","marker":"Richmond et al. 1994"},{"why":"Identifies SN 2011kl and its magnetar-powered interpretation, which TransFit fits.","marker":"Greiner et al. 2015"},{"why":"Provides the SN 2011kl bolometric light-curve data used in the fit.","marker":"Kann et al. 2019"},{"why":"Delivers an earlier SN 1993J model and fit parameters that TransFit is compared against.","marker":"Nagy & Vinkó 2016"}],"fun_headline_variants":["TransFit unifies shock cooling and radioactive peaks in one fast fitter","One diffusion equation now fits supernovae from shock to nickel","Fast transient fits that capture the dark phase and double peaks","Survey-scale light-curve fits that see the shock-cooling dip","From dark phase to magnetar: TransFit's diffusion wave does it all"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes one time-independent gray opacity for the entire ejecta and keeps the diffusion approximation valid down to optical depth near unity; real supernovae have wavelength- and time-dependent line opacity, and gamma-ray leakage matters in the late tail.","fun_headline_variants_meta":{"raw":{"variants":["TransFit unifies shock cooling and radioactive peaks in one fast fitter","One diffusion equation now fits supernovae from shock to nickel","Fast transient fits that capture the dark phase and double peaks","Survey-scale light-curve fits that see the shock-cooling dip","From dark phase to magnetar: TransFit's diffusion wave does it all"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000654,"raw_usage":{"total_tokens":2996,"prompt_tokens":941,"completion_tokens":2055,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":1964}},"tokens_in":557,"tokens_out":2055,"duration_ms":12595,"temperature":1.0,"reasoning_tokens":1964,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:10:40.681523+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit a supernova that has an independent ejecta-mass constraint, for example from nebular spectroscopy, with TransFit using opacities of 0.1 and 0.4 cm$^2$ g$^{-1}$; if the inferred ejecta mass and kinetic energy shift by more than the observational uncertainties, the single-opacity assumption is the practical bottleneck. A second check is to compare the predicted rise time against a frequency-dependent time-dependent radiative-transfer calculation for an event with strong line blanketing.","supporting_citations":[],"review_version":1}