We investigate the inner regions of the Milky Way using data from APOGEE and Gaia EDR3. Our inner Galactic sample has more than 26 500 stars within | X Gal |< 5 kpc, | Y Gal |< 3.5 kpc, | Z Gal |< 1 kpc, and we also carry out the analysis for a foreground-cleaned subsample of 8000 stars that is more representative of the bulge–bar populations. These samples allow us to build chemo-dynamical maps of the stellar populations with vastly improved detail. The inner Galaxy shows an apparent chemical bimodality in key abundance ratios [ α /Fe], [C/N], and [Mn/O], which probe different enrichment timescales, suggesting a star formation gap (quenching) between the high- and low- α populations. Using a joint analysis of the distributions of kinematics, metallicities, mean orbital radius, and chemical abundances, we can characterize the different populations coexisting in the innermost regions of the Galaxy for the first time. The chemo-kinematic data dissected on an eccentricity–| Z | max plane reveal the chemical and kinematic signatures of the bar, the thin inner disc, and an inner thick disc, and a broad metallicity population with large velocity dispersion indicative of a pressure-supported component. The interplay between these different populations is mapped onto the different metallicity distributions seen in the eccentricity–| Z | max diagram consistently with the mean orbital radius and V ϕ distributions. A clear metallicity gradient as a function of | Z | max is also found, which is consistent with the spatial overlapping of different populations. Additionally, we find and chemically and kinematically characterize a group of counter-rotating stars that could be the result of a gas-rich merger event or just the result of clumpy star formation during the earliest phases of the early disc that migrated into the bulge. Finally, based on 6D information, we assign stars a probability value of being on a bar orbit and find that most of the stars with large bar orbit probabilities come from the innermost 3 kpc, with a broad dispersion of metallicity. Even stars with a high probability of belonging to the bar show chemical bimodality in the [ α /Fe] versus [Fe/H] diagram. This suggests bar trapping to be an efficient mechanism, explaining why stars on bar orbits do not show a significant, distinct chemical abundance ratio signature.
more »
« less
Dynamically constraining the length of the Milky way bar
We present a novel method for constraining the length of the Galactic bar using 6D phase-space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named RFreq, based on the maximal extent of trapped bar orbits. We find the RFreq measured from orbits is consistent with the RFreq of the assumed potential only when the length of the bar and pattern speed of said potential is similar to the model from which the initial phase-space coordinates of the orbits are derived. Therefore, one can measure the model’s or the Milky Way’s bar length from 6D phase-space coordinates by determining which assumed potential leads to a self-consistent measured RFreq. When we apply this method to ≈210 000 stars in APOGEE DR17 and Gaia eDR3 data, we find a consistent result only for potential models with a dynamical bar length of ≈3.5 kpc. We find the Milky Way’s trapped bar orbits extend out to only ≈3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential.
more » « less- Award ID(s):
- 1907417
- NSF-PAR ID:
- 10398416
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 520
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 4779-4792
- Size(s):
- ["p. 4779-4792"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract The Galactic bulge is critical to our understanding of the Milky Way. However, due to the lack of reliable stellar distances, the structure and kinematics of the bulge/bar beyond the Galactic center have remained largely unexplored. Here, we present a method to measure distances of luminous red giants using a period–amplitude–luminosity relation anchored to the Large Magellanic Cloud, with random uncertainties of 10%–15% and systematic errors below 1%–2%. We apply this method to data from the Optical Gravitational Lensing Experiment to measure distances to 190,302 stars in the Galactic bulge and beyond out to 20 kpc. Using this sample, we measure a distance to the Galactic center of
R 0= 8108 ± 106stat± 93syspc, consistent with direct measurements of stars orbiting Sgr A*. We cross-match our distance catalog with Gaia DR3 and use the subset of 39,566 overlapping stars to provide the first constraints on the Milky Way’s velocity field (V R ,V ϕ ,V z ) beyond the Galactic center. We show that theV R quadrupole from the bar’s near side is reflected with respect to the Galactic center, indicating that the bar is bisymmetric and aligned with the inner disk. We also find that the vertical heightV Z map has no major structure in the region of the Galactic bulge, which is inconsistent with a current episode of bar buckling. Finally, we demonstrate withN -body simulations that distance uncertainty plays a factor in the alignment of the major and kinematic axes of the bar, necessitating caution when interpreting results for distant stars. -
null (Ed.)ABSTRACT Until the recent advent of Gaia Data Release 2 (DR2) and deep multi-object spectroscopy, it has been difficult to obtain 6D phase space information for large numbers of stars beyond 4 kpc, in particular towards the Galactic Centre, where dust and crowding are significant. We combine line-of-sight velocities from the Abundances and Radial velocity Galactic Origins Survey (ARGOS) with proper motions from Gaia DR2 to obtain a sample of ∼7000 red clump stars with 3D velocities. We perform a large-scale stellar kinematics study of the Milky Way bulge to characterize the bulge velocity ellipsoids in 20 fields. The tilt of the major-axis of the velocity ellipsoid in the radial-longitudinal velocity plane, or vertex deviation, is characteristic of non-axisymmetric systems and a significant tilt is a robust indicator of non-axisymmetry or bar presence. We compare the observations to the predicted kinematics of an N-body boxy-bulge model formed from dynamical instabilities. In the model, the lv values are strongly correlated with the angle (α) between the bulge major-axis and the Sun-Galactic centre line of sight. We use a maximum likelihood method to obtain an independent measurement of α, from bulge stellar kinematics alone, performing a robust error analysis. The most likely value of α given our model is α = (29 ± 3)○, with an additional systematic uncertainty due to comparison with one specific model. In Baade’s window, the metal-rich stars display a larger vertex deviation (lv = −40○) than the metal-poor stars (lv = 10○) but we do not detect significant lv−metallicity trends in the other fields.more » « less
-
null (Ed.)ABSTRACT Stream stars removed by tides from their progenitor satellite galaxy or globular cluster act as a group of test particles on neighbouring orbits, probing the gravitational field of the Milky Way. While constraints from individual streams have been shown to be susceptible to biases, combining several streams from orbits with various distances reduces these biases. We fit a common gravitational potential to multiple stellar streams simultaneously by maximizing the clustering of the stream stars in action space. We apply this technique to members of the GD-1, Palomar 5 (Pal 5), Orphan, and Helmi streams, exploiting both the individual and combined data sets. We describe the Galactic potential with a Stäckel model, and vary up to five parameters simultaneously. We find that we can only constrain the enclosed mass, and that the strongest constraints come from the GD-1, Pal 5, and Orphan streams whose combined data set yields $M(\lt 20\, \mathrm{kpc}) = 2.96^{+0.25}_{-0.26} \times 10^{11} \, \mathrm{ M}_{\odot}$. When including the Helmi stream in the data set, the mass uncertainty increases to $M(\lt 20\, \mathrm{kpc}) = 3.12^{+3.21}_{-0.46} \times 10^{11} \, \mathrm{M}_{\odot}$.more » « less
-
more » « less
ABSTRACT In this paper, we construct the circular velocity curve of the Milky Way out to ∼30 kpc, providing an updated model of the dark matter density profile. We derive precise parallaxes for 120 309 stars with a data-driven model, using APOGEE DR17 spectra combined with GaiaDR3, 2MASS, and WISE photometry. At outer galactic radii up to 30 kpc, we find a significantly faster decline in the circular velocity curve compared to the inner parts. This decline is better fit with a cored Einasto profile with a slope parameter $0.91^{+0.04}_{-0.05}$ than a generalized Navarro–Frenk–White (NFW) profile. The virial mass of the best-fitting dark matter halo profile is only $1.81^{+0.06}_{-0.05}\times 10^{11}$ M⊙, significantly lower than what a generalized NFW profile delivers. We present a study of the potential systematics, affecting mainly large radii. Such a low mass for the Galaxy is driven by the functional forms tested, given that it probes beyond our measurements. It is found to be in tension with mass measurements from globular clusters, dwarf satellites, and streams. Our best-fitting profile also lowers the expected dark matter annihilation signal flux from the galactic centre by more than an order of magnitude, compared to an NFW profile-fit. In future work, we will explore profiles with more flexible functional forms to more fully leverage the circular velocity curve and observationally constrain the properties of the Milky Way’s dark matter halo.
