skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Orbital phase-driven biases in galactic mass constraints from stellar streams
ABSTRACT One of the most promising tracers of the Galactic potential in the halo region is stellar streams. However, individual stream fits can be limited by systematic biases. To study these individual stream systematics, we fit streams in Milky Way-like galaxies from Feedback In Realistic Environments cosmological galaxy formation simulations with an analytical gravitational potential by maximizing the clustering of stream stars in action space. We show that for coherent streams the quality of the constraints depends on the orbital phase of the observed stream stars, despite the fact that the phase information is discarded in action-clustering methods. Streams on intermediate phases give the most accurate results, whereas pericentre streams can be highly biased. This behaviour is tied to the amount of correlation present between positions and momenta in each stream’s data: weak correlation in pericentre streams prohibits efficient differentiation between potentials, while strong correlation in intermediate streams promotes it. Although simultaneous fitting of multiple streams is generally prescribed as the remedy to combat individual stream biases, we find that combining multiple pericentric streams is not enough to yield a bias-free result. We finally show that adopting the two-component Stäckel model does not fundamentally induce a biased mass estimate. With our full data set of two multiwrap streams, we recovered the true rotation curve of the simulated galaxy within $$12{{\ \rm per\ cent}}$$ over the entire range of radii covered by our set of stars (10–176 kpc) and within $$6.5{{\ \rm per\ cent}}$$ between the 5 and 95 percentile distance range (23–109 kpc).  more » « less
Award ID(s):
2007232
PAR ID:
10327059
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
509
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
5365 to 5381
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    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
  2. ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Understanding the evolution of satellite galaxies of the Milky Way (MW) and M31 requires modelling their orbital histories across cosmic time. Many works that model satellite orbits incorrectly assume or approximate that the host halo gravitational potential is fixed in time and is spherically symmetric or axisymmetric. We rigorously benchmark the accuracy of such models against the FIRE-2 cosmological baryonic simulations of MW/M31-mass haloes. When a typical surviving satellite fell in ($$3.4\!-\!9.7\, \rm {Gyr}$$ ago), the host halo mass and radius were typically 26–86 per cent of their values today, respectively. Most of this mass growth of the host occurred at small distances, $$r\lesssim 50\, \rm {kpc}$$, opposite to dark matter only simulations, which experience almost no growth at small radii. We fit a near-exact axisymmetric gravitational potential to each host at z = 0 and backward integrate the orbits of satellites in this static potential, comparing against the true orbit histories in the simulations. Orbital energy and angular momentum are not well conserved throughout an orbital history, varying by 25 per cent from their current values already $$1.6\!-\!4.7\, \rm {Gyr}$$ ago. Most orbital properties are minimally biased, ≲10 per cent, when averaged across the satellite population as a whole. However, for a single satellite, the uncertainties are large: recent orbital properties, like the most recent pericentre distance, typically are ≈20 per cent uncertain, while earlier events, like the minimum pericentre or the infall time, are ≈40–80 per cent uncertain. Furthermore, these biases and uncertainties are lower limits, given that we use near-exact host mass profiles at z = 0. 
    more » « less
  3. ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT The combination of galaxy–galaxy lensing (GGL) and galaxy clustering is a powerful probe of low-redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we use an N-body and halo occupation distribution (HOD) emulator method to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding parameters that describe central galaxy incompleteness, galaxy assembly bias, and a scale-independent multiplicative lensing bias Alens. We use this emulator to forecast cosmological constraints attainable from the GGL surface density profile ΔΣ(rp) and the projected galaxy correlation function wp, gg(rp) in the final (Year 6) DES data set over scales $$r_p=0.3\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$$. For a $$3{{\ \rm per\ cent}}$$ prior on Alens we forecast precisions of $$1.9{{\ \rm per\ cent}}$$, $$2.0{{\ \rm per\ cent}}$$, and $$1.9{{\ \rm per\ cent}}$$ on Ωm, σ8, and $$S_8 \equiv \sigma _8\Omega _m^{0.5}$$, marginalized over all halo occupation distribution (HOD) parameters as well as Alens. Adding scales $$r_p=0.3\!-\!3.0\, h^{-1} \, \mathrm{Mpc}$$ improves the S8 precision by a factor of ∼1.6 relative to a large scale ($$3.0\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$$) analysis, equivalent to increasing the survey area by a factor of ∼2.6. Sharpening the Alens prior to $$1{{\ \rm per\ cent}}$$ further improves the S8 precision to $$1.1{{\ \rm per\ cent}}$$, and it amplifies the gain from including non-linear scales. Our emulator achieves per cent-level accuracy similar to the projected DES statistical uncertainties, demonstrating the feasibility of a fully non-linear analysis. Obtaining precise parameter constraints from multiple galaxy types and from measurements that span linear and non-linear clustering offers many opportunities for internal cross-checks, which can diagnose systematics and demonstrate the robustness of cosmological results. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT The contribution of dissolved globular clusters (GCs) to the stellar content of the Galactic halo is a key constraint on models for GC formation and destruction, and the mass assembly history of the Milky Way. Earlier results from APOGEE pointed to a large contribution of destroyed GCs to the stellar content of the inner halo, by as much as 25 $${{\ \rm per\ cent}}$$, which is an order of magnitude larger than previous estimates for more distant regions of the halo. We set out to measure the ratio between nitrogen-rich (N-rich) and normal halo field stars, as a function of distance, by performing density modelling of halo field populations in APOGEE DR16. Our results show that at 1.5 kpc from the Galactic Centre, N-rich stars contribute a much higher 16.8$$^{+10.0}_{-7.0}\, {{\ \rm per\ cent}}$$ fraction to the total stellar halo mass budget than the 2.7$$^{+1.0}_{-0.8}\, {{\ \rm per\ cent}}$$ ratio contributed at 10 kpc. Under the assumption that N-rich stars are former GC members that now reside in the stellar halo field, and assuming the ratio between first and second population GC stars being 1:2, we estimate a total contribution from disrupted GC stars of the order of 27.5$$^{+15.4}_{-11.5}\, {{\ \rm per\ cent}}$$ at r = 1.5 kpc and 4.2$$^{+1.5}_{-1.3}\, {{\ \rm per\ cent}}$$ at r = 10 kpc. Furthermore, since our methodology requires fitting a density model to the stellar halo, we integrate such density within a spherical shell from 1.5 to 15 kpc in radius, and find a total stellar mass arising from dissolved and/or evaporated GCs of MGC,total = 9.6$$^{+4.0}_{-2.6}\, \times$$ 107 M⊙. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT In this work, we combine spectroscopic information from the SkyMapper survey for Extremely Metal-Poor stars and astrometry from Gaia DR2 to investigate the kinematics of a sample of 475 stars with a metallicity range of $$-6.5 \le \rm [Fe/H] \le -2.05$$ dex. Exploiting the action map, we identify 16 and 40 stars dynamically consistent with the Gaia Sausage and Gaia Sequoia accretion events, respectively. The most metal poor of these candidates have metallicities of $$\rm [Fe/H]=-3.31\, \mathrm{ and }\, -3.74$$, respectively, helping to define the low-metallicity tail of the progenitors involved in the accretion events. We also find, consistent with other studies, that ∼21 per cent of the sample have orbits that remain confined to within 3 kpc of the Galactic plane, that is, |Zmax| ≤ 3 kpc. Of particular interest is a subsample (∼11 per cent of the total) of low |Zmax| stars with low eccentricities and prograde motions. The lowest metallicity of these stars has [Fe/H] = –4.30 and the subsample is best interpreted as the very low-metallicity tail of the metal-weak thick disc population. The low |Zmax|, low eccentricity stars with retrograde orbits are likely accreted, while the low |Zmax|, high eccentricity pro- and retrograde stars are plausibly associated with the Gaia Sausage system. We find that a small fraction of our sample (∼4 per cent of the total) is likely escaping from the Galaxy, and postulate that these stars have gained energy from gravitational interactions that occur when infalling dwarf galaxies are tidally disrupted. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email