{"id":"bbcdd89d-0ca2-438f-8d41-70d993c6825a","arxiv_id":"2607.19540","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The first explicit linear theory of swing amplification in coupled star–gas disks shows transient wave growth of 10–100x in formally stable disks, governed by distance to the axisymmetric stability boundary.","lead":"This paper derives and solves the long-missing equations for how density waves ('swing amplification') grow in a galactic disk containing both stars and gas, treating the stars realistically rather than approximating them as a fluid. It finds waves can be briefly boosted 10–100x in disks that are formally stable, with the boost controlled by how close the disk sits to the instability threshold.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified","rationale":"The reader's verdict ACCEPT is well supported. The derivations in Appendix A are consistent and the code is public. The gas-closure assumption is a genuine limitation for astrophysical application, but it is explicitly stated by the authors and does not undermine the formal linear-theory claims, which are the core contribution. I considered whether the separate time maximization in Eq. (25) constitutes a more serious internal problem: since the numerator and denominator can peak at different initial times t_i, the quoted amplification factors are upper bounds. However, the central '10-100x' result is already demonstrated in Fig. 3 with fixed t_i and raw amplitudes, and the correlation with stability boundary in Fig. 5 is a qualitative pattern likely robust to this definitional choice. Thus I recommend no change to the ACCEPT verdict, with the caveat that the public code be used to verify the quantitative contours under a matched-t_i definition.","tokens_in":16284,"tokens_out":28454,"duration_ms":285006,"concrete_test":"Re-run the public code to produce Fig. 5 using the matched-t_i ratio R(t_i)=max_{t>t_i}|δΣ_s(t;t_i)|/max_{t>t_i}|δΣ_s^0(t;t_i)| and then maximize over t_i. If the 10-100+ contours survive this redefinition, the quantitative claim is robust; if they drop substantially, the formal results still stand but the quoted amplification factors should be re-expressed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"After checking the derivation, I find no load-bearing flaw in the central argument. The master equations (14)-(15) are derived in detail; re-deriving the pure-gas limit from (A20)-(A21) confirms the S^2(t) in (A27) and the absence of an artificial first-derivative term after transforming from surface density to potential. The Volterra kernel (A36-A40) has K(t,t)=0, so the discrete scheme (B53) is well-posed. The main physical claims are supported by the numerical solutions shown, and the paper explicitly acknowledges the idealizations (gas closure, no turbulence/fields, local/linear). The most unconventional step is the separate maximization over t_i in Eq. (25), which technically upper-bounds a matched-t_i amplification factor; this could inflate some quoted numbers but does not threaten the qualitative results or the 10-100x factors already visible in Fig. 3.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops the linear theory of swing amplification (SA) in a two-component shearing sheet consisting of a collisionless stellar component and a gaseous fluid. The authors derive coupled master equations for the gas potential (Eq. 14) and the stellar potential through a Volterra kernel (Eq. 15), provide a finite-thickness extension in §2.2 and Appendix A.3, and solve these equations numerically for impulsively excited swinging waves. Their main claims are: (i) transient non-axisymmetric amplification factors of 10–100+ occur in disks that are stable to axisymmetric perturbations, with the maximum amplification set primarily by the distance to the axisymmetric stability boundary in the (1/Q_s, 1/Q_g) plane (Fig. 5); and (ii) the commonly used 'two-fluid' treatment of stars (Jog 1992) is qualitatively incorrect because stellar phase mixing acts as a sink of perturbation energy, so no neutral non-axisymmetric modes exist in the true collisionless star-gas system. The paper also provides a public Python code for computing SA factors and maximally amplified wavelengths.","tokens_in":16341,"tokens_out":15576,"duration_ms":156748,"significance":"If correct, this paper fills a long-standing gap in the local linear theory of disk instabilities. Previous analytic treatments either considered a single component (Goldreich & Lynden-Bell 1965; Julian & Toomre 1966) or approximated the stars as a fluid (Jog 1992), and the coupled collisionless star-gas equations have not been solved explicitly before. The appendices present a careful derivation of the gas oscillator and the Julian-Toomre-style kernel, and the authors demonstrate recovery of the known single-component limits. The public code and the explicit discussion of caveats (gas closure, no turbulence, no magnetic fields, local approximation) are strengths. The paper also provides a practical result: the axisymmetric criterion (22) is a useful guide for when strong non-axisymmetric structure should appear in high-redshift, gas-rich disks. The main concern is the definition of the 'maximum amplification factor' in Eq. (25), which is used to produce the quantitative headline numbers in the abstract and Fig. 5.","major_comments":[{"comment":"The maximum amplification factor is defined as the ratio of two maxima taken separately over the initial time t_i: [max_{t_i} max_{t>t_i} |δΣ_s(t;t_i)|] / [max_{t_i} max_{t>t_i} |δΣ_s^0(t;t_i)|]. This does not correspond to the amplification factor of any single impulsively excited wave, because the two maxima can occur at different values of t_i. The resulting number is not necessarily an upper or lower bound on the natural matched definition max_{t_i} [max_{t>t_i} |δΣ_s| / max_{t>t_i} |δΣ_s^0|]. Since Fig. 5 and the abstract's 'factors of 10–100+' are based on this metric, the authors should either adopt a matched-t_i definition or demonstrate numerically (e.g., by reporting the argmax t_i for both numerator and denominator) that the separate maximization does not bias the contours or the quoted factors. This is a load-bearing point for the paper's central quantitative claim.","section":"§3.2, Eq. (25)"},{"comment":"The comparison between the true star-gas calculation and the 'two-fluid' calculation (dashed lines) is used to support the claim that two-fluid prescriptions are qualitatively wrong. However, the two-fluid run treats the stellar component as a fluid with sound speed σ while assigning it the same Q_s defined in Eq. (18). A fluid with sound speed σ has Toomre parameter Q_fluid = κσ/(πGΣ_s0), which is (3.36/π)Q_s ≈ 1.07Q_s. Thus the two fluid disks are not matched to the collisionless disks in their distance from the relevant stability boundary. The authors should either run a controlled two-fluid comparison (e.g., adjusting the stellar sound speed or surface density to match Q_fluid to the collisionless Q_s) or explicitly state that the figure is illustrative rather than a quantitative controlled experiment. Without this, the strength of the 'two-fluid is qualitatively wrong' conclusion is","section":"§3.1, Fig. 3"}],"minor_comments":[{"comment":"The text states that in the first row 'the maximum amplification factor in the stars is ~50 in all three panels' and in the second row the gas is amplified 'by up to a factor ~100'. However, panels (c) and (f) have much smaller y-axis ranges (±25 and ±10, respectively). Please qualify these statements or correct them if the quoted values apply only to the leftmost columns.","section":"§3.1, Fig. 3"},{"comment":"The maximum amplification metric (25) uses the stellar surface density response only, while the abstract summarizes the result as 'waves are amplified'. Since gas and star amplification factors can differ by factors of a few (Fig. 3), please specify in the text and abstract that the quoted maximum factors refer to the stellar component, and clarify how the gas response is quantified.","section":"Eq. (25)"},{"comment":"The effective thick-disk master equations (A42)–(A46) are stated without derivation. The presence of the first-derivative term S1(t) in Eq. (A42) is nontrivial. Since finite-thickness behavior is advertised in the abstract, a fuller derivation or at least a more explicit explanation of the approximations entering Eqs. (A42)–(A46) would improve the paper.","section":"§2.2 / Appendix A.3"}],"recommendation":"major_revision","confidential_remarks":"The derivation appears sound and the paper fills a genuine gap. My main concern is the definition of the maximum amplification factor in Eq. (25), which is load-bearing for the headline quantitative claims; the two-fluid comparison also needs to be placed on a more controlled footing. These are fixable with additional analysis, and I would be satisfied with a revised version addressing them."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — this paper does what it says: it writes down the coupled linear swing-amplification equations for a collisionless stellar sheet and a gas sheet, solves them, and ships the code. That pairing has been missing since Julian & Toomre. The derivation in Appendix A is careful; the equations reduce to Goldreich–Lynden-Bell and Julian–Toomre in single-component limits, and the axisymmetric limit reproduces Rafikov and Kim–Ostriker. The validation against Binney’s cloud wake is the kind of check that earns trust.\n\nThe main physical claims are supported: amplification factors of 10–100 occur near the axisymmetric stability boundary, and the contours of maximum amplification track that boundary. The contrast with two-fluid models is the point worth thinking about: the stars phase-mix, so the coupled system does not oscillate forever. That is a real qualitative difference, and a fair warning against treating stars as a fluid.\n\nSoft spots, in proportion. The gas closure is the standard one — constant sound speed, no turbulence, no magnetic fields — and the authors are explicit that real z~4 disks are messier. That is an idealization, not a flaw. The separate maximization over t_i in Eq. (25) technically upper-bounds the matched-t_i amplification factor, so the absolute numbers in Figs. 4–5 should be read as optimistic, but the qualitative pattern and the 10–100 range are already visible in Fig. 3. The footnote on the Goldreich–Tremaine erratum is stated without derivation; I would ask the authors to show the correction. And there is no convergence documentation for the numerical scheme — minor, but easy to add.\n\nOverall, this is a solid, useful paper. It deserves a serious referee. It will be cited as the reference for star-gas swing amplification, and the code will get used. Anyone working on disk instabilities or interpreting high-redshift bars and spirals will want it. I would take it with the gas-closure caveat in mind, but not with suspicion.","headline":"A careful, complete derivation that fills a real gap in swing-amplification theory; the main claims hold up, with the standard gas-closure caveats clearly flagged.","tokens_in":17027,"tokens_out":2373,"would_cite":true,"duration_ms":25765,"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":"Swing amplification in star-gas disks is solved for the first time, showing transient 10–100× amplification in stable disks.","keywords":["swing amplification","shearing sheet","star-gas disks","collisionless stellar dynamics","phase mixing","Toomre Q","axisymmetric stability","high-redshift disk galaxies"],"falsifier":"Run a local shearing-box simulation with collisionless stars and realistic (turbulent, magnetized, multiphase) gas in a disk that is axisymmetrically stable, with parameters near the stability boundary; if maximum non-axisymmetric amplification does not track the distance to that boundary, or if gas perturbations persist as neutral oscillations instead of decaying after the first amplification peak, the paper's central mechanism is refuted.","tokens_in":15983,"feed_emoji":"🌀","tokens_out":4384,"duration_ms":43865,"temperature":0.7,"pith_summary":"The paper writes down and solves the linear equations for non-axisymmetric 'swinging' waves in a two-component galactic disk, with stars treated as a collisionless population and gas as a fluid. Its central finding is that such waves are transiently amplified by factors of ten to a hundred even in disks that are stable to axisymmetric collapse, and that the maximum amplification is set almost entirely by how close the disk sits to the axisymmetric stability boundary in the (1/Qs, 1/Qg) plane. It also finds that the standard two-fluid treatment, which models stars as a gas with an effective sound speed, is qualitatively wrong: collisionless phase mixing makes stars a permanent sink of wave energy, so there are no neutral non-axisymmetric modes. The result matters because high-redshift disks observed by JWST and ALMA are gas rich, and this is the first quantitative baseline for interpreting their spiral and bar features.","feed_headline":"Star-gas swing amplification solved: 10–100× transient boosts","feed_subtitle":"Collisionless stars phase-mix and drain wave energy, so spiral/bar features are transient, not neutral modes.","key_machinery":"The argument is carried by two coupled linear equations. Eq. (14) is a forced harmonic-oscillator equation for the gas potential, with a time-dependent spring constant S^2(t) that encodes rotation, epicyclic motion, gas pressure, and self-gravity. Eq. (15) is a Volterra integral equation for the stellar potential whose kernel is the collisionless stellar response ('JT kernel'): it sums the delayed gravitational influence of all past potential fluctuations on the stellar density. The pair is solved forward in time for impulsive sinusoidal perturbations, and the amplification factor is measured relative to a passive non-self-gravitating response. The same structure carries over to finite-thick","core_discovery":"On the paper's own terms, the discovery is that the coupled star-gas swing-amplification problem admits an explicit linear solution — a gas harmonic-oscillator equation with a time-dependent spring constant coupled through gravity to a stellar response described by the collisionless kernel — and that this solution changes the expected phenomenology. Waves amplify strongly only near the axisymmetric stability boundary; amplification factors can reach 10–100 or more. Finite-thickness disks behave qualitatively the same, with a shifted stability boundary and suppressed small-scale response. The qualitatively new feature with respect to earlier work is that the stellar phase-mixing sink damps th","pith_inferences":["If the phase-mixing sink is the correct picture, swing amplification is a one-shot energy transfer: each episode heats the stellar distribution and therefore reduces the disk's own capacity for the next episode, which may help explain why isolated disk simulations produce repeated but episodic spiral activity.","Because the amplification contours are nearly parallel to the stability boundary, an effective single-Q parameter for two-component disks could be built by inverting the maximum-amplification map; this would extend the paper's two-component result to multi-component modeling without solving the full equations.","A testable extension is to run local shearing-box simulations with collisionless stars plus a multiphase/turbulent gas: if gas perturbations damp and amplification tracks distance to the axisymmetric boundary, the paper's baseline survives; if turbulence re-excites neutral modes, the sink picture needs modification."],"forward_implications":["The axisymmetric stability criterion doubles as a map of expected non-axisymmetric swing amplification: proximity to the boundary in the (1/Qs, 1/Qg) plane predicts where 10–100× amplification occurs.","Spiral and bar features generated by swing amplification in stable star-gas disks are inherently transient; after the first peak, stellar phase mixing damps both components, so persistent features require continuous driving.","Two-fluid idealizations of star-gas dynamics are energetically misleading: they permit neutral oscillations that the true collisionless-star system cannot support, so they should not be used for quantitative predictions.","Finite-thickness disks behave the same qualitatively, so the baseline extends to realistic disks provided the effective thickness parameters are known.","A public code now gives the maximum amplification factor and most-amplified wavelength for given disk parameters, giving observers a direct way to compare high-redshift spiral/bar detections with theory."],"fun_headline_variants":["Swing amplification in star-gas disks now solved explicitly","Star-gas spiral waves amplify 10-100× transiently","Stellar phase mixing quenches swing amplification","New linear solution predicts spiral boost near stability edge","Coupled star-gas swing amplification: transient boosts up to 100×"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The gas is modeled as a single inviscid, isentropic fluid with a constant sound speed in a smooth, homogeneous shearing sheet, so if real turbulent magnetic multiphase gas does not obey that closure, the predicted 10–100× amplification and the permanent damping of gas waves need not transfer to real disks.","fun_headline_variants_meta":{"raw":{"variants":["Swing amplification in star-gas disks now solved explicitly","Star-gas spiral waves amplify 10-100× transiently","Stellar phase mixing quenches swing amplification","New linear solution predicts spiral boost near stability edge","Coupled star-gas swing amplification: transient boosts up to 100×"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000198,"raw_usage":{"total_tokens":1215,"prompt_tokens":763,"completion_tokens":452,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":507,"completion_tokens_details":{"reasoning_tokens":372}},"tokens_in":507,"tokens_out":452,"duration_ms":4936,"temperature":1.0,"reasoning_tokens":372,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T12:29:12.630087+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a local shearing-box simulation with collisionless stars and realistic (turbulent, magnetized, multiphase) gas in a disk that is axisymmetrically stable, with parameters near the stability boundary; if maximum non-axisymmetric amplification does not track the distance to that boundary, or if gas perturbations persist as neutral oscillations instead of decaying after the first amplification peak, the paper's central mechanism is refuted.","supporting_citations":[],"review_version":1}