{"id":"7115394c-79c3-4327-b1e7-d2cffec1ac63","arxiv_id":"2608.13432","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Pre-explosion turbulence in 3D simulations halves electron-capture element yields in Type Ia supernovae, requiring roughly 40% higher central white dwarf densities.","lead":"Type Ia supernovae produce most of the iron in the universe, but the exact explosion mechanism is still debated. This paper's 3D simulations show that small-scale turbulence before the blast cuts the production of electron-capture elements by about half, implying denser white dwarf cores than previously inferred.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Factor-of-two EC reduction and 40% density shift depend on an assumed smoldering-phase turbulence realization (r_D≈40 km, v_rms≈170 km/s); plausible weaker or larger-scale turbulence gives a much smaller effect.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing point: the factor-of-two EC reduction is driven by the specific turbulent initial conditions imported from the authors' own Höflich & Stein (2002) model. My reading of the paper confirms that this is the correct place to focus. The central claim is a quantitative correction to spherical-model density inferences, and every step of that correction flows from f(58Ni)≈0.46 in the 'most realistic' simulation turb4b. The paper's own comparison set shows that changing r_D from 40 km to 80 km (turb1) or reducing the turbulence scale while keeping weaker velocities (wob2) substantially degrades the spherical filling that produces the low f(58Ni). The smoldering-phase turbulence properties are therefore not a minor parameter; plausible alternative simmering models, including Zingale et al. (2011), give materially different initial conditions, and the observational conclusion would collapse if those alternatives are correct. I did not find an internal inconsistency or a reason to reject the paper: the MHD setup is described carefully, the limitations of the early-phase simulation are stated honestly, and the paper is appropriately framed as a first step. But the sensitivity of the headline number to the unresolved turbulence prior justifies the CONDITIONAL verdict the reader gave. No change to that verdict is needed.","tokens_in":25120,"tokens_out":5593,"duration_ms":61055,"concrete_test":"Post-process tracer particles from the already-listed simulations turb2 (r_D=40, v_rms=20 km/s) and turb3 (r_D=40, v_rms=50 km/s) with the same HYDRA EC network used for turb4b, and also run a Zingale-type initial condition with v_rms≈50 km/s distributed over the full WD; compare f(58Ni). If f(58Ni)≳0.7 in either case, the factor-of-two reduction and the 40% density correction are not robust to the stated uncertainty in smoldering-phase turbulence; if f(58Ni)≈0.46, the turbulence-sensitivity concern is weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim requires that pre-existing smoldering-phase turbulence be both small-scale and strong. The adopted values, r_D≈40 km and v_rms≈170 km/s, are not measured; they are imported from Höflich & Stein (2002), a model whose simulated turbulence is confined to the center by a specific C-abundance gradient. Section II itself notes that Zingale et al. (2011), without that gradient, produce v_rms≈50 km/s with convection spread over the whole WD, and the authors flag that this changes the flow. The simulations then show the factor-of-two result is tightly localized: turb1 with r_D=80 km still forms buoyant plumes and incomplete filling (Section III, Figures 5 and 7), and large-scale or absent turbulence (wob2, wob1a) gives much higher f(58Ni). Since f(58Ni) for the 'most realistic' case (turb4b) is the entire basis for the claimed systematic ~40% increase in inferred central density, an unvalidated choice of initial conditions is load-bearing. The concern is not that the simulations are internally inconsistent; it is that the observational conclusion is contingent on one branch of an unsettled smoldering-phase modeling question.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 3D MHD simulations of the early deflagration phase of near-Chandrasekhar-mass white dwarfs with varying initial turbulent velocity fields and magnetic field configurations. It finds that small-scale, strong pre-existing turbulence with a diffusion radius of r_D≈40 km and rms velocity v_rms≈170 km/s, adopted from the authors' model of the smoldering phase, fills the unburned pockets that otherwise plague 3D deflagration models and restores nearly spherical burning. Post-processing tracer particles with a nuclear reaction network yields about a factor-of-two reduction in 58Ni production relative to a spherical delayed-detonation reference model, which the authors translate into a required ~40% increase in the inferred central WD density. The paper also examines the effect of magnetic fields near saturation strength and argues for the need for improved electron-capture rates at low electron fraction. The central observational implication is that near-MCh progenitors may be denser than previously inferred from spherical models.","tokens_in":25317,"tokens_out":5268,"duration_ms":55018,"significance":"If the factor-of-two result holds, the paper would materially revise inferences of central densities and progenitor masses of SNe Ia from late-time spectra, potentially shifting the preferred near-MCh channel toward densities closer to the accretion-induced collapse limit. The work also offers a physically motivated resolution to a long-standing problem in 3D deflagration modeling, namely the formation of large unburned pockets, and it makes a falsifiable prediction that additional low-Ye isotopes should appear in high-density-burning SNe Ia observed with JWST. The simulations are state-of-the-art for this problem: full-star 3D MHD with AMR, tracer-particle nucleosynthesis, and a systematic comparison to a spherical reference model. However, the quantitative claim rests on a single simulation with a specific, and explicitly uncertain, initial turbulence realization, and the paper provides no uncertainty quantification or convergence analysis.","major_comments":[{"comment":"The factor-of-two reduction in f(58Ni) and the consequent 40% density shift are based entirely on simulation turb4b, which uses r_D=40 km and v_rms=170 km/s imported from Höflich & Stein (2002). The manuscript itself acknowledges in Section II that Zingale et al. (2011) obtain v_rms≈50 km/s when the C-abundance gradient is absent, and that turb1 with r_D=80 km still forms buoyant plumes and incomplete filling (Figures 5 and 7). Since no ensemble of turbulence realizations or systematic parameter study is presented, the claimed factor of two should be presented as conditional on a specific and unvalidated smoldering-phase model, or supported by additional simulations spanning the plausible range of r_D and v_rms. As written, the abstract and conclusions state the factor as a general result, which is not justified by the data shown.","section":"Section III B, Table I"},{"comment":"No error bars or convergence tests are reported for the central quantity f(58Ni). The tracer-particle counts are only 7,500–8,500, and the paper does not show that increasing the tracer number or the grid resolution changes the yields by less than the quoted factor. Given that the main conclusion is a ratio of ~0.46, a sensitivity study demonstrating numerical convergence is needed to ensure that the result is not dominated by sampling noise or resolution effects.","section":"Section II (tracer particles); Section III B (f(58Ni))"},{"comment":"The comparison against the spherical reference model is potentially inconsistent in evolution time: the 3D models are truncated in the early deflagration phase, while the spherical model is evolved through the full delayed-detonation explosion. The paper argues that EC production shuts off once expansion begins, but the stopping time itself depends on the burning history; for models with delayed expansion (e.g., wob1a), the tracers may still be at high density when the simulation ends. A demonstration that the EC yields are converged with respect to the evolution time, or a matched-time comparison, would strengthen the central claim.","section":"Section III B"},{"comment":"The prescribed flame speed v_burn=200 km/s is a free parameter that directly sets the burning rate and therefore the onset of expansion, which the paper identifies as the key regulator of EC production. The statement that a 100 km/s run gives 'qualitatively similar behavior' is not backed by a quantitative comparison of f(58Ni). Since the factor-of-two result is sensitive to the timing of expansion, a quantitative test of v_burn sensitivity is needed.","section":"Section II (v_burn); Section III B (expansion timing)"}],"minor_comments":[{"comment":"The text states that the turbulence diffusion radius is set to ≈50 km and v_rms up to ≈200 km/s, but Table I lists r_D=40 km and v_rms=170 km/s for the key simulation turb4b; the numerical values should be consistent between text and table.","section":"Section II"},{"comment":"The caption reads 'simulations 2 and dip2b have very little difference' but the legend and text suggest this should be 'wob2 and dip2b'; this appears to be a typographical error.","section":"Section III B, Figure 8 caption"},{"comment":"The phrase 'the maximally refined innermost (400 km) 3 regions' is garbled; it should likely read '(400 km)^3' to denote a cube of side 400 km.","section":"Section II"},{"comment":"References [52] and [89] both cite the same paper (Höflich & Stein 2002) with identical bibliographic data; one duplicate should be removed.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The core physical idea is interesting and the simulations are competently executed, but the headline numbers (factor-of-two EC reduction and 40% density shift) are contingent on a single, unvalidated choice of initial turbulence parameters from the authors' own prior model. The paper would be much stronger if it included a parameter sweep or at least a quantitative discussion of the sensitivity to r_D and v_rms, as well as convergence tests. There is also a question of scope: the title promises 'Theory vs. Observations,' but the observational comparison is limited to a one-sentence inference about central densities; the paper is more a simulation study with an observational implication. This is acceptable for the journal, but the quantitative claims should be appropriately caveated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is a serious, honest paper with one genuinely new quantitative claim — pre-existing small-scale turbulence from the smoldering phase reduces 58Ni/EC yields by about a factor of two relative to spherical near-MCh models — and one much softer claim built on top of it: that spectral-inferred central densities need to go up by ~40%. The factor of two is a clean ratio from their 3D runs and it is not fitted to observations. The 40% is an interpretation that depends on an assumed rho_c^2 scaling and on which turbulence branch you believe.\n\nWhat's actually good: the paper extends the companion Shiber et al. (2026) result by doing the nucleosynthesis post-processing, which is new. The simulation matrix is well laid out — they vary turbulence scale, strength, and magnetic field morphology, and the effect is physically coherent: small-scale strong turbulence fills the pockets, delays expansion, and cuts EC production; larger-scale turbulence (turb1) or dipole fields don't. The magnetic field comparison is a useful addition. The authors list their limitations plainly at the end (early deflagration only, no code/data deposit, parameter grid to be expanded). That lowers the risk of overclaiming.\n\nSoft spots, in proportion: the factor-two number comes from single runs. No error bars, no tracer convergence test (~7,500–8,500 particles), no ensemble over turbulence realizations. That's a minor-to-moderate issue for a simulation paper at this resolution; the physics is plausible but the quantitative precision is overstated. The bigger issue is that the headline 40% density shift is load-bearing on the assumption that the smoldering phase produces r_D≈40 km, v_rms≈170 km/s turbulence. The paper itself notes Zingale et al. get v_rms≈50 km/s without the C-abundance gradient, and turb1 (r_D=80 km) shows a much weaker effect. So the observational conclusion is contingent on one branch of an unsettled modeling question. The rho_c^2 scaling is asserted without derivation or citation — that should be fixed. These are not fatal; they are 'needs a second round of analysis and a layered conclusions section' issues.\n\nWho this is for: people working on SNe Ia nucleosynthesis, JWST stable-Ni measurements, and near-MCh progenitor channels. It deserves a serious referee. I'd engage with it conditionally — the mechanism is worth taking seriously even if the quantitative density shift shouldn't be quoted yet.\n\nRecommendation: send to peer review. Ask for error estimates or an ensemble, a sensitivity test around the smoldering turbulence parameters, and a derivation or citation for the density scaling before publication.","headline":"A credible 3D result that small-scale smoldering-phase turbulence roughly halves EC yields in near-MCh SNe Ia, but the headline 40% density shift rests on a single turbulence realization and an asserted scaling.","tokens_in":25965,"tokens_out":1802,"would_cite":true,"duration_ms":18025,"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":"Pre-explosion turbulence in near-Chandrasekhar white dwarfs cuts electron-capture element production by about half, so inferred central densities rise by 40 per cent.","keywords":["Type Ia supernovae","electron-capture nucleosynthesis","deflagration-to-detonation transition","smoldering phase turbulence","magnetohydrodynamic simulations","white dwarf progenitors","stable nickel-58","mid-infrared spectroscopy"],"falsifier":"A full-star simulation of the smoldering phase that resolves the convective flow with the actual carbon abundance gradient would settle it: if the turbulence at ignition has a root-mean-square velocity near 50 km/s or eddy scales of 80 km or more, instead of the assumed 170 km/s and 40 km, the factor-of-two reduction and the 40 per cent density correction would not apply to real progenitors.","tokens_in":24867,"feed_emoji":"💥","tokens_out":13617,"duration_ms":124614,"temperature":0.7,"pith_summary":"Type Ia supernovae are the main cosmic source of iron-group elements, and the recent, nearly universal detection of electron-capture elements such as nickel-58 in their mid-infrared spectra points to dense, near-Chandrasekhar-mass white-dwarf explosions. This paper uses three-dimensional magnetohydrodynamic simulations of the early deflagration phase to argue that these explosions must start from a core already churning with the small-scale, fast turbulence left over from the pre-ignition smoldering phase. That turbulence drags the flame into the unburned pockets that otherwise survive between buoyant plumes, making the burning nearly spherical and causing the white dwarf to expand sooner. Because the expansion shuts off electron captures earlier, the production of electron-capture elements drops by about a factor of two relative to spherical models. If this is right, the central densities inferred from spherical models must be raised by about 40 per cent, moving the required progenitors close to the threshold for accretion-induced collapse.","feed_headline":"Turbulence halves electron-capture yields in Type Ia supernovae","feed_subtitle":"Spherical models need 40% denser white-dwarf cores to match observed stable nickel.","key_machinery":"The central object is the pre-existing small-scale turbulent velocity field inherited from the smoldering phase, a Kolmogorov-like spectrum with diffusion radius $r_D \\approx 40\\,\\mathrm{km}$ and rms velocity $v_{\\mathrm{rms}} \\approx 170\\,\\mathrm{km\\,s^{-1}}$, imposed as initial conditions on the three-dimensional magnetohydrodynamic grid. Its job is to drag burned material sideways into the unburned pockets between Rayleigh-Taylor plumes during the first second of deflagration, raising the filling factor toward unity, accelerating the white dwarf's expansion, and thereby shortening the time during which electron captures can run. The electron-capture yields are then recovered by post-processing passive tracer particles through a nuclear reaction network with updated weak rates.","core_discovery":"The paper claims that when a near-Chandrasekhar-mass white dwarf begins its deflagration with the small-scale, fast turbulence expected from the smoldering phase (eddies of about $40\\,\\mathrm{km}$ moving at about $170\\,\\mathrm{km\\,s^{-1}}$), the flame is dragged into the unburned pockets that otherwise remain between buoyant plumes. The burning becomes nearly spherical, the white dwarf expands earlier, and the electron-capture phase is cut short, reducing the production of electron-capture isotopes such as nickel-58 to about 46 per cent of the spherical reference yield. To keep the same nickel-58 mass, the central density inferred from spherical near-Chandrasekhar models must be increased by about 40 per cent, pushing the required progenitors toward densities near the accretion-induced collapse limit. Magnetic fields matter only when they reach roughly one per cent of the saturation field; a large-scale dipole field squeezes the burning and suppresses mixing, while a turbulent field near seven per cent of saturation further reduces the yield to about 42 per cent of the spherical value.","pith_inferences":["If the assumed 40 km, 170 km/s turbulence is close to reality, the 40 per cent density correction goes with it; the paper's own run with 80 km eddies shows the correction shrinks, so pinning down the simmering-phase turbulence is the highest-leverage next step.","A direct observational check would be to map the velocity distribution of stable nickel in many late-time mid-infrared spectra: this model predicts a nearly spherical, low-velocity core of electron-capture material, unlike the plume-dominated structures of turbulence-free three-dimensional deflagrations.","The factor-of-two correction is derived at one central density; because the correction depends on when the white dwarf expands, real events with slightly different ignition conditions should show scatter in nickel-58 yield at fixed brightness, which a dedicated sample of late-time mid-infrared spectra could test.","At the higher central densities implied by the correction, additional neutron-rich isotopes such as calcium-48 and titanium-50 should appear in the most extreme events, which future mid-infrared observations could search for."],"forward_implications":["Observed nickel-58 masses from mid-infrared spectra must be interpreted with a roughly factor-of-two reduction applied to spherical-model expectations, so the central densities inferred for the explosions rise by about 40 per cent.","Near-Chandrasekhar and merger-like scenarios, rather than helium-triggered sub-Chandrasekhar detonations, become the natural explanation for the near-universal detection of nickel-58.","The same turbulence that restores spherical burning also shortens the high-density burning time, so the final electron-capture inventory is set as much by when the white dwarf expands as by its initial central density.","Magnetic fields below about one per cent of the saturation field leave the yields nearly unchanged, while fields near a few per cent of saturation further reduce the yields and redistribute the electron-capture material; large-scale dipole fields, in contrast, suppress mixing and preserve large unburned pockets.","Improved electron-capture rates at low electron fraction are required before the density correction can be made precise, because the yields are set by the duration of high-density burning."],"supporting_citations":[{"why":"Earlier 3D simulations by the same group showing that small-scale pre-existing turbulence fills the unburned pockets; this paper extends them to electron-capture yields.","marker":"[51]"},{"why":"Spherical delayed-detonation model used as the reference for comparing the 3D electron-capture yields and for deriving the density correction.","marker":"[21]"},{"why":"Mid-infrared spectrum of a nearby Type Ia supernova showing stable nickel-58, the key observational anchor for high-density burning.","marker":"[15]"},{"why":"Mid-infrared spectrum of another Type Ia supernova providing an additional electron-capture-element detection that the models must reproduce.","marker":"[16]"},{"why":"Late-time near-infrared survey for stable nickel that defines the observed abundance trend the spherical models are compared against.","marker":"[57]"},{"why":"Source of the smoldering-phase turbulence properties (eddy scale and velocity) imposed as initial conditions in the simulations.","marker":"[52]"},{"why":"Alternative simmering-phase simulations with slower convection, used to bracket the turbulence assumption.","marker":"[53]"},{"why":"Magnetic-field deflagration models that supply the field geometries and strengths, and the positron-trapping interpretation of late-time line widths.","marker":"[54]"},{"why":"Weak-reaction rate tables that set the electron-capture rates used in the post-processing nucleosynthesis network.","marker":"[29]"},{"why":"Turbulence-free 3D deflagration simulations that exhibit the unburned-pocket problem this paper's turbulence is designed to solve.","marker":"[45]"}],"fun_headline_variants":["Turbulence cuts electron-capture yields in half for white-dwarf supernovae","Pre-explosion turbulence forces denser white dwarfs for nickel-58","Electron-capture elements halved by turbulent flame in SNe Ia","White-dwarf turbulence demands 40% more central density for EC yields","Turbulent smoldering phase halves electron-capture isotopes in SNe Ia"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on the assumed pre-explosion turbulence being small-scale and fast, with eddies roughly 40 km wide moving at about 170 km/s; if the real simmering-phase turbulence is weaker or larger-scale, the factor-of-two reduction in electron-capture yields shrinks.","fun_headline_variants_meta":{"raw":{"variants":["Turbulence cuts electron-capture yields in half for white-dwarf supernovae","Pre-explosion turbulence forces denser white dwarfs for nickel-58","Electron-capture elements halved by turbulent flame in SNe Ia","White-dwarf turbulence demands 40% more central density for EC yields","Turbulent smoldering phase halves electron-capture isotopes in SNe Ia"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001005,"raw_usage":{"total_tokens":4315,"prompt_tokens":1075,"completion_tokens":3240,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":691,"completion_tokens_details":{"reasoning_tokens":3139}},"tokens_in":691,"tokens_out":3240,"duration_ms":22909,"temperature":1.0,"reasoning_tokens":3139,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:08:04.068581+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A full-star simulation of the smoldering phase that resolves the convective flow with the actual carbon abundance gradient would settle it: if the turbulence at ignition has a root-mean-square velocity near 50 km/s or eddy scales of 80 km or more, instead of the assumed 170 km/s and 40 km, the factor-of-two reduction and the 40 per cent density correction would not apply to real progenitors.","supporting_citations":[],"review_version":1}