pith. machine review for the scientific record. sign in

arxiv: 2604.17955 · v2 · submitted 2026-04-20 · 🌌 astro-ph.GA

Recognition: unknown

Vaulting the barrier: An intrinsic mechanism to fuel the gas beyond the nuclear ring into the central region of barred galaxies

Kotaro Kobayashi , Naomichi Yutani , Takayuki R. Saitoh , Junichi Baba , Keiichi Wada
Authors on Pith no claims yet

Pith reviewed 2026-05-10 04:29 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords barred galaxiesnuclear ringsgas inflowAGN fuelinghydrodynamic simulationsthree-dimensional dynamicscircumnuclear gasgalactic centers
0
0 comments X

The pith

Three-dimensional gas motions enable bar-driven inflow to vault over nuclear rings and reach galactic centers.

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

The paper demonstrates through three-dimensional simulations that some gas in barred galaxies can bypass the nuclear ring by gaining vertical momentum and crossing over it to deliver material to the central tens of parsecs. This process allows gas from outside the ring to reach the inner region directly. Two-dimensional models miss this entirely and predict gas stagnation at the ring. The mechanism depends on vertical excursions of roughly 100 parsecs and loss of angular momentum, organizing the inflow into specific paths. This helps explain how centers of galaxies can be fueled for starbursts and active nuclei despite the ring acting as a barrier.

Core claim

Using three-dimensional Lagrangian hydrodynamic simulations in a fixed barred potential, gas acquires vertical momentum and vaults across the nuclear ring to reach the inner few tens of parsecs. The circumnuclear material within about 50 pc originates from gas initially located outside the ring at radii greater than or equal to 300 pc. Successful delivery requires a vertical excursion of about 100 pc before encountering the ring and substantial loss of azimuthal angular momentum. The inflow is organized within a spatial region set by the scale height of the ring gas, while most gas still accumulates near the resonance.

What carries the argument

The vaulting pathway in which inflowing gas gains vertical momentum to cross over the nuclear ring in three dimensions.

If this is right

  • Central fueling can occur independently of secondary inflow from the nuclear ring itself.
  • Nuclear rings continue to form and support star formation, but a fraction of gas penetrates inward.
  • The process is organized rather than chaotic, confined to regions determined by the ring's vertical structure.
  • This intrinsic three-dimensional effect can contribute to misalignments observed in nuclear disks.

Where Pith is reading between the lines

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

  • High-resolution kinematic observations could trace whether inner gas shows origins from large radii.
  • Similar vaulting might occur in other dynamical environments with resonances and vertical motions.
  • Two-dimensional models systematically underestimate the rate of central gas delivery in barred systems.

Load-bearing premise

The simulations assume a fixed barred potential and initial gas conditions that permit vertical excursions of order 100 parsecs before the gas encounters the nuclear ring.

What would settle it

Observations of the velocity field or chemical composition around nuclear rings that either confirm or rule out direct trajectories from outside the ring to the central 50 parsecs.

Figures

Figures reproduced from arXiv: 2604.17955 by Junichi Baba, Keiichi Wada, Kotaro Kobayashi, Naomichi Yutani, Takayuki R. Saitoh.

Figure 1
Figure 1. Figure 1: (a) Face-on gas surface-density map at t = 180 Myr with the orbit of the representative particle that reaches the CND (red curve). The x-axis aligns with the bar major axis and the bar rotates counterclockwise. (b) Slice along the line connecting the endpoints of the orange square in panel (a). The horizontal and vertical axes correspond to the direction along the gas flow and perpendicular to the disk, re… view at source ↗
Figure 2
Figure 2. Figure 2: (a) Surface-density maps at t = 200 Myr for the pseudo-2D (left) and 3D, isothermal (right) models. The red circle marks the CND region (r < 50 pc). (b) Radial distributions at t = 150 Myr of gas particles that are in the CND at t = 200 Myr, obtained by tracing those particles back to t = 150 Myr. Blue shows particles already inside the CND at t = 150 Myr, while red shows particles that were outside the CN… view at source ↗
Figure 3
Figure 3. Figure 3: (a) Time evolution of the CND mass (r < 50 pc) in the fiducial calculation that includes heating, cooling, and star formation. (b) CND mass evolution in the isothermal-3D and pseudo-2D models. Self-gravity of the gas is solved in both models. The blue and red curves correspond to the 3D isothermal and pseudo-2D runs, respectively. To enable a direct statistical comparison, we randomly select 8316 of the ri… view at source ↗
Figure 4
Figure 4. Figure 4: (a) Heat maps of the orbital distribution in the Lz −z plane for CND gas (left) and ring gas (right). The color scale shows the particle fraction in each bin (logarithmic scale) and only bins with a fraction > 1×10−4 are plotted. White lines trace representative particle trajectories. Red arrows indicate the direction of particle motion. The red dashed ellipse indicates the nuclear ring region. (b) Left: f… view at source ↗
Figure 5
Figure 5. Figure 5: (a) Histograms of Lz (left) and |z| (right) when CND and ring particles last crossed R = 380 pc. The red dashed lines show the time-averaged half-mass height and Lz inside the ring region. (b) Scatter plot of the same data in the Lz–|z| plane. Red and black points correspond to CND and ring particles, respectively; the dashed lines again show the ring averages. (c) Locations of events where the magnitude o… view at source ↗
read the original abstract

Gas delivery to galactic centers powers nuclear starbursts and active galactic nuclei (AGNs), yet bar-driven inflow is generally expected to stall in a nuclear ring a few hundred parsecs across. Using three-dimensional Lagrangian hydrodynamic simulations in a fixed barred potential, we identify a bypass channel in which a fraction of the inflowing gas acquires vertical momentum, vaults across the ring, and reaches the inner few tens of parsecs. This pathway is absent in two-dimensional calculations, which instead predict long-lived stagnation at the ring. We find that the circumnuclear material within $\sim 50$ pc originates from gas initially located outside the ring ($\gtrsim 300$ pc), rather than from secondary inflow out of the ring itself. Successful delivery requires both a sufficiently large vertical excursion, $|z| \sim 100$ pc before encountering the ring, and substantial loss of azimuthal angular momentum $L_z$. The resulting inflow is organized rather than chaotic: center-reaching trajectories are confined to a limited spatial region set by the scale height of the ring gas. Most bar-driven gas still accumulates near the resonance and fuels star formation in the nuclear ring, but the vaulting stream selects a modest yet sufficient fraction that penetrates to the circumnuclear disk. These results suggest that intrinsically three-dimensional gas motions help link nuclear starbursts, AGN fueling, and the frequent misalignment of nuclear disks with respect to their host galaxies.

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.

Referee Report

2 major / 2 minor

Summary. The paper uses three-dimensional Lagrangian hydrodynamic simulations in a fixed barred galactic potential to identify a 'vaulting' bypass mechanism. Gas parcels acquire vertical momentum, cross over the nuclear ring (typically at a few hundred parsecs), and reach the central ~50 pc region. This channel is absent in equivalent 2D runs, where gas stalls at the ring. Trajectory analysis shows that circumnuclear material originates from gas initially located outside the ring (≳300 pc) rather than from secondary inflow out of the ring, provided |z| excursions reach ~100 pc and substantial Lz is lost. Most gas still accumulates in the ring, but the vaulting stream supplies a modest fraction to the center.

Significance. If the result holds, the work supplies a concrete, intrinsically three-dimensional channel that can link bar-driven inflow to nuclear fueling despite the ring barrier. It offers a natural explanation for observed nuclear starbursts, AGN activity, and frequent misalignments between nuclear disks and host galaxies. The use of Lagrangian trajectory tracking to establish gas origins is a clear strength.

major comments (2)
  1. [Simulation setup] Simulation setup (methods section implied by abstract): the calculations adopt a fixed barred potential with no gas self-gravity. This choice is load-bearing for the vaulting claim, because self-gravity can modify vertical velocity dispersion, ring thickness, and the efficiency of bar-driven vertical resonances or buckling that enable |z| ~ 100 pc excursions before ring encounter. The 2D–3D comparison does not test whether the bypass survives when self-gravity is restored.
  2. [Results on gas origins] Gas-origin claim (abstract and results): the assertion that material within ~50 pc comes from initial radii ≳300 pc rests on specific initial conditions that permit large vertical excursions prior to ring crossing. Without reported resolution studies, convergence tests, or variation of the vertical scale height, it is unclear how robust the fraction of vaulting trajectories is to numerical choices.
minor comments (2)
  1. [Abstract] The abstract is concise but would benefit from a brief statement of the hydro code, particle number, and softening length to allow immediate assessment of the numerical setup.
  2. [Figures] Figure captions and trajectory plots should explicitly label the initial radial range of the tracked particles and the time at which they encounter the ring, to make the origin analysis easier to follow.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of our work and for identifying key points that strengthen the presentation. We respond to each major comment below, indicating revisions where appropriate.

read point-by-point responses

  1. Referee: [Simulation setup] Simulation setup (methods section implied by abstract): the calculations adopt a fixed barred potential with no gas self-gravity. This choice is load-bearing for the vaulting claim, because self-gravity can modify vertical velocity dispersion, ring thickness, and the efficiency of bar-driven vertical resonances or buckling that enable |z| ~ 100 pc excursions before ring encounter. The 2D–3D comparison does not test whether the bypass survives when self-gravity is restored.

    Authors: We agree that the fixed potential without self-gravity is a deliberate simplification chosen to isolate the purely kinematic vaulting mechanism arising from bar-driven vertical motions. Self-gravity would indeed affect vertical dispersion and ring structure, potentially altering the efficiency of |z| excursions. The 2D–3D comparison is intended only to demonstrate that the bypass requires three-dimensional motions, not to claim invariance under self-gravity. In the revised manuscript we will add explicit discussion of this limitation in the Methods and Conclusions sections, noting that the reported channel represents a baseline effect whose robustness with self-gravity remains to be tested in future work. revision: partial

  2. Referee: [Results on gas origins] Gas-origin claim (abstract and results): the assertion that material within ~50 pc comes from initial radii ≳300 pc rests on specific initial conditions that permit large vertical excursions prior to ring crossing. Without reported resolution studies, convergence tests, or variation of the vertical scale height, it is unclear how robust the fraction of vaulting trajectories is to numerical choices.

    Authors: The origin claim is derived directly from Lagrangian trajectory tracking, which shows that particles reaching the central region begin at radii ≳300 pc provided they achieve |z| ~ 100 pc before ring encounter. While the submitted manuscript does not include dedicated resolution or initial-scale-height variation tests, the mechanism is tied to the geometry of the bar potential and the ring scale height. For the revision we will add a new subsection describing the adopted resolution, perform supplementary runs with varied initial vertical scale heights, and report the resulting variation in the vaulting fraction to demonstrate robustness. revision: yes

Circularity Check

0 steps flagged

No circularity: claim is direct outcome of 3D hydro trajectories

full rationale

The paper reports results from Lagrangian hydrodynamic simulations in a fixed barred potential. The central claim—that circumnuclear gas within ~50 pc originates from initial locations ≳300 pc via vertical vaulting—is extracted from particle trajectories and 2D-vs-3D comparisons. No parameter is fitted to the target quantity and then relabeled a prediction; no equation is defined in terms of its own output; no load-bearing uniqueness theorem or ansatz is imported via self-citation. The fixed-potential setup is an explicit modeling choice whose consequences are tested numerically rather than presupposed. The derivation chain is therefore self-contained and non-circular.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review limits visibility into simulation parameters; no explicit free parameters, axioms, or invented entities are stated beyond standard hydrodynamics in a barred potential.

pith-pipeline@v0.9.0 · 5579 in / 1021 out tokens · 49643 ms · 2026-05-10T04:29:18.253427+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 34 canonical work pages

  1. [1]

    1992, MNRAS, 259, 345, doi: 10.1093/mnras/259.2.345

    Athanassoula, E. 1992, MNRAS, 259, 345, doi: 10.1093/mnras/259.2.345

    work page doi:10.1093/mnras/259.2.345 1992

  2. [2]

    1996, ARA&A, 34, 173, doi: 10.1146/annurev.astro.34.1.173

    Buta, R., & Combes, F. 1996, ARA&A, 34, 173, doi: 10.1146/annurev.astro.34.1.173

    work page doi:10.1146/annurev.astro.34.1.173 1996

  3. [3]

    C., & P´ erez-Villegas,´A

    Chaves-Velasquez, L., G´ omez, G. C., & P´ erez-Villegas,´A. 2025, PASA, 42, e014, doi: 10.1017/pasa.2024.130

    work page doi:10.1017/pasa.2024.130 2025

  4. [4]

    2021, Active Galactic Nuclei: Fueling and Feedback (IOP Publishing), doi: 10.1088/2514-3433/ac2a27

    Combes, F. 2021, Active Galactic Nuclei: Fueling and Feedback (IOP Publishing), doi: 10.1088/2514-3433/ac2a27

    work page doi:10.1088/2514-3433/ac2a27 2021

  5. [5]

    2019, Astronomy & Astrophysics (NUGA Survey related)

    Combes, F., et al. 2019, Astronomy & Astrophysics (NUGA Survey related)

    2019

  6. [6]

    , keywords =

    Dehnen, W. 2000, AJ, 119, 800, doi: 10.1086/301226 Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604, doi: 10.1038/nature03335

    work page doi:10.1086/301226 2000

  7. [7]

    2015, Monthly Notices of the Royal Astronomical Society, 446, 2468, doi: 10.1093/mnras/stu2209

    Emsellem, E., Renaud, F., Bournaud, F., et al. 2015, Monthly Notices of the Royal Astronomical Society, 446, 2468, doi: 10.1093/mnras/stu2209

    work page doi:10.1093/mnras/stu2209 2015

  8. [8]

    2000, ApJ, 528, 677, doi: 10.1086/308215

    Englmaier, P., & Shlosman, I. 2000, ApJ, 528, 677, doi: 10.1086/308215

    work page doi:10.1086/308215 2000

  9. [9]

    2000, ApJ, 529, 109, doi: 10.1086/308233 Garc´ ıa-Burillo, S., Alonso-Herrero, A., Ramos Almeida, C., et al

    Fukuda, H., Habe, A., & Wada, K. 2000, ApJ, 529, 109, doi: 10.1086/308233 Garc´ ıa-Burillo, S., Alonso-Herrero, A., Ramos Almeida, C., et al. 2021, Astronomy & Astrophysics, 652, A98 Go’mez, G. C., & Cox, D. P. 2004, ApJ, 615, 744, doi: 10.1086/424499

    work page doi:10.1086/308233 2000

  10. [10]

    1990, ApJ, 356, 359, doi: 10.1086/168845

    Hernquist, L. 1990, ApJ, 356, 359, doi: 10.1086/168845

    work page doi:10.1086/168845 1990

  11. [11]

    G., Tanaka, M., et al

    Hopkins, P. F., & Quataert, E. 2010, MNRAS, 407, 1529, doi: 10.1111/j.1365-2966.2010.17064.x —. 2011, MNRAS, 415, 1027, doi: 10.1111/j.1365-2966.2011.18542.x

    work page doi:10.1111/j.1365-2966.2010.17064.x 2010

  12. [12]

    Yoon, D., & Teuben, P. J. 2012, ApJ, 747, 60, doi: 10.1088/0004-637X/747/1/60

    work page doi:10.1088/0004-637x/747/1/60 2012

  13. [13]

    Kormendy, J., & Ho, L. C. 2013, ARA&A, 51, 511, doi: 10.1146/annurev-astro-082708- 101811

    work page doi:10.1146/annurev-astro-082708- 2013

  14. [14]

    J., Zwaan , M

    Maciejewski, W. 2004, MNRAS, 354, 892, doi: 10.1111/j.1365-2966.2004.08254.x

    work page doi:10.1111/j.1365-2966.2004.08254.x 2004

  15. [15]

    Coleman and Hamilton, Douglas P

    Maciejewski, W., Teuben, P. J., Sparke, L. S., & Stone, J. M. 2002, MNRAS, 329, 502, doi: 10.1046/j.1365-8711.2002.04957.x

    work page doi:10.1046/j.1365-8711.2002.04957.x 2002

  16. [16]

    A., & Cox, D

    Martos, M. A., & Cox, D. P. 1998, ApJ, 509, 703, doi: 10.1086/306514

    work page doi:10.1086/306514 1998

  17. [17]

    1975, PASJ, 27, 533, doi: 10.1093/pasj/27.4.533

    Miyamoto, M., & Nagai, R. 1975, PASJ, 27, 533, doi: 10.1093/pasj/27.4.533

    work page doi:10.1093/pasj/27.4.533 1975

  18. [18]

    F., Frenk, C

    Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563, doi: 10.1086/177173 18 P´ erez, I. 2007, arXiv e-prints, arXiv:0710.2855. https://arxiv.org/abs/0710.2855

    work page doi:10.1086/177173 1996

  19. [19]

    Plummer, H. C. 1911, MNRAS, 71, 460, doi: 10.1093/mnras/71.5.460

    work page doi:10.1093/mnras/71.5.460 1911

  20. [20]

    Rees, M. J. 1984, ARA&A, 22, 471, doi: 10.1146/annurev.aa.22.090184. 002351

    work page doi:10.1146/annurev.aa.22.090184 1984

  21. [21]

    W., & Teuben, P

    Regan, M. W., & Teuben, P. 2003, ApJ, 582, 723, doi: 10.1086/344721

    work page doi:10.1086/344721 2003

  22. [22]

    W., & Teuben, P

    Regan, M. W., & Teuben, P. J. 2004, ApJ, 600, 595, doi: 10.1086/380116

    work page doi:10.1086/380116 2004

  23. [23]

    How Should We Choose Star-Formation Criteria in High-Resolution Simulations of Disk Galaxies?

    Saitoh, T. R., Daisaka, H., Kokubo, E., et al. 2008, PASJ, 60, 667, doi: 10.1093/pasj/60.4.667

    work page doi:10.1093/pasj/60.4.667 2008

  24. [24]

    R., & Makino, J

    Saitoh, T. R., & Makino, J. 2013, ApJ, 768, 44, doi: 10.1088/0004-637X/768/1/44

    work page doi:10.1088/0004-637x/768/1/44 2013

  25. [25]

    H., & Huntley, J

    Sanders, R. H., & Huntley, J. M. 1976, ApJ, 209, 53, doi: 10.1086/154692

    work page doi:10.1086/154692 1976

  26. [26]

    C., & Frank, J

    Shlosman, I., Begelman, M. C., & Frank, J. 1990, Nature, 345, 679, doi: 10.1038/345679a0

    work page doi:10.1038/345679a0 1990

  27. [27]

    Shlosman, I., Frank, J., & Begelman, M. C. 1989, Nature, 338, 45

    1989

  28. [28]

    C., Tress, R

    Sormani, M. C., Tress, R. G., Glover, S. C. O., et al. 2020, Monthly Notices of the Royal Astronomical Society, 497, 5024, doi: 10.1093/mnras/staa1999

    work page doi:10.1093/mnras/staa1999 2020

  29. [29]

    C., Treß, R

    Sormani, M. C., Treß, R. G., Ridley, M., et al. 2018, Monthly Notices of the Royal Astronomical Society, 475, 2383, doi: 10.1093/mnras/stx3258

    work page doi:10.1093/mnras/stx3258 2018

  30. [30]

    C., Treß, R

    Sormani, M. C., Treß, R. G., Glover, S. C. O., et al. 2019, Monthly Notices of the Royal Astronomical Society, 488, 4663, doi: 10.1093/mnras/stz2054

    work page doi:10.1093/mnras/stz2054 2019

  31. [31]

    2025, The Astrophysical Journal, 984, 109

    Su, Y., Zhang, S., Sun, Y., et al. 2025, The Astrophysical Journal, 984, 109

    2025

  32. [32]

    G., Sormani, M

    Tress, R. G., Sormani, M. C., Glover, S. C. O., et al. 2020, Monthly Notices of the Royal Astronomical Society, 499, 4455, doi: 10.1093/mnras/staa3120

    work page doi:10.1093/mnras/staa3120 2020

  33. [33]

    G., Sormani, M

    Tress, R. G., Sormani, M. C., Girichidis, P., et al. 2024, Astronomy & Astrophysics, 691, A303, doi: 10.1051/0004-6361/202450035

    work page doi:10.1051/0004-6361/202450035 2024

  34. [34]

    1992, MNRAS, 258, 82, doi: 10.1093/mnras/258.1.82 —

    Wada, K., & Habe, A. 1992, MNRAS, 258, 82, doi: 10.1093/mnras/258.1.82 —. 1995, MNRAS, 277, 433, doi: 10.1093/mnras/277.2.433

    work page doi:10.1093/mnras/258.1.82 1992

  35. [35]

    P., & Spaans, M

    Wada, K., Papadopoulos, P. P., & Spaans, M. 2009, ApJ, 702, 63, doi: 10.1088/0004-637X/702/1/63

    work page doi:10.1088/0004-637x/702/1/63 2009

  36. [36]

    A., & Cox, D

    Walters, M. A., & Cox, D. P. 2001, The Astrophysical Journal, 549, 353, doi: 10.1086/319057

    work page doi:10.1086/319057 2001

  37. [37]

    M., & Bonner, G

    Widrow, L. M., & Bonner, G. 2015, Monthly Notices of the Royal Astronomical Society, 450, 266, doi: 10.1093/mnras/stv574

    work page doi:10.1093/mnras/stv574 2015