{"id":"51475234-2489-44ed-8cdb-73e24fe4495b","arxiv_id":"2607.24530","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"For liquid Lane–Emden stars, the number of radial growing modes equals the negative Morse index and changes only at mass–radius turning points, with orientation dictating gain or loss.","lead":"A rigorous turning-point law is proved for Newtonian liquid Lane–Emden stars: radial growing modes equal the Morse index of a free-boundary operator and jump only at mass extrema of the mass–radius curve. The result organises liquid-star stability the way the classical static criterion organises neutron-star models, and ties the large-density mode count to a planar spiral/node dichotomy.","discovery_kind":"extension","skeptic_critique":{"model":"moonshotai/kimi-k3","headline":"The local jump law is self-contained and checks out; the global spiral/dichotomy claims rest entirely on the phase portrait of (3.27) imported from the author's prior work [36], exactly as the reader flagged.","rationale":"The reader's weakest_assumption identifies exactly the joint where the strongest claim is least secure: the dichotomy/onset statements inherit their correctness wholesale from the tail analysis of [36]. I verified that the parts proven in this manuscript do not hide an additional soft spot: the marginal-mode identity and its constant cκ recompute correctly; the jump-direction argument is rigorous; the min–max continuity, kernel one-dimensionality (indicial exponents 0 and −d with the Lκ-integrability cutoff), and winding bookkeeping are standard and correctly executed. Corollary 2.10's reduction to [37] is a second external dependence, but it is confined to the non-radial corollary, not the central claim. Because the concern is (a) disclosed by the author, (b) correctly identified by the reader, (c) confined to the global geometric statements, and (d) cheaply testable by independent integration of a planar ODE, it does not warrant downgrading the verdict: ACCEPT stands, with the [36]-dependence understood as the residual risk. If the concrete test revealed a failure in [36]'s phase portrait, parts (ii)(b) and (iii) would need retraction while Theorems 2.7 and 2.8(i) would survive — but there is no internal evidence here suggesting such a failure, and the focus/node discriminant computation (Proposition 3.7, Step 1) and the node-case non-orthogonality (3.37) proven in this paper are correct as far as I can verify.","tokens_in":38962,"tokens_out":4477,"duration_ms":155114,"concrete_test":"Numerically integrate the planar system (3.27) (a two-dimensional autonomous ODE, reproducible in an hour) for a grid of (γ,d) covering γ<γ♯ at d=3,...,12 and γ∈[γ♯,γ*) . For each run check: (i) the orbit from the regular centre converges to v∗ with an exponential rate (validating the Poincaré–Bendixson/Bendixson–Dulac stability claim imported from [36]); (ii) winding of v(τ)−v∗ occurs exactly when D(γ,d)<0 and monotone convergence when D≥0; (iii) support dichotomy at γ=γ♯. If any case with γ<γ♯, D<0 fails to converge to v∗, Theorem 2.8(ii)(b) and (iii)(a) collapse in that regime; if all confirm, the sole external input to the global claims is independently corroborated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I stress-tested the central claim (Theorem 2.8) at its three load-bearing joints. (1) The local jump law (Lemma 3.4) is internally solid: I recomputed the marginal-mode boundary defect (Lemma 3.3, Steps 2–3): from dν(R)+R∂yν(R) = (∂yρ̄(R)/(4πR^{d−1}))∂κM and ∂yρ̄(R)=−M/(γR^{d−1}), the pairing constant is cκ = −M/(4πR^{d−2}) < 0 as stated, and the crossing-direction argument in Step 5 (sign µ = −sign(M′R′) via the limit of µ(κj)/M′(κj) along eigenfunctions converging to tνκ0) is a correct instance of the standard simple-eigenvalue crossing lemma, including the honest treatment of the sign ambiguity t. The eigenvalue-continuity reduction to the fixed pencil (Qκ,ℓκ) and the R′≠0-at-critical-points ODE-uniqueness argument are sound. (2) The parity argument in (iii)(c) and the winding formula (3.25) are bookkeeping-correct. (3) Everything global — the existence of the first mass maximum (ii)(b), the spiral with n_u→∞ (iii)(a), the node convergence (iii)(b) — is transcribed from the phase portrait of (3.27): rest point v∗, its exponential asymptotic stability for γ<γ♯, and the convergence of the gaseous orbit, all explicitly attributed to [36] (\"established in [36] (proven there by a Poincaré–Bendixson and Bendixson–Dulac argument)\"), plus the support dichotomy γ≷γ♯, also from [36]. This paper supplies no independent verification of that input. If [36]'s tail analysis were incomplete for some (γ,d) — e.g., an overlooked periodic orbit in the region where Dulac's criterion is applied, or a failure of exponential stability at some parameter — the onset statement and the automatic D<0 conclusion for d<10 would fail, while the local jump law (Lemma 3.4), the Morse-count identity (Theorem 2.7), and the explicit γ=γ♯ case (Lemma 3.8(ii), which is self-contained) would survive. This is a genuine external dependence, not circularity, and it is disclosed; it is nonetheless the single point where the strongest claim is least secure. Separately, the open node case (d≥10, D≥0) is a承载","agreement_with_reader":"agree"},"referee_report":{"model":"moonshotai/kimi-k3","summary":"The paper studies the one-parameter family (parametrised by central density κ∈(1,∞)) of spherically symmetric steady states of the free-boundary Euler–Poisson system with liquid equation of state p=ρ^γ−1. Theorem 2.7 identifies the number of radial growing modes with the negative Morse index of the linearised radial operator Lκ — the comoving formulation builds in mass conservation, so no constraint or winding-index subtraction is needed — and shows ker Lκ is non-trivial exactly at critical points of κ↦Mκ, spanned by the marginal mode νκ. Theorem 2.8 is a turning point principle: n_u(κ) is locally constant off mass critical points and jumps by ±1 across nondegenerate ones according to the sign change of (∂κM)(∂κR); the paper further proves stability for all κ when γ≥γ*=2(d−1)/d, onset of instability at the first mass maximum for γ<γ* (under γ≥γ♯ or D(γ,d)<0), and a large-κ spiral/node dichotomy governed by the explicit discriminant (1.8), with n_u→∞ along the spiral. Corollaries recover and sharpen the radial (in)stability results of [36] and, via [37], give a non-radial stability criterion.","tokens_in":39348,"tokens_out":7097,"duration_ms":57480,"significance":"This is, to my knowledge, the first complete and rigorous turning point principle for a Newtonian free-boundary stellar family, and a faithful Newtonian analogue of the Hadžić–Lin–Rein relativistic theory, including the mass–radius spiral and n_u→∞. Particular strengths: the local theory (Lemmas 3.1–3.5) is fully self-contained and variational, with the marginal-mode identity (2.19) carrying the explicit nonzero constant cκ=−Mκ/(4πR^{d−2}); the comoving reduction yields a plain Morse-index count with no winding subtraction, a genuine bookkeeping simplification; the spiral/node dichotomy is decided by an explicit parameter-free discriminant with a clean dimension-ten threshold; the borderline case γ=γ♯ is exactly solvable (κ1 given in closed form, stabilised count 1), supplying falsifiable predictions; and the open node case d≥10 is honestly delineated (Remark 2.8.3) rather than glossed over. The result also sharpens prior work (large-κ instability for d≥10, γ∈[γ♯,γ*) is new relative to [36]).","major_comments":[],"minor_comments":[{"comment":"Step 5, displayed limit for µ(κj)/M′(κj): the factor ℓκ0[˜νκ0,˜νκ0] appears in the numerator, but dividing (3.22) by M′ and passing to the limit gives µ/M′ → cκ0 R′(κ0)/(R^{d+3}_{κ0} ℓκ0[˜νκ0,˜νκ0]) — the ℓ-factor should be in the denominator. The conclusion is unaffected (ℓ>0, so the limit is still a nonzero number of sign opposite to R′(κ0)), but the displayed formula should be corrected.","section":"Lemma 3.4"},{"comment":"Step 1 imports the phase portrait of (3.27) from [36] (rest point v*, exponential asymptotic stability for γ<γ♯, orbit convergence, support dichotomy; also the explicit γ=γ♯ profile used in Lemma 3.8(ii)). Since Theorem 2.8(ii)(b)–(iii) hinge entirely on these inputs, please give theorem-level references within [36] rather than a global citation, so the dependence is directly checkable.","section":"Proposition 3.7"},{"comment":"The proof of Corollary 2.10 relies on the spherical-harmonic block decomposition and non-negativity of the l≥1 blocks from [37], which at present is a preprint by the same author. Please state precisely which statements of [37] are used and note its publication status, since this corollary is one of the headline applications.","section":"Corollary 2.10"},{"comment":"The phrase 'creating extrema of the mass when γ<2(d−1)/d' (also in §1.3) is unconditional, but per Remark 2.8.3 the existence of a mass maximum is not guaranteed in the node case d≥10, D(γ,d)≥0. A brief qualification would avoid overstating the scope.","section":"Abstract"},{"comment":"The extension to general liquid equations of state is asserted ('can be used to prove', tails 'would match') rather than shown. Consider softening the language or labelling it as a programme, since the tail analysis for asymptotically polytropic P is not carried out here.","section":"Remark 2.8.4"},{"comment":"Real-analyticity at the centre is justified via a 'classical majorant argument' with a citation to [28], which gives the series expansions; a reference to a convergence result (e.g., Hukuhara-type theory for the regular singular point) would make this step fully self-supporting. Also, notation 'f0' vs 'f0,κ' is used somewhat loosely in §2.1.","section":"Lemma 3.2"}],"recommendation":"accept","confidential_remarks":"The manuscript depends on the author's own prior work at two key joints: [36] (published, Quart. Appl. Math. 2024) supplies the entire phase-portrait input for the global statements, and [37] (currently a preprint, arXiv:2603.03548) supplies the non-radial block analysis behind Corollary 2.10. This is legitimate and transparently disclosed, and the main theorems (2.7, 2.8) are independent of [37]; still, the editor may wish to confirm the status of [37] since one of the advertised corollaries inherits it. The acknowledgements disclose AI assistance with language; the journal may wish to verify this conforms to its policy. On the science: the local turning point law is proven in full and I verified the crossing-direction argument independently; the paper is careful, complete, and unusually honest about its one open case. I recommend acceptance."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real news is a clean radial turning-point principle for liquid (stiffened-gas) Lane–Emden stars: nu equals the negative Morse index of the comoving operator, stays constant off mass extrema, and jumps by ±1 with the bend of the mass–radius curve. That is the Zel’dovich–Wheeler static criterion made rigorous in a free-boundary Newtonian setting where the surface term actually creates the extrema.\n\nWhat is new is not the slogan but the free-boundary bookkeeping. The liquid surface produces a Robin condition and a boundary term in the energy form; the comoving mass-preserving formulation then lets the growing-mode count be plain n−(Lκ) with no mean-zero constraint and no winding-index subtraction. Lemma 3.3 (marginal mode as kernel exactly at mass critical points) and Lemma 3.4 (simple eigenvalue crossing with sign µ = −sign(M′R′)) are the load-bearing pieces, and they check out by standard variational/ODE arguments. The winding form and the explicit γ = γ♯ case are clean. Combined with the author’s non-radial work, you also get a d=3 criterion that all growing modes are radial—useful and not available relativistically.\n\nSoft spots in proportion: everything about onset, the spiral with nu → ∞, and the automatic D < 0 claim for d < 10 is transcribed from the phase portrait of the gaseous-tail system in [36]. This paper does not re-prove Poincaré–Bendixson/Dulac or exponential attraction to v∗. That is disclosed external dependence, not circularity; the local jump law and Morse identification survive even if some tail detail in [36] were incomplete. The node regime d ≥ 10, D ≥ 0 is left open, again honestly. Citation pattern is appropriate: Lin–Zeng, Hadžić–Lin–Rein, and the author’s own profile/non-radial papers are the right inputs.\n\nThis is for people who work on free-boundary Euler–Poisson or stellar spectral stability. Math is standard and written out; no data, no free parameters. I would send it to referees without hesitation. Engage if you care about turning-point organisation or liquid stars; skip only if you are purely gaseous/relativistic and already saturated on the analogy.","headline":"Solid turning-point theorem for liquid Lane–Emden stars; local jump law is self-contained, global spiral claims ride on the author’s prior planar-tail analysis.","tokens_in":39417,"tokens_out":574,"would_cite":true,"duration_ms":13937,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q35","85A15","76E20","35B35"],"pacs":[],"model":"grok-4.5","headline":"For liquid Lane–Emden stars, the number of radial growing modes equals the negative Morse index and jumps only at mass extrema of the mass–radius curve.","keywords":["Lane-Emden stars","Euler-Poisson system","turning point principle","mass-radius curve","liquid free boundary","radial stability","growing modes","Morse index"],"falsifier":"For a fixed γ below the mass-critical index and dimension three, numerically continue the mass–radius curve through the first mass maximum and compute the lowest eigenvalue of L_κ on either side: the eigenvalue must cross from positive to negative exactly there, and further crossings must track later mass extrema (or accumulate if the curve spirals).","tokens_in":38963,"feed_emoji":"⭐","tokens_out":992,"duration_ms":29048,"temperature":0.7,"pith_summary":"This paper proves a turning point principle for spherically symmetric liquid stars modeled by the free-boundary Euler–Poisson system with stiffened polytropic pressure. Fixing the adiabatic index, the steady states form a one-parameter family labeled by central density, and the liquid surface breaks the self-similarity that makes gaseous stars scale-invariant. The result is that the count of radially unstable modes is locally constant along the family and can change only where the mass has an extremum; the orientation of the bend on the mass–radius curve decides whether a growing mode is gained or lost. At large central density a planar dynamical system for the gaseous tail decides whether that curve spirals—driving the mode count to infinity—or settles with finitely many turns. The law recovers and sharpens known radial stability thresholds and, with existing non-radial work, yields a mass–radius criterion for full linear stability.","feed_headline":"Liquid stars gain unstable modes only at mass turning points","feed_subtitle":"Bends in the mass–radius curve count radial growing modes; spirals drive the count to infinity","key_machinery":"The reduced radial operator L_κ in comoving coordinates, whose energy form includes a liquid-surface boundary term; the growing-mode count is its negative Morse index, and the marginal mode along the family detects mass turning points via a boundary identity linking L_κ ν_κ to ∂_κ M_κ.","core_discovery":"Along the liquid Lane–Emden family, the number of radial growing modes equals the negative Morse index of the linearised radial operator, stays constant between mass critical points, and changes by exactly the jump of a turning index at each nondegenerate mass extremum, according to the sign change of the product of the mass and radius derivatives. In the infinite-support spiral regime the mode count tends to infinity; otherwise it stabilises to a finite value, odd and at least one when the curve returns to the origin.","pith_inferences":["The open node regime in dimensions ten and higher is the natural next target: a single sign computation along the attracting eigendirection would close the large-density stability question there.","The same comoving mass-preserving reduction that eliminates the winding index may simplify turning-point arguments for other free-boundary stellar models with a density jump.","Because the liquid surface term is what creates the turning points, any equation of state that retains a fixed surface density and an asymptotically polytropic tail should inherit an analogous mass–radius counting law."],"forward_implications":["When γ is at least the mass-critical index, every liquid Lane–Emden star is radially linearly stable, with no mass turning points.","When γ is below that index (and the discriminant condition holds, automatic in dimensions below ten), stars are stable up to the first mass maximum and unstable immediately beyond it and at large central density.","In the spiral regime the number of radial growing modes tends to infinity with central density.","In three dimensions every exponentially growing mode is radial, so the same mass–radius turning-point rule governs full linear stability.","The precise mode count at any non-critical central density equals the net number of counter-clockwise horizontal crossings of the mass–radius tangent up to that point."],"fun_headline_variants":["Liquid stars gain radial modes only at mass turning points","Mass-radius bends set the count of growing radial modes","Mode number jumps solely at mass extrema of liquid stars","Spiral mass-radius curves drive unbound growth of modes","Turning index at mass peaks dictates mode gain or loss"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"The large-density spiral-versus-node picture and the onset of the first mass maximum rest on the phase portrait of a planar system for the gaseous tail already established in earlier work on the same family.","fun_headline_variants_meta":{"raw":{"variants":["Liquid stars gain radial modes only at mass turning points","Mass-radius bends set the count of growing radial modes","Mode number jumps solely at mass extrema of liquid stars","Spiral mass-radius curves drive unbound growth of modes","Turning index at mass peaks dictates mode gain or loss"]},"model":"grok-4.5","effort":"low","cost_usd":0.004161,"raw_usage":{"total_tokens":1334,"prompt_tokens":852,"num_sources_used":0,"completion_tokens":61,"cost_in_usd_ticks":41608000,"prompt_tokens_details":{"text_tokens":852,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":421,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":852,"tokens_out":61,"duration_ms":6916,"temperature":1.0,"reasoning_tokens":421,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T12:18:07.674738+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"For a fixed γ below the mass-critical index and dimension three, numerically continue the mass–radius curve through the first mass maximum and compute the lowest eigenvalue of L_κ on either side: the eigenvalue must cross from positive to negative exactly there, and further crossings must track later mass extrema (or accumulate if the curve spirals).","supporting_citations":[],"review_version":1}