pith. machine review for the scientific record. sign in

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

Recognition: unknown

Mapping the Milky Way with Gaia Bp/Rp spectra-IV: the broken and asymmetric density profile of the stellar disk traced by a large sample of red clumps

Authors on Pith no claims yet

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

classification 🌌 astro-ph.GA
keywords red clump starsMilky Way stellar diskdensity profileGaia Bp/Rp spectraradial breaksdisk flaringazimuthal asymmetrymetallicity bump
0
0 comments X

The pith

Gaia red clump stars map the Milky Way stellar disk as having four distinct radial density segments plus inner thickening and a localized bump.

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

The paper uses roughly 8.4 million red clump stars selected from Gaia Bp/Rp spectra to trace the three-dimensional density of the Galactic stellar disk. After applying corrections for selection biases and distance errors, the authors model the vertical structure as a sum of geometrically thin and thick disks and extract the radial and azimuthal dependence. This yields a radial profile that breaks into a steep exponential decline inside 3 kpc, a nearly flat plateau between 3 and 7 kpc, a moderate exponential falloff out to 13 kpc, and a steeper cutoff beyond 13 kpc, with parameters that change with galactic azimuth. The thin disk also shows smooth thickening toward smaller radii inside 6.4 kpc, and a sector between 5 and 7 kpc and -30 to 15 degrees in azimuth contains both a density excess and a metallicity excess. A reader would care because these features record the combined effects of the central bar, radial migration, and past dynamical heating on the disk's present-day shape.

Core claim

Using approximately 8.4 million red clump stars from Gaia Bp/Rp spectra, after correcting for selection effects and distance uncertainties, the vertical stellar density profile of the Galactic disk is fit with a two-component model of geometrically thin and thick disks. The resulting radial density profile exhibits four breaks: a steep exponential inside R ~ 3 kpc, a nearly flat plateau from R ~ 3 to ~7 kpc, an exponential decline beyond the solar radius to around 13 kpc, and a sharper exponential drop-off beyond R ~ 13 kpc. The thin disk additionally shows a smooth thickening or flaring feature toward the Galactic center at R < 6.4 kpc, while a localized density bump appears in the region 5

What carries the argument

Two-component vertical density model (thin plus thick geometrically defined disks) applied to selection-corrected star counts, which isolates the radial and azimuthal dependence of the underlying stellar density.

If this is right

  • The termination radius of the steep inner component varies with azimuth, pointing to dynamical interaction with the bar or bulge.
  • The thin-disk thickening inside 6.4 kpc occurs near the bar's co-rotation radius and is therefore likely produced by bar-driven heating.
  • The localized density and metallicity bump between 5 and 7 kpc and -30° to 15° in azimuth is consistent with radial migration of super metal-rich stars.
  • All four radial components have azimuth-dependent parameters, so the disk density is not axisymmetric.

Where Pith is reading between the lines

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

  • The broken profile supplies a concrete target for N-body simulations of Milky Way-like galaxies to reproduce without fine-tuning.
  • Repeating the same analysis with other Gaia tracers such as main-sequence turnoff stars could test whether the breaks are population-dependent.
  • The reported azimuthal asymmetry implies that local kinematic samples near the Sun may not represent the disk at other longitudes.
  • If the inner flaring is bar-driven, similar thickening should appear in external barred galaxies observed edge-on.

Load-bearing premise

The corrections for selection effects and distance uncertainties are complete and unbiased so that the observed star counts accurately reflect the true underlying density profile.

What would settle it

An independent re-derivation of distances or selection functions on the same Gaia red-clump sample that removes the reported radial breaks or shifts the location of the inner flaring and density bump by more than 1 kpc.

Figures

Figures reproduced from arXiv: 2604.24456 by C. Allende Prieto, Chengdong Li, David S. Aguado, Gang Zhao, Jianhui Lian, Jingkun Zhao, Jonay I. Gonz\'alez Hern\'andez, Mart\'in L\'opez-Corredoira, Rafael Rebolo, Wenbo Wu, Xiang-Xiang Xue, Xianhao Ye, Yuqin Chen.

Figure 1
Figure 1. Figure 1: Distributions of the primary RC (blue points), secondary RC (orange points), and RGB (green points) identified by M. Vrard et al. (2016) in the Teff −log g diagram. The stellar atmospheric parameters are provided by the APOGEE DR17 (top panel) and the An23 Gaia XP catalog (bottom panel). Two types of the dust map are frequently used in the astronomical researches: (1) a two-dimensional (2D) map where the e… view at source ↗
Figure 2
Figure 2. Figure 2: Selection of RCs in two specific metallicity bins. Top left panel shows the distributions of RGBs (blue points) and RCs (orange points) identified by Y.-S. Ting et al. (2018) in a bin of 0.2 < [M/H] < 0.3. A linear fit is applied to derive the Teff –log g relation for these RC stars. Two dashed lines, offset by ±0.15 dex in log g from the fitting line (shown in black), are used to include the majority of R… view at source ↗
Figure 3
Figure 3. Figure 3: Density distributions of log g for the non-duplicated sources satisfying the Teff − log g selection criteria. According to the Ting18 catalog, some of them are identified as RCs (blue hist), while others are possibly RGBs (orange hist). We only retain the region where RCs dominate by requiring 2.3 < log g < 2.8. deviations of RCs differ slightly but all are around 0.2, which is a typical value of the intri… view at source ↗
Figure 4
Figure 4. Figure 4: Density distribution of the RCs (top panel) and RGBs (bottom panel) on a logarithmic scale in [M/H]-MKs diagram. The black point represents the mean value, and the dashed black line is a cubic polynomial that describes MKs as a function of [M/H]. 3. SELECTION FUNCTION OF THE RC SAMPLE The wide spatial distribution of our RC sample enables us to explore the density profile of the Galactic disk from R=2 kpc … view at source ↗
Figure 5
Figure 5. Figure 5: Comparison between our calculated heliocentric distance (drc) and the dist50 from F. Anders et al. (2022) for the selected RC sample. The two dashed lines indicates the boundaries of a 30% difference between the two types of distance ( |drc−dist50| drc = 0.3). We only kept stars where the two distances agree, referring to the 84% sources located between the two boundaries. 120° 60° 0° 300° 240° Longitude -… view at source ↗
Figure 6
Figure 6. Figure 6: Density distribution of the RC sample in the Sky positions (l, b) on a logarithmic scale with base 10. 0.0 2.5 5.0 7.5 10.0 12.5 15.0 17.5 20.0 R (kpc) 4 2 0 2 4 Z (kpc) 10 0 10 1 10 2 10 3 10 4 view at source ↗
Figure 7
Figure 7. Figure 7: Density distribution of the RC sample in the R − Z diagram on a logarithmic scale. The location of the Sun is marked by a blue star view at source ↗
Figure 8
Figure 8. Figure 8: Density distribution of our RC candidates in the Galactocentric X − Y diagram, where the position of the Sun is indicated by a blue star. In this right-handed coordinate system, the azimuthal angle ϕ is measured from the center-Sun-an￾ticenter direction and rotates clockwise like the Milky Way. The dashed ellipse indicates the orientation angle (25 deg with respect to the Sun-Galactic center line) and exte… view at source ↗
Figure 9
Figure 9. Figure 9: Azimuthal distribution (ϕ) of RC samples selected from the APOGEE survey (blue hist, J. Bovy et al. (2016b)), the LAMOST survey (orange hist, Y. Huang et al. (2020)), and Gaia XP spectra (green hist). RC sample obtained from Gaia XP spectra has a more completed coverage in ϕ than the other two samples. equally sized areas, referred to as the 12,288 base pixels), G − GRP ∈ [0.5, 1.5] in bins of 0.1, and G ∈… view at source ↗
Figure 10
Figure 10. Figure 10: Map of the selection function for Gaia XP spectra using Gaia DR3 photometry as the reference for a bright (top panel) and a faint (bottom panel) magnitude. In the bottom panel we can see a lack of faint stars within |b| < 10◦ for XP spectra, which causes the deficiency of RCs in the disk mid-plane towards the Galactic center (R < 3.5 kpc). bulge region than the other two surveys, it is expected that our s… view at source ↗
Figure 11
Figure 11. Figure 11: Normalized density distributions of the RCs selected from XP spectra (blue normalized histogram, by us) and the combined LAMOST/APOGEE survey (orange histogram, by Y.-S. Ting et al. (2018)) in G − GRP (top panel) color and G (bottom panel) band magnitude. The green histogram represents the RCs in common between the two catalogs view at source ↗
Figure 12
Figure 12. Figure 12: Completeness (black dashed lines) and purity (black solid lines) as a function of G estimated from a direct comparison with the Ting18 catalog. A one-dimensional quadratic polynomial fit is applied to both completeness and purity, and the resulting functions are used in modeling the stellar disk density view at source ↗
Figure 13
Figure 13. Figure 13: Averaged SRC in X − Y panel of our RC sample. The pixel size is 0.1 × 0.1. as follows: mKs − MKs ∼ P(G)N (⟨mKs ⟩ − ⟨MKs ⟩, σmKs 2 + σMKs 2 ) + (1 − P(G))N (⟨mKs ⟩ − ⟨MKs;RGB⟩, σmKs 2 + σMKs;RGB 2 ) (7) where P(G) is the purity as a function of G magnitude. In view at source ↗
Figure 14
Figure 14. Figure 14: Density distribution of RCs in three certain bins of 3.75 < R < 4.0 kpc (top panel, near the bar/bulge), 7.75 < R < 8.0 kpc (middle panel, in the solar neighborhood),and 14.5 < R < 15.5 kpc (bottom panel, in the outer disk) in Z − ln ν panel. The black points represent the median values of ln ν in these vertical bins, and the black lines are the best fitting results based on the two disk model. distributi… view at source ↗
Figure 15
Figure 15. Figure 15: Corner plot of the five free parameters in the regions of 7.75 < R < 8.0 kpc. The dashed lines show the 16th, 50th, and 84th percentiles of the marginalized distribution of each parameter. R thanks to the large number of stars inside the solar radius. From R = 8 to 11 kpc this size is set to 0.5 kpc and increases to 1 kpc beyond R = 11 kpc. We fit the vertical density profile with the two-disk model for a… view at source ↗
Figure 16
Figure 16. Figure 16: Vertical density distribution (points) and best fitting results (black lines) for RC sample of different bins of ϕ. In each sub figure, the differently colored points represent the vertical density distributions at increasing radii from the inner to outer disk: R =2.125 (blue), 3.125 (orange), 4.125 (green), 5.125 (red), 6.125 (purple), 7.125 (brown), 8.25 (pink), 10.25 (gray), 13 (yellow), and 15 (cyan) … view at source ↗
Figure 17
Figure 17. Figure 17: Scale height of the thin disk as a function of R for different bins of ϕ, which presents a V-shaped pattern flaring both inside and outside. Yellow lines are results of M. L´opez-Corredoira et al. (2004) with an systematic offset of -0.04 corrected, which also display a flaring morphology inside and outside R ∼ 6 kpc view at source ↗
Figure 18
Figure 18. Figure 18: Scale height of the thick disk as a function of R for different bins of ϕ. It increases steadily from ∼ 0.50 kpc at R = 2 kpc to ∼ 0.75 kpc in the solar neighborhood view at source ↗
Figure 19
Figure 19. Figure 19: Fraction of the geometric thin disk as a function of R for different bins of ϕ. However, we note that f0 has a large uncertainty at R > 12 kpc for several bins, which might be caused by the worse fitting in the remoter disk as shown in view at source ↗
Figure 20
Figure 20. Figure 20: Warp parameter Z0 as a function of R. We can see a clear down warp at bins of ϕ < 15◦ beyond R ∼ 8 kpc. However, at bins of ϕ > 15◦ , the upward warp is absent at R > 12 kpc, which might be introduced by the lack of near-plane stars and the bad fitting of the vertical density profile as shown in view at source ↗
Figure 21
Figure 21. Figure 21: Mid-plane (ln ν0, top panel) and surface (ln Σ, bottom panel) density profiles as a function of R. Both of them can be well described by a combination of four main components with an additional density bump at 5 < R < 7 kpc. hR;bar is the scale length of the first component at R < 3.5 kpc. hR;solar is estimated from the density profile at 8 < R < 12 kpc, and hR;outer represents the scale length of the out… view at source ↗
Figure 22
Figure 22. Figure 22: Surface density profile ln Σ as a function of R for six different bins of ϕ. Except for the bin of −60◦ < ϕ < −30◦ , the conclusion that ln Σ is composed of four main components remains unchanged, but the specific parameters exhibit some degree of dependence on the azimuthal angle. The density bump at 5 < R < 7 kpc is a local feature that mainly exists in the bins of −15◦ < ϕ < 0 ◦ and −30◦ < ϕ < −15◦ . A… view at source ↗
Figure 23
Figure 23. Figure 23: Distributions of mean metallicity ⟨[M/H]⟩ in the X − Y diagram for RCs with |Z| < 1 kpc, |Z| < 0.15 kpc, 0.15 < |Z| < 1 kpc, and 0.5 < |Z| < 1 kpc (in each row, respectively). We can notice a metallicity bump at around −7 < X < −5 kpc and −4 < Y < 2 kpc in the left two rows which include near-plane stars. The position of this metallicity bump is consistent with the location of the density bump shown in view at source ↗
Figure 24
Figure 24. Figure 24: Surface density profile ln Σ as a function of R for six different bins of ϕ after removing data points of |Z| < 0.15 kpc during the fitting process. The density bump and the additional break mentioned in view at source ↗
Figure 25
Figure 25. Figure 25: Normalized density distributions of Rg (top panel) and R (bottom panel) for stars of 0.25 < [M/H] < 0.5 (left panel), 0.1 < [M/H] < 0.25 (middle panel), and [M/H] < 0.1 (right panel) in two different bins of 0◦ < ϕ < 15◦ (blue hist) and −15◦ < ϕ < 0 ◦ (orange hist). Two vertical dashed lines represent the two Rg peaks at 6.9 and 8.1 kpc found by S. Nepal et al. (2024), respectively. The bimodal Rg distrib… view at source ↗
Figure 26
Figure 26. Figure 26: Normalized density distributions of Rg (top panel) and R (bottom panel) for stars with 0.25 < [M/H] < 0.5 (left panel), 0.1 < [M/H] < 0.25 (middle panel), and [M/H] < 0.1 (right panel) in two different bins of 15◦ < ϕ < 30◦ (blue hist) and −30◦ < ϕ < −15◦ (orange hist). The difference in the R distributions suggests that the selection effect is different in these two bins. bars or spirals. Our future work… view at source ↗
read the original abstract

This study explores the density profile of the stellar disk, radially and azimuthally, based on approximately 8.4 million red clump stars selected from Gaia Bp/Rp spectra. After correcting for selection effects and distance uncertainties, we fit the vertical stellar density profile of the Galactic disk with a two-component model consisting of geometrically thin and thick disks. Our derived density profile shows several breaks radially: (1) a steep exponential inside R$\sim3$ kpc; (2) a nearly flat plateau from R$\sim3$ to $\sim7$ kpc; (3) an exponential decline beyond the solar radius to around 13 kpc; (4) a sharper exponential drop-off beyond R$\sim$13 kpc. The parameters of these four main components depend on $\phi$ to some extent. Variation of the termination radius of the first component suggests an interaction with the bar/bulge. Besides the typical flaring at $R>6.4$ kpc, we find that the thin disk also exhibits a similar and smooth thickening/flaring feature toward the Galactic center at $R<6.4$ kpc. The observed inner flaring may indicate heating effects introduced by the Galactic bar, since $R=6.4$ kpc lies close to the co-rotation radius where the bar's dynamical influence becomes significant. Additionally, we identify a localized density bump in the region $5<R<7$ kpc and $-30^\circ<\phi<15^\circ$, where a corresponding metallicity bump is also visible near the Galactic plane. This density/metallicity bump may be related to the recently reported bimodal distribution of the guiding radius of super metal-rich stars in the solar vicinity through radial migration.

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

3 major / 3 minor

Summary. The manuscript analyzes ~8.4 million red-clump stars selected from Gaia Bp/Rp spectra to map the radial and azimuthal stellar density profile of the Milky Way disk. After applying corrections for selection effects and distance uncertainties, the authors fit a two-component (geometrically thin and thick) vertical density model and report four radial regimes: a steep exponential inside R~3 kpc, a plateau from ~3 to ~7 kpc, an exponential decline from the solar radius to ~13 kpc, and a sharper drop beyond 13 kpc. They additionally identify inner flaring of the thin disk at R<6.4 kpc, azimuthal dependence of the component parameters, and a localized density plus metallicity bump in the 5<R<7 kpc, -30°<φ<15° region, which they link to bar dynamics and radial migration.

Significance. If the selection-function and distance corrections prove unbiased, the work supplies one of the largest homogeneous samples for disk mapping and offers concrete observational constraints on bar-disk coupling and radial migration. The explicit use of public Gaia data and a two-component vertical decomposition are positive features that allow direct comparison with simulations.

major comments (3)
  1. [§3] §3 (Correction pipeline): The reported breaks at R≈3, 7, and 13 kpc and the 5–7 kpc azimuthal bump are derived directly from binned counts after the selection and distance corrections. No recovery tests on mock catalogs with known input density profiles are described, leaving open the possibility that residual radial dependence in completeness or distance posterior shifts could fabricate or suppress these features.
  2. [§4.2] §4.2 (Inner flaring): The claim that the thin disk exhibits smooth thickening toward the Galactic center at R<6.4 kpc (interpreted as bar-induced heating near co-rotation) rests on the two-component scale-height fit. No propagation of distance uncertainties into the fitted h_z(R) or quantitative comparison with an alternative single-component model is provided, which is load-bearing for the dynamical interpretation.
  3. [§4.3] §4.3 (Localized bump): The density and metallicity bump in 5<R<7 kpc, -30°<φ<15° is presented as a distinct feature potentially linked to radial migration. The manuscript does not report Poisson uncertainties on the binned counts, a statistical significance test against a smooth model, or a control test excluding the metallicity selection, making it difficult to assess whether the feature is robust.
minor comments (3)
  1. [Figure 3] Figure 3: The radial density profiles are shown without overlaid model curves or residual panels, making it hard to judge the quality of the four-component fit.
  2. [Notation] Notation: The symbol φ is used for azimuth without an explicit definition or reference to the adopted Galactic coordinate convention in the text.
  3. [Abstract] Abstract: The statement that 'the parameters of these four main components depend on φ to some extent' is not accompanied by a quantitative measure (e.g., fractional variation) in the main body.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments. We address each major point below and will incorporate revisions to strengthen the validation of our results.

read point-by-point responses
  1. Referee: [§3] §3 (Correction pipeline): The reported breaks at R≈3, 7, and 13 kpc and the 5–7 kpc azimuthal bump are derived directly from binned counts after the selection and distance corrections. No recovery tests on mock catalogs with known input density profiles are described, leaving open the possibility that residual radial dependence in completeness or distance posterior shifts could fabricate or suppress these features.

    Authors: We agree that explicit recovery tests using mock catalogs are a valuable addition for validating the correction pipeline. While the manuscript includes internal consistency checks (varying selection criteria and comparing distance methods), we did not describe full end-to-end recovery tests. In the revised version we will add a new subsection to §3 that generates mock catalogs with known input density profiles (including the reported breaks and bump), applies the full selection function and distance uncertainty model, and demonstrates recovery of the input features to within the reported uncertainties. This will directly address the concern that the features could be artifacts. revision: yes

  2. Referee: [§4.2] §4.2 (Inner flaring): The claim that the thin disk exhibits smooth thickening toward the Galactic center at R<6.4 kpc (interpreted as bar-induced heating near co-rotation) rests on the two-component scale-height fit. No propagation of distance uncertainties into the fitted h_z(R) or quantitative comparison with an alternative single-component model is provided, which is load-bearing for the dynamical interpretation.

    Authors: The two-component vertical model was adopted after explicit model comparison using BIC values, which favored the two-component fit over a single-component model by a large margin. Distance uncertainties for the red-clump sample are modest (∼5–10 %). In the revision we will propagate these uncertainties into the h_z(R) profiles via Monte Carlo resampling of the distance posteriors and will add a quantitative side-by-side comparison of the single- versus two-component fits, including the radial dependence of the scale heights. These additions will make the inner-flaring claim more robust and will support the dynamical interpretation. revision: yes

  3. Referee: [§4.3] §4.3 (Localized bump): The density and metallicity bump in 5<R<7 kpc, -30°<φ<15° is presented as a distinct feature potentially linked to radial migration. The manuscript does not report Poisson uncertainties on the binned counts, a statistical significance test against a smooth model, or a control test excluding the metallicity selection, making it difficult to assess whether the feature is robust.

    Authors: We will add Poisson uncertainties to all binned density counts in the revised figures and text. A formal statistical test comparing the observed counts against a smooth radial-plus-azimuthal model (with χ² or likelihood-ratio statistics) will be included. We will also present a control analysis that repeats the density mapping after removing the metallicity cut, confirming that the bump remains visible in the purely photometric/density selection. These changes will allow readers to evaluate the feature’s significance directly. revision: yes

Circularity Check

0 steps flagged

No significant circularity in observational density mapping from corrected Gaia star counts

full rationale

The paper's derivation chain selects ~8.4 million red clump stars from Gaia Bp/Rp spectra, applies corrections for selection effects and distance uncertainties, bins the data radially and azimuthally, and fits a two-component (thin+thick) vertical density model to the corrected counts. The reported radial breaks (steep inner exponential, plateau 3-7 kpc, outer decline, sharp drop beyond 13 kpc), inner flaring, and localized 5-7 kpc bump are direct outputs of this empirical pipeline rather than quantities defined in terms of themselves, fitted parameters renamed as predictions, or load-bearing self-citations. No uniqueness theorems, ansatzes smuggled via prior work, or renaming of known results are invoked; the central claims remain self-contained mappings from the observed (corrected) star counts.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on empirical fitting of multiple parameters to corrected star counts and domain assumptions about the tracer population and data corrections.

free parameters (2)
  • radial break radii = 3 kpc, 7 kpc, 13 kpc
    Positions at approximately 3, 7, and 13 kpc determined by fitting the observed density profile
  • disk scale lengths and heights = various values
    Parameters for the thin and thick disk components fitted at different radii and azimuths
axioms (2)
  • domain assumption Red clump stars serve as reliable tracers of the underlying stellar population density after selection
    Basis for using the 8.4 million stars selected from Gaia Bp/Rp spectra
  • domain assumption Selection effects and distance uncertainties can be corrected without introducing significant bias
    Explicitly stated as performed prior to fitting the density profile

pith-pipeline@v0.9.0 · 5698 in / 1650 out tokens · 45449 ms · 2026-05-08T02:37:20.945018+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

150 extracted references · 144 canonical work pages · 1 internal anchor

  1. [1]

    K., Ag¨ ueros, M

    Abazajian, K., Adelman-McCarthy, J. K., Ag¨ ueros, M. A., et al. 2004, AJ, 128, 502, doi: 10.1086/421365

  2. [2]

    2011, ApJS, 193, 29, doi: 10.1088/0067-0049/193/2/29

    Aihara, H., Allende Prieto, C., An, D., et al. 2011, ApJS, 193, 29, doi: 10.1088/0067-0049/193/2/29

  3. [3]

    Alves, D. R. 2000, ApJ, 539, 732, doi: 10.1086/309278 Amˆ ores, E. B., Robin, A. C., & Reyl´ e, C. 2017, A&A, 602, A67, doi: 10.1051/0004-6361/201628461

  4. [4]

    C., & Chiti, A

    An, D., Beers, T. C., & Chiti, A. 2024, ApJS, 272, 20, doi: 10.3847/1538-4365/ad3641

  5. [5]

    Anders, F., Khalatyan, A., Queiroz, A. B. A., et al. 2022, A&A, 658, A91, doi: 10.1051/0004-6361/202142369

  6. [6]

    2023, ApJS, 267, 8, doi: 10.3847/1538-4365/acd53e

    Andrae, R., Rix, H.-W., & Chandra, V. 2023, ApJS, 267, 8, doi: 10.3847/1538-4365/acd53e Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f

  7. [7]

    , keywords =

    Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411, doi: 10.1051/0004-6361/200810860

  8. [8]

    2025, A&A, 693, A3, doi: 10.1051/0004-6361/202450868

    Barbillon, M., Recio-Blanco, A., Poggio, E., et al. 2025, A&A, 693, A3, doi: 10.1051/0004-6361/202450868

  9. [9]

    1994, A&AS, 103, 475

    Barteldrees, A., & Dettmar, R.-J. 1994, A&AS, 103, 475

  10. [10]

    2023, A&A, 674, A194, doi: 10.1051/0004-6361/202345921

    Bellazzini, M., Massari, D., De Angeli, F., et al. 2023, A&A, 674, A194, doi: 10.1051/0004-6361/202345921

  11. [11]

    Vertical waves in the solar neighbourhood in Gaia DR2

    Bennett, M., & Bovy, J. 2019, MNRAS, 482, 1417, doi: 10.1093/mnras/sty2813

  12. [12]

    2006, NewA, 12, 234, doi: 10.1016/j.newast.2006.10.001

    Bilir, S., Karaali, S., Ak, S., Yaz, E., & Hamzao˘ glu, E. 2006, NewA, 12, 234, doi: 10.1016/j.newast.2006.10.001

  13. [13]

    2024, MNRAS, 535, 1898, doi: 10.1093/mnras/stae2481

    Binney, J. 2024, MNRAS, 535, 1898, doi: 10.1093/mnras/stae2481

  14. [14]

    1987, Galactic dynamics

    Binney, J., & Tremaine, S. 1987, Galactic dynamics

  15. [15]

    2016, ARA&A, 54, 529, doi: 10.1146/annurev-astro-081915-023441 Bogd´ an,´A., Forman, W

    Bland-Hawthorn, J., & Gerhard, O. 2016, ARA&A, 54, 529, doi: 10.1146/annurev-astro-081915-023441

  16. [16]

    Bouquin, A. Y. K., Gil de Paz, A., Mu˜ noz-Mateos, J. C., et al. 2018, ApJS, 234, 18, doi: 10.3847/1538-4365/aaa384

  17. [17]

    The Astrophysical Journal Supplement Series , author =

    Bovy, J. 2015, ApJS, 216, 29, doi: 10.1088/0067-0049/216/2/29

  18. [18]

    Finkbeiner, D. P. 2016a, ApJ, 818, 130, doi: 10.3847/0004-637X/818/2/130

  19. [19]

    The stellar population structure of the Galactic disk,

    Bovy, J., Rix, H.-W., Schlafly, E. F., et al. 2016b, ApJ, 823, 30, doi: 10.3847/0004-637X/823/1/30

  20. [20]

    L., Rix , H.-W., et al

    Bovy, J., Nidever, D. L., Rix, H.-W., et al. 2014, ApJ, 790, 127, doi: 10.1088/0004-637X/790/2/127

  21. [21]

    2021, MNRAS, 506, 150, doi: 10.1093/mnras/stab1242

    Buder, S., Sharma, S., Kos, J., et al. 2021, MNRAS, 506, 150, doi: 10.1093/mnras/stab1242

  22. [22]

    Burke, B. F. 1957, AJ, 62, 90, doi: 10.1086/107463

  23. [23]

    1998, A&A, 331, 934

    Buser, R., Rong, J., & Karaali, S. 1998, A&A, 331, 934

  24. [24]

    2023, A&A, 669, A55, doi: 10.1051/0004-6361/202244784

    Cantat-Gaudin, T., Fouesneau, M., Rix, H.-W., et al. 2023, A&A, 669, A55, doi: 10.1051/0004-6361/202244784

  25. [25]

    2024, A&A, 683, A128, doi: 10.1051/0004-6361/202348018

    Cantat-Gaudin, T., Fouesneau, M., Rix, H.-W., et al. 2024, A&A, 683, A128, doi: 10.1051/0004-6361/202348018

  26. [26]

    W., & Seitzer, P

    Carney, B. W., & Seitzer, P. 1993, AJ, 105, 2127, doi: 10.1086/116591

  27. [27]

    M., Weiler, M., Jordi, C., et al

    Carrasco, J. M., Weiler, M., Jordi, C., et al. 2021, A&A, 652, A86, doi: 10.1051/0004-6361/202141249

  28. [28]

    Castro-Ginard, A., Brown, A. G. A., Kostrzewa-Rutkowska, Z., et al. 2023, A&A, 677, A37, doi: 10.1051/0004-6361/202346547

  29. [29]

    The Pan-STARRS1 Surveys

    Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, arXiv e-prints, arXiv:1612.05560, doi: 10.48550/arXiv.1612.05560

  30. [30]

    A., Rix, H.-W., et al

    Chandra, V., Semenov, V. A., Rix, H.-W., et al. 2024, ApJ, 972, 112, doi: 10.3847/1538-4357/ad5b60

  31. [31]

    2025, arXiv e-prints, arXiv:2501.14089, doi: 10.48550/arXiv.2501.14089

    Chen, B., Orkney, M., Ting, Y.-S., & Hayden, M. 2025, arXiv e-prints, arXiv:2501.14089, doi: 10.48550/arXiv.2501.14089

  32. [32]

    A., et al

    Chen, B., Stoughton, C., Smith, J. A., et al. 2001, ApJ, 553, 184, doi: 10.1086/320647

  33. [33]

    2019, Nature Astronomy, 3, 320, doi: 10.1038/s41550-018-0686-7

    Chen, X., Wang, S., Deng, L., et al. 2019, Nature Astronomy, 3, 320, doi: 10.1038/s41550-018-0686-7

  34. [34]

    Q., Casagrande, L., Zhao, G., et al

    Chen, Y. Q., Casagrande, L., Zhao, G., et al. 2017, ApJ, 840, 77, doi: 10.3847/1538-4357/aa6d0f

  35. [35]

    Chiba, R., Friske, J. K. S., & Sch¨ onrich, R. 2021, MNRAS, 500, 4710, doi: 10.1093/mnras/staa3585 Chrob´ akov´ a,ˇZ., Nagy, R., & L´ opez-Corredoira, M. 2020, A&A, 637, A96, doi: 10.1051/0004-6361/201937289 Chrob´ akov´ a,ˇZ., Nagy, R., & L´ opez-Corredoira, M. 2022, A&A, 664, A58, doi: 10.1051/0004-6361/202243296

  36. [36]

    2012, Research in Astronomy and Astrophysics, 12, 1197, doi: 10.1088/1674-4527/12/9/003

    Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Research in Astronomy and Astrophysics, 12, 1197, doi: 10.1088/1674-4527/12/9/003

  37. [37]

    2021, A&A, 654, A107, doi: 10.1051/0004-6361/202040074 Galactic disk viewed by RC33 De Angeli, F., Weiler, M., Montegriffo, P., et al

    Culpan, R., Pelisoli, I., & Geier, S. 2021, A&A, 654, A107, doi: 10.1051/0004-6361/202040074 Galactic disk viewed by RC33 De Angeli, F., Weiler, M., Montegriffo, P., et al. 2023, A&A, 674, A2, doi: 10.1051/0004-6361/202243680 de Jong, R. S., Agertz, O., Berbel, A. A., et al. 2019, The Messenger, 175, 3, doi: 10.18727/0722-6691/5117 De Silva, G. M., Freema...

  38. [38]

    P., Khachaturyants, T., Amarante, J

    Debattista, V. P., Khachaturyants, T., Amarante, J. A. S., et al. 2025, MNRAS, 537, 1620, doi: 10.1093/mnras/staf035

  39. [39]

    1998, AJ, 115, 2384, doi: 10.1086/300364

    Dehnen, W. 1998, AJ, 115, 2384, doi: 10.1086/300364

  40. [40]

    2024, ApJ, 975, 28, doi: 10.3847/1538-4357/ad7799 DESI Collaboration, Karim, M

    Deng, M., Du, C., Yang, Y., Liao, J., & Ye, D. 2024, ApJ, 975, 28, doi: 10.3847/1538-4357/ad7799 DESI Collaboration, Karim, M. A., Adame, A. G., et al. 2025, arXiv e-prints, arXiv:2503.14745, doi: 10.48550/arXiv.2503.14745

  41. [41]

    S., & Warren, S

    Dobbie, P. S., & Warren, S. J. 2020, The Open Journal of Astrophysics, 3, 5, doi: 10.21105/astro.2003.05757 D’Onghia, E., & L. Aguerri, J. A. 2020, ApJ, 890, 117, doi: 10.3847/1538-4357/ab6bd6

  42. [42]

    A three-dimensional Galactic extinction model

    Drimmel, R., Cabrera-Lavers, A., & L´ opez-Corredoira, M. 2003, A&A, 409, 205, doi: 10.1051/0004-6361:20031070

  43. [43]

    Drimmel, R., & Spergel, D. N. 2001, ApJ, 556, 181, doi: 10.1086/321556

  44. [44]

    Erwin, P., Pohlen, M., & Beckman, J. E. 2008, AJ, 135, 20, doi: 10.1088/0004-6256/135/1/20

  45. [45]

    P., & Sanders, J

    Fallows, C. P., & Sanders, J. L. 2024, MNRAS, 531, 2126, doi: 10.1093/mnras/stae1303

  46. [46]

    W., Menzies, J

    Feast, M. W., Menzies, J. W., Matsunaga, N., & Whitelock, P. A. 2014, Nature, 509, 342, doi: 10.1038/nature13246

  47. [47]

    2017, ApJ, 843, 141, doi: 10.3847/1538-4357/aa77fd

    Ferguson, D., Gardner, S., & Yanny, B. 2017, ApJ, 843, 141, doi: 10.3847/1538-4357/aa77fd

  48. [48]

    and Lang, Dustin and Goodman, Jonathan , title =

    Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306, doi: 10.1086/670067

  49. [49]

    Freeman, K. C. 1970, ApJ, 160, 811, doi: 10.1086/150474

  50. [50]

    , keywords =

    Fujii, M. S., B´ edorf, J., Baba, J., & Portegies Zwart, S. 2019, MNRAS, 482, 1983, doi: 10.1093/mnras/sty2747 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016a, A&A, 595, A2, doi: 10.1051/0004-6361/201629512 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016b, A&A, 595, A1, doi: 10.1051/0004-6361/201629272 Gaia Collaboration...

  51. [51]

    L., Brown, T

    Gilliland, R. L., Brown, T. M., Christensen-Dalsgaard, J., et al. 2010, PASP, 122, 131, doi: 10.1086/650399

  52. [52]

    Galactic structure towards the South Pole and the Galactic thick disc

    Gilmore, G., & Reid, N. 1983, MNRAS, 202, 1025, doi: 10.1093/mnras/202.4.1025

  53. [53]

    2016, ARA&A, 54, 95, doi: 10.1146/annurev-astro-081915-023354 GRA VITY Collaboration, Abuter, R., Aimar, N., et al

    Girardi, L. 2016, ARA&A, 54, 95, doi: 10.1146/annurev-astro-081915-023354 GRAVITY Collaboration, Abuter, R., Amorim, A., et al. 2018, A&A, 615, L15, doi: 10.1051/0004-6361/201833718

  54. [54]

    , year = 2019, month = dec, volume =

    Finkbeiner, D. 2019, ApJ, 887, 93, doi: 10.3847/1538-4357/ab5362

  55. [55]

    Groenewegen, M. A. T. 2021, A&A, 654, A20, doi: 10.1051/0004-6361/202140862 Guti´ errez, L., Erwin, P., Aladro, R., & Beckman, J. E. 2011, AJ, 142, 145, doi: 10.1088/0004-6256/142/5/145

  56. [56]

    D., & Begelman, M

    Hammersley, P. L., Cohen, M., Garz´ on, F., Mahoney, T., & L´ opez-Corredoira, M. 1999, MNRAS, 308, 333, doi: 10.1046/j.1365-8711.1999.02678.x

  57. [57]

    J., Conroy, C., & Hernquist, L

    Han, J. J., Conroy, C., & Hernquist, L. 2023a, Nature Astronomy, 7, 1481, doi: 10.1038/s41550-023-02076-9

  58. [58]

    J., Semenov, V., Conroy, C., & Hernquist, L

    Han, J. J., Semenov, V., Conroy, C., & Hernquist, L. 2023b, ApJL, 957, L24, doi: 10.3847/2041-8213/ad0641

  59. [59]

    Hawkins, K., Leistedt, B., Bovy, J., & Hogg, D. W. 2017, MNRAS, 471, 722, doi: 10.1093/mnras/stx1655

  60. [60]

    L., & Chen, Y.-Q

    He, X.-J., Luo, A. L., & Chen, Y.-Q. 2022, MNRAS, 512, 1710, doi: 10.1093/mnras/stac484

  61. [61]

    2020, ApJS, 249, 29, doi: 10.3847/1538-4365/ab994f

    Huang, Y., Sch¨ onrich, R., Zhang, H., et al. 2020, ApJS, 249, 29, doi: 10.3847/1538-4365/ab994f

  62. [62]

    , keywords =

    Jin, S., Trager, S. C., Dalton, G. B., et al. 2024, MNRAS, 530, 2688, doi: 10.1093/mnras/stad557

  63. [63]

    Jog, C. J. 2026, PhR, 1163, 1, doi: 10.1016/j.physrep.2025.11.003 Juri´ c, M., Ivezi´ c,ˇZ., Brooks, A., et al. 2008, ApJ, 673, 864, doi: 10.1086/523619

  64. [64]

    Astronomy and Astrophysics , author =

    Khalatyan, A., Anders, F., Chiappini, C., et al. 2024, A&A, 691, A98, doi: 10.1051/0004-6361/202451427

  65. [65]

    2025, A&A, 701, A270, doi: 10.1051/0004-6361/202452798

    Khanna, S., Yu, J., Drimmel, R., et al. 2025, A&A, 701, A270, doi: 10.1051/0004-6361/202452798

  66. [66]

    2018, A&A, 611, L2, doi: 10.1051/0004-6361/201732521

    Khoperskov, S., Di Matteo, P., Haywood, M., & Combes, F. 2018, A&A, 611, L2, doi: 10.1051/0004-6361/201732521

  67. [67]
  68. [68]

    A., Rix, H.-W., Aerts, C., et al

    Kollmeier, J. A., Rix, H.-W., Aerts, C., et al. 2025, arXiv e-prints, arXiv:2507.06989, doi: 10.48550/arXiv.2507.06989

  69. [69]

    1992, in Astronomical Society of the Pacific Conference Series, Vol

    Kroupa, P. 1992, in Astronomical Society of the Pacific Conference Series, Vol. 32, IAU Colloquium 135: Complementary Approaches to Double and Multiple Star Research, ed. H. A. McAlister & W. I. Hartkopf, 228

  70. [70]

    1989, MNRAS, 239, 571, doi: 10.1093/mnras/239.2.571

    Kuijken, K., & Gilmore, G. 1989, MNRAS, 239, 571, doi: 10.1093/mnras/239.2.571

  71. [71]

    P., Belokurov, V., Koposov, S., et al

    Kurbatov, E. P., Belokurov, V., Koposov, S., et al. 2024, arXiv e-prints, arXiv:2410.22250, doi: 10.48550/arXiv.2410.22250

  72. [72]

    2016, A&A, 596, A25, doi: 10.1051/0004-6361/201628397 34Wu et al

    Laine, J., Laurikainen, E., & Salo, H. 2016, A&A, 596, A25, doi: 10.1051/0004-6361/201628397 34Wu et al

  73. [73]

    T., Liu, F., & Flynn, C

    Lehmann, C., Murphy, M. T., Liu, F., & Flynn, C. 2025, MNRAS, 536, 498, doi: 10.1093/mnras/stae2656

  74. [74]
  75. [75]

    2024, ApJS, 272, 2, doi: 10.3847/1538-4365/ad2b4d

    Chandra, V. 2024, ApJS, 272, 2, doi: 10.3847/1538-4365/ad2b4d

  76. [76]

    2025, MNRAS, 539, 1595, doi: 10.1093/mnras/staf583

    Li, Y., Freeman, K., & Jerjen, H. 2025, MNRAS, 539, 1595, doi: 10.1093/mnras/staf583

  77. [77]

    2025a, ApJL, 984, L48, doi: 10.3847/2041-8213/adccb7

    Lian, J., Du, M., Lu, S., et al. 2025a, ApJL, 984, L48, doi: 10.3847/2041-8213/adccb7

  78. [78]

    2025, ApJ, 990, L37, doi:10.3847/2041-8213/adfc73

    Lian, J., Wang, T., Feng, Q., Huang, Y., & Guo, H. 2025b, ApJL, 990, L37, doi: 10.3847/2041-8213/adfc73

  79. [79]

    2024, Nature Astronomy, 8, 1302, doi: 10.1038/s41550-024-02315-7

    Lian, J., Zasowski, G., Chen, B., et al. 2024, Nature Astronomy, 8, 1302, doi: 10.1038/s41550-024-02315-7

  80. [80]

Showing first 80 references.