Invariant manifolds in barred galaxy simulations. I. Material density waves

Merc\`e Romero-G\'omez; Santi Roca-F\`abrega; Toni Soler-Terricabras

arxiv: 2604.26742 · v1 · submitted 2026-04-29 · 🌌 astro-ph.GA

Invariant manifolds in barred galaxy simulations. I. Material density waves

Toni Soler-Terricabras , Merc\`e Romero-G\'omez , Santi Roca-F\`abrega This is my paper

Pith reviewed 2026-05-07 10:36 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords barred galaxiesspiral armsinvariant manifoldsdensity wavesN-body simulationsJacobi energygalactic dynamicsradial migration

0 comments

The pith

In barred galaxy simulations, spiral arms are material structures built by particles on invariant-manifold orbits while the surrounding disc responds as density waves.

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

The paper tests whether invariant manifold theory can separate the particles that actually build and sustain spiral arms in an N-body model of an isolated barred galaxy from those that merely respond to them. By sorting disc particles according to their Jacobi energy relative to the bar's equilibrium points, the authors isolate a manifold-compatible population that starts near the bar and follows transit orbits along the arms. These particles alone produce the observed outward-migrating ridge in the radius-velocity plane and the characteristic streaming motions. The remaining low-energy particles stay on nearly circular orbits but feel the spiral's self-gravity and add to the density contrast in the manner of classical density waves. The combined picture is that the arms themselves are material features kept in place by manifold-guided transport, while the bulk disc behaves as a system of material density waves.

Core claim

Barred spiral arms emerge as material structures sustained by manifold-guided transport, with the surrounding disc behaving as a system of material density waves. The Jacobi constant cleanly separates three kinematic populations: low-energy particles on nearly circular orbits that dominate the disc, high-energy particles on banana orbits, and manifold-compatible particles on transit orbits that originate near the bar and trace the arms. Only the last group generates the prominent outward-migrating ridge in the R-v_phi plane and reproduces the spiral streaming pattern, while the low-energy population enhances the density contrast through small perturbations consistent with traditional density

What carries the argument

Jacobi-energy classification of disc particles relative to the energies of the bar's equilibrium points, which isolates the subset of orbits compatible with invariant manifolds and therefore able to follow the spiral arms.

If this is right

Only the manifold-compatible population produces the outward-migrating ridge observed in the R-v_phi plane.
These particles alone reproduce the characteristic spiral streaming motions.
The low-energy population shows global quasi-circular motion with small perturbations induced by the spiral's self-gravity.
The arms are material structures kept in place by manifold-guided radial transport rather than pure wave patterns.

Where Pith is reading between the lines

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

Kinematic maps of real galaxies could be searched for the same Jacobi-energy-separated populations to test whether observed spirals show the predicted manifold-driven signatures.
The framework suggests that radial migration in barred spirals is concentrated in the manifold-compatible orbits rather than distributed across the whole disc.
If the separation holds in other simulations with different bar strengths or gas content, it would provide a dynamical criterion for distinguishing material arms from density-wave arms in observations.

Load-bearing premise

The Jacobi-energy classification relative to the bar's equilibrium points reliably isolates orbits compatible with invariant manifolds inside a fully self-gravitating, time-evolving N-body simulation.

What would settle it

Re-run the same simulation after removing or re-assigning the manifold-compatible particles and check whether the spiral arms and their streaming ridge in the R-v_phi plane disappear or weaken.

Figures

Figures reproduced from arXiv: 2604.26742 by Merc\`e Romero-G\'omez, Santi Roca-F\`abrega, Toni Soler-Terricabras.

**Figure 1.** Figure 1: Time evolution of the bar and spiral arms diagnostics for the B1 model, following the method of Dehnen et al. (2023). Time is measured from the start of the simulation. (a): Fourier m = 2 amplitudes, A2, for the bar (dotted red; averaged over [R0, R1] kpc, as defined in Appendix B of Dehnen et al. (2023)) and for the spiral arms (solid red; averaged from R1 to 10 kpc). The dotted and solid horizontal lines… view at source ↗

**Figure 2.** Figure 2: b, aligned with the bar’s semi-minor axis. As expected, L3 corresponds to a local minimum of the effective potential and is linearly stable, whereas L4 and L5 are located at local maxima and are also stable throughout the analysed time interval. Therefore, in our model the dominant contribution to the spirals is directly provided by the manifolds associated with the Lyapunov orbits around L1 and L2 view at source ↗

**Figure 3.** Figure 3: Temporal evolution of the Jacobi energies at the unstable equilibrium points. (a): Jacobi energy at the Lagrange points L1–L2 (red) and L4–L5 (blue) as a function of time. Shaded regions indicate the dispersion associated with the determination of the Jacobi energies. (b): Energy separation ∆EJ = EL4,5 − EL1,2 as a function of time. The grey band marks the corresponding dispersion. The maximum value of t… view at source ↗

**Figure 4.** Figure 4: Dynamical diagnostics for the different kinematic populations at t = 1.141 Gyr after the beginning of the simulation. Left column: lowenergy population. Middle column: manifold-compatible population; Right column: high-energy population. First row: Face-on particle distribution colour-coded by their surface density. The dashed circles indicate the extent of the bar. Second row: Tangential rotation curve, … view at source ↗

**Figure 5.** Figure 5: Temporal evolution of the relative contribution of each dynamical population, colour-coded as follows: low-energy (pink), manifoldcompatible (blue), high-energy (orange). Article number, page 7 of 11 view at source ↗

read the original abstract

We investigate the dynamical origin and kinematic signatures of spiral structure in an N-body simulation of an isolated barred galaxy, assessing whether invariant manifold theory provides a consistent dynamical framework to disentangle the disc particle populations and to identify those that genuinely build, trace, and sustain the spiral arms. We compute the Jacobi energy of disc particles and classify them relative to the energies of the equilibrium points, thereby isolating manifold-compatible orbits. We analyse their spatial distribution and velocity structure to characterise spiral-related streaming motions. The Jacobi constant provides a physically motivated dynamical separator that reveals three distinct kinematic populations: (i) low-energy particles on nearly circular orbits populating most of the disc, (ii) high-energy particles associated with banana orbits, and (iii) manifold-compatible particles originating near the bar and following transit orbits along the spiral arms. Only the manifold-compatible population generates the prominent outward-migrating ridge observed in the R - v_phi plane and reproduces the characteristic spiral streaming pattern. In contrast, the low-energy population exhibits a global quasi-circular motion with small perturbations induced by the self-gravity of the spiral structure. Our results demonstrate that the spiral arms are dynamically traced by the manifold-compatible population, which forms the backbone of the structure and drives effective radial transport. The bulk of low-energy disc particles responds to the spiral perturbation similarly to the traditional density wave picture, enhancing the density contrast caused by the invariant-manifold compatible particles. In this framework, barred spiral arms emerge as material structures sustained by manifold-guided transport, with the surrounding disc behaving as a system of material density waves.

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

Manifold theory applied to this N-body run gives a plausible split for the spiral ridge, but the time-dependent potential leaves the classification on shaky ground.

read the letter

Your colleague should know that this work takes invariant manifold theory and applies it to classify particles in a self-consistent N-body simulation of a barred galaxy. They compute Jacobi energies relative to the bar's L1 and L2 points and separate the disc into low-energy nearly circular orbits, high-energy banana orbits, and a manifold-compatible transit population. The key result is that only the manifold group shows the outward-migrating ridge in the R-v_phi plane and the spiral streaming motions, while the low-energy particles have small perturbations that boost the density contrast in a density-wave manner. This leads to their picture of barred spirals as material structures sustained by manifold transport with the rest of the disc responding passively. What the paper does well is extract specific kinematic signatures from the simulation and tie them back to the dynamical classification. The separation reveals distinct behaviors that match some observed features in galaxies, and they present it as a hybrid between manifold and density wave ideas. The abstract is clear on the populations and their roles. The main soft spot is the use of Jacobi energy in a time-evolving N-body potential. The classification assumes something close to a fixed potential for the energy to isolate manifold orbits reliably, but here the potential fluctuates as the bar and spirals grow and change. Without reported checks on how long the selected particles stay in their groups or whether they actually follow the manifold paths over several bar periods, the separation could be mixing things. The abstract also skips details on the exact manifold computation method or any robustness tests against potential variations. This is for galactic dynamicists who follow bar and spiral structure work. Someone looking for simulation-based tests of analytic orbit theory would find the kinematic analysis useful. It deserves peer review because the approach is worth exploring further, even with the current limitations on the evidence. I recommend sending it out for refereeing, with requests for the missing stability tests and more on how the manifolds are identified in the live run.

Referee Report

2 major / 1 minor

Summary. The manuscript analyzes an N-body simulation of an isolated barred galaxy and classifies disc particles by their Jacobi energy relative to the bar's L1/L2 equilibrium points. This isolates three populations: low-energy particles on nearly circular orbits, high-energy particles on banana orbits, and manifold-compatible particles on transit orbits. The central claim is that only the manifold-compatible population traces and sustains the spiral arms as material structures, producing the observed outward-migrating ridge in the R-v_φ plane and characteristic streaming motions, while the low-energy population responds to the spiral perturbation in a manner consistent with traditional density-wave theory and merely enhances the density contrast.

Significance. If the classification and dynamical separation hold, the work supplies a physically motivated framework that reconciles material and wave interpretations of barred spirals, with direct implications for radial transport and kinematic signatures. The analysis relies on direct numerical computation of Jacobi energies and orbit classification from the simulation output, drawing on established dynamical-systems theory without introducing fitted parameters or ad-hoc reductions.

major comments (2)

[Methods (classification procedure)] The classification of manifold-compatible orbits rests on computing Jacobi energies relative to the bar's equilibrium points in simulation snapshots (as described in the abstract and the methods section on orbit classification). However, the underlying N-body simulation is fully self-gravitating and time-evolving, so the potential is not autonomous and Jacobi energy is not conserved; no explicit test of classification stability across multiple bar rotations or verification that selected particles follow manifold geometry is provided, which directly undermines the isolation of the transit-orbit population.
[Results (R - v_φ analysis)] In the results on the R - v_φ plane and spiral streaming (corresponding to the claims about the outward-migrating ridge and kinematic populations), the assertion that only the manifold-compatible particles generate the prominent ridge and drive effective radial transport is presented qualitatively. No quantitative decomposition is given showing the fractional contribution of each population to the ridge amplitude or density contrast, leaving the 'backbone' claim without a direct metric of dominance.

minor comments (1)

[Abstract] The abstract introduces 'banana orbits' and 'transit orbits' without a brief definition or reference to their standard properties in the barred-galaxy context, which would aid readability for a general audience.

Simulated Author's Rebuttal

2 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 incorporate to strengthen the presentation.

read point-by-point responses

Referee: The classification of manifold-compatible orbits rests on computing Jacobi energies relative to the bar's equilibrium points in simulation snapshots (as described in the abstract and the methods section on orbit classification). However, the underlying N-body simulation is fully self-gravitating and time-evolving, so the potential is not autonomous and Jacobi energy is not conserved; no explicit test of classification stability across multiple bar rotations or verification that selected particles follow manifold geometry is provided, which directly undermines the isolation of the transit-orbit population.

Authors: We agree that the fully self-gravitating and time-evolving nature of the N-body simulation implies that the potential is not autonomous and Jacobi energy is not exactly conserved. Our classification is performed snapshot-by-snapshot using the instantaneous bar potential and equilibrium points at each output time. To address the concern, the revised manuscript will include an explicit test tracking the Jacobi energies of the classified particles over several bar rotation periods, demonstrating the stability of the manifold-compatible assignment with only minor variations. We will additionally show that the selected particles follow trajectories consistent with manifold geometry when integrated in the time-averaged potential. revision: yes
Referee: In the results on the R - v_φ plane and spiral streaming (corresponding to the claims about the outward-migrating ridge and kinematic populations), the assertion that only the manifold-compatible particles generate the prominent ridge and drive effective radial transport is presented qualitatively. No quantitative decomposition is given showing the fractional contribution of each population to the ridge amplitude or density contrast, leaving the 'backbone' claim without a direct metric of dominance.

Authors: The manuscript figures show that the outward-migrating ridge appears exclusively among the manifold-compatible particles, while the low-energy population exhibits only small perturbations around circular motion and the high-energy population remains on banana orbits. We acknowledge that a quantitative metric would provide stronger support for the dominance claim. In the revised version we will add a decomposition of the ridge region, reporting the fractional contribution of each population to both the particle count in the ridge and the local density contrast, thereby supplying a direct numerical measure of the manifold-compatible population's role. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper computes Jacobi energies directly from N-body simulation snapshots, classifies particles relative to bar equilibrium points, and reports observed spatial/kinematic differences among the resulting populations. This is post-processing of simulation output using standard dynamical-systems separators; no parameter is fitted to a data subset and then renamed as a prediction, no central premise reduces to a self-citation chain, and no ansatz or uniqueness theorem is smuggled in. The claims about which population traces the arms follow from the classified particles' measured distributions rather than from definitional equivalence to the inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Central claim depends on applicability of invariant manifold theory to the N-body system and validity of energy-based orbit classification in a self-gravitating disc.

axioms (2)

domain assumption The rotating frame tied to the bar allows computation of a conserved Jacobi energy and identification of equilibrium points whose manifolds govern particle motion.
Invoked when classifying particles relative to the energies of the equilibrium points.
domain assumption Invariant manifold structures identified in the bar potential remain meaningful guides for particle orbits even when the full disc self-gravity is included.
Required to treat manifold-compatible particles as the backbone of the spiral arms.

pith-pipeline@v0.9.0 · 9923 in / 1283 out tokens · 85613 ms · 2026-05-07T10:36:03.295200+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages

[1]

, " * write output.state after.block = add.period write newline

ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

work page
[2]

write newline

" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

work page
[3]

@esa (Ref

\@ifclassloaded agu2001 natbib The agu2001 class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later remove the command natbib from the document \@ifclassloaded aguplus natbib The aguplus class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later rem...

work page
[4]

@stdbsttrue NAT@ctr \@lbibitem[ NAT@ctr ] \@lbibitem[#1]#2 \@extra@b@citeb \@ifundefined br@#2\@extra@b@citeb \@namedef br@#2 \@nameuse br@#2\@extra@b@citeb \@ifundefined b@#2\@extra@b@citeb @num @parse #2 [ @natanchorstart #2\@extra@b@citeb \@biblabel @num @natanchorend] @ifcmd#1(@)(@)\@nil #2 @lbibitem\@undefined @lbibitem\@lbibitem \@lbibitem[#1]#2 @lb...

work page
[5]

@open @close @open @close and [1] URL: #1 \@ifundefined chapter * \@mkboth \@ifundefined NAT@sectionbib * \@mkboth * \@mkboth\@gobbletwo \@ifclassloaded amsart * \@ifclassloaded amsbook * \@ifundefined bib@heading @heading NAT@ctr thebibliography [1] @ \@biblabel NAT@ctr \@bibsetup #1 NAT@ctr 0 @openbib .11em \@plus.33em \@minus.07em 4000 4000 `\.=1000 \@...

work page

[1] [1]

, " * write output.state after.block = add.period write newline

ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

work page

[2] [2]

write newline

" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

work page

[3] [3]

@esa (Ref

\@ifclassloaded agu2001 natbib The agu2001 class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later remove the command natbib from the document \@ifclassloaded aguplus natbib The aguplus class already includes natbib coding, so you should not add it explicitly Type <Return> for now, but then later rem...

work page

[4] [4]

@stdbsttrue NAT@ctr \@lbibitem[ NAT@ctr ] \@lbibitem[#1]#2 \@extra@b@citeb \@ifundefined br@#2\@extra@b@citeb \@namedef br@#2 \@nameuse br@#2\@extra@b@citeb \@ifundefined b@#2\@extra@b@citeb @num @parse #2 [ @natanchorstart #2\@extra@b@citeb \@biblabel @num @natanchorend] @ifcmd#1(@)(@)\@nil #2 @lbibitem\@undefined @lbibitem\@lbibitem \@lbibitem[#1]#2 @lb...

work page

[5] [5]

@open @close @open @close and [1] URL: #1 \@ifundefined chapter * \@mkboth \@ifundefined NAT@sectionbib * \@mkboth * \@mkboth\@gobbletwo \@ifclassloaded amsart * \@ifclassloaded amsbook * \@ifundefined bib@heading @heading NAT@ctr thebibliography [1] @ \@biblabel NAT@ctr \@bibsetup #1 NAT@ctr 0 @openbib .11em \@plus.33em \@minus.07em 4000 4000 `\.=1000 \@...

work page