No Arabic abstract
We introduce a new technique to constrain the line-of-sight integrated electron density of our Galactic halo $text{DM}_text{MW,halo}$ through analysis of the observed dispersion measure distributions of pulsars $text{DM}_text{pulsar}$ and fast radio bursts $text{DM}_text{FRB}$. We model these distributions, correcting for the Galactic interstellar medium, with kernel density estimation---well-suited to the small data regime---to find lower/upper bounds to the corrected $text{DM}_text{pulsar}$/$text{DM}_text{FRB}$ distributions: $max[text{DM}_text{pulsar}] approx 7pm2 text{ (stat)} pm 9 text{ (sys) pc cm}^{-3}$ and $min[text{DM}_text{FRB}] approx 63^{+27}_{-21} text{ (stat)} pm 9 text{ (sys) pc cm}^{-3}$. Using bootstrap resampling to estimate uncertainties, we set conservative limits on the Galactic halo dispersion measure $-2 < text{DM}_text{MW,halo} < 123 text{pc cm}^{-3}$ (95% c.l.). The upper limit is especially conservative because it may include a non-negligible contribution from the FRB host galaxies and a non-zero contribution from the cosmic web. It strongly disfavors models where the Galaxy has retained the majority of its baryons with a density profile tracking the presumed dark matter density profile. Last, we perform Monte Carlo simulations of larger FRB samples to validate our technique and assess the sensitivity of ongoing and future surveys. We recover bounds of several tens $text{pc cm}^{-3}$ which may be sufficient to test whether the Galaxy has retained a majority of its baryonic mass. We estimate that a sample of several thousand FRBs will significantly tighten constraints on $text{DM}_text{MW,halo}$ and offer a valuable complement to other analyses.
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 use a distribution function analysis to estimate the mass of the Milky Way out to 100 kpc using a large sample of halo stars. These stars are compiled from the literature, and the vast majority (~98%) have 6D phase-space information. We pay particular attention to systematic effects, such as the dynamical influence of the Large Magellanic Cloud (LMC), and the effect of unrelaxed substructure. The LMC biases the (pre-LMC infall) halo mass estimates towards higher values, while realistic stellar halos from cosmological simulations tend to underestimate the true halo mass. After applying our method to the Milky Way data we find a mass within 100 kpc of M(< 100 kpc) = 6.07 +/- 0.29 (stat.) +/- 1.21 (sys.) x 10^11 M_Sun. For this estimate, we have approximately corrected for the reflex motion induced by the LMC using the Erkal et al. model, which assumes a rigid potential for the LMC and MW. Furthermore, stars that likely belong to the Sagittarius stream are removed, and we include a 5% systematic bias, and a 20% systematic uncertainty based on our tests with cosmological simulations. Assuming the mass-concentration relation for Navarro-Frenk-White haloes, our mass estimate favours a total (pre-LMC infall) Milky Way mass of M_200c = 1.01 +/- 0.24 x 10^12 M_Sun, or (post-LMC infall) mass of M_200c = 1.16 +/- 0.24 x 10^12 M_Sun when a 1.5 x 10^11 M_Sun mass of a rigid LMC is included.
Theoretical and observational arguments suggest that there is a large amount of hot ($sim 10^6$ K), diffuse gas residing in the Milky Ways halo, while its total mass and spatial distribution are still unclear. In this work, we present a general model for the gas density distribution in the Galactic halo, and investigate the gas evolution under radiative cooling with a series of 2D hydrodynamic simulations. We find that the mass inflow rate in the developed cooling flow increases with gas metallicity and the total gas mass in the halo. For a fixed halo gas mass, the spatial gas distribution affects the onset time of the cooling catastrophe, which starts earlier when the gas distribution is more centrally-peaked, but does not substantially affect the final mass inflow rate. The gravity from the Galactic bulge and disk affects gas properties in inner regions, but has little effect on the final inflow rate either. We confirm our results by investigating cooling flows in several density models adopted from the literature, including the Navarro-Frenk-White (NFW) model, the cored-NFW model, the Maller & Bullock model, and the $beta$ model. Typical mass inflow rates in our simulations range from $sim 5 M_{odot}$ yr$^{-1}$ to $sim 60 M_{odot}$ yr$^{-1}$, and are much higher than the observed star formation rate in our Galaxy, suggesting that stellar and active galactic nucleus feedback processes may play important roles in the evolution of the Milky Way (MW) and MW-type galaxies.
The halo of the Milky Way provides a laboratory to study the properties of the shocked hot gas that is predicted by models of galaxy formation. There is observational evidence of energy injection into the halo from past activity in the nucleus of the Milky Way; however, the origin of this energy (star formation or supermassive-black-hole activity) is uncertain, and the causal connection between nuclear structures and large-scale features has not been established unequivocally. Here we report soft-X-ray-emitting bubbles that extend approximately 14 kiloparsecs above and below the Galactic centre and include a structure in the southern sky analogous to the North Polar Spur. The sharp boundaries of these bubbles trace collisionless and non-radiative shocks, and corroborate the idea that the bubbles are not a remnant of a local supernova but part of a vast Galaxy-scale structure closely related to features seen in gamma-rays. Large energy injections from the Galactic centre are the most likely cause of both the {gamma}-ray and X-ray bubbles. The latter have an estimated energy of around 10$^{56}$ erg, which is sufficient to perturb the structure, energy content and chemical enrichment of the circumgalactic medium of the Milky Way.
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.