{"id":"0e1b8647-f000-46dd-9705-ae37d38b120e","arxiv_id":"2608.10264","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A new energy-based classifier separates proto-neutron star oscillation modes into four families and identifies the dominant high-frequency gravitational-wave feature as the PNS fundamental mode.","lead":"Oscillation modes of newborn neutron stars can be sorted into four families by measuring which restoring force, pressure or buoyancy, dominates the mode's energy. Applied to 28 supernova simulations, this new scheme identifies the dominant gravitational-wave signal as the fundamental oscillation of the proto-neutron star surface.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The HFF identification as the PNS f-mode rests on a frequency-ordering convention rather than a direct computation of the free-surface energy; this should be tested.","rationale":"I read the paper in good faith and credit the genuine independent support for the HFF identification: the parameter-free analytic estimate of Eq. (41) reproduces the HFF track, and the spectrogram overlay in Fig. 2 shows the dominant emission following this mode. However, the physical nature of that mode as an f-mode is not established by the energy diagnostic because the free-surface energy Mf is never computed. The paper itself flags this in Sec. IV A, noting that Mf+p is computed as M2 - Mg to avoid second-derivative inaccuracies, and in Sec. IV B, where the f-mode label is assigned by frequency ordering and domain-dependence arguments rather than by a direct measurement of the restoring force. This is precisely the reader's weakest assumption. A direct evaluation of Mf is the decisive test: if the free-surface energy is not dominant, the central claim's wording (HFF equals PNS f-mode) would need revision, even if the mode-tracking formalism remains useful. The concern does not invalidate the paper's methodological advance or the practical identification of the HFF track, so the CONDITIONAL verdict remains appropriate; no change to the reader's verdict is needed.","tokens_in":40114,"tokens_out":3458,"duration_ms":34337,"concrete_test":"Directly evaluate Mf from Eq. (32) for the lowest f+p-PNS mode of the reference model at several post-bounce times (e.g., 0.5, 1.0, 1.4 s), using the GREAT eigenfunctions on a high-resolution (1600-point) grid to control numerical differentiation error, and compare Mf, Mp, M-alpha, and Mg. If Mf is not the dominant restoring-force energy (or at least comparable to Mp), the f-mode label for the HFF is unsupported; if Mf dominates, the central claim survives. As a secondary check, recompute F including the psi^2 term from Eq. (B35) to confirm the classification threshold is stable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the dominant high-frequency GW feature is the PNS f-mode (Sec. VI). The physical basis for this label is the energy classification, but the f-mode is not actually identified by its defining restoring force. In Sec. IV A the authors explicitly avoid computing Mf (Eq. 32) because of second-derivative numerical inaccuracies, defining instead Mf+p = M2 - Mg (Eq. 37). In Sec. IV B the lowest-frequency f+p-PNS mode is classified as the f-mode by convention: it is 'the lowest frequency mode of the family' appearing only when the PNS surface is included. The automatic classification (Sec. V A) similarly orders f+p-PNS modes by increasing frequency and calls the first one the f-mode. Thus the HFF mode is labeled f-mode without an independent calculation showing that the free-surface restoring force dominates its energy. If a direct evaluation of Mf showed that Mp is larger for this mode, the mode would be p1-PNS, and the central claim as stated would be incorrect even though the frequency tracking and spectrogram match might remain. The approximations in the energy functional (dropped psi^2 terms in Eq. 23, approximate force density in Eq. 28) are secondary here because the HFF mode likely has F well above 0.5, but they further weaken the formal support for the label.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a physically motivated classification of linear oscillation modes of proto-neutron stars surrounded by a stalled accretion shock, based on the contribution of individual restoring forces (compression, free surface, gravitational potential, buoyancy) to the mode energy. The formalism is derived from a perturbation of the ADM mass and applied to eigenmodes computed with GREAT for 28 one-, two-, and three-dimensional core-collapse supernova simulations using two codes and six equations of state. The authors identify four mode families: PNS f+p modes, shock f+p modes, core g-modes, and surface g-modes, and they propose an automatic classification procedure based on an energy fraction F=(M2-Mg)/M2 and on the spatial localization of the mode energy. The central claim is that the dominant high-frequency GW feature (HFF) is the f-mode of the PNS, associated with the strong density gradient at the PNS surface, and that a prominent avoided crossing around 0.4 s post-bounce involves this f-mode and the first core g-mode.","tokens_in":40353,"tokens_out":5106,"duration_ms":58524,"significance":"If the central claim is correct, this is a valuable contribution to PNS asteroseismology: it offers a physically grounded resolution of a long-standing controversy about the nature of the HFF, and the automatic classification procedure is a step toward systematic universal-relation studies on large simulation ensembles. The analytic energy decomposition in Appendices B-C is carefully presented, the application to a diverse 28-model, two-code, multi-EoS dataset is a strength, and the paper is honest about its approximations (omitted psi^2 terms, approximate force density, empirical thresholds). The inclusion of independent checks via propagation diagrams and simple analytic estimates is commendable, although, as detailed below, the identification of the HFF mode as the PNS f-mode still rests partly on a frequency-ordering convention rather than on a direct computation of the free-surface energy contribution.","major_comments":[{"comment":"The label 'f-mode' for the HFF is assigned by convention rather than by direct evaluation of the defining restoring force. Section IV B states that the lowest-frequency f+p-PNS mode is classified as the PNS f-mode because it is the lowest-frequency mode of the family and appears only when the PNS surface is included, and Section V A applies the same frequency-ordering convention automatically. However, Mf (Eq. 32) is never computed separately; the quantity actually used, Mf+p = M2 - Mg (Eq. 37), includes Mp and Malpha, so the mode labeled f-PNS could in principle be p1-PNS. The independent Kelvin-model estimate in Eq. (41) provides supporting evidence, but it is an idealized scaling relation and is not a direct computation of Mf. Since the central claim of the paper (abstract and Sec. VI) is that the HFF is the PNS f-mode, I request a direct evaluation of Mf (e.g., with a regularized or alternative treatment of the second derivatives in Eq. 32) for the HFF mode, or an equivalent diagnostic that explicitly shows that the free-surface restoring force dominates the energy of this mode.","section":"Sec. IV A/B and Sec. V A"},{"comment":"The energy fraction F, which is the basis of the automatic classification, inherits the approximations in the energy functional: Eq. (23) omits psi^2 contributions and Eq. (28) uses an approximate force density. The paper acknowledges these omissions, but it does not quantify their impact. This matters because the classification threshold F=0.5 is used globally, and several modes have values of F near 0.5, for example the surface g-modes in restricted domains (Sec. IV B, Figs. 6-7) and modes near avoided crossings (Fig. 3). A sensitivity estimate of the neglected terms for a representative mode, especially those near the threshold, would establish that the classification is robust rather than dependent on the approximations.","section":"Sec. II E/F and Sec. V A"}],"minor_comments":[{"comment":"The sentence after Eq. (11) repeats the phrase 'and the normalization of the eigenfunction' twice; please remove the duplicate.","section":"Sec. II C"},{"comment":"The underbrace labels for the restoring-force terms (compression, free surface, gravitational potential, buoyancy) are typeset after Eqs. (13)-(14) in a way that makes it unclear which terms they refer to; please reposition the braces so each label is directly under its corresponding term.","section":"Sec. II D"},{"comment":"The comparison between analytic models and numerical eigenfrequencies in Fig. 11 uses frequency-rescaling factors of order unity for each family; please report the actual values of these rescaling factors, since without them the 'excellent matching' is partly by construction.","section":"Sec. IV E"},{"comment":"The automatic classification relies on several empirical thresholds (e.g., the density threshold rho_thr=1e14 g/cm^3, the N^2 threshold 7e5 s^-2, and the velocity criterion 8e7 cm/s), but only the velocity versus density surface definitions are compared; a brief sensitivity test for rho_thr and the F=0.5 threshold would strengthen the claim of robustness across the 28 simulations.","section":"Sec. V A and Appendix D"},{"comment":"The caption lists symbols for 1D, 2D, and 3D simulations but the description of the code symbols is a bit compressed; please clarify the notation in the caption or in the legend.","section":"Fig. 13 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is a serious and potentially high-impact contribution, and the authors are well positioned to address the main gap. The referee report focuses on one load-bearing point: the f-mode label of the HFF is assigned by frequency ordering within the f+p family rather than by a direct computation of the free-surface energy contribution. I would encourage the editor to request this additional computation or an equivalent diagnostic, together with a quantitative sensitivity estimate for the approximations in the energy functional, before final acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere's my read of Tseneklidou et al. The genuinely new thing is the energy-based classification scheme: instead of counting nodes, they compute the work done by compression, buoyancy, free-surface, and metric terms (Eqs. 31-34) and define F = (M2 - Mg)/M2. That lets them sort modes into four families—f+p-PNS, f+p-shock, g-core, g-surface—and do it automatically across 28 simulations with two codes and six EoS. That is a real advance over the node-counting and eigenvalue-matching approaches that have disagreed on the nature of the high-frequency feature.\n\nThe paper is mostly careful. The derivations in Appendices B and C are detailed, the analytic estimates (Kelvin's incompressible f-mode, p-mode scaling, g-mode asymptotic periods) match the numerical trends, and the propagation diagrams and radial energy distributions independently support the four-family picture. The explanation of the avoided crossing at ~0.4 s between the PNS f-mode and the first core g-mode is convincing and nicely ties the frequency break to a physical interaction.\n\nThe soft spot is exactly where the stress-test note lands. The paper never computes the free-surface energy Mf directly; it defines Mf+p = M2 - Mg and then labels the lowest-frequency f+p-PNS mode as the f-mode by convention. So the central claim—that the HFF is the PNS f-mode—does not follow from a direct measurement of the dominant restoring force. If that lowest mode were actually p-dominated, the identification would be wrong even though the frequency tracking and spectrogram match would still look fine. The authors are transparent about this, and they add propagation diagrams and a parameter-free Kelvin estimate as independent checks, which helps. But the direct computation of Mf is not that hard in principle, and it would close the gap. The neglected psi^2 terms and approximate force density are secondary, since the HFF mode likely has F well above 0.5, but they add noise to the formal argument.\n\nAlso note the classifier and its validation share the same energy functional; the BhC matching uses the same integrand. That is a mild circularity, not fatal, but worth an explicit test with an independent mode diagnostic.\n\nOverall: this is a solid, useful paper that deserves a serious referee. The recommendation would be: engage with it, and ask for a direct Mf computation or an equivalent demonstration that the free-surface force dominates for the HFF mode. If that check passes, the central claim is likely right and the classification framework is a genuine step forward for asteroseismology.","headline":"Useful energy-based mode classification for PNS oscillations; central f-mode claim is credible but rests on a frequency-ordering convention rather than a direct free-surface energy measurement.","tokens_in":40966,"tokens_out":2179,"would_cite":true,"duration_ms":21015,"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":"By classifying modes by the energy of their restoring forces, this paper identifies four families of oscillations in proto-neutron stars and shows that the main high-frequency gravitational-wave emission is the PNS f-mode.","keywords":["proto-neutron stars","core-collapse supernovae","gravitational wave asteroseismology","stellar oscillation modes","mode classification","restoring-force energy","high-frequency feature","buoyancy frequency"],"falsifier":"Recompute $F$ for the reference model at several post-bounce times using the full second-order metric perturbation $\\psi^2$ and the potential-energy form of $M_2$ rather than the kinetic proxy; if the mode that follows the high-frequency gravitational-wave track then has $F<0.5$, or if the lowest-frequency f+p-PNS mode is shown by an independent eigenfunction-tracking analysis to be a p-mode, the central claim fails.","tokens_in":39871,"feed_emoji":"🌟","tokens_out":10994,"duration_ms":96086,"temperature":0.7,"pith_summary":"This paper aims to settle what actually oscillates when a newborn proto-neutron star rings behind a stalled supernova shock. It argues that the spectrum contains four coexisting families of modes—PNS f/p modes, shock f/p modes, core g-modes, and surface g-modes—and that the dominant high-frequency gravitational-wave emission is the PNS f-mode, an interface mode tied to the steep density gradient at the star's surface. To reach this, the authors replace node-counting with an energy-based classifier that measures the restoring forces doing the work of each mode. A sympathetic reader would care because, if the paper is right, the main gravitational-wave track becomes a direct probe of the proto-neutron star's mass and radius, turning a disputed feature into a tool for asteroseismology.","feed_headline":"The dominant high-frequency ring is the proto-neutron star's f-mode","feed_subtitle":"An energy-based classifier on 28 simulations ties the main gravitational-wave track to the star's surface f-mode.","key_machinery":"The central object is the energy fraction $F=(M_2-M_g)/M_2$ (equivalently $F=(M_f+M_p+M_\\alpha)/M_2$), computed by splitting the mode energy into work done by compression, a free surface, the gravitational potential, and buoyancy. It classifies a mode as f+p when $F>0.5$ and as a g-mode otherwise, and spatial integration over the core, PNS, and post-shock regions assigns each mode to one of four families. The total $M_2$ is derived from the perturbed ADM mass using a standing-wave displacement, with the $\\psi^2$ metric contribution dropped.","core_discovery":"On the paper's own terms, the central discovery is that the mode spectrum of a proto-neutron star with a stalled accretion shock is not a single sequence but four distinct families, each living in a different region: f+p modes inside the PNS, f+p modes in the cavity between the PNS surface and the shock, g-modes in the core stable layer, and g-modes in the surface stable layer. The energy fraction $F=(M_2-M_g)/M_2$ separates f+p from g character, and the lowest-frequency f+p-PNS mode—present only when the PNS surface is included in the domain—matches the rising high-frequency track in the gravitational-wave spectrogram. The paper concludes that this track is the PNS f-mode, an interface mode produced by the strong density gradient at the PNS surface, not a g-mode.","pith_inferences":["If the HFF really is an interface mode, its frequency should respond sensitively to the accretion rate and the surface density scale height; a test would compare models with different accretion histories at fixed PNS mass.","The classifier's reliance on an approximate total energy means a natural next step is to include the full $\\psi^2$ metric terms for a subset of snapshots; the paper's own numbers suggest the f-mode is exactly where this approximation is least safe.","The unresolved power gap—the narrow band where gravitational-wave emission is suppressed—probably needs interference or cancellation effects beyond the avoided crossing; the four mode families give a concrete setting to test that.","The same restoring-force energy split could be adapted to other multi-cavity oscillators, such as neutron-star merger remnants or accreting white dwarfs, wherever the total mode energy can be computed."],"forward_implications":["If the HFF is the PNS f-mode, its rising track tracks the contracting PNS's surface density gradient, giving a direct asteroseismic handle on PNS mass and radius.","The automatic classifier makes it practical to run the same mode identification over hundreds of simulations, replacing hand-tracked node counts with a consistent label for each mode.","PNS-only eigenvalue calculations misclassify surface g-modes and cannot capture shock f/p modes; calculations must extend to the shock to see the complete spectrum.","The strong avoided crossing between the PNS f-mode and the first core g-mode explains the break in the HFF frequency evolution around 0.4 s after bounce seen in several simulations."],"supporting_citations":[{"why":"supplies the general-relativistic linear perturbation eigenvalue problem and the shooting-method code used to compute the eigenmodes","marker":"[34]"},{"why":"provides the framework neglecting spacetime perturbations and the limiting g-mode and p-mode frequency estimates on which the energy decomposition builds","marker":"[31]"},{"why":"is the earlier study identifying the dominant emission with the fundamental mode, the claim this paper's energy classification confirms and extends","marker":"[33]"},{"why":"established the f-mode frequency scaling with PNS mass and radius that the PNS f-mode identification must match","marker":"[36]"},{"why":"built universal relations from 25 simulations, the systematic large-sample analysis this automatic classifier is designed to enable","marker":"[41]"},{"why":"identified core g-modes and the entropy and lepton gradients shaping the two stably stratified regions that the two g-mode families inhabit","marker":"[46]"},{"why":"describes the PNS surface as an interface where accreted matter excites the mode this paper identifies as the PNS f-mode","marker":"[17]"},{"why":"reports the empirical high-frequency feature in gravitational-wave spectrograms whose physical nature this paper aims to settle","marker":"[19]"}],"fun_headline_variants":["PNS f-mode dominates gravitational-wave ring","Four-mode families from core and shock","Energy-based classifier splits PNS oscillation modes","Neutron star f-mode: the dominant supernova ring","Supernova modes: four families, one dominant f-mode"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the approximate energy fraction $F=(M_2-M_g)/M_2$ (Eqs. 37–38), built from a total energy that drops the $\\psi^2$ metric terms (Eq. 23) and from an approximate force density (Eq. 28), still places the mode that tracks the high-frequency emission on the f+p side of the 0.5 threshold. The paper states that the $\\psi^2$ contribution is small for most modes but is largest for the f-mode, which is exactly the mode the conclusion names.","fun_headline_variants_meta":{"raw":{"variants":["PNS f-mode dominates gravitational-wave ring","Four-mode families from core and shock","Energy-based classifier splits PNS oscillation modes","Neutron star f-mode: the dominant supernova ring","Supernova modes: four families, one dominant f-mode"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000451,"raw_usage":{"total_tokens":2332,"prompt_tokens":1067,"completion_tokens":1265,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":683,"completion_tokens_details":{"reasoning_tokens":1193}},"tokens_in":683,"tokens_out":1265,"duration_ms":12185,"temperature":1.0,"reasoning_tokens":1193,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:11:31.232519+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute $F$ for the reference model at several post-bounce times using the full second-order metric perturbation $\\psi^2$ and the potential-energy form of $M_2$ rather than the kinetic proxy; if the mode that follows the high-frequency gravitational-wave track then has $F<0.5$, or if the lowest-frequency f+p-PNS mode is shown by an independent eigenfunction-tracking analysis to be a p-mode, the central claim fails.","supporting_citations":[],"review_version":1}