We use recent proper motion measurements of the tangential velocity of M31, along with its radial velocity and distance, to derive the likelihood of the sum of halo masses of the Milky Way and M31. This is done using a sample halo pairs in the Bolsho
i cosmological simulation of $Lambda$CDM cosmology selected to match properties and environment of the Local Group. The resulting likelihood gives estimate of the sum of masses of $M_{rm MW,200}+M_{rm M31,200}=$ $2.40_{-1.05}^{+1.95}times10^{12},M_{odot}$ ($90%$ confidence interval). This estimate is consistent with individual mass estimates for the Milky Way and M31 and is consistent, albeit somewhat on the low side, with the mass estimated using the timing argument. We show that although the timing argument is unbiased on average for all pairs, for pairs constrained to have radial and tangential velocities similar to that of the Local Group the argument overestimates the sum of masses by a factor of $1.6$. Using similar technique we estimate the total dark matter mass enclosed within $1$ Mpc from the Local Group barycenter to be $M_{rm LG}(r<1, {rm Mpc})=4.2_{-2.0}^{+3.4}times10^{12},M_{odot}$ ($90%$ confidence interval).
Previous studies showed that an estimate of the likelihood distribution of the Milky Way halo mass can be derived using the properties of the satellites similar to the Large and Small Magellanic Clouds (LMC and SMC). However, it would be straightforw
ard to interpret such an estimate only if the properties of the Magellanic Clouds (MCs) are fairly typical and are not biased by the environment. In this study we explore whether the environment of the Milky Way affects the properties of the SMC and LMC such as their velocities. To test for the effect of the environment, we compare velocity distributions for MC-sized subhalos around Milky Way hosts in a sample selected simply by mass and in the second sample of such halos selected with additional restrictions on the distance to the nearest cluster and the local galaxy density, designed to mimic the environment of the Local Group (LG). We find that satellites in halos in the LG-like environments do have somewhat larger velocities, as compared to the halos of similar mass in the sample without environmental constraints. We derive the host halo likelihood distribution for the samples in the LG-like envirionment and in the control sample and find that the environment does not significantly affect the derived likelihood. We use the updated properties of the SMC and LMC to derive the constraint on the MW halo mass $log{({rm M}_{200} /msol)}=12.06^{+0.31}_{-0.19}$ (90% confidence interval). We also explore the incidence of close pairs with relative velocities and separations similar to those of the LMC and SMC and find that such pairs are quite rare among $Lambda$CDM halos. Taking into account the close separation of the MCs in the Busha et al. 2011 method results in the shift of the MW halo mass estimate to smaller masses, with the peak shifting approximately by a factor of two.[Abridged]
The tight relation of star formation with molecular gas indicated by observations and assumed in recent models implies that the efficiency with which galaxies convert their gas into stars depends on gas metallicity. This is because the abundance of m
olecular hydrogen is sensitive to the abundance of dust, which catalyzes the formation of H_2 and helps to shield it from dissociating radiation. In this study we point out that in the absence of significant pre-enrichment by Population III stars forming out of zero metallicity gas, such H_2-based star formation is expected to leave an imprint in the form of bi-modality in the metallicity distribution among dwarf galaxies and in the metallicity distribution of stars within individual galaxies. The bi-modality arises because when gas metallicity (and dust abundance) is low, formation of molecular gas is inefficient, the gas consumption time scale is long, and star formation and metal enrichment proceed slowly. When metallicity reaches a critical threshold value star formation and enrichment accelerate, which leads to rapid increase in both stellar mass and metallicity of galaxies. We demonstrate this process both using a simple analytical model and full cosmological simulations. In contrast, observed metallicity distributions of dwarf galaxies or stars within them are not bi-modal. We argue that this discrepancy points to substantial early stochastic pre-enrichment by population III stars to levels Z ~ 0.01 Z_sun in dense, star forming regions of early galaxies.
We calculate the probability distribution function (PDF) of the expected annihilation luminosities of dark matter subhalos as a function of subhalo mass and distance from the Galactic center using a semi-analytical model of halo evolution. We find th
at the PDF of luminosities is relatively broad, exhibiting a spread of as much as an order of magnitude at fixed subhalo mass and halo-centric distance. The luminosity PDF allows for simple construction of mock samples of gamma-ray luminous subhalos and assessment of the variance in among predicted gamma-ray signals from dark matter annihilation. Other applications include quantifying the variance among the expected luminosities of dwarf spheroidal galaxies, assessing the level at which dark matter annihilation can be a contaminant in the expected gamma-ray signal from other astrophysical sources, as well as estimating the level at which nearby subhalos can contribute to the antimatter flux.
Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate
that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.
(Abridged) We perform dissipationless N-body simulations to elucidate the dynamical response of thin disks to bombardment by cold dark matter (CDM) substructure. Our method combines (1) cosmological simulations of the formation of Milky Way (MW)-size
d CDM halos to derive the properties of substructure and (2) controlled numerical experiments of consecutive subhalo impacts onto an initially-thin, fully-formed MW type disk galaxy. The present study is the first to account for the evolution of satellite populations over cosmic time in such an investigation of disk structure. We find that accretions of massive subhalos onto the central regions of host halos, where the galactic disks reside, since z~1 should be common. One host halo accretion history is used to initialize the controlled simulations of satellite-disk encounters. We show that these accretion events severely perturb the thin galactic disk and produce a wealth of distinctive dynamical signatures on its structure and kinematics. These include (1) considerable thickening and heating at all radii, with the disk thickness and velocity ellipsoid nearly doubling at the solar radius; (2) prominent flaring associated with an increase in disk thickness greater than a factor of 4 in the disk outskirts; (3) surface density excesses at large radii, beyond ~5 disk scale lengths, resembling those of observed antitruncated disks; (4) lopsidedness at levels similar to those measured in observational samples of disk galaxies; and (5) substantial tilting. The interaction with the most massive subhalo drives the disk response while subsequent bombardment is much less efficient at disturbing the disk. We conclude that substructure-disk encounters of the kind expected in the LCDM paradigm play a significant role in setting the structure of disk galaxies and driving galaxy evolution.
