Skip to main content
Sandbox Physics
Stellar streams, encounters and observed projections

Stellar streams · L4

Read the 2025 GD-1 thin-stream and cocoon paper

Sandbox Physics ·

GD-1 contains more than a thin line. A fixed 2025 paper identifies a wider candidate population with a broader velocity distribution around the thin stream: the cocoon. The interesting questions are how that population was selected, how its inferred width changes with background assumptions, and which histories remain possible.

Start with Figure 2: which selection are you seeing?

Figure 2 connects stream longitude, both proper motions and stream latitude. It is not an unbiased all-sky census: colour, magnitude, parallax and proper-motion choices already shape the structure. Figure 7 adds uneven spectroscopic footprints, so the cocoon analysis uses a better-covered interval, roughly minus thirty to plus two degrees in longitude.

Only after recognising those selections, examine Figure 8. Residual latitude and radial velocity show a narrow concentration plus a broader distribution. The 715 Gaia projection points and 261 DESI region targets are different products; subsequent member counts depend on the particular selection.

Reading: Valluri et al. (2025), GD-1 Stellar Stream and Cocoon in the DESI Early Data Release, arXiv:2407.06336v2

Figures 10–11 and Table 3: quantify narrow and broad

Section 4 describes the thin stream and cocoon with Gaussian mixtures. A star can plausibly belong to either component; membership is probabilistic rather than fixed by two sharp sky boundaries. The velocity comparison must also include individual measurement uncertainties.

N\mathcal N denotes a bell-shaped probability density with unit area, σ\sigma its width, and QQ the narrow-component weight in the currently selected sample, between zero and one.

p(x)=Q N(x∣xˉ1,σ12)+(1−Q) N(x∣xˉ2,σ22)p(x)=Q\,\mathcal N(x\mid\bar x_1,\sigma_1^2)+(1-Q)\,\mathcal N(x\mid\bar x_2,\sigma_2^2)
The variable is a latitude offset; each component has its own mean and width. Source note: equation 5’s printed expansion has ambiguous typesetting. We define normalized Gaussian densities here instead of transcribing that expansion; this says nothing about the authors’ software or the origin of our orbit residuals.

Reading: Valluri et al. (2025), GD-1 Stellar Stream and Cocoon in the DESI Early Data Release, arXiv:2407.06336v2

Worked example: turn an angular width into a length

Table 3’s bivariate fit gives latitude dispersions of about 0.149∘0.149^\circ and 1.252∘1.252^\circ. A Gaussian full width at half maximum is about 2.355 times its dispersion: approximately 0.351 and 2.948 degrees.

Full width at half maximum does not enclose all members. Dispersion, half-maximum width, a hard selection boundary and full three-dimensional thickness are different quantities.

w≃d Δϕ,9 kpc×0.351π180≃55 pcw\simeq d\,\Delta\phi,\qquad9\,\mathrm{kpc}\times\frac{0.351\pi}{180}\simeq55\,\mathrm{pc}
Convert degrees to radians first. At the same assumed nine-kiloparsec distance, the cocoon FWHM is about 463 parsecs. These are projected widths at the assumed distance; line-of-sight thickness needs other information.

Reading: Valluri et al. (2025), GD-1 Stellar Stream and Cocoon in the DESI Early Data Release, arXiv:2407.06336v2

Change the background, change the inferred velocity dispersion

The Table-3 bivariate result gives about 3.09 kilometres per second for the thin stream and 7.81 for the cocoon. An alternative contamination treatment in section 4 gives about 2.85 and 4.84. The paper’s approximate five-to-eight cocoon range reflects that model dependence, not two supposedly identical precise measurements.

Figure 12 concerns metallicity, not another dynamical simulation. The paper notes a systematic precision scale near 0.3 dex. Similar metallicity distributions do not make the observed scatter a clean intrinsic chemical spread, and a non-rejecting distribution test does not prove a common birthplace.

Reading: Valluri et al. (2025), GD-1 Stellar Stream and Cocoon in the DESI Early Data Release, arXiv:2407.06336v2

One broad component can remember several histories

Section 5 considers a cluster stripped within a dwarf galaxy before accretion, debris of the dwarf itself, dark-substructure heating and encounters with a massive satellite such as Sagittarius. These processes can affect widths, dispersions, folds or off-track structure, and need not be mutually exclusive.

More generally, bars, spirals, giant molecular clouds and selection also deserve consideration. Rotating or changing external fields can transfer energy and angular momentum, but this paper does not complete a GD-1 comparison of every such process. For GD-1 specifically, its relatively large pericentre makes inner-disc cloud encounters less likely in the cited argument; a list of ordinary perturbers is not a list of equally plausible explanations.

Our experiment integrates only a spherical host and a prescribed Plummer passage. It includes no bar, spiral, tidal release or selection function. It supplies the response to one prescribed encounter; distinguishing histories requires more complete joint predictions.

Reading: Valluri et al. (2025), GD-1 Stellar Stream and Cocoon in the DESI Early Data Release, arXiv:2407.06336v2 · Erkal & Belokurov (2015), Forensics of Subhalo-Stream Encounters, arXiv:1412.6035v2

Why inferring the Galactic potential is also degenerate

Gap growth depends on more than the kick. The host determines how orbital frequencies change after a perturbation, so identical kicks accumulate different phase offsets in different hosts. Erkal section 6 isolates this dependence; Figure 8 supplies a harmonic-host counterexample with no secular gap growth.

In the cold circular-stream limit, the secular expansion coefficient includes (4−γ2)/γ2(4-\gamma^2)/\gamma^2, where γ=κ/Ω\gamma=\kappa/\Omega. That host factor trades against mass, encounter speed and elapsed time. Potential inference needs joint position and velocity predictions with release, perturbation and selection models. The fixed-host orbit comparison does not solve that inverse problem.

A bar or spiral introduces another distinction: the force pattern itself rotates. A rigid perturbation with fixed amplitude and constant pattern speed obeys the relation below, and can leave correlated energy and angular-momentum changes. Departures from one static orbit are therefore not automatically evidence for a localized dark subhalo.

dEdt=ΩpdLzdt\frac{dE}{dt}=\Omega_p\frac{dL_z}{dt}
Here Ωp\Omega_p is the pattern angular speed. Changing amplitude adds explicit-time work. Changing only the pattern speed preserves this instantaneous relation, but energy minus pattern speed times angular momentum gains a term from the pattern acceleration and is not conserved. No real GD-1 fit is made with it.

Reading: Erkal & Belokurov (2015), Forensics of Subhalo-Stream Encounters, arXiv:1412.6035v2 · Valluri et al. (2025), GD-1 Stellar Stream and Cocoon in the DESI Early Data Release, arXiv:2407.06336v2 · Sandbox Physics standalone stellar-stream scientific reader v1

Finish with a prediction that could distinguish histories

If a former dwarf environment produced the broad component, how should chemistry, velocity and longitude trends correlate? If a localized encounter did it, can one parameter set explain the gap, velocity changes and off-track structure together? If a feature appears only in counts, does it survive a different completeness or selection model?

These questions probe mechanisms more directly than asking whether a broad patch is visible. Erkal’s derivation already shows mass, speed, impact and elapsed-time degeneracies in a single gap profile. More dimensions add constraints without automatically removing host and formation-history dependence.

Finally inspect Appendix A. Figure 13 shifts one distance-proxy curve for display. Do not reconstruct solar reflex by reading that shifted ordinate. The runnable method uses the stated formula and fixed products, while preserving its own poor conditional real-data orbit comparison.

Reading: Valluri et al. (2025), GD-1 Stellar Stream and Cocoon in the DESI Early Data Release, arXiv:2407.06336v2 · Erkal & Belokurov (2015), Forensics of Subhalo-Stream Encounters, arXiv:1412.6035v2 · Sandbox Physics standalone stellar-stream scientific reader v1

Exercises and explanations

  1. Does a cocoon FWHM near three degrees mean every star within three degrees belongs to it?

    Read the explanation

    No. Gaussian components overlap and extend beyond their half-maximum widths. That width is a shape measure, not a membership boundary; probabilities also depend on other measurements, background and selection.

  2. Does a larger velocity dispersion uniquely establish dark-subhalo heating?

    Read the explanation

    No. A prior dwarf environment, satellite encounter, contamination and selection can affect it. Compare joint predictions for the same observations instead of deciding history from one dispersion number.

Try the prediction

Compare structures before choosing a history

  1. Compare compact and diffuse flybys without assigning a dark-matter identity to the imposed perturber.
  2. Keep finite simulated counts separate from observed mixture widths and posterior samples.
  3. Recompute the real projection exercise and see why adding one orbit does not automatically explain a cocoon.

Standalone scientific materials

Download the scientific ZIP · Read the setup and calculation scope

The ZIP freshly integrates 25 finite encounter cases and independent limits without network. Figures consume new saved records. The explicit optional real-data exercise fetches fixed inputs privately and preserves its poor conditional fit.

Sources and calculation scope

  1. Erkal & Belokurov (2015), Forensics of Subhalo-Stream Encounters, arXiv:1412.6035v2

    All 14 pages actually read. Sections 2–3 and equations 1–23 derive the short-encounter circular-stream model; sections 4–6 compare simulations and physical degeneracies. Printed pages 3/9 actually viewed. The imposed warm ensemble is not the paper’s tidal-release calculation.

  2. Valluri et al. (2025), GD-1 Stellar Stream and Cocoon in the DESI Early Data Release, arXiv:2407.06336v2

    Fixed 30 January 2025 v2, 25 pages fully read; original pages 7/12/14/21 actually viewed. Sections 3–5, Figures 2/7–12, Table 3 and Appendix A are distinguished from our conditional orbit exercise. DOI 10.3847/1538-4357/ada690; no original paper, figure or observed row redistributed.

  3. Sandbox Physics standalone stellar-stream scientific reader v1

    Own Python integrates 25 finite cases anew, checks an independent high-order reference, half steps, analytic limits and saved-record counts. Default zero network. Optional exact-source galpy exercise retains poor fit/nonconvergence/no-coverage failures. No App/Engine, completed tracks, originals or observed-derived figures in the ZIP.

FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy