{"id":"c15f3307-e5a9-4812-b3fe-e53b7c3d39ac","arxiv_id":"1908.06990","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The luminosity shell in a supernova's planar phase is not fixed at the breakout shell; its mass grows logarithmically with time, changing predicted temperatures by two orders of magnitude for compact progenitors.","lead":"This paper derives a self-similar solution for radiation diffusion in a supernova envelope and shows that, during the planar phase, the source of the observed light moves inward in mass, growing logarithmically with time. The correction leaves the bolometric luminosity almost unchanged but shifts the predicted observed temperature for blue supergiant and Wolf-Rayet explosions down by about two orders of magnitude.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (23) fixes an undetermined proportionality constant in Eq. (22) to unity without derivation; this constant controls the claimed factor ~10 in the luminosity-shell mass.","rationale":"The reader's weakest-assumption call identified the single power-law density profile and the gravity-free Sakurai solution as the main risk. My reading of the manuscript points to a more immediate, internal gap: the claimed logarithmic law and the factor ~10 depend on an undetermined constant in Eq. (22), which the paper does not evaluate. The paper is otherwise careful: it states its approximations, gives the self-similar ODE, and acknowledges the 50 keV validity limit. The qualitative statement that the luminosity shell recedes logarithmically in mass is likely robust, since it follows from the heuristic leakage argument in Eq. (17) as well as from the self-similar ansatz. But the specific quantitative claim that the shell reaches about ten times the breakout-shell mass, and the corresponding statement that the radiation originates in regions about ten times denser, is not pinned down unless the missing constant is computed. Because a re-derivation can settle this directly, the appropriate verdict is conditional acceptance rather than rejection. I disagree with the reader's identification of the weakest assumption because the same concern applies even for a perfect power-law envelope and an exact Sakurai profile.","tokens_in":28003,"tokens_out":26495,"duration_ms":276852,"concrete_test":"Re-derive the self-similar ODE by substituting the full ansatz (21)/(24) into Eq. (19), keeping the proportionality constant in Eq. (22) explicit as F'(t)F(t)^{-μn/(n+1)} = C/t. Solve for C from the requirement that all explicit t-dependence cancels in the transformed diffusion equation, then compute A = C(1 - μn/(n+1)) mbo^{-μn/(n+1)} and evaluate mls(ts)/mbo = [1 + A ln(ts/tbo)]^{(n+1)/(n+1-μn)}. If A differs from 1 by more than ~30%, the headline 'factor ~10' and the corresponding density and thermalization enhancement should be revised. An independent check would be a 1D radiation-hydrodynamics simulation of a single power-law envelope, tracking the Lagrangian mass coordinate of the escaping luminosity during the planar phase, but the analytic re-derivation is the direct and decisive test.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative claim is Eq. (23), mls(t) = mbo[1 + ln(t/tbo)]^{(n+1)/(n+1-μn)}, which at the end of the planar phase gives mls ≈ 10 mbo. The derivation of this equation rests on Eq. (22), F'(t) F(t)^{-μn/(n+1)} ∝ t^{-1}, together with the initial condition F(tbo) = mbo. Integrating this first-order equation gives F(t) = mbo [1 + A ln(t/tbo)]^{(n+1)/(n+1-μn)}, where A is a combination of the undisplayed proportionality constant and the normalization mbo^{μn/(n+1)}. The paper states Eq. (23) with A = 1, but it never derives A by substituting the full ansatz (24) into Eq. (19) and enforcing that the transformed equation is time-independent. The text's later observation that F = 0 at t = tbo/e is a consequence of A = 1, not an independent physical condition, since the power-law profile is explicitly acknowledged to be unphysical as m → 0. The endpoint factor is sensitive to A: for ts/tbo = 10^3 and n = 3/2 (exponent q ≈ 1.13), A = 1 gives mls(ts)/mbo ≈ 10, A = 0.1 gives ≈ 1.8, and A = 10 gives ≈ 120. Because the abstract's 'density ~10 times higher' and the thermalization enhancement scale as powers of this factor, the headline result is controlled by an unstated constant. This is an internal gap in the self-similar derivation, not a question of realistic progenitor structure.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a self-similar solution for the specific energy of radiation diffusing through a supernova envelope during the planar post-breakout phase, building on the Sakurai (1960) power-law hydrodynamic profiles. Its central claim is that the luminosity shell is not fixed at the breakout shell: the Lagrangian mass of the shell from which the observed radiation escapes grows logarithmically with time, m_ls(t) = m_bo [1 + ln(t/t_bo)]^{(n+1)/(n+1-\\mu n)} (Eq. 23), reaching roughly 10 m_bo by the end of the planar phase. The authors argue that this logarithmic recession exposes regions about ten times denser than the breakout shell, enhancing free-free photon production and accelerating thermalization. They apply the result to RSG, BSG, and WR progenitors, and derive revised bolometric luminosities, observed temperatures, and the planar-to-spherical transition.","tokens_in":28369,"tokens_out":7143,"duration_ms":82345,"significance":"If the derivation is completed, the paper corrects a longstanding assumption in early supernova emission models: the planar-phase luminosity shell is not identical to the breakout shell. The claimed effect on the observed temperature, especially for BSG and WR explosions, is observationally consequential and would shift predicted early-time spectra from X-rays toward the UV/soft X-ray bands. The paper's self-similar framework is transparent, the exponents in the luminosity and temperature scalings are internally consistent, and the authors compare with previous analytic and numerically calibrated models rather than fitting their own output. The main quantitative claim, however, rests on an undetermined normalization in the derivation of Eq. (23), which must be addressed before the factor-of-ten result can be considered established.","major_comments":[{"comment":"The derivation of Eq. (23) from Eq. (22) is incomplete and load-bearing. Integrating F'(t) F(t)^{-\\mu n/(n+1)} ∝ t^{-1} with F(t_bo)=m_bo yields F(t) = m_bo [1 + A ln(t/t_bo)]^{(n+1)/(n+1-\\mu n)}, where A is an undetermined combination of the proportionality constant in Eq. (22), the normalization m_bo^{\\mu n/(n+1)}, and the exponent. The paper simply states Eq. (23) with A=1, and Eq. (24) builds that choice into the self-similar variables. The subsequent central estimate m_ls(t_s)≈10 m_bo and the associated 'density ~10 times higher' claim are controlled by A; for t_s/t_bo=10^3 and n=3/2, A=0.1 gives m_ls/m_bo≈1.8 while A=10 gives ≈120. The observation in the paragraph following Eq. (23) that the first leaking shell sits at m=0 when t=t_bo/e is a consequence of the A=1 choice, not an independent physical boundary condition, because the power-law profile is explicitly unrealistic as m→0. Please derive A by substituting a general ansatz for F(t) into Eq. (19) and enforcing that the transformed equation is time-independent, or equivalently show that the proportionality in Eq. (22) carries the required coefficient.","section":"§4.2, Eqs. (22)–(24)"},{"comment":"The quantitative application assumes a single power-law density profile ρ ∝ (R-r)^n over the entire radial range probed by the receding luminosity shell. The authors acknowledge this limitation in Section 8, but because the luminosity shell reaches ~10 m_bo, this range is not confined to the outermost envelope where the single power law is best calibrated. I request a quantitative sensitivity test, for example evaluating Eq. (23) for a broken power-law profile or against one of the MESA models mentioned in Section 8.1, to show that the factor ~10 is not an artifact of extrapolating the power law to larger depth.","section":"§8 and Eq. (23)"}],"minor_comments":[{"comment":"Please specify that log denotes the natural logarithm and show the intermediate integration step between Eq. (22) and Eq. (23), including the proportionality constant explicitly.","section":"§4.2"},{"comment":"Table 1 is captioned for γ=5/3 while Table 2 is captioned for γ=4/3; please clarify which adiabatic index is used in the S60 shock solution and which is used for the radiation-dominated internal energy, and whether the table entries are consistent with the equations that follow.","section":"Tables 1 and 2"},{"comment":"The order-of-magnitude estimate in Eq. (18) gives m_ls ≈ m_bo [2 log(t/...)]^{(n+1)/(n+1-\\mu n)}, while the self-similar solution in Eq. (23) gives m_bo [1 + log(t/t_bo)]^{...}; a short paragraph reconciling the factor of two and the additive constant would help readers understand which prefactor is physical.","section":"§4.1 and §4.2"},{"comment":"The abstract states that the observed temperature will decrease by two orders of magnitude, but the quantitative comparison with the uncorrected models in Figures 7 and 8 appears closer to one order of magnitude at the end of the planar phase; please state the comparison point explicitly so the claim is unambiguous.","section":"Abstract and §8"}],"recommendation":"major_revision","confidential_remarks":"The paper is interesting and likely correct in its qualitative picture, but the undetermined constant A in Eqs. (22)-(23) is currently not derived; since the factor-of-ten and density-enhancement claims scale directly with A, this is a blocking technical point. I think it is fixable within the manuscript's scope: the authors need to complete the similarity reduction and either derive A or show that it is forced to unity. The second major comment about profile robustness is less severe and could be addressed by a short sensitivity test."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper finds a genuinely new piece of physics—the luminosity shell during the planar phase is not the breakout shell, it recedes logarithmically in mass—and the self-similar machinery is mostly sound. But the headline numbers, m_ls ~ 10 m_bo and the temperature drop, sit on a constant that the paper never derives.\n\nThe heuristic argument in Sec. 4.1 is plausible, and the self-similar solution is a real step beyond Nakar & Sari (2010) and Piro et al. (2010). The authors are right that a fixed mass coordinate was an assumption, not a result. The luminosity itself is only weakly affected, which they say honestly. The temperature effect is the interesting part, and the improvement in thermal coupling from moving to higher-density shells is physical.\n\nThe soft spot is Eqs. (22)–(23). Integrating F' F^{-μn/(n+1)} ∝ t^{-1} with F(tbo)=mbo gives F(t) = mbo [1 + A ln(t/tbo)]^q, with A an undetermined constant of proportionality. The paper states A=1. That constant controls everything quantitative: mpl ≈ 10 mbo at ts, the 'density ~10 times higher' in the abstract, and the T_obs ∝ [1+ln(t/tbo)]^{-2.6} drop. If A were 0.1, the factor is ~1.8 and the temperature drop is a factor of a few, not two orders of magnitude. The authors could derive A by substituting the full ansatz (24) into Eq. (19) and requiring time-independence, or fix it by comparison to a numerical calculation. As written, the central quantitative claim is underdetermined.\n\nOther caveats are milder: the single power-law density profile is flagged in Sec. 8, and the 50 keV validity bound is stated. There is no independent numerical validation of the temperature evolution, which would help given the A issue. The self-citation to NS10 is not a problem—the paper is transparently correcting an earlier assumption that Sari co-authored.\n\nWho is this for: anyone modeling early shock-cooling emission from blue supergiants, Wolf-Rayet stars, and compact progenitors. It deserves a serious referee, but the referee should ask for the missing derivation before the quantitative predictions can be trusted.","headline":"The luminosity-shell recession is real, but the claimed factor ~10 and the two-order-of-magnitude temperature drop rest on an undisplayed constant in Eq. (23).","tokens_in":28882,"tokens_out":5486,"would_cite":true,"duration_ms":56037,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["97.60.Bw","95.30.Lz"],"model":"deepseek-v4-flash","headline":"During the planar phase, the shell that supplies a supernova's observed radiation recedes inward logarithmically in mass, reaching about ten times the breakout mass.","keywords":["supernova shock breakout","planar phase","photon diffusion","self-similar solution","luminosity shell","thermalization","free-free emission","early supernova emission"],"falsifier":"Compute a numerical radiation-hydrodynamics model of a realistic RSG, BSG, or WR progenitor from breakout through the planar phase and track, in Lagrangian mass, the shell whose optical depth satisfies $\\tau=c/v$ and the shell from which the escaping luminosity originates. If the emitting shell's mass stays within a few percent of $m_{\\rm bo}$ rather than following $[1+\\ln(t/t_{\\rm bo})]^{(n+1)/(n+1-\\mu n)}$ up to $\\sim 10\\,m_{\\rm bo}$, the paper's central claim is wrong. A cheaper observational check is to measure the early temperature decline of a nearby type Ib/c supernova: the predicted steepening produces a drop to roughly 100 eV at the end of the planar phase, whereas the fixed-shell model keeps the radiation in hard X-rays.","tokens_in":27830,"feed_emoji":"💥","tokens_out":5835,"duration_ms":58302,"temperature":0.7,"pith_summary":"During the planar phase of a supernova, before the ejecta radius doubles, the radiation that reaches the observer was previously thought to leak from a single fixed mass shell, the breakout shell. This paper derives a self-similar diffusion solution and shows instead that the luminosity shell recedes logarithmically into the envelope in Lagrangian mass, $m_{\\rm ls}(t)=m_{\\rm bo}[1+\\ln(t/t_{\\rm bo})]^{(n+1)/(n+1-\\mu n)}$. By the end of the planar phase the emitting mass is about ten times the breakout-shell mass, so the radiation that escapes comes from regions roughly ten times denser. That shift barely changes the bolometric luminosity, but it strongly accelerates thermalization and lowers the predicted observed temperature, in some cases by two orders of magnitude, for blue supergiant and Wolf-Rayet explosions.","feed_headline":"Early supernova light leaks from a shell that sinks 10× deeper","feed_subtitle":"A logarithmic correction pushes the diffusion wave into denser layers, cutting predicted temperatures by up to two orders of magnitude.","key_machinery":"The engine of the argument is a two-stage self-similar construction. First the post-shock planar hydrodynamics is taken from the Sakurai shock solution, expressed in Lagrangian mass $m$, with density $\\rho\\propto (R-r)^n$ and shock velocity $v_{\\rm sh}\\propto (R-r)^{-\\mu n}$. Then the energy-diffusion equation $\\partial u/\\partial t=-\\partial(4\\pi r^2 F)/\\partial m-u/(3t)$ is transformed by normalizing mass and specific energy to the luminosity shell $m_{\\rm ls}(t)$; requiring the transformed equation to be an ordinary differential equation fixes $m_{\\rm ls}(t)$ to the logarithmic growth law. The same solution yields the outer diffusive energy profile $u\\propto m^{-\\mu n/(n+1)}$, whose flux sets the bolometric luminosity, and the thermal-coupling coefficient $\\eta$ (free-free photon production versus the photon number needed for a blackbody) that sets the observed temperature.","core_discovery":"The paper's central claim is that the Lagrangian mass coordinate from which the observed post-breakout radiation escapes is not constant. Because the planar phase contains many dynamical times, photons from shells deeper than the breakout shell have time to diffuse outward and refill the outer diffusive profile; the shell from which luminosity effectively emerges therefore satisfies a growing condition rather than the naive $\\tau=c/v$ condition. The self-similar solution forces this mass to grow as $m_{\\rm ls}(t)=m_{\\rm bo}[1+\\ln(t/t_{\\rm bo})]^{(n+1)/(n+1-\\mu n)}$, reaching $m_{\\rm pl}\\approx 10\\,m_{\\rm bo}$ at the planar-to-spherical transition. This means the radiation observed during the planar phase originates in layers about ten times denser and with an order of magnitude more mass than the breakout shell. As a direct consequence the free-free photon production rate rises roughly a hundredfold, thermal coupling improves with time, and the observed temperature declines faster than adiabatic cooling alone would predict, although bolometric luminosity is almost unchanged because the internal luminosity profile is flat in mass.","pith_inferences":["If the correction holds, fast-cadence UV and soft-X-ray observations of nearby type Ib/c supernovae could measure the inward recession of the luminosity shell directly by fitting the predicted $[1+\\ln(t/t_{\\rm bo})]$ temperature slope, without resolving the breakout itself.","The same Lagrangian-recession mechanism should appear in other radiation-dominated diffusion problems with power-law envelopes and many dynamical times, such as failed supernovae or neutron-star merger outflows, with the logarithmic exponent set by the local density index.","Because the predicted breakout temperatures for Wolf-Rayet stars exceed the pair-production threshold near 50 keV, an extension that includes pair production could cap $T_{\\rm obs}$ while preserving the density-enhancement effect on thermalization.","Since the bolometric luminosity is almost insensitive to the correction, observers should prioritize spectral and temperature measurements in early-time surveys rather than relying on light-curve shape alone."],"forward_implications":["The bolometric luminosity stays close to $L\\propto t^{-4/3}$, so light-curve observations alone will not reveal the correction.","For blue supergiant and Wolf-Rayet explosions, radiation remains out of thermal equilibrium through the planar phase, but the observed peak shifts from X-rays toward the UV or soft X-ray band, with $T_{\\rm obs}$ dropping by two orders of magnitude by the end of the phase.","For red supergiant explosions, the radiation reaches thermal equilibrium during the planar phase, whereas previous fixed-shell models predicted a longer nonthermal phase.","At the planar-to-spherical transition, the luminosity shell has grown to roughly $10m_{\\rm bo}$ and the transition lasts about $3t_s$, after which the planar correction leaves no trace in the observed properties.","Light-travel-time effects smear the early emission over a timescale $R/c$ and broaden the spectrum into a power law whose slope is set by the temperature decay index."],"supporting_citations":[{"why":"Supplies the previously accepted fixed-shell planar-phase model and the thermal-coupling method that the paper corrects and extends.","marker":"Nakar & Sari (2010)"},{"why":"Provides the self-similar shock and planar expansion profiles in power-law envelopes that serve as initial conditions for the diffusion solution.","marker":"Sakurai (1960)"},{"why":"Used to estimate the normalization constants $C_1$ and $C_2$ that connect the analytic profiles to progenitor mass, radius, and explosion energy.","marker":"Matzner & McKee (1999)"},{"why":"Another earlier planar-phase treatment assuming a fixed luminosity shell, which the paper's logarithmic correction revises.","marker":"Piro et al. (2010)"},{"why":"Provides the spherical-phase solution that the ejecta rejoins after the transition phase.","marker":"Rabinak & Waxman (2011)"},{"why":"Gives numerically calibrated luminosities used as a comparison for the analytic planar-phase light curves.","marker":"Shussman et al. (2016)"}],"fun_headline_variants":["Supernova light escapes from a shell that sinks deeper","Logarithmic diffusion pushes supernova emission to denser layers","Early supernova emission traced to 10× denser shell","Supernova breakout light emerges from deeper, denser matter","Planar-phase supernova light comes from sinking shell"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the pre-explosion envelope is a single power-law density profile $\\rho\\propto (R-r)^n$ and that the Sakurai shock solution holds all the way to the stellar edge; if a realistic progenitor's density deviates over the shells swept by the receding luminosity shell, the logarithmic exponent and the factor of ten would change.","fun_headline_variants_meta":{"raw":{"variants":["Supernova light escapes from a shell that sinks deeper","Logarithmic diffusion pushes supernova emission to denser layers","Early supernova emission traced to 10× denser shell","Supernova breakout light emerges from deeper, denser matter","Planar-phase supernova light comes from sinking shell"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001011,"raw_usage":{"total_tokens":4304,"prompt_tokens":1013,"completion_tokens":3291,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":629,"completion_tokens_details":{"reasoning_tokens":3210}},"tokens_in":629,"tokens_out":3291,"duration_ms":24078,"temperature":1.0,"reasoning_tokens":3210,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:30:01.454167+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute a numerical radiation-hydrodynamics model of a realistic RSG, BSG, or WR progenitor from breakout through the planar phase and track, in Lagrangian mass, the shell whose optical depth satisfies $\\tau=c/v$ and the shell from which the escaping luminosity originates. If the emitting shell's mass stays within a few percent of $m_{\\rm bo}$ rather than following $[1+\\ln(t/t_{\\rm bo})]^{(n+1)/(n+1-\\mu n)}$ up to $\\sim 10\\,m_{\\rm bo}$, the paper's central claim is wrong. A cheaper observational check is to measure the early temperature decline of a nearby type Ib/c supernova: the predicted steepening produces a drop to roughly 100 eV at the end of the planar phase, whereas the fixed-shell model keeps the radiation in hard X-rays.","supporting_citations":[{"cited_title":"2010, ApJ, 725, 904 —","cited_arxiv_id":null,"evidence_quote":"Supplies the previously accepted fixed-shell planar-phase model and the thermal-coupling method that the paper corrects and extends."},{"cited_title":"1960, Comm","cited_arxiv_id":null,"evidence_quote":"Provides the self-similar shock and planar expansion profiles in power-law envelopes that serve as initial conditions for the diffusion solution."},{"cited_title":"D., & McKee, C","cited_arxiv_id":null,"evidence_quote":"Used to estimate the normalization constants $C_1$ and $C_2$ that connect the analytic profiles to progenitor mass, radius, and explosion energy."},{"cited_title":"L., Chang, P., & Weinberg, N","cited_arxiv_id":null,"evidence_quote":"Another earlier planar-phase treatment assuming a fixed luminosity shell, which the paper's logarithmic correction revises."},{"cited_title":"2011, ApJ, 728, 63","cited_arxiv_id":null,"evidence_quote":"Provides the spherical-phase solution that the ejecta rejoins after the transition phase."},{"cited_title":"Type II supernovae Early Light Curves","cited_arxiv_id":"1610.05323","evidence_quote":"Gives numerically calibrated luminosities used as a comparison for the analytic planar-phase light curves."}],"review_version":1}