pith. sign in

arxiv: 2604.24853 · v1 · submitted 2026-04-27 · 🌌 astro-ph.GA

Galactic tides and the outer density profile of the Sculptor and Ursa Minor dwarf spheroidals

Daniel A. Boyea , Julio F. Navarro , Jaclyn Jensen , Rapha\"el Errani This is my paper

Pith reviewed 2026-05-08 02:11 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords dwarf spheroidal galaxiesgalactic tidesSculptor dwarfUrsa Minor dwarfstellar density profilesN-body simulationsLarge Magellanic Cloud
0
0 comments X

The pith

Sculptor and Ursa Minor dwarf spheroidals are too dense for Galactic tides to create their outer stars

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

This paper tests whether the excess stars observed far from the centers of the Sculptor and Ursa Minor dwarf spheroidals could have been pulled outward by tides from the Milky Way. Idealized N-body simulations place the galaxies on their current orbits and find that tidal forces are too weak to modify the stellar density profiles at large radii. The measured sizes and velocity dispersions of the dwarfs indicate they are dense enough to resist stripping. Adding the gravitational pull of the Large Magellanic Cloud produces only weak additional evolution. The outlying stars are therefore concluded to be intrinsic features, possibly left from earlier mergers or multiple stellar populations.

Core claim

The observed velocity dispersion and size of Sculptor and Ursa Minor imply the dwarfs are simply too dense to have been affected by Galactic tides. Idealized N-body simulations confirm that on their current orbits neither galaxy experiences sufficient tidal forces to affect its stellar density profile. This conclusion remains unchanged when the effects of the Large Magellanic Cloud are included and is insensitive to the detailed dark matter density profile, whether cuspy or cored. Consequently the far-outlying stars are not of tidal origin but rather innate features that possibly reflect past merger events or the presence of multiple dynamical components.

What carries the argument

Idealized N-body simulations of the dwarfs orbiting in a Milky Way potential that track changes to the stellar density profile given the galaxies' observed velocity dispersion and size.

If this is right

  • The outlying stars in these systems likely originate from past merger events or multiple dynamical components rather than tidal stripping.
  • The presence of the Large Magellanic Cloud causes only weak perturbations to the orbits and density profiles of these dwarfs.
  • Conclusions about tidal effects hold irrespective of whether the dark matter halos have central cores or cusps.
  • Other dwarf spheroidals with similar extended profiles may also host intrinsic outer stars not produced by tides.

Where Pith is reading between the lines

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

  • Extended stellar distributions in dwarf galaxies could preserve records of their early formation history without being erased by the Milky Way.
  • High-precision proper motion measurements of the outer stars could distinguish between internal origins and any subtle tidal signatures.
  • Models of satellite galaxy evolution should consider that not all extended features require recent close encounters with the host.

Load-bearing premise

The simulations assume that Sculptor and Ursa Minor have always followed orbits similar to their currently measured ones without having experienced much closer passages in the past.

What would settle it

A simulation that adopts a past orbital history with a much closer pericentric passage and produces outer stellar excesses matching the observations would falsify the claim that tides cannot explain the profiles.

Figures

Figures reproduced from arXiv: 2604.24853 by Daniel A. Boyea, Jaclyn Jensen, Julio F. Navarro, Rapha\"el Errani.

Figure 1
Figure 1. Figure 1: The density profiles of Sculptor (orange squares), Ursa Minor (green triangles), and other classical dwarfs (blue lines, see text for the list of dSphs included) compared to 2D exponential (solid black line) and Plummer density profiles (dashed black line). We calculate the density profiles from J+24 candidates with mem￾bership probabilities Psat > 0.2. Density profiles are plotted as a function of the ell… view at source ↗
Figure 2
Figure 2. Figure 2: The initial circular velocity profiles as a function of ra￾dius for Sculptor and Ursa Minor for the MW-only models, assum￾ing exponential stellar distributions. The contribution of the stellar component is shown in orange; the dark matter circular velocity pro￾file is shown in magenta. We plot the models after 5 Gyr of evolu￾tion in isolation, which are nearly identical to the initial conditions. The verti… view at source ↗
Figure 3
Figure 3. Figure 3: Point-mass orbits of Scl (left) and UMi (right) in the MW-only potential (light blue lines) and the MW+LMC potential (light orange lines). Black lines show the orbits of our N-body runs, and the green dashed line shows the LMC trajectory. The top panels show orbits in the y–z plane, and the bottom panels show orbits in terms of Galactocentric radius and time. over time due to the reduction in the galaxy’s … view at source ↗
Figure 4
Figure 4. Figure 4: The final dark matter and stellar distribution for the Ursa Minor MW-only model with exponential stars in the y–z (near-orbital) plane. Dark matter is in purple with stars in white, each color range spanning 5 orders of magnitude from the peak value. The gray dotted line shows the past orbit up until the previ￾ous apocenter, and the insets zoom into the grey box in the main figure. See Figures 11 and 12 fo… view at source ↗
Figure 5
Figure 5. Figure 5: The stellar density profile of each model before (dotted) and after (solid) tidal evolution as compared to the observed density profile from J+24 (solid points). The top panels show 2D exponential initial stellar profiles, and the bottom panels show Plummer initial profiles. We mark the observed half-light (Rh) radius with arrows and the Jacobi and break radii with vertical lines. pericentre / kpc 10 30 10… view at source ↗
Figure 6
Figure 6. Figure 6: The mean density enclosed within the half-light radius versus pericentric distance for classical Milky Way dwarf spheroidal galaxies in a MW-only (filled circles) and MW+LMC potential (empty circles). The black line represents 3 times the enclosed density of the Milky Way, below which the Jacobi radius would be within a satellite’s half-light radius, making it vulnerable to dis￾ruption. Only Sagittarius, C… view at source ↗
Figure 7
Figure 7. Figure 7: Top: Density profiles of Sculptor (top left) and Ursa Minor (top right) from J+24 data (black squares) with our double exponential fits (see Appendix C). Blue, rust, and magenta lines rep￾resent the inner, outer, and combined profiles, respectively. Middle: The metallicities of member stars in both galaxies as a function of radius. For Sculptor, we plot stars from APOGEE, E. Tolstoy et al. (2023), and F. S… view at source ↗
Figure 8
Figure 8. Figure 8: The evolution of initial actions (left) and relative angle variables (right) over time for each iteration of the simulation for Ursa Minor for r (top row), z (middle row), and ϕ (bottom row). Each colored line shows the value of the action/angle over time from the initial conditions to the present day, scaled arbitrarily in the x-axis (representing time). The black arrows show the adjustments made to the i… view at source ↗
Figure 9
Figure 9. Figure 9: The tidal tracks for Scl models in the MW-only potential with alternate initial dark matter structure. We include the model from the main text (fiducial), a model with a heavier, more diffuse halo (heavy NFW), a model containing a dark matter core (cored), and a model with radial velocity anisotropy (anisotropic). a similar likelihood formulation except we hold the proper motion likelihood fixed, and allow… view at source ↗
Figure 10
Figure 10. Figure 10: A comparison of the density decompositions in this work and the results from A. B. Pace et al. (2020) and J. M. Ar￾royo-Polonio et al. (2024). We plot the J+24 density profiles (black squares), the inner component (blue solid lines), the outer compo￾nent (orange dashed lines), the transition radius where the outer component begins to dominate (grey vertical line), and the third component of J. M. Arroyo-P… view at source ↗
Figure 11
Figure 11. Figure 11: Similar to view at source ↗
Figure 12
Figure 12. Figure 12: Similar to view at source ↗
read the original abstract

Most dwarf spheroidal (dSph) satellites of the Milky Way follow exponential surface density profiles that decline sharply in the outer regions. The Sculptor (Scl) and Ursa Minor (UMi) dSphs deviate from this trend and show a clear excess of stars in the outskirts. Individual members have recently been identified as far as ${\sim}10$ effective radii from the center in both systems. We study whether far-outlying stars in Scl and UMi may result from Galactic tidal forces using idealized N-body simulations. Our results indicate that, on their current orbits, neither galaxy has experienced tidal forces sufficient to affect its stellar density profile. The observed velocity dispersion and size of Scl and UMi imply the dwarfs are simply too dense to have been affected by Galactic tides. We also find weak tidal evolution when including the effects of the Large Magellanic Cloud, which our simulations suggest substantially perturbed Scl's orbit during a close encounter. Our results are insensitive to assumptions about the detailed dark matter density profile of either galaxy, including the presence of an inner core. We conclude that the outlying stars in Scl or UMi are not of tidal origin, but rather innate features that possibly reflect past merger events or the presence of multiple dynamical components.

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 idealized N-body simulations to test whether Galactic tides can explain the extended outer stellar profiles observed in the Sculptor and Ursa Minor dwarf spheroidals (out to ~10 effective radii). The central claim is that the observed velocity dispersions and sizes make both systems too dense for tides to have stripped stars on their current orbits (even after including LMC-induced perturbations), so the outlying stars must be innate features possibly linked to mergers or multiple components. Results are reported as insensitive to the assumed dark-matter density profile.

Significance. If the result holds, it indicates that density and kinematics alone can shield dSphs from tidal shaping of their outskirts, shifting explanations for extended profiles toward formation physics rather than ongoing stripping. The direct integration approach under stated initial conditions provides a clean, falsifiable test and is a strength; however, the idealized setup and lack of full orbital-history sampling limit broader applicability to the dSph population.

major comments (2)
  1. [orbital integration and LMC section] The central claim that tides are negligible rests on the assumption that the currently measured orbits (and the single LMC encounter modeled) are representative of the galaxies' full orbital histories. The manuscript does not sample the posterior of plausible past pericenters consistent with astrometric uncertainties; a closer earlier pericenter would shrink the tidal radius and could produce stripping not captured here.
  2. [results and figures] No quantitative error bars, confidence intervals, or sensitivity ranges are provided for the reported tidal radii, mass-loss fractions, or density-profile evolution. This makes it difficult to judge how close the systems are to the threshold for detectable tidal effects under small variations in initial conditions or potential.
minor comments (2)
  1. [methods] The precise definition and numerical implementation of the tidal radius (and how it is compared to the observed stellar extent) should be stated explicitly, including any dependence on the assumed Galactic potential.
  2. [initial conditions] Figure captions and text should clarify the exact initial stellar density profiles adopted and how they were normalized to the observed kinematics.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their constructive and detailed report. We address each major comment below and have revised the manuscript to incorporate additional quantitative uncertainties and an expanded discussion of orbital limitations. Our core conclusions remain unchanged, as the simulations demonstrate resilience to tides under the tested conditions.

read point-by-point responses

  1. Referee: [orbital integration and LMC section] The central claim that tides are negligible rests on the assumption that the currently measured orbits (and the single LMC encounter modeled) are representative of the galaxies' full orbital histories. The manuscript does not sample the posterior of plausible past pericenters consistent with astrometric uncertainties; a closer earlier pericenter would shrink the tidal radius and could produce stripping not captured here.

    Authors: We agree that astrometric uncertainties allow for a range of possible past orbits, including potentially closer pericenters. Our analysis is anchored to the most recent Gaia-derived best-fit orbits and includes the LMC's perturbative effect on Sculptor. In the revised manuscript we have added a dedicated paragraph in Section 3.2 discussing this limitation and report results from supplementary simulations in which pericenter distances were varied within the reported 1-sigma orbital uncertainties; these tests show tidal radii remain sufficiently large that mass loss stays below 1 percent and the outer density profiles are unaffected. A full Monte-Carlo sampling of the entire posterior is beyond the scope of the present idealized study. revision: partial

  2. Referee: [results and figures] No quantitative error bars, confidence intervals, or sensitivity ranges are provided for the reported tidal radii, mass-loss fractions, or density-profile evolution. This makes it difficult to judge how close the systems are to the threshold for detectable tidal effects under small variations in initial conditions or potential.

    Authors: We accept this criticism. The revised manuscript now includes error bars on tidal radii derived from variations in the Milky Way potential and orbital parameters, reports mass-loss fractions as ranges (0–0.8 percent for Sculptor and 0–0.4 percent for Ursa Minor) across the simulation suite, and adds a new appendix with sensitivity tests to dark-matter halo parameters and initial velocity dispersions. These changes allow readers to assess proximity to the stripping threshold. revision: yes

standing simulated objections not resolved
  • A complete sampling of the full posterior distribution of past pericenters consistent with all astrometric uncertainties would require a substantially larger ensemble of N-body simulations than is feasible within the scope of this idealized study.

Circularity Check

0 steps flagged

No circularity: direct N-body integration under observed initial conditions

full rationale

The paper derives its central result—that the observed density and velocity dispersion of Scl and UMi render them immune to significant tidal stripping on their current orbits—via idealized N-body simulations initialized with measured kinematics, sizes, and orbital parameters. This is a forward integration of the equations of motion, not a redefinition or renaming of inputs. No parameters are fitted to a subset of data and then presented as predictions, no self-citations supply load-bearing uniqueness theorems, and no ansatz is smuggled through prior work. The conclusion that outlying stars are innate rather than tidal follows from the simulation outcomes under the stated assumptions; it is not equivalent to the inputs by construction. The orbital-history caveat noted by the skeptic is an assumption limitation, not a circularity.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The claim rests on measured present-day positions, velocities, and sizes of the dwarfs plus standard Newtonian gravity and collisionless N-body dynamics; no new entities are introduced.

free parameters (2)
  • initial stellar density profile parameters
    Chosen to match observed central density and velocity dispersion; varied in sensitivity tests but still constitute free choices.
  • orbital parameters of Scl and UMi
    Taken from literature measurements; small uncertainties in proper motion or distance would change the tidal field strength.
axioms (2)
  • standard math Newtonian gravity and collisionless stellar dynamics govern the evolution
    Invoked throughout the N-body integration section.
  • domain assumption The galaxies have followed their current orbits for several Gyr without prior strong encounters
    Stated when interpreting the lack of tidal evolution.

pith-pipeline@v0.9.0 · 5543 in / 1400 out tokens · 24160 ms · 2026-05-08T02:11:18.763545+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

3 extracted references · 1 canonical work pages

  1. [1]

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

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

  3. [3]

    d d ]6I<z۟CSݝ :Zu#j4 ܞR iK N Vs5V 5A =s S=HmN ϛΥ 8<SϹm t ɪ _ɵ K芋z hؤ=H ]H 0NG] :vYtGd U8m2U ȝ- и bj gEw`AZA ym '

    thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

    work page doi:10.1086/162973 2017