{"id":"64527e4f-6c83-4c9a-9412-805397f9d12a","arxiv_id":"2607.00081","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Hierarchical Bayesian power-law fits to early ZTF light curves of 972 SNe Ia yield population parameters for rise time, index, and color evolution, revealing a bifurcation with SALT2 stretch.","lead":"Researchers fit power-law models to the early light curves of 972 Type Ia supernovae from ZTF using a hierarchical Bayesian approach to derive population-level rise parameters. These measurements can refine models of supernova explosions and their application in cosmology.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"Power-law truncation at 30% peak flux lacks explicit validation that it avoids bias in population parameters despite noted sensitivity beyond ~40%.","rationale":"The reader's weakest assumption correctly isolates the single methodological choice whose validity is least secured by the information given. The abstract's convergence statement is reassuring but insufficient without quantitative checks across cutoffs; confirming stability would strengthen the claim, while instability would require re-deriving the population parameters and bifurcation under a validated truncation.","tokens_in":2014,"tokens_out":370,"duration_ms":13853,"concrete_test":"Re-run the hierarchical Bayesian fit on the same 972-event sample with truncation thresholds of 20%, 25%, and 35% of peak flux; if the population means for t_rise or α shift by more than 1σ relative to the 30% values, or if the x1 bifurcation disappears or changes regime boundaries, the 30% choice is not bias-free.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim of robust population constraints on t_rise, α, and the x1 bifurcation rests on the power-law model being an adequate description when data are cut at 30% of peak. The abstract states that fits are sensitive if data beyond ~40% are included but converge at earlier epochs; however, this does not demonstrate that 30% itself is free of systematic offset in the hierarchical posterior means or that the bifurcation is insensitive to the exact cutoff. If the true early-time flux has curvature, nickel mixing signatures, or color evolution that the single power-law cannot capture before 30%, the reported μ=18.55 d, σ=1.42 d, and the high-stretch vs low-stretch regime split could be artifacts of the truncation choice rather than intrinsic population properties.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper claims to derive robust population parameters for the early rise of 972 Type Ia supernovae using power-law fits to ZTF light curves truncated at 30% of peak flux within a hierarchical Bayesian model. Key results include t_rise with mean 18.55 days and dispersion 1.42 days, rise index α mean 2.10 and dispersion 0.48 in r-band, color evolution difference mean 0.20, and a bifurcation in the rise morphology versus SALT2 x1 relation, with implications for 56Ni mixing in high-stretch events.","tokens_in":2224,"tokens_out":478,"duration_ms":33164,"significance":"The study benefits from an order-of-magnitude larger sample than previous work and employs a hierarchical framework that properly accounts for uncertainties, providing empirical benchmarks for SN Ia models. The identification of distinct regimes in rise behavior and potential flux excesses offers testable predictions for explosion physics. However, the significance is tempered by the need to confirm that the truncation choice does not introduce systematic biases in the reported means and dispersions.","major_comments":[{"comment":"The claim of robust constraints on t_rise, α, and the x1 bifurcation depends on the power-law model being adequate at the 30% peak flux truncation. The abstract acknowledges sensitivity when including data beyond ~40% but does not quantify whether the population hyperparameters (μ, σ) shift significantly when the cutoff is varied between 20-35%, which is necessary to support the 'robust' qualifier.","section":"Abstract"},{"comment":"The reported bifurcation where high-stretch SNe show trends with x1 but low-stretch do not is central to the interpretation of widespread outward 56Ni mixing. Without an explicit test of how this split behaves under different truncation epochs or alternative early-time models, it is unclear if the dichotomy is physical or influenced by the modeling cutoff.","section":"Results on bifurcation"}],"minor_comments":[{"comment":"The rise index is denoted α in r-band and the color difference as α_g - α_r; a clear definition early in the text would aid readability.","section":"Notation"},{"comment":"Clarify how the 972 events were selected from the one thousand ZTF supernovae mentioned in the title to ensure reproducibility.","section":"Sample description"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading and constructive comments, which identify key areas where additional quantification would strengthen the manuscript's claims of robustness. We address each major comment below and will revise the paper to incorporate the suggested tests.","responses":[{"response":"We agree that the abstract's use of 'robust' would be better supported by an explicit demonstration that the population hyperparameters remain stable when the truncation epoch is varied between 20% and 35% of peak flux. The manuscript notes general convergence for earlier epochs but does not report the specific shifts in μ and σ across this range. In the revised version we will add a supplementary analysis (or figure) repeating the hierarchical fits at 20%, 25%, 30%, and 35% truncations and tabulating the resulting changes in the reported means and dispersions for t_rise, α, and color evolution. This will either confirm stability or allow us to qualify the robustness statement appropriately.","revision_made":"yes","referee_comment":"[Abstract] The claim of robust constraints on t_rise, α, and the x1 bifurcation depends on the power-law model being adequate at the 30% peak flux truncation. The abstract acknowledges sensitivity when including data beyond ~40% but does not quantify whether the population hyperparameters (μ, σ) shift significantly when the cutoff is varied between 20-35%, which is necessary to support the 'robust' qualifier."},{"response":"We acknowledge that the bifurcation result is central to the physical interpretation and that its stability under changes in truncation has not been explicitly demonstrated. While the manuscript already flags sensitivity beyond ~40%, we have not tested the high- versus low-stretch split at other cutoffs within 20-35%. In revision we will repeat the bifurcation analysis at 25% and 35% truncations and report whether the two-regime behavior persists. We note that a full exploration of alternative early-time models (e.g., broken power laws or physically motivated templates) lies beyond the scope of the present power-law study but will be flagged as important future work; the truncation tests alone will clarify whether the reported dichotomy is robust to the specific modeling cutoff chosen.","revision_made":"yes","referee_comment":"[Results on bifurcation] The reported bifurcation where high-stretch SNe show trends with x1 but low-stretch do not is central to the interpretation of widespread outward 56Ni mixing. Without an explicit test of how this split behaves under different truncation epochs or alternative early-time models, it is unclear if the dichotomy is physical or influenced by the modeling cutoff."}],"tokens_in":1614,"tokens_out":553,"duration_ms":18759,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main thing here is that 972 ZTF Type Ia events give population constraints on early rise time (mean 18.55 days, dispersion 1.42), power-law index (mean 2.10, dispersion 0.48 in r), and g-r color difference, plus a reported split where high-stretch events show x1 trends in rise and color that low-stretch events lack.\n\nThe sample size is the clear step forward. Earlier datasets were an order of magnitude smaller, so this volume-complete set under a hierarchical Bayesian fit produces concrete means and sigmas with uncertainties. That scale also lets them identify the bifurcation in rise morphology versus SALT2 x1, which smaller samples did not catch, and to note that long-duration flux excesses are common in the high-stretch group while tying some anomalies to possible nickel mixing.\n\nThe soft spot is the power-law truncation at 30% of peak flux. The abstract states that fits converge at early epochs but shift if data past ~40% are included, yet it does not test whether the 30% choice itself introduces offset in the population parameters or the bifurcation due to unmodeled curvature or early color evolution. The nickel-mixing and anomaly interpretations rest on the same model, so they carry the same caveat.\n\nThis is for the SN Ia early-time and progenitor community. Someone comparing explosion models to data or needing empirical rise statistics would get direct use from the numbers and the split. The work engages the literature with a straightforward empirical approach and a large dataset, so it deserves a serious referee even with the truncation question to resolve.","headline":"Large ZTF sample pins down SN Ia rise population numbers and flags a stretch bifurcation, though the 30% flux cutoff needs checks for bias.","tokens_in":2749,"tokens_out":398,"would_cite":true,"duration_ms":24117,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Type Ia supernovae exhibit a bifurcated rise morphology versus stretch relation, with high-stretch events showing shallower rises and persistent early blue colors.","keywords":["Type Ia supernovae","early light curves","rise time","power-law fits","ZTF survey","population parameters","nickel mixing","SALT2 stretch"],"falsifier":"Re-fitting the same 972 events with data included up to 40 percent of peak flux and obtaining substantially different population means or scatters for rise time or index would show the truncation choice biases the results.","tokens_in":2938,"feed_emoji":"","tokens_out":796,"duration_ms":23858,"temperature":0.7,"pith_summary":"This paper measures the early rise of 972 Type Ia supernovae from a volume-complete ZTF sample using power-law fits restricted to data before 30 percent of peak flux. It reports population averages for rise time of 18.55 days, rise index of 2.10 in r-band, and g-r color difference of 0.20, all with quantified scatter. The central finding is that the link between rise shape and overall light-curve width splits into two regimes: high-stretch supernovae display clear correlations with shallower rises and longer-lasting blue colors, while low-stretch events show no such trends. A reader would care because these measurements directly constrain the distribution of nickel in the outer ejecta and therefore the explosion mechanism.","feed_headline":"ZTF sample pins Type Ia rise at 18.55 days average","feed_subtitle":"972 events show rise-stretch relation splits into two regimes, with high-stretch events tied to outward nickel mixing.","key_machinery":"Hierarchical Bayesian power-law fits applied to light curves truncated at 30 percent of peak flux across a volume-complete sample of 972 events.","core_discovery":"Power-law modeling of early light curves for 972 ZTF SNe Ia yields population parameters t_rise of 18.55 plus or minus 0.08 days with sigma 1.42 days, alpha of 2.10 plus or minus 0.04 with sigma 0.48 in ZTF r, and alpha_g minus alpha_r of 0.20 plus or minus 0.02 with sigma 0.17. The relation between rise morphology and SALT2 x1 stretch bifurcates into two regimes where high-stretch events correlate higher x1 with shallower rises and more persistent blue colors while low-stretch events lack these trends; long-duration flux excesses around 5 days are common in the high-stretch population and tied to those features, indicating widespread outward 56Ni mixing.","pith_inferences":["The two-regime behavior could serve as an observable to separate different progenitor channels for Type Ia supernovae.","Early-time photometry from future wide-field surveys could be used to assign individual events to one regime or the other before peak.","If the regimes also differ in peak luminosity or color at maximum, the bifurcation would affect how these events are standardized for cosmology."],"forward_implications":["High-stretch SNe Ia commonly display long-duration flux excesses linked to shallow rises and early blue colors.","The dichotomy between high- and low-stretch populations in rise morphology must be reproduced by multi-dimensional explosion models that include realistic progenitor setups.","Several normal SNe Ia display unusually long rise times that may include short-duration flux excesses over a smooth rise.","Rise times correlate positively with x1 overall, but this relation flattens within the high-stretch population."],"fun_headline_variants":["18.55 days average Type Ia rise from ZTF 972 sample","Two rise regimes by stretch in 972 Type Ia supernovae","High stretch events show shallow rises and blue colors","5 day flux excesses tied to high stretch SNe Ia","Rise index 2.10 with color evolution 0.20 in ZTF r"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The power-law model remains an adequate description of the light curves when the fit is restricted to data up to 30 percent of peak flux.","fun_headline_variants_meta":{"raw":{"variants":["18.55 days average Type Ia rise from ZTF 972 sample","Two rise regimes by stretch in 972 Type Ia supernovae","High stretch events show shallow rises and blue colors","5 day flux excesses tied to high stretch SNe Ia","Rise index 2.10 with color evolution 0.20 in ZTF r"]},"model":"grok-4.3","cost_usd":0.007482,"raw_usage":{"total_tokens":3575,"prompt_tokens":949,"num_sources_used":0,"completion_tokens":88,"cost_in_usd_ticks":74824500,"prompt_tokens_details":{"text_tokens":949,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2538,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":949,"tokens_out":88,"duration_ms":17683,"temperature":1.0,"reasoning_tokens":2538,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-02T17:53:26.079900+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Re-fitting the same 972 events with data included up to 40 percent of peak flux and obtaining substantially different population means or scatters for rise time or index would show the truncation choice biases the results.","supporting_citations":[],"review_version":1}