{"id":"74d73461-a837-4db2-8b84-530e7be75371","arxiv_id":"2507.13429","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In a suite of 12 supernova models, electron-neutrino lepton number crossings appear at larger radii when proto-neutron star convection is included and at smaller radii when muon production is included.","lead":"This paper simulates neutrino angular distributions in 12 core-collapse supernova models that vary the nuclear equation of state, muon production, and proto-neutron star convection. It finds that convection shifts electron-neutrino lepton number crossings outward, while muon creation shifts them inward.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The known νe/¯νe density flip in the post-processing solver at tpb≥0.5 s can manufacture the very ELN-crossing shifts attributed to convection and muons; the paper's relative-trend claim lacks a test that crossings lie outside the flip region.","rationale":"The reader's weakest assumption—that the simplified Boltzmann solver reliably captures relative differences across models—is precisely the load-bearing point. I agree and sharpen it: the solver's known νe/¯νe number-density flip at tpb≥0.5 s directly alters the sign of the ELN, so any crossing location computed in that regime is suspect. The paper's own Appendix A admits the flip 'might affect the formation of ELN crossings,' yet the conclusions for convection and muons rely on crossing positions at exactly these times. The fact that the radial domain is model-dependent (Table II) means the artifact could produce the observed ordering without physics. This is not an accusation; it is an untested alternative explanation. The proposed test would separate the artifact from the physical deleptonization-dip effect. If the test fails, the abstract's central claim is an artifact of the transport scheme; if it passes, the paper's conclusions are robust. Therefore I recommend keeping the reader's conditional verdict, with the added condition that the flip artifact be ruled out.","tokens_in":17466,"tokens_out":6575,"duration_ms":73220,"concrete_test":"For each model and tpb≥0.5 s, compute the flip radius r_flip where the Boltzmann n_νe and n_¯νe profiles cross. Overlay r_flip on Fig. 6 and test whether the reported ELN crossing intervals lie at/beyond r_flip and whether the convection/muon-induced shift in crossing radii correlates with the shift in r_flip. As a stronger check, rerun the solver for the benchmark and +c/+m models at tpb=0.5 and 1.0 s with the full energy range to 380 MeV and inelastic scattering, or with the νe/¯νe normalization forced to match VERTEX (n_νe > n_¯νe at all r), and compare the ELN crossing radial ordering. If crossings persist outside the flip region with unchanged relative ordering, the claim survives; if not, the central trend is an artifact.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Appendix A documents that, for tpb ≳ 0.5 s, the Boltzmann solver produces n_νe < n_¯νe at large radii (the 'flip'), whereas VERTEX has n_νe > n_¯νe everywhere. Since ELN crossings are defined by sign changes of n_νe − n_¯νe, this artifact can directly create or move crossings. The central claims—convection shifts crossings to larger radii and muons to smaller radii—are derived from crossing locations at tpb = 0.5, 0.75, 1.0, and 3.0 s (Fig. 6), precisely where the flip is present. Moreover, the radial domain [rmin, rmax] is shifted outward for convection models (Table II), so the flip footprint may shift with the domain and mimic the reported radial ordering. The paper asserts that relative changes are 'sufficiently reliable' but provides no analysis separating physical crossing locations from flip-induced ones. If the reported radial shifts track the flip radius rather than the deleptonization dip, the abstract's central conclusion would be an artifact.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper investigates how electron-neutrino lepton number (ELN) angular crossings in core-collapse supernovae depend on nuclear equation of state, muon production, and proto-neutron-star (PNS) convection. The authors use a suite of 12 spherically symmetric neutrino-hydrodynamics simulations of an 18.6 M_sun progenitor, extract static fluid profiles at six post-bounce times, and solve the Boltzmann equation in post-processing to obtain neutrino angular distributions. ELN crossings are identified with criteria from the authors' earlier work. The central claim is that PNS convection shifts the PNS radius outward and therefore favors ELN crossings at larger radii, while muon creation contracts the PNS and favors crossings at smaller radii, with only mild dependence on the nuclear equation of state.","tokens_in":17732,"tokens_out":8225,"duration_ms":97481,"significance":"If established, the result would be a useful, falsifiable prediction: the location of fast-flavor-conversion conditions in core-collapse supernovae would depend measurably on PNS convection and on muonic microphysics. The paper is commendable for using external hydrodynamic simulations, for solving the Boltzmann equation without fitting any target quantity, and for honestly quantifying discrepancies with VERTEX, including O(10%) number-density differences, up to about 50% flux differences, and the acknowledged νe/anti-νe number-density flip at large radii for tpb ≳ 0.5 s. However, the main claims are drawn precisely from the time window and radial region in which the post-processing solver is known to disagree with VERTEX in a way that can affect ELN crossings. The significance therefore currently rests on an unresolved artifact that must be isolated before the central conclusion can be accepted.","major_comments":[{"comment":"The central radial-ordering claim is made at tpb = 0.5, 0.75, 1.0, and 3.0 s, precisely the post-bounce times for which the Boltzmann solver produces the acknowledged νe/anti-νe number-density flip at large radii (Appendix A), while VERTEX has νe > anti-νe everywhere. Since ELN crossings are defined by sign changes of the νe − anti-νe angular distribution, this flip can create or move crossings. The paper states that the flip 'might affect the formation of ELN crossings' (Appendix A) but does not test whether the crossing radii reported in Fig. 6 lie inside or outside the flip region, nor whether the convection/muon ordering survives when the affected radii are excluded or when VERTEX angular distributions for the benchmark model are used. The assertion in Section V that relative changes are 'sufficiently reliable' is therefore not yet supported by the evidence presented.","section":"Appendix A; Sec. IV, Fig. 6"},{"comment":"For the convection models the radial integration domain starts at systematically larger rmin (e.g., 18 vs 15 km for LS220 at 0.5 s, 16 vs 13 km at 0.75 s, and 14 vs 12 km at 1.0 s; Table II). Because crossings below rmin are not computed, the reported shift of convection crossings to larger radii could partly reflect the shifted domain rather than a physical effect. The analysis should either use a common radial domain for paired runs or explicitly verify that the no-convection models have no crossings inside the convection-model domain.","section":"Sec. IV, Fig. 6; Table II"},{"comment":"Beyond the flip, the Boltzmann solver differs from VERTEX in its collisional kernel, energy range (1–100 MeV vs up to 380 MeV), and boundary conditions, with number densities differing by O(10%) and fluxes by up to about 50% (Section V). The paper's conclusions concern the locations of angular crossings, a quantity that is sensitive to the shape of the angular distributions and is not among the validated moments. The authors should demonstrate stability of the model-to-model crossing locations under these known systematic differences, for example by varying the boundary prescription, by checking sensitivity to the energy cutoff near the decoupling region, or by comparing the benchmark model against VERTEX's stored angular distributions. Without such a test, the claim that relative changes are sufficiently reliable is not established.","section":"Sec. V; Appendix A"}],"minor_comments":[{"comment":"The collision operator for antineutrinos is written with a bar in the second line of Eq. (1), but the text introduces only C; please define the barred operator explicitly.","section":"Sec. III.A, Eq. (1)"},{"comment":"The statement that ELN crossings appear at all post-bounce times except tpb = 3 s refers to the benchmark model only; Section IV and Fig. 5 report crossings at 3 s for some muon models. Please state the model dependence explicitly in the caption.","section":"Sec. III.B and Fig. 4 caption"},{"comment":"The caption lists specific radii (133.2 km, 20.8 km, 21.8 km, 22.8 km) but does not state which model each value belongs to; specify whether these are the extraction radii for the benchmark model or give the full set.","section":"Fig. 5 caption"},{"comment":"The vertical lines in Fig. 6 are described as the radial range where ELN crossings are found, but the paper does not define 'crossing radius' precisely; since crossings occur in (r, cos θ) space, please define what it means for a crossing to be located at a given radius.","section":"Fig. 6 and Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the journal's scope and the authors are transparent about the limitations of their post-processing solver. My main reservation is that the headline conclusion is taken from the regime where the solver is known to produce a νe/anti-νe reversal not present in VERTEX, so a quantitative robustness test is needed before the central claim can be trusted. This is fixable in revision, which is why I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the systematic sweep: 12 spherically symmetric 18.6 M⊙ models, with three EoSs, muons on/off, and PNS convection on/off, all post-processed through a Boltzmann solver to give ELN crossing locations across six post-bounce times. That is a useful matrix. The qualitative trends—convection pushes crossings outward, muons pull them inward—are concrete and physically motivated through the PNS radius and deleptonization dip. The paper also deserves credit for its honesty: Appendix A quantifies the disagreement with VERTEX (O(10%) number densities, up to 50% fluxes), tests energy resolution, and explicitly admits the νe/ν̄e flip at large radii for tpb ≥ 0.5 s.\n\nBut that flip is the soft spot, and it is a load-bearing one. The central claims in the abstract and Fig. 6 are derived from exactly those late times—0.5, 0.75, 1.0, and 3.0 s—when the flip is present. The flip is a sign change in n_νe − n_ν̄e, and ELN crossings are defined by sign changes of that same quantity. So the artifact can directly create or move crossings. The paper states the flip \"might affect the formation of ELN crossings,\" but then asserts that relative changes are \"sufficiently reliable\" without showing that the crossings sit outside the flip region, or that the trends survive once flip-affected radii are excluded. The radial domain itself is shifted outward for convection models (Table II), so the flip footprint could shift with the domain and mimic the reported radial ordering. The authors need to demonstrate, not assert, that the crossings they count are physical.\n\nMinor issues: no code or data release, which makes it hard to check the crossing selection criteria against the raw distributions. The comparison to Ref. [23] is a side point, and the discussion of multidimensional effects is appropriately hedged.\n\nTo be clear, the paper is not sloppy and the transport method is not novel in itself—it is inherited from prior work. The value is in the systematic comparison. But as it stands, the abstract's central conclusion is vulnerable to a documented artifact in the method. A serious referee should ask for a quantitative separation of flip-induced and physical crossings before the radial-shift claims are accepted. If the trends survive that test, this becomes a solid within-subfield reference. As written, it is a conditional accept at best.\n\nI would send it to peer review—the question matters and the suite is genuinely useful—but I would not yet cite the directional trends as established.","headline":"Systematic 12-model post-processing study of ELN crossings, but the headline radial-shift claims sit exactly where the solver's documented νe/ν̄e number-density flip could manufacture the crossings.","tokens_in":18253,"tokens_out":1836,"would_cite":false,"duration_ms":24191,"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":"This paper shows that proto-neutron-star convection shifts electron-neutrino lepton number crossings to larger radii, while muon creation moves them to smaller radii in a suite of 12 supernova models.","keywords":["electron-neutrino lepton number crossings","fast neutrino flavor conversion","core-collapse supernovae","proto-neutron star convection","muon creation","neutrino transport","Boltzmann equation","nuclear equation of state"],"falsifier":"A direct calculation that would settle the claim: rerun the same twelve static fluid profiles through a Boltzmann solver with the full collision kernel and boundary conditions matched to the simulations, storing the angular distributions, and check whether convection still moves ELN crossings to larger radii and muons to smaller radii; if the relative ordering by model changes or disappears, the central claim fails.","tokens_in":17327,"feed_emoji":"🔄","tokens_out":7207,"duration_ms":81590,"temperature":0.7,"pith_summary":"This paper aims to establish that the microphysics and hydrodynamics of the supernova core leave a measurable imprint on the angular distributions of electron neutrinos and antineutrinos, and therefore on where conditions for fast neutrino flavor conversion arise. The authors solve the Boltzmann equations on static fluid profiles from 12 spherically symmetric simulations of an 18.6 solar-mass star, varying the nuclear equation of state, the presence of muons, and a mixing-length treatment of proto-neutron-star convection. They find that convection shifts the proto-neutron-star radius outward, pushing electron-lepton-number (ELN) crossings to larger radii, while muon creation contracts the star and pulls crossings inward. These effects are mild across equations of state. A reader should care because ELN crossings are the standard diagnostic that fast flavor conversion is possible, and such conversion may influence the explosion, neutrino signals, and nucleosynthesis.","feed_headline":"Convection and muons shift where supernova neutrinos flip flavor","feed_subtitle":"Whether supernova neutrinos can undergo fast flavor conversion depends on core convection and muon creation","key_machinery":"The central object is the electron-neutrino lepton number (ELN) angular distribution, $d(n_{\\nu_e} - n_{\\bar{\\nu}_e})/d\\cos\\theta$ as a function of propagation angle; an ELN crossing is a sign change of this distribution in $\\cos\\theta$, which signals favorable conditions for fast flavor conversion. The machinery is a steady-state Boltzmann transport solver in spherical symmetry, applied to static radial profiles of density, temperature, chemical potentials, and lepton fractions extracted from each supernova simulation at six post-bounce times. The solver uses 100 energy bins, 300 angular bins, and 150 radial bins, and a collision kernel inspired by an open-source neutrino radiation hydrodynamics code; the solution is compared across models with and without convection and muons. The comparison of crossing radii across models is the argument: convection changes the PNS radius and the location of the deleptonization dip, and muons change the contraction rate.","core_discovery":"The central discovery is a systematic correlation between core physics and the location of ELN crossings in neutrino angular distributions. Across all 12 models and six post-bounce times, ELN crossings for forward directions appear in almost all snapshots except late cooling phases, but their radial locations and shapes change in a pattern: adding a mixing-length treatment of PNS convection moves the deleptonization dip in the electron-fraction profile outward and thereby shifts crossings to larger radii; including muon production softens the equation of state, speeds up PNS contraction, and makes crossings appear at smaller radii. The nuclear equation of state has only mild influence. The paper further argues that the crossings appear after neutrino decoupling in all its spherical models, contrary to an earlier conclusion based on different spherical models, and attributes the difference to the treatment of neutrino transport.","pith_inferences":["This is an editorial extension: if the radial shift is as systematic as reported, fast flavor conversion could ignite at different depths in supernovae with muons, potentially changing the neutrino spectra at Earth in a way the paper does not calculate.","A natural extension is to repeat the twelve-profile comparison with the full simulation collision kernel and time-dependent boundaries; the central claim would be strengthened if the relative ordering of crossing radii by model survives that test.","Another editorial inference: the paper's late-time result suggests the Kelvin-Helmholtz cooling phase may be a less promising epoch for fast flavor conversion, but the paper does not simulate the flavor conversion itself, so the observable neutrino signal could differ."],"forward_implications":["ELN crossings, and hence favorable conditions for fast flavor conversion, appear at almost all post-bounce times in all twelve models, so the phenomenon is not limited to a special progenitor or equation of state.","Models that include proto-neutron-star convection will have broader forward-peaked neutrino angular distributions and crossings at larger radii than otherwise identical models without convection.","Models that include muon production will have crossings at smaller radii, because faster PNS contraction shifts neutrino decoupling inward.","The paper supports the view that moment-based closure schemes alone are insufficient to reliably infer ELN crossings, and that locating them reliably requires solving the Boltzmann equation for the angular distributions.","The same microphysical effects should also affect fast-flavor conditions in neutron-star merger remnants, where similar neutrino decoupling physics operates."],"supporting_citations":[{"why":"Supplies the earlier spherical-symmetry Boltzmann result reporting no ELN crossings, which this paper directly addresses and contradicts with its own solver.","marker":"[23]"},{"why":"Provides the method-comparison framework and selection criteria for robust ELN crossings, and defines three of the models used here as benchmarks.","marker":"[37]"},{"why":"Introduces the post-processing approach of solving Boltzmann equations on static supernova profiles, which this paper follows.","marker":"[38]"},{"why":"Further develops the post-processing Boltzmann method and is cited as the basis for the transport setup.","marker":"[39]"},{"why":"Gives the one-dimensional supernova models with and without muon creation whose ELN features motivate the muon comparison.","marker":"[32]"},{"why":"Reports PNS convection effects on ELN crossings in three-dimensional models, providing the prior finding that this paper extends.","marker":"[35]"},{"why":"Establishes the physical effect of muon creation in supernova matter that softens the equation of state and accelerates PNS contraction.","marker":"[47]"},{"why":"Supplies the mixing-length treatment of proto-neutron-star convection used in half of the models.","marker":"[48]"},{"why":"Provides the collision-kernel modeling that the Boltzmann solver's interactions are inspired by.","marker":"[54]"}],"fun_headline_variants":["Muons and convection steer supernova neutrino crossing sites","Core convection and muon creation shift neutrino flavor crossings","Supernova core physics reshapes neutrino flavor crossing zones","Neutrino flavor crossings follow core convection and muons","PNS convection and muons relocate neutrino flavor conversion triggers"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the post-processing Boltzmann solver, with its simplified collision kernel and static boundary conditions, captures the relative differences in neutrino angular distributions across models well enough that the predicted radial shifts in ELN crossings are real, despite known discrepancies with the full hydrodynamics code (number densities agree to about 10%, fluxes differ by up to 50%, and the solver can flip the $\\nu_e$/ $\\bar{\\nu}_e$ number-density ordering at large radii for $t_{\\rm pb} \\ge 0.5$ s).","fun_headline_variants_meta":{"raw":{"variants":["Muons and convection steer supernova neutrino crossing sites","Core convection and muon creation shift neutrino flavor crossings","Supernova core physics reshapes neutrino flavor crossing zones","Neutrino flavor crossings follow core convection and muons","PNS convection and muons relocate neutrino flavor conversion triggers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000227,"raw_usage":{"total_tokens":1487,"prompt_tokens":975,"completion_tokens":512,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":434}},"tokens_in":591,"tokens_out":512,"duration_ms":5839,"temperature":1.0,"reasoning_tokens":434,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:23:50.960487+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct calculation that would settle the claim: rerun the same twelve static fluid profiles through a Boltzmann solver with the full collision kernel and boundary conditions matched to the simulations, storing the angular distributions, and check whether convection still moves ELN crossings to larger radii and muons to smaller radii; if the relative ordering by model changes or disappears, the central claim fails.","supporting_citations":[],"review_version":1}