{"id":"2a0f21ef-34fb-4f08-934b-41d45b7729b0","arxiv_id":"2502.05304","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A new multi-component stability parameter Q_delta = 5/3 maps disc swing amplification more consistently than older Q definitions and requires gas Q to be rescaled close to Q^(5/3).","lead":"The paper argues that the standard Toomre Q parameter, built for axisymmetric stability, is misused as a measure of spiral activity, and proposes a redefinition in which gas disc Q is rescaled close to its square. For the Solar Neighbourhood this gives Q = 1.58, matching simulation expectations better than previous multi-component definitions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Q_δ's one-to-one map to swing amplification is calibrated and validated with the same GLB model, which the paper itself calls only 'indicative'; no independent test links it to real spiral activity.","rationale":"The paper's central contribution is the claim that Q_δ=5/3 gives a close to one-to-one map between Q and spiral activity (maximum swing amplification) across discs with different gas fractions, and should therefore replace Q_Rk as the standard. For this to be true, the quantity used to define 'spiral activity' must be a reliable predictor of real spiral growth. The GLB maximum amplification is a local, linear, shearing-sheet quantity whose own inventors and the present authors describe as only 'indicative'. The paper adduces one cross-check (Fig. 2) between GLB and the more complete JT formalism, but only for single-component stellar discs; the multi-component extension (Eq. 39) is made by inserting the Rafikov WKB dispersion relation (Eq. 10), which is formally valid only for tightly wound waves, into a regime where the wave passes through nearly unwound configurations. No independent calculation or simulation verifies that the GLB maximum for gas-star mixtures tracks the amplitude of spirals in real discs. Thus the calibration of δ=5/3 is a fit to an unvalidated model. This does not mean the paper is wrong: the basic point that gas Q scales roughly quadratically may be robust, and the authors are transparent about the approximations. But the specific recommendation to adopt Q_δ as the new standard, and the quantitative predictions (Q≈1.58 for the Solar Neighbourhood, gas-driven m≳10 patterns in the Milky Way and M74), hang on this premise. The proposed simulation test would settle whether the composition-induced scatter in spiral activity at fixed Q_δ is truly smaller than at fixed Q_Rk. If it is, the CONDITIONAL verdict can be upgraded; if not, the central claim reduces to an internal reparametrization. The reader's weakest assumption identifies exactly this premise, so I agree; the conditional verdict is the appropriate outcome pending such validation.","tokens_in":26726,"tokens_out":14070,"duration_ms":135769,"concrete_test":"Run a suite of local shearing-box simulations (with collisionless star particles plus isothermal gas) spanning the same grid of gas fractions and velocity ratios as Fig. 10, holding Q_δ fixed (e.g., 1.5) by adjusting surface densities and velocity dispersions; measure the resulting spiral amplitude (e.g., maximum Fourier density contrast or m-integrated power). If the composition-induced scatter in spiral amplitude at fixed Q_δ is comparable to that at fixed Q_Rk, the one-to-one premise fails. As a secondary check, compute the linear maximum amplification for the same two-component discs with a full local-shearing-sheet treatment (JT-type for stars, fluid equations for gas) without the WKB quasi-static approximation for gas-dominated cases, and compare to the GLB curves; systematic deviations would indicate that the extension in Eq. 39 is the weak link.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The single load-bearing premise is that the maximum swing amplification factor from the GLB local shearing-sheet calculation, extended to multi-component discs via the Rafikov (2001) WKB dispersion relation (Eq. 10), is the correct physical invariant for spiral activity. This premise is not externally validated. The paper's own Section 3 states that the GLB calculation 'cannot be used to describe the evolution of real spiral patterns in terms of swing amplification', and Section 1 calls it only a 'useful indicative measure'. The extension to multiple components (Eqs. 38–39) applies the WKB dispersion relation where waves are not tightly wound (near γ=0), and the only cross-check against the more complete JT formalism (Fig. 2) is for single-component stellar discs, not for the gas-rich discs where the proposed redefinition changes Q most. Moreover, the calibration of δ=5/3 (Fig. 10) minimizes ΔGLB, a quantity defined using the same GLB model, so the resulting one-to-one map is an internal consistency property rather than a demonstrated match to real spiral growth. If the GLB maximum does not track real spiral amplitude, Q_δ is a reparametrization of an approximate model, and the recommendation to replace Q_Rk as the standard collapses. The same premise also underlies the application sections, where the GLB amplification-vs-m curves (Figs. 14, 17, 18) are used to predict dominant spiral multiplicity for the Milky Way and M74.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that Toomre's Q conflates two distinct notions: the axisymmetric stability boundary at Q=1 and the susceptibility to spiral activity for Q>1. It proposes redefinitions that preserve Q=1 but give a more nearly one-to-one map between Q and the maximum swing-amplification factor computed from the Goldreich-Lynden-Bell (GLB) shearing-sheet formalism. The first proposal, Q_dr (Eq. 52), equals the square of the classical gas Q and is close to the classical stellar Q. A refined hybrid parameter Q_delta (Eq. 64) is then introduced, with an exponent delta=5/3 chosen by minimizing the composition scatter of the GLB-based measure Delta_GLB (Eq. 46). The paper provides an iterative formula and code, applies the new parameter to the Solar Neighbourhood (Q≈1.58) and to M74, and predicts that high-m, gas-dominated instabilities produce flocculent ISM structure and Milky Way spurs.","tokens_in":27015,"tokens_out":5948,"duration_ms":64212,"significance":"If the proposed redefinition is correct, it would resolve a known inconsistency in the use of Q for star-gas discs and would provide a practical replacement for Q_Rk, with quantitative consequences for Milky Way and M74 spiral structure. The paper contains a clean analytic identity (Eqs. 50-52) showing Q_dr=Q_g^2 for a single gas disc, a transparent numerical comparison of several Q definitions, and publicly available code for the iterative computation. The authors are also candid about the limitations of the GLB formalism and about the lack of error propagation in the applications. However, the headline recommendation Q_delta=5/3 is calibrated and validated with the same approximate GLB model that defines the criterion, and no independent test connects the resulting parameter to real spiral amplitudes in simulations or observations. The significance is therefore real but conditional on an external validation that the current manuscript does not provide.","major_comments":[{"comment":"The exponent delta=5/3 is selected by minimizing the spread of Delta_GLB (Eq. 46), and the same Delta_GLB is then used to demonstrate the improvement of Q_delta over previous definitions. This is an internal-consistency calibration, not an independent validation. Since Sections 1 and 3 state that the GLB calculation 'cannot be used to describe the evolution of real spiral patterns' and is only a 'useful indicative measure', the central claim that Q_delta=5/3 is the 'correct' definition for spiral activity requires an external test, e.g. against N-body/gas simulations or against the Julian-Toomre formalism extended to gas-rich discs. Without such a test, the quantitative recommendation rests entirely on the GLB model the paper itself calls approximate.","section":"Sec. 4.2, Eq. (64), Fig. 10"},{"comment":"Requirement (ii) demands a decreasing one-to-one map between Q and maximum swing amplification, but Fig. 9 shows that the maximum amplification flattens for Q_dr≳3, and the text admits that flat gradients allow discs of different composition with very different Q to share the same maximum amplification. Thus the proposed definition fails its own one-to-one criterion in the high-Q regime. Moreover, the recommended Q_delta is never subjected to the Delta_JT test that is used for Q_dr in Fig. 8, so it is unknown whether Q_delta satisfies requirement (ii) for multi-component stellar discs. The paper should either demonstrate one-to-one behaviour for Q_delta over the range it recommends, or explicitly restrict the claimed validity to Q≲3.","section":"Sec. 4.1, Figs. 8-9; requirement (ii) in Sec. 2.2"},{"comment":"The comparison with the Aumer et al. (2016) simulation value Q≈1.5 is not a strong validation because the authors state 'Without full error propagation' and vary only the gas surface density, ignoring the quoted uncertainties in kappa (41±2.5 km/s/kpc) and stellar surface density (38±4 M_sun/pc^2). The resulting range Q_delta=1.06-2.13 completely brackets both the new and old central values, so the claim that Q≈1.58 is 'closer to results from simulations' is not robust. A full error estimate should be provided, or the claim should be weakened accordingly.","section":"Sec. 5.1, Eqs. (90)-(93)"}],"minor_comments":[{"comment":"The abstract and conclusion state that the gas Q should be redefined to 'close to the square' of the traditional definition, but the recommended Q_delta=5/3 (Eq. 64) reduces to Q_g^{5/3} for a pure gas disc, not Q_g^2; the Q^2 relation holds only for the intermediate Q_dr. The wording should be adjusted to match the final recommended definition.","section":"Abstract and Sec. 6"},{"comment":"The text 'epsilon between 11 6 and 2' contains a typographical artifact ('11 6' should be '11/6'). Since the recommended epsilon is stated as approximately 2, this is a presentation issue, but the broken fraction should be fixed.","section":"Sec. 4.3, around Eq. (79)"},{"comment":"The choice of kappa=85 km/s/kpc to obtain Q_delta=1.03, after noting that the inferred value 81.6 gives Q_delta=0.99, is presented without uncertainty quantification. Since the galaxy is face-on and the rotation curve is poorly constrained, the statement that the disc is 'marginally unstable' or 'marginally stable' needs an error bar or a more careful caveat.","section":"Sec. 5.2, after Table 3"},{"comment":"The definition of X=k_crit/k_y and the sampling choices (250 values of X, 50 initial phases, 1000 timesteps) are given in the text, but there is no convergence test for these numerical parameters; a brief convergence statement would strengthen confidence in the reported amplification factors.","section":"Sec. 3, Eq. (45)"},{"comment":"The statement that code 'will be shared on reasonable request' is followed by a GitHub link; the link should be described as the canonical source and, ideally, the notebook should be archived with a version identifier to support reproducibility.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its approximations and provides useful analytic insight, but the central recommendation Q_delta=5/3 is calibrated and validated within the same GLB framework that it aims to replace. I would advise the editor that this is a correctness-risk issue rather than a presentation issue: the authors should add an external validation (e.g. comparing predicted m-dependence or Q values against simulations) or substantially soften the claim that Q_delta should become the new standard. The paper's scope is appropriate for the journal, and the material is potentially publishable after such an addition."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper makes one genuinely new claim: the Toomre Q for gas should be redefined close to its square when used as a swing-amplification measure, and a multi-component Q_δ with δ≈5/3 gives a nearly one-to-one map to maximum GLB amplification. The analytic relation Q_dr = Q^2 for gas is clean, and the authors are transparent that the GLB formalism is approximate and that their thickness treatment is a 'curiosity.' The code is a plus.\n\nThe soft spots are real but not fatal. The δ=5/3 exponent is calibrated by minimizing Δ_GLB, and the improvement is then demonstrated using the same Δ_GLB. That is circular in part, and the paper doesn't hide it. The cross-check with the JT formalism covers only single-component stellar discs, so it doesn't validate the gas-rich regime where the redefinition bites. The applications to the Milky Way and M74 are qualitative; the Solar Neighbourhood Q change from 1.45 to 1.58 sits inside the error bars they quote.\n\nI'd still send it to review. The gas Q^2 point is an important correction that will affect how people interpret gas-rich discs, and the paper is honest about what it can't do. A referee should push for one of two things: an independent test (simulation or observation) of the amplified-m parity, or at least an explicit statement that Q_δ is designed to be a good single-number summary of GLB amplification, not a certified measure of real spiral strength. The current text goes some way toward the second, but the abstract and conclusions oversell it as the new standard.\n\nWho gets value: anyone measuring or simulating disc stability in gas-rich galaxies. It's a serious paper, not a quick desk-reject.","headline":"Gas Q should indeed be roughly squared for swing amplification, but the calibration and validation both lean on the same GLB model, so treat δ=5/3 as a useful internal consistency result until independent tests appear.","tokens_in":27586,"tokens_out":1962,"would_cite":true,"duration_ms":18667,"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":"A redefined disc-stability Q maps one-to-one to spiral amplification, and raises the Solar Neighbourhood value to 1.58.","keywords":["Toomre Q","swing amplification","disc stability","multi-component discs","interstellar gas","spiral structure","Milky Way","M74"],"falsifier":"Run a hydro/N-body simulation that holds $Q_{\\delta=5/3}$ fixed while changing the gas fraction from about 0.25 to 0.75; if the measured spiral amplitude differs substantially between the runs, the claimed one-to-one map between $Q$ and swing amplification fails. A second test is observational: resolved cold-gas maps of the Milky Way should show recurrent $m \\approx 10$–$25$ spirals independent of the stellar arms if the predicted gas instability is real.","tokens_in":26485,"feed_emoji":"🌀","tokens_out":12809,"duration_ms":103849,"temperature":0.7,"pith_summary":"Toomre's Q is the standard stability parameter for galactic discs, but it was built to mark the boundary of axisymmetric instability at $Q=1$, not to measure spiral activity above that boundary. The paper argues that using $Q$ at $Q>1$ as a gauge of spiral strength is nevertheless common practice and, for gas-rich discs, numerically wrong. Working from the Goldreich-Lynden-Bell swing-amplification formalism, the authors show that the gas $Q$ should be redefined close to the square of its traditional value, and they construct a multi-component parameter $Q_{\\delta=5/3}$ that gives a close one-to-one map between $Q$ and the maximum swing amplification factor while preserving $Q=1$ as the axisymmetric stability boundary. Applied to the Solar Neighbourhood the new definition raises $Q$ from about 1.44 to 1.58, closer to simulation values, and it predicts that high-$m$ gas-dominated instabilities, not stellar spiral arms, produce the flocculent gas structure seen in the Milky Way and M74.","feed_headline":"Galaxy stability Q redefined to match spiral activity directly","feed_subtitle":"New multi-component Q keeps Q=1 as the stability edge and raises the Solar Neighbourhood estimate to 1.58.","key_machinery":"The load-bearing machinery is the Goldreich-Lynden-Bell (GLB) swing-amplification calculation, extended to discs with multiple stellar and gas components through Rafikov's WKB dispersion relation. The GLB equation gives the reduced spring rate of a shearing density wave as a function of wave inclination; when this rate is negative the wave grows, and integrating the growth over the swing from leading to trailing gives a maximum amplification factor. The paper uses this factor as the target: a useful $Q$ must make maximum amplification a decreasing, one-to-one function of $Q$ (requirement (ii) in Section 2.2). The new parameter $Q_\\delta$ (equation 64) is defined by a fixed-point condition on the dispersion relation—multiply all surface densities by $Q_\\delta$, scale gas sound speeds by $Q_\\delta^{1-1/\\delta}$, and require marginal axisymmetric stability—so that the map between $Q$ and amplification is as composition-independent as possible. The choice $\\delta=5/3$ is the value that minimizes the spread $\\Delta_{\\mathrm{GLB}}$ across gas fractions and sound-speed ratios.","core_discovery":"On the paper's own terms, the discovery is that $Q$ as traditionally defined conflates two different instabilities: the axisymmetric stability boundary ($Q=1$) and the susceptibility to non-axisymmetric spiral growth via swing amplification (the regime $Q>1$). While this conflation is harmless for a single stellar population, it becomes a real numerical error when gas is present. Using the Goldreich-Lynden-Bell shearing-sheet calculation extended to multi-component discs through the Rafikov WKB dispersion relation, the authors find that the standard gas $Q$ must be redefined close to its square to put gas and stellar discs on the same amplification scale. They then define a new multi-component parameter $Q_\\delta$ (equation 64, optimum $\\delta = 5/3$) by demanding that multiplying all surface densities by $Q_\\delta$ while scaling the gas sound speed by $Q_\\delta^{1-1/\\delta}$ brings the disc to marginal axisymmetric stability. This parameter preserves $Q=1$ as the stability boundary, reduces to the standard stellar $Q$ for a pure stellar disc, and in the authors' tests gives a close one-to-one relationship between $Q$ and maximum swing amplification across gas fractions. The paper concludes that $Q_\\delta$ should replace $Q_{\\mathrm{Rk}}$ as the standard stability measure for thin multi-component discs.","pith_inferences":["Beyond the paper: if the gas $Q$ is indeed squared, published comparisons of stability between gas-rich and gas-poor galaxies using the classical gas $Q$ may have systematically underestimated the stability of gas-rich systems; re-analysis of existing $Q$ catalogues would reveal how large the correction is in practice.","Beyond the paper: a direct controlled N-body/hydro experiment that fixes $Q_{\\delta=5/3}$ and varies gas fraction would test the one-to-one claim more cleanly than the paper's existing simulation comparison, which is not designed for this purpose.","Beyond the paper: the same fixed-point construction that defines $Q_\\delta$ could be applied to multi-age stellar populations, providing a composition-independent stability parameter for discs with arbitrary stellar age distributions, an extension the paper only sketches.","Beyond the paper: if the gas-instability link is right, high-resolution cold-gas surveys of the Milky Way should reveal ubiquitous $m \\approx 10$–$25$ spiral structure independent of stellar arms, which can be checked against existing molecular-line observations."],"forward_implications":["For any disc containing gas, the traditional gas $Q$ must be redefined close to its square, so gas-rich discs are more stable at a given classical $Q$ than previously assumed, changing the relative weight of stellar and gas contributions to spiral activity.","The new parameter $Q_{\\delta=5/3}$ (equation 64) is recommended as the standard stability measure for thin multi-component discs, replacing $Q_{\\mathrm{Rk}}$; an iterative formula (equation 65) makes it straightforward to compute.","The Solar Neighbourhood's $Q$ rises from about 1.44 under the Rafikov definition to 1.58 under $Q_{\\delta=5/3}$, bringing the local disc estimate into line with self-regulated N-body simulations that settle near $Q \\approx 1.5$.","Swing amplification peaks broadly in azimuthal wavenumber $m$: stellar-dominated patterns favour low $m$, while gas-dominated instability appears at $m \\gtrsim 10$, which the paper links to the flocculent gas structure in M74 and local spurs in the Milky Way.","For the Milky Way, $m=2$ waves have maximum amplification below 1.25, too weak to sustain a two-armed pattern through feedback, whereas $m=4$ waves are viable; M74's $m=2$ grand design is close to its amplification peak."],"supporting_citations":[{"why":"Supplies the shearing-sheet equation of motion (the GLB equation) that the paper integrates to obtain maximum swing amplification.","marker":"Goldreich & Lynden-Bell (1965)"},{"why":"Adapted the GLB formalism to stellar discs and defined the swing amplification calculation that the paper takes as its target measure of spiral activity.","marker":"Toomre (1981)"},{"why":"Provides the WKB dispersion relation for multi-component discs of gas and stellar components, the basis for the paper's $Q_\\delta$ fixed-point definition.","marker":"Rafikov (2001)"},{"why":"Gives the approximate multi-component Q (Q_RW) that the paper extends into the iterative approximation $Q_{\\delta,\\mathrm{RW}}$.","marker":"Romeo & Wiegert (2011)"},{"why":"Surveys existing multi-component Q definitions and supplies the error scaling used to justify avoiding them for many components.","marker":"Romeo & Falstad (2013)"},{"why":"Provides a Julian-Toomre implementation of swing amplification used to cross-check the GLB computations for stellar discs.","marker":"Binney (2020)"},{"why":"N-body simulations of star-forming discs that settle near Q ≈ 1.5, the simulation comparison point for the Solar Neighbourhood value.","marker":"Aumer et al. (2016)"},{"why":"Supplies the multi-phase ISM surface densities and sound speeds used in the Solar Neighbourhood and M74 applications.","marker":"Holmberg & Flynn (2000)"}],"fun_headline_variants":["Galaxy stability Q redefined to fix gas disc mismatch","New Q for gas discs: stability and swing amplification unified","Q=1 stays stability edge, but gas Q now squares for swing","Redefined Q keeps axisymmetric edge, links to swing amplification","Gas disc Q redefined: now matches swing amplification scale"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the maximum swing-amplification factor from the local Goldreich-Lynden-Bell shearing-sheet model, extended to stars and gas through Rafikov's dispersion relation, is the right measure of spiral activity, so that minimizing the composition-dependent scatter in that factor is the correct design criterion for a new Q; if that model does not track real spiral growth, the proposed Q is merely a re-parametrisation of an approximate calculation.","fun_headline_variants_meta":{"raw":{"variants":["Galaxy stability Q redefined to fix gas disc mismatch","New Q for gas discs: stability and swing amplification unified","Q=1 stays stability edge, but gas Q now squares for swing","Redefined Q keeps axisymmetric edge, links to swing amplification","Gas disc Q redefined: now matches swing amplification scale"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000634,"raw_usage":{"total_tokens":3010,"prompt_tokens":1115,"completion_tokens":1895,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":731,"completion_tokens_details":{"reasoning_tokens":1810}},"tokens_in":731,"tokens_out":1895,"duration_ms":13677,"temperature":1.0,"reasoning_tokens":1810,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T19:51:44.413907+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a hydro/N-body simulation that holds $Q_{\\delta=5/3}$ fixed while changing the gas fraction from about 0.25 to 0.75; if the measured spiral amplitude differs substantially between the runs, the claimed one-to-one map between $Q$ and swing amplification fails. A second test is observational: resolved cold-gas maps of the Milky Way should show recurrent $m \\approx 10$–$25$ spirals independent of the stellar arms if the predicted gas instability is real.","supporting_citations":[{"cited_title":"R., 2001, @doi [ ] 10.1046/j.1365-8711.2001.04201.x , https://ui.adsabs.harvard.edu/abs/2001MNRAS.323..445R 323, 445","cited_arxiv_id":null,"evidence_quote":"Provides the WKB dispersion relation for multi-component discs of gas and stellar components, the basis for the paper's $Q_\\delta$ fixed-point definition."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides a Julian-Toomre implementation of swing amplification used to cross-check the GLB computations for stellar discs."}],"review_version":1}