{"id":"10d31594-3574-48e0-8d26-423701f5f6b3","arxiv_id":"2412.15033","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Linear response theory is extended to rotating oblate Stäckel clusters, reproducing N-body bending and bar mode growth rates and mapping stability in the flattening-rotation plane.","lead":"This paper computes the instabilities of rotating, flattened stellar clusters using linear response theory, and shows that the predicted growth rates match N-body simulations to about 10%. It maps how bending and bar modes depend on flattening and rotation, a step toward systematic stability analyses of galaxy models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No basis-truncation convergence study is shown; the flattest KK models in Figs. 2 and 4 may not have converged modes, so the claimed 10% agreement with SV97 is not yet established.","rationale":"We focused on the numerical convergence of the response matrix because it is the most direct internal requirement for the paper's headline validation claim. The reader's weakest assumption concerns the match between SV97 equilibria and the Kuzmin-Kutuzov DF; that is external and may be satisfied (the paper's authors likely selected SV97 for that reason), but the paper does not demonstrate it. The basis-truncation issue is unambiguously internal: the paper states that flat clusters need more ℓ harmonics, yet reports no convergence test in ℓmax, nmax, or kmax. Since the paper's own Fig. F1 only covers action-space sampling, a reader cannot tell whether the steeply rising growth rates at small c/a are physical or truncation artifacts. The concrete test with higher truncations would settle this with modest cost, given the public code. We also note that Eq. (13) appears inconsistent with Appendix F.4: the product form omits the sum-of-separable-terms structure and the 1/pu and 1/pv factors. This should be corrected or clarified, but the code check will reveal which expression is implemented. Our recommendation therefore leaves the reader's conditional acceptance unchanged, with the condition sharpened to explicitly require basis-convergence tests for the reported modes.","tokens_in":24457,"tokens_out":15505,"duration_ms":101436,"concrete_test":"Recompute the fastest-growing m=2 mode for the flattest points used in the SV97 comparison (e.g., c/a=0.11, αr=0 and αr=1; and c/a=0.136, αr=0) with successive truncations (ℓmax,nmax,kmax) = (30,10,10), (40,20,20), and (50,30,30), using the public SPOCK code. If the growth rate and pattern speed change by more than ~5–10% between the largest two truncations for any of these cases, the presented results are not converged and the 10% agreement with SV97 is not a robust validation; if they remain stable, the concern is resolved. The same run should be cross-checked against the Appendix F.4 expression for W rather than the printed Eq. (13), which appears inconsistent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim is the ~10% agreement of the computed m=2 growth rates, pattern speeds, and shapes with SV97. This requires the response-matrix modes to be numerically converged. The paper presents only one convergence test, Figure F1, which varies the action-space sampling (nJ, nLz). No test is shown for the truncation parameters of the biorthogonal basis (ℓmax, nmax) or the Fourier harmonic sum (kmax), despite the statement in Section 4.2 that 'as increasingly flattened clusters are considered, a greater number of ℓ harmonics are required to achieve convergence.' The quoted production run uses (ℓmax,nmax,kmax)=(30,10,10). For the flattest cases in the SV97 comparison (c/a≲0.15), the growth rate rises steeply with decreasing c/a; if the basis is not saturated at these truncations, the zero of det(I−M) will shift systematically. A truncated Gram-Schmidt basis generally underestimates the response until saturated, so the published growth rates may be lower limits. Without a per-case convergence check against higher truncations, the 10% agreement—especially at small c/a—is not established. This is an internal, resolvable gap, not a challenge to the SV97 identification.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the linear response (matrix) method to rotating and flattened oblate stellar clusters, using the integrable Kuzmin–Kutuzov Stäckel family as equilibria. The authors construct angle-action variables, build a biorthogonal basis in prolate spheroidal coordinates, compute the response matrix for m=2 modes, and compare their predicted growth rates, pattern speeds, and mode shapes with the N-body simulations of Sellwood & Valluri (1997, SV97). They report ~10% agreement for non-rotating and rotating models, recover the non-rotating stability threshold at c/a=0.208, identify bending and bar-growing instabilities, and provide a public implementation. The paper also includes an analytic-continuation technique to probe damped modes and a discussion of the impact of the sign-function discontinuity in the Lynden-Bell rotation prescription.","tokens_in":24713,"tokens_out":4445,"duration_ms":27588,"significance":"If the numerical convergence and equilibrium-matching issues are resolved, this would be a valuable and timely contribution: it is the first implementation of the linear response matrix for rotating, flattened Stäckel clusters beyond the shell-orbit approximation of Robijn (1995), it provides a systematic stability map in the (flattening, rotation) plane, and it makes the code public, facilitating future applications to more realistic multi-component models. The validation against SV97 is a strong point of the paper: the comparison is an external benchmark, and the growth rates, pattern speeds, and mode shapes are all predicted without fitting to the simulation data. The identification of two distinct instability families (saddle-shaped bending modes and spiral bar-growing modes) is physically interesting and supported by the eigenvector shapes.","major_comments":[{"comment":"The only numerical convergence test shown (Fig. F1) varies the action-space sampling (nJ, nLz) for a single case (a=0.85, ω=0.01i). No test is presented for the basis truncation parameters (ℓmax, nmax) or the harmonic sum kmax, even though Section 4.2 states that more ℓ harmonics are required for increasingly flattened clusters. For the flattest cases in Figs. 2 and 4 (c/a ≲ 0.15), a truncated basis will underestimate the response until saturation, so the reported growth rates may be lower limits. The claimed ~10% agreement with SV97 is therefore not established for those cases unless a per-case convergence study against higher (ℓmax, nmax, kmax) is provided.","section":"Section 4.2 and Appendix F.5, Fig. F1"},{"comment":"The validation against SV97 assumes that the N-body equilibria are well represented by the two-integral Kuzmin–Kutuzov DF of Eq. (27) with the Lynden-Bell rotation prescription of Eq. (29), which fixes σ_R=σ_z by construction. The paper does not compare the phase-space structure (e.g., the anisotropy parameter or velocity moments) of the SV97 models with this DF. Without such a comparison, the ~10% agreement could be coincidental. The authors should compute, for instance, the projected velocity dispersions of their DF and compare them with the SV97 initial conditions, or state explicitly that the comparison tests only the response matrix in a KK-like equilibrium.","section":"Section 4.2, Eqs. (27) and (29)"},{"comment":"The '10% agreement' claim is made without a quantitative definition or error bars on either the linear-theory mode frequencies or the SV97 measurements. Since the mode frequencies are obtained by locating zeros of det(I−M) on a discrete complex-frequency grid, the grid resolution and the zero-finding/interpolation procedure should be stated, and at least approximate uncertainties on ω0 and γ should be given. As written, the claim is not falsifiable from the figures alone, particularly for the pattern speeds in the bottom panel of Fig. 4.","section":"Figs. 2 and 4 and Section 6"}],"minor_comments":[{"comment":"Typos: 'useful per say' should be 'useful per se' (Introduction); 'let use define' should be 'let us define' (beginning of Section 4.1).","section":"Section 1 and Section 4.1"},{"comment":"'bashed red' should be 'dashed red'; also the caption gives 'c/a = 0.2' while the text repeatedly states the threshold as c/a = 0.208, so the caption should be made consistent.","section":"Figure 6 caption"},{"comment":"Captions referring to 'a = 0.88' and 'a = 0.9' are ambiguous: a is a potential parameter and the flattening ratio is c/a. Since c=1−a in the adopted units, please state the corresponding c/a value explicitly in each caption.","section":"Figure captions (Figs. 1, 3, 7, H1)"},{"comment":"The text states that λ_r is proportional to α_r for fixed c/a, but the proportionality constant is not given. For reproducibility, provide the conversion factor (or a small table) between α_r and λ_r for the models shown in Figs. 5 and 6.","section":"Section 5.2, Eq. (31)"},{"comment":"The analytic continuation via Padé approximants (Eq. I75) is used to locate damped modes in Fig. 8, but the sensitivity of the mode location to the number of sampling points and the approximation order is not discussed. A short convergence test for the specific case in Fig. 8 would strengthen the claim of a damped-to-unstable transition.","section":"Appendix I and Section 5.4"}],"recommendation":"major_revision","confidential_remarks":"This is a technically ambitious paper that addresses a long-standing gap (linear response of rotating, flattened Stäckel clusters) and validates against an external N-body benchmark. The main weakness is the absence of basis-truncation convergence tests, which is directly load-bearing for the central 10% agreement claim. The issue is resolvable within the manuscript's scope: the authors should add per-case convergence checks for (ℓmax, nmax, kmax) and, ideally, a comparison of the KK DF with the SV97 equilibria. I also note that the paper's wording in the conclusions ('converged m=2 modes') is stronger than what the single convergence test in Fig. F1 supports; the revised version should temper or substantiate that claim. The code release is commendable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Kerwann, here's my read on arXiv:2412.15033.\n\nThe headline is that this is the first full three-dimensional linear response computation for rotating flattened stellar clusters, and it lands on something real: the m=2 stability map for the Kuzmin-Kutuzov family, including the bending-to-bar transition and bar modes that persist in very round rotators. The derivation is careful, the appendices are substantive, and the code is public. The validation against SV97 is a good instinct, and where it works, the 10% agreement is meaningful. The authors also handle the discontinuous rotation DF honestly, showing in Appendix K that the modes are not edge modes.\n\nNow the soft spots. The convergence evidence is incomplete. Figure F1 only varies the action-space sampling; there is no test that the basis truncation (ℓmax, nmax) or the harmonic sum (kmax) are saturated for the flattest cases in Figures 2 and 4. The paper itself says flatter clusters need more ℓ harmonics, so this is not a nitpick. A truncated Gram-Schmidt basis will underestimate the response until it is saturated, which means the published growth rates for c/a≲0.15 could be systematically low. If so, the claimed 10% agreement with SV97 at small flattening is not yet established. That's a fixable gap, but it is a load-bearing one for the central validation claim. A per-case convergence check against higher truncations would settle it.\n\nAlso, the equivalence between the SV97 equilibria and this two-integral Kuzmin-Kutuzov DF is assumed, not demonstrated. No phase-space comparison. That could make the 10% agreement partly coincidental. Minor, but worth stating.\n\nNeither of these flaws kills the paper. The method is sound, the new results are genuinely new, and the code is a public good. The right audience is anyone doing stability analysis of flattened stellar systems or linear response with Stackel potentials. I'd send this to a serious referee, with a request to address the convergence gap before publication.\n\nIn short: worth refereeing, needs a revision, not a desk reject.","headline":"A serious and mostly solid 3D linear response computation for rotating flattened stellar clusters; the missing basis-convergence checks make the 10% SV97 agreement provisional, not established.","tokens_in":25214,"tokens_out":3271,"would_cite":true,"duration_ms":29205,"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":"Linear theory reproduces N-body galaxy instabilities to 10 percent.","keywords":["linear response theory","Stäckel potentials","Kuzmin-Kutuzov cluster","galaxy stability","bar instability","bending instability","rotating stellar systems","action-angle variables"],"falsifier":"Construct N-body equilibria from the exact Kuzmin-Kutuzov distribution function, run them, and measure the fastest m=2 mode; if its growth rate and pattern speed differ from the matrix-method prediction by much more than the quoted 10 percent, the agreement with SV97 was an artifact of the equilibrium assumption rather than a vindication of the linear response calculation.","tokens_in":24282,"feed_emoji":"🌌","tokens_out":7804,"duration_ms":50033,"temperature":0.7,"pith_summary":"This paper extends linear response theory to flattened, rotating stellar systems by working with the integrable Kuzmin-Kutuzov Stäckel family, which spans shapes from a spherical isochrone to a cold Toomre disk. It aims to show that a matrix calculation in angle-action variables predicts the same growing bi-symmetric modes that N-body simulations reveal, with growth rates, pattern speeds, and mode shapes matching to about 10 percent. The result matters because stability thresholds and mode shapes determine whether a galaxy forms bars or buckles, and such questions have mostly been studied one N-body run at a time. If the method is sound, the stability of an entire family of oblate rotating clusters can be computed systematically, not simulated case by case.","feed_headline":"Bar and bending modes of star clusters reproduced to 10%","feed_subtitle":"A matrix method maps where flattened, rotating stellar systems go unstable, matching more costly N-body runs.","key_machinery":"The central object is the response matrix M(ω) in angle-action space, whose determinant zeros give the system's linear normal modes. The calculation uses prolate spheroidal coordinates for the oblate Stäckel potential, the analytically known two-integral Kuzmin-Kutuzov distribution function, and a bi-orthogonal basis of spheroidal harmonics; rotation is added through the Lynden-Bell daemon prescription, which multiplies the distribution function by 1 + α_r sgn(L_z). The matrix elements are built from Fourier-transformed basis functions and three-dimensional action-space integrals, and analytic continuation lets the authors follow modes across the stability transition into damped territory.","core_discovery":"The paper's central claim is that the converged m=2 eigenmodes of flattened and rotating Kuzmin-Kutuzov spheroids, computed from the linear response matrix, match the growing modes measured by Sellwood & Valluri (1997) in N-body simulations: growth rates, pattern speeds, and mode shapes agree at the roughly 10 percent level, and the non-rotating stability threshold at c/a=0.208 is reproduced. Rotation changes the character of the dominant instability, producing two families: saddle-shaped bending modes at low spin and spiral-shaped bar-growing modes at high spin. The bar-growing modes persist to much rounder shapes than the bending modes, with maximally rotating clusters stabilizing only at the spherical limit. The paper also maps the stability boundaries in the flattening-rotation plane and shows that the least unstable rotating configurations become flatter as rotation increases.","pith_inferences":["A direct resimulation of the exact Kuzmin-Kutuzov distribution function would decouple the linear-response check from the fidelity of the SV97 equilibria; if the mode match survives, the matrix method stands on its own for axisymmetric systems.","The bending-to-bar transition suggests that the dominant m=2 instability may serve as a dynamical discriminant between slow and fast rotators, a connection the paper does not pursue observationally.","Extending the same machinery to m=0 and m=1 modes, once inertial pseudo-forces are included, would complete the stability portrait and may reveal additional axisymmetric or lop-sided instabilities in these clusters.","The quoted computational cost, about 38 hours on 512 cores per cluster, means systematic parameter surveys are feasible; one could map the full flattening-spin plane at higher resolution, including the multiple bar modes the paper identifies."],"forward_implications":["For any model in the Kuzmin-Kutuzov family, the growth rate, pattern speed, and shape of the dominant m=2 instability can now be computed directly, without running an N-body simulation.","The stability map in the flattening-spin plane separates a bending-mode regime from a bar-growing-mode regime, with the bar mode persisting to rounder shapes up to the sphere at maximal rotation.","The least unstable rotating models become flatter as spin increases, giving a concrete, testable trend for observed elliptical and spheroidal galaxies.","Because the formalism is built on Stäckel separability, it can be extended to multi-component Stäckel discs and halos, moving linear stability analysis closer to realistic galaxy models.","The public code allows independent groups to reproduce the stability maps and to compute modes for their own Stäckel-based equilibria."],"supporting_citations":[{"why":"It supplies the N-body growth rates, pattern speeds, and mode shapes that the linear calculation reproduces.","marker":"SV97"},{"why":"It provides the two-integral Kuzmin-Kutuzov distribution function used in equation (27) and its connecting limits to the isochrone and Toomre disk.","marker":"DZ88"},{"why":"It derives angle-action coordinates for oblate Stäckel potentials, the framework on which the whole calculation rests.","marker":"de Zeeuw et al. 1986"},{"why":"It gives the prolate spheroidal bi-orthogonal basis elements used to expand the potential and density perturbations.","marker":"Robijn & Earn 1996"},{"why":"It introduced the matrix method for the linear response of stellar systems, which this paper extends to rotating oblate clusters.","marker":"Kalnajs 1977"},{"why":"It supplies the daemon rotation prescription F(1 + α_r sgn L_z) used to introduce rotation into the distribution function.","marker":"Lynden-Bell 1960"},{"why":"It provides the angle-variable construction used here to compute orbital angles for Stäckel potentials.","marker":"Binney 2012"},{"why":"It contributes the analytic-continuation technique used to identify damped bar modes across the stability transition.","marker":"Fouvry & Prunet 2022"}],"fun_headline_variants":["Star cluster instabilities matched to 10% with matrix method","Two instabilities shape rotating star clusters: bending and bar modes","Rotation flattens star clusters, new linear response maps it","Bar and bending modes in stellar clusters: linear theory matches N-body","Flattened star clusters: linear response pinpoints instabilities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison assumes that the SV97 N-body clusters are faithfully described by the two-integral Kuzmin-Kutuzov distribution function of equation (27) with the Lynden-Bell rotation prescription of equation (29), but the paper never checks the phase-space structure of those simulations against this distribution function.","fun_headline_variants_meta":{"raw":{"variants":["Star cluster instabilities matched to 10% with matrix method","Two instabilities shape rotating star clusters: bending and bar modes","Rotation flattens star clusters, new linear response maps it","Bar and bending modes in stellar clusters: linear theory matches N-body","Flattened star clusters: linear response pinpoints instabilities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000171,"raw_usage":{"total_tokens":1241,"prompt_tokens":883,"completion_tokens":358,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":499,"completion_tokens_details":{"reasoning_tokens":271}},"tokens_in":499,"tokens_out":358,"duration_ms":4095,"temperature":1.0,"reasoning_tokens":271,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:41:26.331752+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Construct N-body equilibria from the exact Kuzmin-Kutuzov distribution function, run them, and measure the fastest m=2 mode; if its growth rate and pattern speed differ from the matrix-method prediction by much more than the quoted 10 percent, the agreement with SV97 was an artifact of the equilibrium assumption rather than a vindication of the linear response calculation.","supporting_citations":[{"cited_title":"1986, MNRAS, 221, 1001, 10.1093/mnras/221.4.1001","cited_arxiv_id":null,"evidence_quote":"It derives angle-action coordinates for oblate Stäckel potentials, the framework on which the whole calculation rests."},{"cited_title":"1960, , 120, 204, 10.1093/mnras/120.3.204","cited_arxiv_id":null,"evidence_quote":"It supplies the daemon rotation prescription F(1 + α_r sgn L_z) used to introduce rotation into the distribution function."},{"cited_title":"2022, , 509, 2443, 10.1093/mnras/stab3020","cited_arxiv_id":null,"evidence_quote":"It contributes the analytic-continuation technique used to identify damped bar modes across the stability transition."}],"review_version":1}