Pith. sign in

REVIEW 2 major objections 1 minor 4 references

The Milky Way bar's pattern speed can be measured to roughly 1 km/s/kpc precision from the shapes of tidal tails around a few well-chosen nearby open clusters.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.3

2026-06-28 21:57 UTC pith:OXTBEZF3

load-bearing objection Simulations show bar pattern speed shapes cluster tails more than spirals or satellites, but the 1 km/s/kpc precision claim rests on untested separation from stronger GMC effects. the 2 major comments →

arxiv 2605.31439 v1 pith:OXTBEZF3 submitted 2026-05-29 astro-ph.GA

Dynamics of tidal tails of open clusters: I. effects of bar, spiral arms and giant molecular clouds

Janez Kos , Jovana Risojevi\'c , Samo Ilc This is my paper
classification astro-ph.GA
keywords tidal tailsopen clustersGalactic barpattern speedMilky Way potentialn-body simulationsstellar streamsgiant molecular clouds
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Open clusters dissolve and leave behind elongated tidal tails whose shapes trace the gravitational forces along their orbits. Simulations show that the inner Galaxy's bar, especially how fast it rotates, produces clear changes in tail orientation and structure, while spiral arms cause smaller shifts and distant satellites leave nearby tails untouched. Giant molecular clouds create stronger random distortions than the bar for older, disk-plane clusters. Statistical comparisons of tail shapes across many potential models identify which nearby clusters carry the cleanest signal of the bar's rotation rate. If the metrics isolate that signal in real data, observations of only a handful of clusters would fix the pattern speed independently of the bar's length or tilt.

Core claim

N-body simulations of synthetic and real open clusters demonstrate that the Galactic bar's pattern speed exerts the dominant influence on the shapes and orientations of their tidal tails; statistical metrics applied to the stellar distributions allow recovery of this pattern speed to a precision of order 1 km s^{-1} kpc^{-1} from a small set of sensitive nearby clusters, with limited confounding from spiral arms and negligible effects from satellites, though giant molecular clouds introduce stronger perturbations for old in-plane clusters.

What carries the argument

Non-parametric statistical metrics that quantify differences between tail morphologies across large grids of gravitational potential models.

Load-bearing premise

The statistical metrics will separate the bar pattern speed signal from the larger distortions produced by giant molecular clouds when the same metrics are applied to actual observations.

What would settle it

Observed tail morphologies for the selected clusters show no systematic match to the bar pattern speed values that the simulations predict, or the scatter from unmodeled GMCs exceeds the differences produced by changing the pattern speed.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The bar pattern speed dominates tail morphology for nearby clusters over spiral arms or satellites.
  • A small number of targeted observations can constrain the pattern speed without depending on bar length or orientation.
  • GMC perturbations are stronger than bar effects for old clusters lying in the plane.
  • Spiral arms produce only limited changes in tail shape compared with the bar.

Where Pith is reading between the lines

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

  • The same tail-shape metrics could be applied to other dissolving structures to test consistency of the bar speed measurement.
  • Clusters at different ages or heights above the plane could be used to separate bar-driven signals from GMC noise in follow-up work.
  • If the method works, repeated observations over years might reveal whether the bar pattern speed itself changes slowly.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 1 minor

Summary. The manuscript uses n-body simulations of open clusters in a Milky Way potential including bar, spiral arms, GMCs and satellites. Non-parametric statistical metrics are applied to quantify differences in tidal tail morphologies across varied potential parameters. The central claim is that the bar (especially its pattern speed) strongly affects tail shapes, while GMCs produce stronger distortions for old/in-plane clusters; observations of tidal tails from a handful of well-selected nearby clusters should allow measurement of the bar pattern speed to ~1 km s^{-1} kpc^{-1} precision, independent of bar length and orientation. Observability is deferred to Paper II.

Significance. If the metrics can isolate the bar pattern-speed signal, the work would supply an independent dynamical probe of the inner Milky Way using dissolving clusters. The forward n-body approach with explicit parameter variation and non-parametric metrics is a clear strength, enabling broad exploration of the potential without assuming specific functional forms. This could complement gas-dynamical or orbit-based bar constraints.

major comments (2)
  1. [Abstract, final paragraph] Abstract, final paragraph: The headline claim of ~1 km s^{-1} kpc^{-1} precision (independent of bar length/orientation) is load-bearing. The same paragraph states that GMC perturbations produce stronger tail distortions than the bar for old and in-plane clusters. No test is shown that the non-parametric metrics recover the bar signal when realistic GMC amplitudes and observational errors are present.
  2. [Abstract] Abstract: The quantitative precision of order 1 km s^{-1} kpc^{-1} is asserted without reference to the specific simulation ensemble, the exact statistical metrics, the error budget, or the procedure used to derive the figure. This absence prevents assessment of whether the claimed sensitivity is robust.
minor comments (1)
  1. A dedicated methods subsection defining the non-parametric statistical metrics (including any distance measures or summary statistics) would improve reproducibility.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive and detailed report. Below we respond point-by-point to the major comments, all of which concern the abstract. We agree that the abstract's claims require clarification and will revise accordingly.

read point-by-point responses

  1. Referee: [Abstract, final paragraph] Abstract, final paragraph: The headline claim of ~1 km s^{-1} kpc^{-1} precision (independent of bar length/orientation) is load-bearing. The same paragraph states that GMC perturbations produce stronger tail distortions than the bar for old and in-plane clusters. No test is shown that the non-parametric metrics recover the bar signal when realistic GMC amplitudes and observational errors are present.

    Authors: We agree that the precision claim is prominent and that the abstract does not demonstrate recovery of the bar signal under combined GMC and error perturbations. The quoted precision is an estimate based on the spread of the non-parametric metrics across our bar-pattern-speed ensemble (Sections 4–5). Because observability, error budgets and combined perturbations are deferred to Paper II, no such joint test appears here. We will revise the abstract to qualify the precision statement and note the conditions under which it is expected. revision: yes

  2. Referee: [Abstract] Abstract: The quantitative precision of order 1 km s^{-1} kpc^{-1} is asserted without reference to the specific simulation ensemble, the exact statistical metrics, the error budget, or the procedure used to derive the figure. This absence prevents assessment of whether the claimed sensitivity is robust.

    Authors: The abstract is a concise summary; the ensemble is defined in Section 2, the metrics in Section 3, and the procedure yielding the ~1 km s^{-1} kpc^{-1} figure is shown in Section 5. We will revise the abstract to state explicitly that the quoted precision follows from the simulation results presented in the paper. revision: partial

standing simulated objections not resolved
  • Explicit demonstration that the non-parametric metrics recover the bar pattern-speed signal in the joint presence of realistic GMC amplitudes and observational errors (deferred to Paper II).

Circularity Check

0 steps flagged

No circularity; forward simulations yield independent sensitivity claims

full rationale

The paper conducts forward n-body simulations that vary gravitational potential parameters (bar pattern speed, spiral arms, GMCs) and quantifies resulting tail morphologies with non-parametric statistical metrics. No step reduces a claimed prediction to a fitted input, self-definition, or self-citation chain; the 1 km s^{-1} kpc^{-1} precision statement is an extrapolation from simulation outcomes to future observations (Paper II), not a quantity forced by the paper's own equations. The derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The paper relies on standard n-body integration methods and multi-component Galactic potential models drawn from prior literature; the abstract does not list explicit free parameters or new entities but implies variation of bar pattern speed, spiral arm strength, and GMC encounter rates.

free parameters (2)
  • bar pattern speed
    Explicitly varied across simulation runs to quantify its effect on tail shape and orientation.
  • GMC perturbation parameters
    Included as variable perturbers whose strength affects older in-plane clusters.
axioms (1)
  • domain assumption Orbits of open clusters in the Solar neighbourhood are sensitive to the gravitational potential of the inner Galaxy, which is dominated by the Galactic bar.
    Stated directly in the abstract as the motivation for focusing on bar effects.

pith-pipeline@v0.9.1-grok · 5873 in / 1275 out tokens · 32144 ms · 2026-06-28T21:57:00.122077+00:00 · methodology

0 comments
read the original abstract

Open clusters gradually dissolve, and their stars disperse into the Galactic field. Lost stars form tidal tails-elongated streams that trace the cluster orbit ahead of and behind its core. From the shape and orientation of the tidal tails, it is possible to infer the shape of the gravitational potential governing the cluster's motion. The orbits of open clusters, including those in the Solar neighbourhood, are sensitive to the gravitational potential of the inner Galaxy, which is dominated by the Galactic bar. Using n-body simulations of synthetic and real open clusters, we investigate how sensitive the shapes and orientations of tidal tails are to variations of the gravitational potential of the Milky Way. We consider the effects of the bar as well as spiral arms, giant molecular clouds (GMCs) and satellite galaxies. We analyse the stellar distributions within tidal tails using statistical metrics that quantify the differences between tail morphologies. Such non-parametric approach enables us to efficiently explore tidal tails across a large parameter space of gravitational potential models. We find that the Galactic bar-particularly its pattern speed-has a strong influence on the orbits of open clusters and the shape of their tails. Spiral arms have a limited effect, and satellite galaxies do not disturb the tidal tails of nearby open clusters. Perturbations by GMCs affect most clusters, with distortions stronger than those by the bar observed in old and in-plane clusters. We identify nearby open clusters that are most sensitive to the pattern speed of the bar. By observing the tidal tails of a handful of well-selected nearby clusters, we should be able to measure the pattern speed of the bar with a precision in the order of $1\ \mathrm{km\,s^{-1}\,kpc^{-1}}$ independently from length and orientation of the bar. We will present the observability of tidal tails in paper II.

Figures

Figures reproduced from arXiv: 2605.31439 by Janez Kos, Jovana Risojevi\'c, Samo Ilc.

Figure 1
Figure 1. Figure 1: Rotational curve of the common gravitational potential (blue line). Rotational curves of individual components that constitute the common potential are also shown. Rotational curves from Sofue (2020) and Ou et al. (2024) are shown in black. Poggio 2018). The rotational velocity of our model and its com￾ponents is shown in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Potentials used in this work. Left: the complete potential including spiral arms and GMCs. Top panels show the potential in the Z = 0 plane, and in the Y = 0 plane. Below, we plot the density corresponding to the above potential. Following from left to right: potential and density of disks (thin and thick), the bar, the spiral arms, and the GMCs. Spiral arms are implemented as a perturbation of the common … view at source ↗
Figure 3
Figure 3. Figure 3: Coordinate system used in this work. Z axis (not shown) points out of the drawing. PA marks the position angle of the bar. In this work, we use a Cartesian galactocentric coordi￾nate system in which the Sun is at a position (X, Y, Z) = (−8122.0, 0.0, 20.8) pc (GRAVITY Collaboration et al. 2018; Bennett & Bovy 2019). The local standard of rest (LSR) at X = −8122.0 pc has a velocity (U, V, W) = (0.0, 233.4, … view at source ↗
Figure 4
Figure 4. Figure 4: Left: Similarities between all pairs of simulations 1.5 Gyr long with different pattern speed of spiral arms (Ωs) and cluster’s current vR (panels left to right) and vT (panels top to bottom). In each panel, we plot the KLD and MMD distances between pairs of simulations with different Ωs . Above the diagonals, we plot the MMD distance; below, the KLD distance. We note that the numerical values for MMD and … view at source ↗
Figure 5
Figure 5. Figure 5: Left: similarities between all pairs of simulations with a varied bar pattern speed Ωb and cluster’s current vR (panels left to right) and vT (panels top to bottom). In each panel, we plot the KLD and MMD distances between pairs of simulations with different Ωb. Above the diagonals, we plot the MMD distance; below, the KLD distance. We note that the numerical values for MMD and KLD distances are not on the… view at source ↗
Figure 6
Figure 6. Figure 6: Left: similarities between the tidal tails developed in a simulation with vR = 0.0 km s−1 , vT = −30.0 km s−1 , Ωb = 39.0 km s−1 kpc−1 , and PA = 0.41 and sim￾ulations with the same vR and vT but with 0.21 ≤ PA ≤ 0.77 and 29.0 km s−1 kpc−1 ≤ Ωb ≤ 39.0 km s−1 kpc−1 . KLD distances and MMD distances are plotted in the top and bottom panels, respectively. Black lines mark slopes of 0.114 rad/km s−1 kpc−1 = 11… view at source ↗
Figure 8
Figure 8. Figure 8: Similarities between the tidal tails developed in a simulation with vR = 0.0 km s−1 , vT = −30.0 km s−1 , Ωb = 39.0 km s−1 kpc−1 , and bar’s semi-major axis ab = 4.5 kpc and simulations with 2.0 kpc ≤ ab ≤ 6.5 kpc and 28.0 km s−1 kpc−1 ≤ Ωb ≤ 51.0 km s−1 kpc−1 . KLD distances and MMD distances are plotted in the top and bottom panels, respectively. apocenter. Although the differences are detectable with st… view at source ↗
Figure 7
Figure 7. Figure 7: Orbits of clusters in four simulations marked a, b, i, and f in [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: Left: Stellar densities along the principal axis of 25 simulated clusters. Plots for individual clusters are shifted vertically for clarity and share the same scale as shown for the bottom-most plot. D is the position along the principal axis of the cluster. Right: six clusters marked a to f in the left plot are shown in the XY plane. They all share the same scale marked on the bottom-right. Images of clus… view at source ↗
Figure 10
Figure 10. Figure 10: Similarities between all pairs of simulations with no satellite galaxies, with LMC, SMC, and SDG. 25 panels show a combination of present day velocities vR (left to right) and vT (top to bottom). Above the diagonals in each panel, we plot the MMD distance; below the diag￾onals, the KLD distance. We note that the numerical values for MMD and KLD distances are not on the same scale, hence two colour bars. 5… view at source ↗
Figure 11
Figure 11. Figure 11: Simulations of tidal tails for selected real open clusters. In each panel, we plot a cluster’s tidal tails in the Galactic plane, for 7 different values of bar pattern speed Ωb. Stars in each panel are plotted in random order. Simulation results for all simulated clusters can be found in Appendix E. 28.0 32.0 36.0 40.0 44.0 48.0 Alessi_191 ASCC_123 HSC_1074 HSC_1724 NGC_752 28.0 32.0 36.0 40.0 44.0 48.0 2… view at source ↗
Figure 12
Figure 12. Figure 12: Similarities between all pairs of simulations with varied bar pattern speed Ωb for clusters shown in [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Simulation results for clusters whose tidal tails have been studied in the literature. In each panel, we plot a cluster’s tidal tails in the Galactic plane, for 7 different values of bar pattern speed Ωb. Stars in each panel are plotted in random order. Simulation results for all simulated clusters can be found in Appendix E (see Figures E.1–E.5). Way. We built the gravitational potential with a combinati… view at source ↗
Figure 14
Figure 14. Figure 14: Sensitivity of tidal tails to the variations of the bar pattern speed when stars far from the cluster are discarded. The plot shows the case for a synthetic cluster in a CR with the bar at vR = 0 km s−1 , vT = −30 km s−1 , and Ωb = 34.0 km s−1 kpc−1 . To compute the similari￾ties between this cluster and a cluster where Ωb is 1, 3, or 5 km s−1 kpc−1 different, we only used stars up to the maximum distance… view at source ↗
Figure 15
Figure 15. Figure 15: Distribution of CR and OLR for open clusters in the Solar neighbourhood, assuming a bar pattern speed given on the horizontal axis. before the current global bar pattern speed is actually observed via the shape of the tidal tails, the observability of any other complication is rather speculative. We show that at least half of the real open clusters studied in this work are inherently unaffected by the var… view at source ↗

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

4 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    Training generative neural networks via Maximum Mean Discrepancy optimization

    Antoja, T., Helmi, A., Dehnen, W., et al. 2014, A&A, 563, A60 Bekki, K. & Stanimirovi´c, S. 2009, MNRAS, 395, 342 Bennett, M. & Bovy, J. 2019, MNRAS, 482, 1417 Bhattacharya, S., Rao, K. K., Agarwal, M., Balan, S., & Vaidya, K. 2022, MN- RAS, 517, 3525 Binney, J. & Tremaine, S. 2008, Galactic Dynamics: Second Edition Bland-Hawthorn, J. & Gerhard, O. 2016, ...

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  2. [2]

    In this work, we set the scale radius of the halo toR=12.0 kpc

    with large variations in the literature. In this work, we set the scale radius of the halo toR=12.0 kpc. Virial radius isR vir =320 kpc. To model the nucleus, we used a Plummer potential withR=300 pc and set its mass such that the rotational curve in the Galactic nucleus matches the observations. Here, we provide Python code that reproduces the common gra...

    2012

  3. [3]

    Here we test the performance of the following kernels

    Appendix B.1.3: Kernels The core of the MMD algorithm is a kernel mean embedding, so for MMD we also compare different kernels. Here we test the performance of the following kernels. Gaussian kernelhas a form K(d)=exp − d2 2σ2 ! , (B.7) whereσcontrols the bandwidth. Polynomial kernelhas a form K(x,y)= x⊤y σ2 +c !d , (B.8) wherepandqare vectors consisting ...

    1998

  4. [4]

    We also marked angular velocity, pericentre, apocentre and eccentricity of each orbit in the grid

    Each panel shows an orbit integrated back in time for 500 Myr from the initial conditions matching the position of the Sun and with the velocityv z =0, and vR andv T marked on the edge of the grid. We also marked angular velocity, pericentre, apocentre and eccentricity of each orbit in the grid. Orbits in perturbed potentials can differ from those shown i...

    2024