arxiv: 2604.06862 · v1 · submitted 2026-04-08 · 🌌 astro-ph.GA

Recognition: no theorem link

Phase spirals induced by the gas warp

Shuyu Wang , Anthony G. A. Brown , Victor P. Debattista , Tigran Khachaturyants

Authors on Pith no claims yet

Pith reviewed 2026-05-10 18:26 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords phase spiralsgas warpMilky Way discbending wavesN-body SPH simulationsGaia observationsstellar kinematicsgalactic dynamics

0 comments

The pith

Irregular gas accretion along the galactic warp induces global, long-lived phase spirals in the stellar disc.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tests whether irregular gas inflow along a galactic warp can generate phase spirals in the stellar disc, like those detected by Gaia in the Milky Way's solar neighbourhood. In N-body plus SPH simulations of an isolated galaxy, one run includes an imposed warp with gas accretion while a control run does not; the warped case produces strong one-armed spirals that appear across the disc and persist for roughly 10 billion years. Their emergence rate tracks the gas inflow and bending-wave strength, implying the spirals act as a record of warp-driven perturbations. This mechanism offers an internal source for the observed spirals without needing external satellite encounters. Readers would care because it alters how we reconstruct the Milky Way's recent dynamical history from stellar kinematics.

Core claim

Using two N-body + SPH simulations of an isolated galaxy, one with a warp and gas inflow and one without, we detect strong one-armed phase spirals in the warped simulation that are prevalent across the disc and persist over ~10 Gyr. The morphology of these spirals varies with location and evolves with time, with their emergence rate tied to the gas inflow at the outer disc and the bending-wave amplitude. These spirals reach amplitudes comparable to those in Gaia DR3 data. Only weak and stochastically distributed spirals appear in the unwarped control simulation. We conclude that phase spirals can be induced by the irregular gas accretion along the warp and that these structures occur both on

What carries the argument

Warp-induced bending waves from irregular gas accretion, tracked via Fourier decomposition of stellar phase-space structure in the simulations.

If this is right

Phase spirals can form globally throughout the disc rather than remaining localized features.
Spiral properties evolve directly with bending-wave amplitude, allowing them to serve as a historical record of warp activity.
Spiral amplitudes can reach observed Gaia levels solely from warp-driven gas inflow.
Unwarped discs produce only weak, randomly distributed spirals, highlighting the warp's role.
The spirals remain detectable for ~10 Gyr, implying they preserve information about ancient gas accretion events.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

This internal mechanism could account for phase spirals in galaxies lacking recent mergers or satellite encounters.
Gas-dynamical processes may shape large-scale stellar phase-space features more than previously emphasized in merger-dominated models.
Observers could test the scenario by searching for correlations between gas warp strength and phase-spiral visibility in external galaxies.
The long lifetime raises the possibility of using phase spirals to date past episodes of irregular gas accretion in disc galaxies.

Load-bearing premise

The isolated-galaxy N-body+SPH model with an imposed warp and gas inflow sufficiently captures the relevant physics and timescales of the real Milky Way disc so that the detected spirals are not simulation artifacts.

What would settle it

High-resolution mapping of bending-wave amplitudes across the Milky Way disc that shows no spatial or temporal correlation with the strength or presence of phase spirals.

Figures

Figures reproduced from arXiv: 2604.06862 by Anthony G. A. Brown, Shuyu Wang, Tigran Khachaturyants, Victor P. Debattista.

**Figure 3.** Figure 3: The phase spiral in region 4 of the warped galaxy at 8 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

**Figure 5.** Figure 5: The phase spirals in the Rg-θϕ plane for star particles with JR < median(JR) in the warped galaxy at 8.7 Gyr, color-coded by density contrast. Each normalized phase space coordinate covers the range [−3, 3]. expect from an idealized model with a single impact time and no self-gravity (Widmark et al. 2025). Indeed, phase spirals appear to be most tightly wound at Rg ∈ [6.9, 9.1] kpc. The detected phase spir… view at source ↗

**Figure 6.** Figure 6: Incidence map of the m = 1 phase spirals weighted by residual Vϕ at the 12 regions at R ∈ [6.277, 8.277] kpc ( [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗

**Figure 7.** Figure 7: The time evolution of the phase spiral incidence rate at R ∈ [6.277, 8.277] kpc from 1.91 to 11.61 Gyr, in the warped galaxy (dashed black line) versus the gas inflow through a spherical shell at 15 kpc (dashed blue line) and the amplitude of the bending waves in the same region (dashed red line). The Butterworth filter implemented in scipy with a period of 0.175 Gyr is applied to smooth the signals, where… view at source ↗

**Figure 9.** Figure 9: The time evolution of the phase spiral incidence rate at [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗

**Figure 8.** Figure 8: Incidence map of the m = 1 phase spirals weighted by ∆Vϕ at 12 regions in the unwarped galaxy, from 10.75 to 11.71 Gyr. no strong phase spirals have been detected. The phase spirals in the unwarped simulation may be induced by stochastic heating (Tremaine et al. 2023; Tepper-García et al. 2025), [PITH_FULL_IMAGE:figures/full_fig_p007_8.png] view at source ↗

**Figure 10.** Figure 10: The phase spiral in the Gaia DR3 data compared to that in an equivalent region in the warped galaxy at 8.7 Gyr, color-coded by density contrast. The observational sample contains 55 880 stars. The simulation sample contains 35 989 star particles. Each axis shows the normalized phase space coordinate, covering the range [−3, 3]. The vertical position and velocity dispersions of each sample are indicated… view at source ↗

read the original abstract

The discovery of the phase space spirals in the Solar neighborhood in Gaia Data Release 2 has prompted various attempts to understand their origin. A source of bending waves, which has been neglected as a cause of the phase spiral, is irregular gas inflow along the warp. We aim to study whether perturbations by the gas warp could induce phase spirals. Accounting for this additional formation scenario for phase spirals could improve our current understanding of the perturbation history of the Milky Way disc. We use two N-body + SPH (Smooth Particle Hydrodynamics) simulations of an isolated galaxy to search for, and study, warp-induced phase spirals. We study the emergence and propagation of the detected phase spirals using Fourier decomposition. We detect strong one-armed phase spirals in the warped simulation. These phase spirals are prevalent and persist over ~10 Gyr. The morphology of these phase spirals varies with location and evolves with time. In particular, the emergence rate of the phase spiral evolves with the gas inflow at the outer disc and the bending wave amplitude, indicating that these phase spirals are a record of warp-induced bending waves. We find that these phase spirals can reach amplitudes comparable to those in the Gaia DR3. We only detect weak and stochastically distributed phase spirals in an unwarped control simulation. We conclude that phase spirals can be induced by the irregular gas accretion along the warp. These phase spirals occur globally and are long-lived.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The warped simulation produces strong global phase spirals tied to gas inflow over 10 Gyr while the control does not, but the imposed warp leaves open whether this is physical or setup-driven.

read the letter

The main thing to know is that these two N-body+SPH runs show irregular gas accretion along a warp can generate strong, one-armed phase spirals that last ~10 Gyr, reach Gaia-like amplitudes, and appear globally, while the unwarped control only has weak stochastic ones. The authors link the spiral emergence rate to the inflow and bending-wave amplitude via Fourier analysis, which gives a direct handle on the mechanism. This is new because most prior work on phase spirals has focused on satellite passages or bars, and the abstract treats gas-warp inflow as a neglected channel. The controlled contrast is the clearest part of the evidence and avoids obvious circularity. The isolated-galaxy setup with SPH gas is a reasonable first step for testing the idea. That said, the warp is imposed and the abstract gives no resolution tests, convergence checks, or quantitative error bars on the amplitudes or inflow correlation. Without those, it is possible the spirals partly reflect how the warp is maintained rather than pure accretion physics. The model also omits external torques or satellite effects that operate in the real Milky Way, so the longevity and strength may not translate directly. This paper is for people who model Milky Way disc dynamics or interpret Gaia phase-space data. Anyone already running warp or bending-wave simulations would get value from the comparison. It is worth sending to peer review because the simulation contrast supplies a concrete, falsifiable claim even if the methods section needs tightening and more diagnostics.

Referee Report

3 major / 2 minor

Summary. The paper uses two controlled N-body+SPH simulations of an isolated galaxy (one with an imposed warp and irregular gas inflow, one unwarped control) to show that strong, global, one-armed phase spirals can be induced by warp-driven gas accretion. Fourier decomposition reveals these spirals emerge globally, persist for ~10 Gyr, evolve with inflow rate and bending-wave amplitude, reach Gaia-comparable strengths, and are absent or weak in the control run. The authors conclude that such spirals record the perturbation history from gas warps.

Significance. If validated, the result supplies a new, physically motivated channel for phase-spiral formation that links directly to observed Milky Way warps and accretion. The controlled simulation design (warped vs. unwarped) and use of Fourier decomposition to track emergence, propagation, and correlation with inflow provide a clean baseline and falsifiable signature. This strengthens the case that phase spirals can be long-lived global features rather than short-lived local responses.

major comments (3)

[Simulation methods] Simulation methods section: no resolution or convergence tests are reported for the N-body+SPH runs. Because the central claim rests on the detection of coherent, long-lived phase spirals whose amplitudes and morphology could be sensitive to particle number, softening, or SPH kernel, the absence of such tests leaves open the possibility that the reported spirals are numerically seeded rather than physically driven by the warp.
[Results (Fourier decomposition)] Results on Fourier analysis: the persistence over ~10 Gyr and the correlation with gas inflow rate are stated qualitatively, but no quantitative error analysis, bootstrap uncertainties, or statistical significance thresholds for the detected m=1 Fourier modes are provided. This is load-bearing for the claim that the spirals are a robust record of warp-induced bending waves rather than stochastic fluctuations.
[Results (amplitude comparison)] Comparison with Gaia: the assertion that spiral amplitudes are comparable to Gaia DR3 is made without explicit metrics, error bars, or a table/figure showing measured amplitudes, wavelengths, or phase offsets in both the simulation and the data. Without these, the quantitative match cannot be assessed and the claim that the mechanism operates at observed strengths remains unverified.

minor comments (2)

[Methods] The exact variables and radial bins used in the Fourier decomposition should be stated explicitly in the methods so that the one-armed spiral detection can be reproduced.
[Figures] Figure captions for the time-evolution panels would benefit from indicating the precise simulation times shown and the color scale normalization.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments on our manuscript. We address each major comment point by point below, indicating the revisions we will undertake.

read point-by-point responses

Referee: Simulation methods section: no resolution or convergence tests are reported for the N-body+SPH runs. Because the central claim rests on the detection of coherent, long-lived phase spirals whose amplitudes and morphology could be sensitive to particle number, softening, or SPH kernel, the absence of such tests leaves open the possibility that the reported spirals are numerically seeded rather than physically driven by the warp.

Authors: We agree that explicit resolution and convergence tests are not presented and that this represents a limitation. The primary evidence that the phase spirals are physically driven rather than numerically seeded is the direct comparison with the unwarped control simulation, which uses identical numerical parameters, particle numbers, softening lengths, and SPH settings but produces only weak, stochastic spirals. In the revised manuscript we will expand the methods section to report the numerical parameters employed and to discuss the robustness of the results based on this controlled differential experiment. We will also note the absence of full multi-resolution convergence tests as a limitation arising from computational cost. revision: partial
Referee: Results on Fourier analysis: the persistence over ~10 Gyr and the correlation with gas inflow rate are stated qualitatively, but no quantitative error analysis, bootstrap uncertainties, or statistical significance thresholds for the detected m=1 Fourier modes are provided. This is load-bearing for the claim that the spirals are a robust record of warp-induced bending waves rather than stochastic fluctuations.

Authors: The m=1 Fourier amplitudes are extracted from azimuthal decompositions of the stellar density and vertical velocity fields and exhibit clear, sustained evolution that tracks the gas inflow rate and bending-wave amplitude, remaining well above the levels seen in the control run. We acknowledge that formal uncertainties were not included. In the revision we will add bootstrap-derived error estimates on the m=1 amplitudes (obtained by resampling particles within each radial bin) and will define a significance threshold relative to the control-run fluctuations. These quantitative measures will be incorporated into the relevant time-evolution figures and text. revision: yes
Referee: Comparison with Gaia: the assertion that spiral amplitudes are comparable to Gaia DR3 is made without explicit metrics, error bars, or a table/figure showing measured amplitudes, wavelengths, or phase offsets in both the simulation and the data. Without these, the quantitative match cannot be assessed and the claim that the mechanism operates at observed strengths remains unverified.

Authors: We accept that the Gaia comparison was stated qualitatively without supporting metrics. In the revised manuscript we will add a new figure (or table) that reports the measured m=1 phase-spiral amplitudes, characteristic wavelengths, and phase offsets extracted from the simulation at solar-neighbourhood radii and late times, placed alongside the corresponding quantities measured from Gaia DR3 using an identical Fourier procedure. Bootstrap uncertainties will be shown on the simulation values to enable a direct, quantitative assessment of comparability. revision: yes

Circularity Check

0 steps flagged

No circularity: results follow from controlled N-body+SPH simulation outputs

full rationale

The paper's central claim rests on direct comparison of two isolated-galaxy N-body+SPH runs (warped with irregular gas inflow versus unwarped control). Phase spirals are identified via Fourier decomposition of the simulation particle data, with their amplitude, prevalence, and time evolution reported as emergent outputs. No equations, fitted parameters, or self-citations are invoked to derive the spirals; the unwarped run supplies an independent baseline that is not constructed from the same inputs. The derivation chain is therefore self-contained against external benchmarks and does not reduce any reported result to its own setup by definition.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Information is limited to the abstract; the simulations rely on standard assumptions of N-body+SPH galactic modeling whose detailed parameter choices are not stated here.

axioms (1)

domain assumption The N-body + SPH simulation of an isolated galaxy with imposed warp accurately represents the dynamical response of a real disc galaxy to irregular gas accretion.
Invoked when using the simulation results to infer real Milky Way behavior.

pith-pipeline@v0.9.0 · 5561 in / 1290 out tokens · 51129 ms · 2026-05-10T18:26:18.111385+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

49 extracted references · 49 canonical work pages

[1]

J., & Bensby, T

Alinder, S., McMillan, P. J., & Bensby, T. 2023, A&A, 678, A46

work page 2023
[2]

J., & Bensby, T

Alinder, S., McMillan, P. J., & Bensby, T. 2024, A&A, 690, A15

work page 2024
[3]

2018, Nature, 561, 360

Antoja, T., Helmi, A., Romero-Gómez, M., et al. 2018, Nature, 561, 360

work page 2018
[4]

2023, A&A, 673, A115

Antoja, T., Ramos, P., García-Conde, B., et al. 2023, A&A, 673, A115

work page 2023
[5]

S., Baba, J., Portegies Zwart, S., & Bédorf, J

Asano, T., Fujii, M. S., Baba, J., Portegies Zwart, S., & Bédorf, J. 2025, A&A, 700, A109 Astropy Collaboration, Price-Whelan, A. M., Sip˝ocz, B. M., et al. 2018, AJ, 156, 123 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33

work page 2025
[6]

C., & Weinberg, M

Banik, U., van den Bosch, F. C., & Weinberg, M. D. 2023, ApJ, 952, 65

work page 2023
[7]

2012, MNRAS, 426, 1324

Binney, J. 2012, MNRAS, 426, 1324

work page 2012
[8]

& Schönrich, R

Binney, J. & Schönrich, R. 2018, MNRAS, 481, 1501

work page 2018
[9]

2019, MNRAS, 486, 1167

Bland-Hawthorn, J., Sharma, S., Tepper-Garcia, T., et al. 2019, MNRAS, 486, 1167

work page 2019
[10]

& Tepper-García, T

Bland-Hawthorn, J. & Tepper-García, T. 2021, MNRAS, 504, 3168

work page 2021
[11]

S., Dekel, A., Kolatt, T

Bullock, J. S., Dekel, A., Kolatt, T. S., et al. 2001, ApJ, 555, 240

work page 2001
[12]

Chequers, M. H. & Widrow, L. M. 2017, MNRAS, 472, 2751

work page 2017
[13]

W., & Tremaine, S

Chiba, R., Ding, J., Hamilton, C., Kunz, M. W., & Tremaine, S. 2025, MNRAS, 543, 190

work page 2025
[14]

2025, MNRAS, 543, 2159

Chiba, R., Frankel, N., & Hamilton, C. 2025, MNRAS, 543, 2159

work page 2025
[15]

& Widrow, L

Darling, K. & Widrow, L. M. 2018, MNRAS, 484, 1050

work page 2018
[16]

Darragh-Ford, E., Hunt, J. A. S., Price-Whelan, A. M., & Johnston, K. V . 2023, ApJ, 955, 74

work page 2023
[17]

P., Khachaturyants, T., Amarante, J

Debattista, V . P., Khachaturyants, T., Amarante, J. A. S., et al. 2025, MNRAS, 537, 1620

work page 2025
[18]

P., van den Bosch, F

Debattista, V . P., van den Bosch, F. C., Roškar, R., et al. 2015, MNRAS, 452, 4094

work page 2015
[19]

Fiteni, K., Caruana, J., Amarante, J. A. S., Debattista, V . P., & Beraldo e Silva, L. 2021, MNRAS, 503, 1418

work page 2021
[20]

Frankel, N., Bovy, J., Tremaine, S., & Hogg, D. W. 2023, MNRAS, 521, 5917 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1 Gaia Collaboration, Recio-Blanco, A., Kordopatis, G., et al. 2023a, A&A, 674, A38 Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023b, A&A, 674, A1 García-Conde, B., Antoja, T., Roca-Fàbrega, ...

work page 2023
[21]

P., & Khachaturyants, T

Ghosh, S., Debattista, V . P., & Khachaturyants, T. 2022, MNRAS, 511, 784

work page 2022
[22]

2025, ApJ, 980, 24

Gilman, D., Bovy, J., Frankel, N., & Benson, A. 2025, ApJ, 980, 24

work page 2025
[23]

Grand, R. J. J., Pakmor, R., Fragkoudi, F., et al. 2023, MNRAS, 524, 801 GRA VITY Collaboration, Abuter, R., Aimar, N., et al. 2022, A&A, 657, L12

work page 2023
[24]

2024, ApJ, 960, 133

Guo, R., Li, Z.-Y ., Shen, J., Mao, S., & Liu, C. 2024, ApJ, 960, 133

work page 2024
[25]

R., Millman, K

Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357

work page 2020
[26]

Hunt, J. A. S., Price-Whelan, A. M., Johnston, K. V ., & Darragh-Ford, E. 2022, MNRAS, 516, L7

work page 2022
[27]

Hunt, J. A. S., Price-Whelan, A. M., Johnston, K. V ., et al. 2024, MNRAS, 527, 11393

work page 2024
[28]

Hunt, J. A. S., Stelea, I. A., Johnston, K. V ., et al. 2021, MNRAS, 508, 1459

work page 2021
[29]

Hunt, J. A. S. & Vasiliev, E. 2025, New A Rev., 100, 101721

work page 2025
[30]

Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90

work page 2007
[31]

P., & Daniel, K

Khachaturyants, T., Beraldo e Silva, L., Debattista, V . P., & Daniel, K. J. 2022, MNRAS, 512, 3500

work page 2022
[32]

P., Ghosh, S., Beraldo e Silva, L., & Daniel, K

Khachaturyants, T., Debattista, V . P., Ghosh, S., Beraldo e Silva, L., & Daniel, K. J. 2022, MNRAS, 517, L55

work page 2022
[33]

2019, A&A, 622, L6

Khoperskov, S., Di Matteo, P., Gerhard, O., et al. 2019, A&A, 622, L6

work page 2019
[34]

Laporte, C. F. P., Minchev, I., Johnston, K. V ., & Gómez, F. A. 2019, MNRAS, 485, 3134

work page 2019
[35]

2023, MNRAS, 524, 6331

Li, C., Siebert, A., Monari, G., Famaey, B., & Rozier, S. 2023, MNRAS, 524, 6331

work page 2023
[36]

2025, ApJ, 988, 254

Lin, J., Li, Z.-Y ., Guo, R., et al. 2025, ApJ, 988, 254

work page 2025
[37]

2016, A&A, 595, A4

Lindegren, L., Lammers, U., Bastian, U., et al. 2016, A&A, 595, A4

work page 2016
[38]

Miller, G. E. & Scalo, J. M. 1979, ApJS, 41, 513

work page 1979
[39]

F., Frenk, C

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563

work page 1996
[40]

S., et al

Pontzen, A., Roškar, R., Stinson, G. S., et al. 2013, pynbody: Astrophysics Sim- ulation Analysis for Python, astrophysics Source Code Library, ascl:1305.002

work page 2013
[41]

2006, MNRAS, 373, 1074 Tepper-García, T., Bland-Hawthorn, J., Bedding, T

Stinson, G., Seth, A., Katz, N., et al. 2006, MNRAS, 373, 1074 Tepper-García, T., Bland-Hawthorn, J., Bedding, T. R., Federrath, C., & Agertz, O. 2025, MNRAS, 542, 1987

work page 2006
[42]

2023, MNRAS, 521, 114

Tremaine, S., Frankel, N., & Bovy, J. 2023, MNRAS, 521, 114

work page 2023
[43]

2019, MNRAS, 482, 1525

Vasiliev, E. 2019, MNRAS, 482, 1525

work page 2019
[44]

& Belokurov, V

Vasiliev, E. & Belokurov, V . 2020, MNRAS, 497, 4162

work page 2020
[45]

E., et al

Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261

work page 2020
[46]

W., Stadel, J., & Quinn, T

Wadsley, J. W., Stadel, J., & Quinn, T. 2004, New A, 9, 137

work page 2004
[47]

Widmark, A., Laporte, C., & de Salas, P. F. 2021, A&A, 650, A124

work page 2021
[48]

V ., & Hunt, J

Widmark, A., Tavangar, K., Kalish, J., Johnston, K. V ., & Hunt, J. A. S. 2025, arXiv e-prints, arXiv:2507.19579

work page arXiv 2025
[49]

Widrow, L. M. 2023, MNRAS, 522, 477 Article number, page 9 A&A proofs:manuscript no. aa58812-25 Appendix A: Supplementary figures A.1. Star formation history Figure A.1 shows the global star formation history (SFH) of the warped simulation. The global star formation rate (SFR) displays an overall trend of decline with time. The SFR in Fig. A.1 ap- pears s...

work page 2023