No Arabic abstract
Dynamical models allow us to connect the motion of a set of tracers to the underlying gravitational potential, and thus to the total (luminous and dark) matter distribution. They are particularly useful for understanding the mass and spatial distribution of dark matter (DM) in a galaxy. Globular clusters (GCs) are an ideal tracer population in dynamical models, since they are bright and can be found far out into the halo of galaxies. We aim to test how well Jeans-Anisotropic-MGE (JAM) models using GCs (positions and line-of-sight velocities) as tracers can constrain the mass and radial distribution of DM halos. For this, we use the E-MOSAICS suite of 25 zoom-in simulations of L* galaxies. We find that the DM halo properties are reasonably well recovered by the JAM models. There is, however, a strong correlation between how well we recover the mass and the radial distribution of the DM and the number of GCs in the galaxy: the constraints get exponentially worse with fewer GCs, and at least 150 GCs are needed in order to guarantee that the JAM model will perform well. We find that while the data quality (uncertainty on the radial velocities) can be important, the number of GCs is the dominant factor in terms of the accuracy and precision of the measurements. This work shows promising results for these models to be used in extragalactic systems with a sample of more than 150 GCs.
We measure the total stellar halo luminosity using red giant branch (RGB) stars selected from Gaia data release 2. Using slices in magnitude, colour and location on the sky, we decompose RGB stars belonging to the disc and halo by fitting 2-dimensional Gaussians to the Galactic proper motion distributions. The number counts of RGB stars are converted to total stellar halo luminosity using a suite of isochrones weighted by age and metallicity, and by applying a volume correction based on the stellar halo density profile. Our method is tested and calibrated using Galaxia and N-body models. We find a total luminosity (out to 100 kpc) of L_halo = 7.9 +/- 2.0 x 10^8 L_Sun excluding Sgr, and L_halo = 9.4 +/- 2.4 x 10^8 L_Sun including Sgr. These values are appropriate for our adopted stellar halo density profile and metallicity distribution, but additional systematics related to these assumptions are quantified and discussed. Assuming a stellar mass-to-light ratio appropriate for a Kroupa initial mass function (M*/L = 1.5), we estimate a stellar halo mass of M*_halo = 1.4 +/- 0.4 x 10^9 M_Sun. This mass is larger than previous estimates in the literature, but is in good agreement with the emerging picture that the (inner) stellar halo is dominated by one massive dwarf progenitor. Finally, we argue that the combination of a ~10^9 M_Sun mass and an average metallicity of <[Fe/H]> ~ -1.5 for the Galactic halo points to an ancient (~10 Gyr) merger event.
The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase-space distribution of dynamical tracers. We test this approach using realistic mock stellar halos constructed from the Aquarius N-body simulations of dark matter halos in the $Lambda$CDM cosmology. We extend the standard treatment to include a Navarro-Frenk-White (NFW) potential and use a maximum likelihood method to recover the parameters describing the simulated halos from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is highly correlated with the estimate of halo concentration. The best-fit halo masses within the virial radius, $R_{200}$, are biased, ranging from a 40% underestimate to a 5% overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass within the half-mass radius of tracers even when restricted to tracers within 60kpc. The recovered velocity anisotropies of tracers, $beta$, are biased systematically, but this does not affect other parameters if tangential velocity data are used as constraints.
New mass estimates and cumulative mass profiles with Bayesian credible regions (c.r.) for the Milky Way (MW) are found using the Galactic Mass Estimator (GME) code and dwarf galaxy (DG) kinematic data from multiple sources. GME takes a hierarchical Bayesian approach to simultaneously estimate the true positions and velocities of the DGs, their velocity anisotropy, and the model parameters for the Galaxys total gravitational potential. In this study, we incorporate meaningful prior information from past studies and simulations. The prior distributions for the physical model are informed by the results of Eadie & Juric (2019), which used globular clusters instead of DGs, as well as by the subhalo distributions of the Ananke Gaia-like surveys from Feedback In Realistic Environments-2 (Fire-2) cosmological simulations (see Sanderson et al. 2020). Using DGs beyond 45 kpc, we report median and 95% c.r estimates for $r_{200}$ = 212.8 (191.12,238.44) kpc, and for the total enclosed mass $M_{200}$ = 1.19 (0.87,1.68)$times10^{12}M_{odot}$ (adopting $Delta_c=200$). Median mass estimates at specific radii are also reported (e.g., $M(<50text{ kpc})=0.52times10^{12}M_{odot}$ and $M(100text{ kpc})=0.78times10^{12}M_{odot}$). Estimates are comparable to other recent studies using GAIA DR2 and DGs, but notably different from the estimates of Eadie & Juric (2019). We perform a sensitivity analysis to investigate whether individual DGs and/or a more massive Large Magellanic Cloud (LMC) on the order of $10^{11}M_{odot}$ may be affecting our mass estimates. We find possible supporting evidence for the idea that some DGs are affected by a massive LMC and are not in equilibrium with the MW.
We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{rm GC} approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing line-of-sight velocities. We confirm the association of Arp 2, Pal 12, Terzan 7, and Terzan 8 with the Sagittarius (Sgr) stream. These clusters and NGC 6101 have tangential velocity $V_{rm tan} > 290$ km s$^{-1}$, whereas all other clusters have $V_{rm tan} < 200$ km s$^{-1}$. NGC 2419, the most distant GC in our sample, is also likely associated with the Sgr stream, whereas NGC 4147, NGC 5024, and NGC 5053 definitely are not. We use the distribution of orbital parameters derived using the 3D velocities to separate halo GCs that either formed within the MW or were accreted. We also assess the specific formation history of e.g. Pyxis and Terzan 8. We constrain the MW mass via an estimator that considers the full 6D phase-space information for 16 of the GCs from $R_{rm GC} = 10$ to 40 kpc. The velocity dispersion anisotropy parameter $beta = 0.609^{+0.130}_{-0.229}$. The enclosed mass $M (<39.5 rm{kpc}) = 0.61^{+0.18}_{-0.12} times 10^{12}$ M$_{odot}$, and the virial mass $M_rm{vir} = 2.05^{+0.97}_{-0.79} times 10^{12}$ M$_{odot}$, are consistent with, but on the high side among recent mass estimates in the literature.
We report on the extent of the effects of the Milky Ways gravitational field in shaping the structural parameters and internal dynamics of its globular cluster population. We make use of a homogeneous, up-to-date data set with kinematics, structural properties, current and initial masses of 156 globular clusters. In general, cluster radii increase as the Milky Way potential weakens; with the core and Jacobi radii being those which increase at the slowest and fastest rate respectively. We interpret this result as the innermost regions of globular clusters being less sensitive to changes in the tidal forces with the Galactocentric distance. The Milky Ways gravitational field also seems to have differentially accelerated the internal dynamical evolution of individual clusters, with those toward the bulge appearing dynamically older. Finally we find a sub-population consisting of both compact and extended globular clusters (as defined by their rh/rJ ratio) beyond 8 kpc that appear to have lost a large fraction of their initial mass lost via disruption. Moreover, we identify a third group with rh/rJ > 0.4, which have lost an even larger fraction of their initial mass by disruption. In both cases the high fraction of mass lost is likely due to their large orbital eccentricities and inclination angles, which lead to them experiencing more tidal shocks at perigalacticon and during disc crossings. Comparing the structural and orbital parameters of individual clusters allows for constraints to be placed on whether or not their evolution was relaxation or tidally dominated.