No Arabic abstract
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.
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.
The total number and luminosity function of the population of dwarf galaxies of the Milky Way (MW) provide important constraints on the nature of the dark matter and on the astrophysics of galaxy formation at low masses. However, only a partial census of this population exists because of the flux limits and restricted sky coverage of existing Galactic surveys. We combine the sample of satellites recently discovered by the Dark Energy Survey (DES) with the satellites found in Sloan Digital Sky Survey (SDSS) Data Release 9 (together these surveys cover nearly half the sky) to estimate the total luminosity function of satellites down to $M_{rm V}=0$. We apply a new Bayesian inference method in which we assume that the radial distribution of satellites independently of absolute magnitude follows that of subhaloes selected according to their peak maximum circular velocity. We find that there should be at least $124^{+40}_{-27}$ (68 per cent CL, statistical error) satellites brighter than $M_{rm V}=0$ within $300$ kpc of the Sun. As a result of our use of new data and better simulations, and a more robust statistical method, we infer a much smaller population of satellites than reported in previous studies using earlier SDSS data only; we also address an underestimation of the uncertainties in earlier work by accounting for stochastic effects. We find that the inferred number of faint satellites depends only weakly on the assumed mass of the MW halo and we provide scaling relations to extend our results to different assumed halo masses and outer radii. We predict that half of our estimated total satellite population of the MW should be detected by the Large Synoptic Survey Telescope. The code implementing our estimation method is available online.
We propose a novel method to constrain the Milky Way (MW) mass $M_{rm vir}$ with its corona temperature observations. For a given corona density profile, one can derive its temperature distribution assuming a generalized equilibrium model with non-thermal pressure support. While the derived temperature profile decreases substantially with radius, the X-ray-emission-weighted average temperature, which depends most sensitively on $M_{rm vir}$, is quite uniform toward different sight lines, consistent with X-ray observations. For an Navarro-Frenk-White (NFW) total matter distribution, the corona density profile should be cored, and we constrain $M_{rm vir}=(1.19$ - $2.95) times 10^{12} M_{rm sun}$. For a total matter distribution contributed by an NFW dark matter profile and central baryons, the corona density profile should be cuspy and $M_{rm vir,dm}=(1.34$ - $5.44) times 10^{12} M_{rm sun}$. Non-thermal pressure support leads to even higher values of $M_{rm vir}$, while a lower MW mass may be possible if the corona is accelerating outward. This method is independent of the total corona mass, its metallicity, and temperature at very large radii.
We introduce the ARTEMIS simulations, a new set of 42 zoomed-in, high-resolution (baryon particle mass of ~ 2x10^4 Msun/h), hydrodynamical simulations of galaxies residing in haloes of Milky Way mass, simulated with the EAGLE galaxy formation code with re-calibrated stellar feedback. In this study, we analyse the structure of stellar haloes, specifically the mass density, surface brightness, metallicity, colour and age radial profiles, finding generally very good agreement with recent observations of local galaxies. The stellar density profiles are well fitted by broken power laws, with inner slopes of ~ -3, outer slopes of ~ -4 and break radii that are typically ~ 20-40 kpc. The break radii generally mark the transition between in situ formation and accretion-driven formation of the halo. The metallicity, colour and age profiles show mild large-scale gradients, particularly when spherically-averaged or viewed along the major axes. Along the minor axes, however, the profiles are nearly flat, in agreement with observations. Overall, the structural properties can be understood by two factors: that in situ stars dominate the inner regions and that they reside in a spatially-flattened distribution that is aligned with the disc. Observations targeting both the major and minor axes of galaxies are thus required to obtain a complete picture of stellar haloes.