pith. sign in

arxiv: 2605.19882 · v1 · pith:WWVWBHNSnew · submitted 2026-05-19 · 🌌 astro-ph.GA · math.DS

Arm morphology in off-centre barred galaxies

P. S\'anchez-Mart\'in , M. L\'opez-Vilamaj\'o , M. Romero-G\'omez , J. J. Masdemont This is my paper

Pith reviewed 2026-05-20 03:38 UTC · model grok-4.3

classification 🌌 astro-ph.GA math.DS
keywords barred galaxiesoff-centre barsinvariant manifoldsLagrangian pointsarm morphologygalactic dynamicsasymmetric galaxiesLarge Magellanic Cloud
0
0 comments X

The pith

Off-centre asymmetric bars distort invariant manifolds to produce arms of unequal density and shape, or a single arm after a critical offset triggers a bifurcation.

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

The paper builds a model of a barred galaxy with an internally lopsided and displaced bar to see how these features change the points of equilibrium and the paths that stars follow to form spiral arms. Small displacements that keep the system's centre of mass inside the bar keep the usual five equilibrium points but warp the associated paths so that the two arms differ in density and shape. When the displacement grows large enough for the centre of mass to leave the bar, a bifurcation removes two of the unstable points and leaves only one arm-forming structure. This supplies a dynamical account for the lopsided arms observed in galaxies such as the Large Magellanic Cloud and connects bar position directly to arm asymmetry. A reader would care because the same basic orbital structures appear to control arm appearance across a range of realistic bar misalignments.

Core claim

In a barred galaxy model that includes an off-centred and asymmetric bar, the positions, stability, and bifurcations of the Lagrangian equilibrium points are tracked as functions of the displacement of the asymmetric mass along the bar and the offset between the bar and the system's centre of mass. For modest offsets that keep the centre of mass inside the bar, the classical five-equilibrium-point configuration is preserved, yet the invariant manifolds of planar Lyapunov orbits around the unstable points become strongly distorted, producing two arms with different densities and shapes. At the threshold where the galactic centre of mass exits the bar ellipsoid, a pitchfork bifurcation removes

What carries the argument

Invariant manifolds of planar Lyapunov orbits around unstable Lagrangian points in an off-centred asymmetric barred potential, which trace the global paths that form spiral arms.

If this is right

  • Modest internal bar lopsidedness and bar-disc offsets that keep the centre of mass inside the bar preserve five equilibrium points but distort manifolds to yield two arms of unequal density and shape.
  • When the offset reaches the point at which the centre of mass exits the bar ellipsoid, a pitchfork bifurcation removes the collinear unstable points.
  • The resulting three-equilibrium-point configuration is supported by a single unstable point and its manifold, which forms one arm.
  • The framework is compatible with the observed correlation between off-centre bars and photometric lopsidedness and explains the strongly asymmetric arm morphology of galaxies such as the LMC.

Where Pith is reading between the lines

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

  • Arm asymmetry in other observed galaxies could serve as a photometric indicator of bar-disc misalignment even when the centre of mass position is not directly measurable.
  • Self-consistent N-body simulations that include the vertical structure of real galaxies could test whether the planar manifold predictions remain dominant or are altered by out-of-plane motions.
  • The critical offset threshold might be inverted from observed arm properties to estimate the degree of bar-disc misalignment in large photometric surveys.

Load-bearing premise

The invariant manifolds of planar Lyapunov orbits around the unstable points continue to organise arm morphology after the bar is displaced and made asymmetric, without major effects from vertical motions or time-dependent changes in the potential.

What would settle it

A numerical simulation or dynamical model of a barred galaxy with the centre of mass placed outside the bar ellipsoid that still develops two distinct arms matching the modest-offset case, or that retains five equilibrium points, would contradict the predicted bifurcation and single-arm transition.

Figures

Figures reproduced from arXiv: 2605.19882 by J. J. Masdemont, M. L\'opez-Vilamaj\'o, M. Romero-G\'omez, P. S\'anchez-Mart\'in.

Figure 1
Figure 1. Figure 1: Schematic of the lopsided galaxy configurations. Top: Model A, where the bar major axis is aligned with the system’s centre of mass (CM). Bottom: Models B and C, where the bar is offset from the CM, with Λ representing the distance from the CM to the bar major axis. The bar and the asymmetric mass component are outlined by green dashed lines; the disc contour is shown as a black dash-dotted line; the black… view at source ↗
Figure 2
Figure 2. Figure 2: Equilibrium points in the symmetric reference model. Bar and the asymmetric mass component outlined by green dashed lines. Five Lagrangian points shown as magenta dots. Contour levels of the effective potential highlight the dynamics near the unstable points L1 and L2. collinear points L1 and L2 located at the bar ends, the central point L3, and the triangular points L4 and L5. Under standard parameter val… view at source ↗
Figure 5
Figure 5. Figure 5: Equilibrium point continuation for Model C (Λ = 2.42 kpc) as a function of δ. Only three equilibrium points exist: L4, L3, and L5 (shown from top to bottom in spatial projection). Top panels: 3D curves and δ–x projection showing the high sensitivity of L4 to displacement and the change in stability of L4 relative to Models A and B. Bottom panels: δ–y projection where L4 approaches the y-coordinate of L3, a… view at source ↗
Figure 6
Figure 6. Figure 6: Equilibrium point continuation as Λ varies with δ = 0. The pitchfork bifurcation at Λ ≈ 1.5 kpc (bar semi-minor axis) marks the structural transition from Models B to C. Top panels: 3D paths in (x, y,Λ) space showing the clear coalescence, and Λ–x projection displaying how L1 and L2 approach and merge with L4 at the critical value. Bottom panels: Λ–y projection showing remaining equilibrium points primaril… view at source ↗
Figure 8
Figure 8. Figure 8: As the asymmetric mass component displaces along the major axis of the bar (increasing δ), the density distribution develops systematic asymmetry, particularly in the right region and around the unstable points. This asymmetric potential structure directly influences the manifold geometry, as discussed below. -10 0 10 x -10 -5 0 5 10 y -10 0 10 x -10 -5 0 5 10 y -10 0 10 x -10 -5 0 5 10 y -10 0 10 x -10 -5… view at source ↗
Figure 9
Figure 9. Figure 9: Model A isopotential contours (Λ = 0). Equilibrium points marked by magenta dots. From top left to bottom right: δ = 2.5,3.5,4.5,5.5 kpc. The level structure of the potential deforms with increasing δ, reflecting the asymmetric distribution of the displaced mass component. The isopotential contours ( [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Model A unstable invariant manifolds associated with Lyapunov periodic orbits of L1 and L2 (Λ = 0). Equilibrium points marked by crosses. Bar and the asymmetric mass component outlined by dotted black curves. From top left to bottom right: δ = 2.5,3.5,4.5,5.5 kpc. The arm structure retains a two-armed pattern with asymmetric arm opening, becoming more pronounced as δ increases. The unstable invariant mani… view at source ↗
Figure 11
Figure 11. Figure 11: Model A (y,y˙) projection of the intersection of plane S ≡ {x = 0} with the stable manifold associated to the Lyapunov family around L2 (top row) and L1 (bottom row). From left to right: δ = 2.5,3.5,4.5,5.5 kpc. The curves show the narrowing and stretching of the tube as the asymmetric mass component moves away from L1, marking the constriction of transport channels and differences in escape orbit populat… view at source ↗
Figure 12
Figure 12. Figure 12: Model A transit orbits resulting from integration of initial conditions from [PITH_FULL_IMAGE:figures/full_fig_p009_12.png] view at source ↗
Figure 16
Figure 16. Figure 16: Model C isodensity contours (Λ = 2.42 kpc, large offset with CM exterior to bar). Equilibrium points marked by black crosses. From top left to bottom right: δ = 2.5,3.5,4.5,5.5 kpc. The density distribution exhibits strong asymmetry. -5 0 5 10 x -5 0 5 10 y -5 0 5 10 x -5 0 5 10 y -5 0 5 10 x -5 0 5 10 y -5 0 5 10 x -5 0 5 10 y [PITH_FULL_IMAGE:figures/full_fig_p010_16.png] view at source ↗
Figure 15
Figure 15. Figure 15: Model B unstable invariant manifolds (Λ = 1.21 kpc). Bar and the asymmetric mass component outlined by dotted black curves. From top left to bottom right: δ = 2.5,3.5,4.5,5.5 kpc. The two-armed pattern is maintained, but the manifold shapes differ significantly from Model A: the manifold related to L2 acquires a distinctly different pattern compared to that of L1. by Models C, in which the centre of mass … view at source ↗
Figure 18
Figure 18. Figure 18: Model C unstable invariant manifolds associated with the Lyapunov periodic orbits of the remaining unstable point L4 (Λ = 2.42 kpc). Bar and the asymmetric mass component outlined by dotted black curves. From top left to bottom right: δ = 2.5,3.5,4.5,5.5 kpc. Only one arm is dynamically supported by the invariant manifold. of the offset between bar and disc could therefore act as an indicator of the under… view at source ↗
read the original abstract

Many barred galaxies, including the Large Magellanic Cloud (LMC), display strong lopsidedness and off-centre bars. The dynamical connection between bar-disc misalignments, internal mass asymmetries, and arm morphology is not yet fully characterised. We investigate how internal mass imbalances within the bar and global offsets between the bar and the centre of mass of the system modify the equilibrium-point structure and the invariant manifolds that organise arms. We construct a barred galaxy model which includes an off-centred and asymmetric in shape bar. Using numerical continuation, we track the position, stability, and bifurcations of the Lagrangian equilibrium points as functions of the displacement of the asymmetric mass component along the bar and of the offset between the bar and the system's centre of mass. For representative configurations we compute the invariant manifolds of planar Lyapunov orbits around unstable points and analyse the resulting arm structures. Internal bar lopsidedness and modest bar-disc offsets that keep the centre of mass inside the bar preserve the classical configuration with five equilibrium points, but strongly distort the associated invariant manifolds, producing two arms with different densities and shapes. The bar-disc offset reaches a threshold at the point at which the galactic centre of mass exits the bar ellipsoid, in which a pitchfork bifurcation removes the collinear unstable points and the system transitions to a three-equilibrium-point configuration in which a single unstable point and its associated manifold supports one arm. This framework is compatible with the observed correlation between off-centre bars and photometric lopsidedness, and it provides a dynamical explanation for the strongly asymmetric arm morphology of galaxies such as the LMC.

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

1 major / 0 minor

Summary. The paper constructs a barred galaxy potential with an off-centred and internally asymmetric bar, then uses numerical continuation to follow the positions and stabilities of the Lagrangian equilibria as functions of bar displacement and centre-of-mass offset. For representative parameter values it computes the invariant manifolds of planar Lyapunov orbits around the remaining unstable points and shows that modest offsets produce two arms of unequal density and shape while a threshold offset triggers a pitchfork bifurcation that eliminates the collinear unstable points, leaving a single unstable equilibrium whose manifold supports one arm. The results are presented as a dynamical explanation for the observed correlation between off-centre bars and lopsided arm morphology, including in the LMC.

Significance. If the planar-manifold organisation of arms survives the symmetry breaking, the work supplies a concrete dynamical pathway from bar-disc misalignment to asymmetric spiral structure and thereby links photometric lopsidedness to the equilibrium-point topology. The numerical-continuation approach is parameter-free once the potential form is fixed, and the reported bifurcation threshold is a falsifiable prediction that could be tested against observed bar offsets.

major comments (1)
  1. The central claim that the invariant manifolds of planar Lyapunov orbits continue to organise global arm density and shape after the bar is displaced and made lopsided rests on an untested extrapolation. No 3-D orbit sampling, vertical-frequency analysis, or comparison with a live N-body disc is described that would confirm the planar structures remain dominant once symmetry is broken; without such a check the reported transition to a single-arm regime cannot be regarded as robust.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for identifying a key limitation in the scope of the analysis. We respond to the major comment below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: The central claim that the invariant manifolds of planar Lyapunov orbits continue to organise global arm density and shape after the bar is displaced and made lopsided rests on an untested extrapolation. No 3-D orbit sampling, vertical-frequency analysis, or comparison with a live N-body disc is described that would confirm the planar structures remain dominant once symmetry is broken; without such a check the reported transition to a single-arm regime cannot be regarded as robust.

    Authors: We agree that the present study is restricted to the planar potential and does not contain three-dimensional orbit sampling, vertical-frequency analysis, or live N-body comparisons. The work deliberately isolates the effect of bar displacement and internal asymmetry on the planar equilibrium points and their invariant manifolds using numerical continuation. This approach yields a concrete, parameter-controlled demonstration of how modest offsets distort the two-arm manifold structure and how a threshold offset triggers a pitchfork bifurcation to a single-arm configuration. We acknowledge that vertical motions and self-gravity could alter the dominance of these planar structures in a fully three-dimensional, live disc. In the revised manuscript we will add an explicit limitations subsection that states the planar assumption, notes the absence of 3-D validation, and outlines how the reported bifurcation threshold could be tested in future N-body or 3-D orbit studies. We maintain that the planar mechanism remains a useful and falsifiable dynamical pathway even if it is later modified by three-dimensional effects. revision: partial

Circularity Check

0 steps flagged

Numerical continuation of equilibria and manifolds yields self-contained results with no circular reduction.

full rationale

The derivation proceeds by constructing an explicit parameterized potential for an off-centre asymmetric bar, then applying numerical continuation to locate and track the Lagrangian points and their bifurcations as functions of the displacement and offset parameters. Invariant manifolds are computed directly from planar Lyapunov orbits around the resulting unstable equilibria for representative cases. Arm morphologies emerge as direct outputs of these manifold integrations rather than as quantities fitted to or defined by the same data. No self-citations, ansatzes smuggled via prior work, or renaming of known results appear in the load-bearing steps; the framework remains independent of external fitted inputs and is therefore self-contained.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on standard Newtonian gravity and the assumption that a rigid, time-independent potential with an offset asymmetric bar component is sufficient to capture the equilibrium structure; no new physical entities are introduced.

free parameters (2)
  • bar displacement parameter
    The offset between the bar and the system center of mass is varied continuously as a control parameter in the numerical continuation.
  • asymmetric mass component displacement
    The position of the lopsided mass distribution inside the bar is treated as an independent control parameter.
axioms (2)
  • domain assumption The gravitational potential is time-independent and planar motion dominates arm organization
    Invoked when restricting analysis to planar Lyapunov orbits and their invariant manifolds.
  • standard math Numerical continuation accurately tracks stability changes and bifurcations without missing branches
    Implicit in the use of continuation methods to follow equilibrium points as parameters vary.

pith-pipeline@v0.9.0 · 5839 in / 1612 out tokens · 28552 ms · 2026-05-20T03:38:27.735911+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    Using numerical continuation, we track the position, stability, and bifurcations of the Lagrangian equilibrium points as functions of the displacement of the asymmetric mass component along the bar and of the offset between the bar and the system's centre of mass. For representative configurations we compute the invariant manifolds of planar Lyapunov orbits around unstable points

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages

  1. [1]

    2012, MNRAS: Letters, 426, L46

    Athanassoula, E. 2012, MNRAS: Letters, 426, L46

    work page 2012

  2. [2]

    1983, A&A, 127, 349

    Athanassoula, E., Bienayme, O., Martinet, L., & Pfenniger, D. 1983, A&A, 127, 349

    work page 1983

  3. [3]

    1997, MNRAS, 286, 284

    Athanassoula, E., Puerari, I., & Bosma, A. 1997, MNRAS, 286, 284

    work page 1997

  4. [4]

    & Tremaine, S

    Binney, J. & Tremaine, S. 2011, Galactic Dynamics, 2nd edn. (Princeton, NJ: Princeton University Press)

    work page 2011

  5. [5]

    & Athanassoula, E

    Colin, J. & Athanassoula, E. 1989, A&A, 214, 99

    work page 1989

  6. [6]

    & Schilder, F

    Dankowicz, H. & Schilder, F. 2013, Recipes for continuation (SIAM) de Swardt, B., Sheth, K., Kim, T., et al. 2015, ApJ, 808, 90 De Vaucouleurs, G. & Freeman, K. 1972, Vistas in Astronomy, 14, 163 de Vaucouleurs, G. & Freeman, K. C. 1970, in IAU Symposium, V ol. 38, The Spiral Structure of our Galaxy, ed. W. Becker & G. I. Kontopoulos, 356

    work page 2013

  7. [7]

    2023, MNRAS, 518, 2712

    Dehnen, W., Semczuk, M., & Schönrich, R. 2023, MNRAS, 518, 2712

    work page 2023

  8. [8]

    2010, The European Physical Journal Special Topics, 186, 91

    Efthymiopoulos, C. 2010, The European Physical Journal Special Topics, 186, 91

    work page 2010

  9. [9]

    2019, MNRAS, 484, 1487 El Youssoufi, D., Cioni, M.-R

    Efthymiopoulos, C., Kyziropoulos, P., Páez, R., Zouloumi, K., & Gravvanis, G. 2019, MNRAS, 484, 1487 El Youssoufi, D., Cioni, M.-R. L., Bell, C. P., et al. 2019, MNRAS, 490, 1076

    work page 2019

  10. [10]

    B., Frogel, J

    Eskridge, P. B., Frogel, J. A., Pogge, R. W., et al. 2002, ApJS, 143, 73

    work page 2002

  11. [11]

    & Masdemont, J

    Gidea, M. & Masdemont, J. J. 2007, International journal of bifurcation and chaos, 17, 1151

    work page 2007

  12. [12]

    Grady, J., Belokurov, V ., & Evans, N. W. 2021, ApJ, 909, 150

    work page 2021

  13. [13]

    2016, MNRAS, 459, 3419

    Harsoula, M., Efthymiopoulos, C., & Contopoulos, G. 2016, MNRAS, 459, 3419

    work page 2016

  14. [14]

    M., Skowron, D

    Jacyszyn-Dobrzeniecka, A. M., Skowron, D. M., Mróz, P., et al. 2017, Acta Astron., 67, 1 Jiménez-Arranz, Ó., Horta, D., van der Marel, R., et al. 2025, A&A, 698, A88 Article number, page 11 of 13 A&A proofs:manuscript no. Off_centre_galaxies Jiménez-Arranz, Ó., Roca-Fàbrega, S., Romero-Gómez, M., et al. 2024, A&A, 688, A51 Jiménez-Arranz, Ó., Romero-Gómez...

    work page 2017

  15. [15]

    Jog, C. J. 1997, ApJ, 488, 642

    work page 1997

  16. [16]

    Jog, C. J. & Combes, F. 2009, Physics Reports, 471, 75

    work page 2009

  17. [17]

    J., Lintott, C

    Kruk, S. J., Lintott, C. J., Simmons, B. D., et al. 2017, MNRAS, 469, 3363

    work page 2017

  18. [18]

    Levine, S. E. & Sparke, L. S. 1998, ApJ, 496, L13 Łokas, E. L. 2021, A&A, 655, A97 Martínez-García, E. E. 2011, ApJ, 744, 92

    work page 1998

  19. [19]

    & Nagai, R

    Miyamoto, M. & Nagai, R. 1975, PASJ, 27, 533

    work page 1975

  20. [20]

    A., D’Onghia, E., Athanassoula, E., Wilcots, E

    Pardy, S. A., D’Onghia, E., Athanassoula, E., Wilcots, E. M., & Sheth, K. 2016, ApJ, 827, 149

    work page 2016

  21. [21]

    Patra, N. N. & Jog, C. J. 2019, MNRAS, 488, 4942

    work page 2019

  22. [22]

    2023, A&A, 673, A36

    Pfenniger, D., Saha, K., & Wu, Y .-T. 2023, A&A, 673, A36

    work page 2023

  23. [23]

    Plummer, H. C. 1911, MNRAS, 71, 460

    work page 1911

  24. [24]

    2022, MNRAS, 512, 563 Romero-Gómez, M., Athanassoula, E., Masdemont, J

    Ripepi, V ., Chemin, L., Molinaro, R., et al. 2022, MNRAS, 512, 563 Romero-Gómez, M., Athanassoula, E., Masdemont, J. J., & García-Gómez, C. 2007, A&A, 472, 63 Romero-Gómez, M., Masdemont, J. J., Athanassoula, E., & García-Gómez, C. 2006, A&A, 453, 39 Romero-Gómez, M., Masdemont, J. J., García-Gómez, C., & Athanassoula, E. 2009, Comm. Nonlinear Science Nu...

    work page 2022

  25. [25]

    2018, MNRAS, 475, 676

    Springel, V ., Pakmor, R., Pillepich, A., et al. 2018, MNRAS, 475, 676

    work page 2018

  26. [26]

    & Patsis, P

    Tsigaridi, L. & Patsis, P. 2013, MNRAS, 434, 2922

    work page 2013

  27. [27]

    2008, MNRAS, 387, 1264 van der Marel, R

    Tsoutsis, P., Efthymiopoulos, C., & V oglis, N. 2008, MNRAS, 387, 1264 van der Marel, R. P. 2001, AJ, 122, 1827 van der Marel, R. P. & Kallivayalil, N. 2014, ApJ, 781, 121

    work page 2008

  28. [28]

    W., Lintott, C

    Willett, K. W., Lintott, C. J., Bamford, S. P., et al. 2013, MNRAS, 435, 2835

    work page 2013

  29. [29]

    G., Adelman, J., Anderson, Jr., J

    York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 120, 1579

    work page 2000

  30. [30]

    2013, ApJ, 772, 135 Article number, page 12 of 13 P

    Zaritsky, D., Salo, H., Laurikainen, E., et al. 2013, ApJ, 772, 135 Article number, page 12 of 13 P. Sánchez-Martín et al.: Arm morphology in off-centre barred galaxies Appendix A: Appendix A: Influence of the disc centre location To ensure that our results are robust with respect to the choice of placing the geometric centre of the disc at the coordinate...

    work page 2013