No Arabic abstract
Halo stars orbit within the potential of the Milky Way and hence their kinematics can be used to understand the underlying mass distribution. However, the inferred mass distribution depends sensitively upon assumptions made on the density and the velocity anisotropy profiles of the tracers. Also, there is a degeneracy between the parameters of the halo and that of the disk or bulge. Here, we decompose the Galaxy into bulge, disk and dark matter halo and then model the kinematic data of the halo BHB and K-giants from the SEGUE. Additionally, we use the gas terminal velocity curve and the Sgr A$^*$ proper motion. With $R_odot = 8.5$kpc, our study reveals that the density of the stellar halo has a break at $17.2^{+1.1}_{-1.0}$ kpc, and an exponential cut-off in the outer parts starting at $97.7^{+15.6}_{-15.8}$kpc. Also, we find the velocity anisotropy is radially biased with $beta_s= 0.4pm{0.2}$ in the outer halo. We measure halo virial mass $M_{text{vir}} = 0.80^{+0.31}_{-0.16} times 10^{12} M_{odot}$, concentration $c=21.1^{+14.8}_{-8.3}$, disk mass of $0.95^{+0.24}_{-0.30}times10^{11} M_{odot}$, disk scale length of $4.9^{+0.4}_{-0.4}$ kpc and bulge mass of $0.91^{+0.31}_{-0.38} times 10^{10} M_{odot}$. The mass of halo is found to be small and this has important consequences. The giant stars reveal that the outermost halo stars have low velocity dispersion interestingly suggesting a truncation of the stellar halo density rather than a small overall mass of the Galaxy. Our estimates of local escape velocity $v_{rm esc} = 550.9^{+32.4}_{-22.1}$ kms$^{-1}$ and dark matter density $rho^{rm DM}_{odot} = 0.0088^{+0.0024}_{-0.0018} M_{odot} {rm pc^{-3}} $ ($0.35^{+0.08}_{-0.07}$ GeV cm$^{-3}$) are in good agreement with recent estimates. Some of the above estimates are depended on the adopted value of $R_odot$ and of outer power-law index of the tracer number density.
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.
In the $Gaia$ era stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern datasets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ``action diamond) is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially-anisotropic $Gaia$-Enceladus merger remnant, which we find achieves a sample purity of 82 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic datasets can noticeably impact the $E-L_{z}$ distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo.
The first and second moments of stellar velocities encode important information about the formation history of the Galactic halo. However, due to the lack of tangential motion and inaccurate distances of the halo stars, the velocity moments in the Galactic halo have largely remained known unknowns. Fortunately, our off-centric position within the Galaxy allows us to estimate these moments in the galacto-centric frame using the observed radial velocities of the stars alone. We use these velocities coupled with the Hierarchical Bayesian scheme, which allows easy marginalisation over the missing data (the proper-motion, and uncertainty-free distance and line-of-sight velocity), to measure the velocity dispersions, orbital anisotropy ($beta$) and streaming motion ($v_{rm rot}$) of the halo main-sequence turn-off (MSTO) and K-giant (KG) stars in the inner stellar halo (r $lesssim 15$ kpc). We study the metallicity bias in kinematics of the halo stars and observe that the comparatively metal-rich ([Fe/H]$>-1.4$) and the metal-poor ([Fe/H]$leq - 1.4$) MSTO samples show a clear systematic difference in $v_{rm rot} sim 20-40$ km s$^{-1}$, depending on how restrictive the spatial cuts to cull the disk contamination are. The bias is also detected in KG samples but with less certainty. Both MSTO and KG populations suggest that the inner stellar halo of the Galaxy is radially biased i.e. $sigma_r>sigma_theta$ or $sigma_phi$ and $beta simeq 0.5$. The apparent metallicity contrariety in the rotation velocity among the halo sub-populations supports the co-existence of multiple populations in the galactic halo that may have formed through distinct formation scenarios, i.e. in-situ versus accretion.
We analyse systems analogous to the Milky Way (MW) in the EAGLE cosmological hydrodynamics simulation in order to deduce the likely structure of the MWs dark matter halo. We identify MW-mass haloes in the simulation whose satellite galaxies have similar kinematics and spatial distribution to those of the bright satellites of the MW, specifically systems in which the majority of the satellites (8 out of 11) have nearly co-planar orbits that are also perpendicular to the central stellar disc. We find that the normal to the common orbital plane of the co-planar satellites is well aligned with the minor axis of the host dark matter halo, with a median misalignment angle of only $17.3^circ$. Based on this result, we infer that the minor axis of the Galactic dark matter halo points towards $(l,b)=(182^circ,-2^circ)$, with an angular uncertainty at the 68 and 95 percentile confidence levels of 22$^circ$ and 43$^circ$ respectively. Thus, the inferred minor axis of the MW halo lies in the plane of the stellar disc. The halo, however, is not homologous and its flattening and orientation vary with radius. The inner parts of the halo are rounder than the outer parts and well-aligned with the stellar disc (that is the minor axis of the halo is perpendicular to the disc). Further out, the halo twists and the minor axis changes direction by $90^circ$. This twist occurs over a very narrow radial range and reflects variations in the filamentary network along which mass was accreted into the MW.
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.