Search

We investigate the presence and importance of dark matter discs in a sample of 24 simulated Milky Way galaxies in the APOSTLE project, part of the EAGLE programme of hydrodynamic simulations in Lambda-CDM cosmology. It has been suggested that a dark disc in the Milky Way may boost the dark matter density and modify the velocity modulus relative to a smooth halo at the position of the Sun, with ramifications for direct detection experiments. From a kinematic decomposition of the dark matter and a real space analysis of all 24 halos, we find that only one of the simulated Milky Way analogues has a detectable dark disc component. This unique event was caused by a merger at late time with an LMC-mass satellite at very low grazing angle. Considering that even this rare scenario only enhances the dark matter density at the solar radius by 35% and affects the high energy tail of the dark matter velocity distribution by less than 1%, we conclude that the presence of a dark disc in the Milky Way is unlikely, and is very unlikely to have a significant effect on direct detection experiments., Comment: 6 pages, 2 figures. Accepted for publication in MNRAS. Data available on request

We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $\alpha$-low (`thin disc') and $\alpha$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $\alpha$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrm{\phi}}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrm{\phi}}$/$\partial$[Fe/H] for the $\alpha$-high disc, consistent with results from high-resolution surveys. We also find differences between the $\alpha$-low and $\alpha$-high discs in all three components of velocity dispersion. We discuss the implications of an $\alpha$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$\alpha$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components., Comment: 11 pages, 7 figures, accepted for publication in MNRAS

The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius, and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we make use of the APOSTLE cosmological hydrodynamical simulations of the Local Group to test the validity and accuracy of such mass estimators for both dispersion and rotation-supported galaxies, for field and satellite galaxies, and for galaxies of varying masses, shapes, and velocity dispersion anisotropies. We find that the mass estimators of Walker et al. and Wolf et al. are able to recover the masses of dispersion-dominated systems with little systematic bias, but with a 1-sigma scatter of 25 and 23 percent, respectively. The error on the estimated mass is dominated by the impact of the 3D shape of the stellar mass distribution, which is difficult to constrain observationally. This intrinsic scatter becomes the dominant source of uncertainty in the masses estimated for galaxies like the dwarf spheroidal (dSph) satellites of the Milky Way, where the observational errors in their sizes and velocity dispersions are small. Such scatter may also affect the inner density slopes of dSphs derived from multiple stellar populations, relaxing the significance with which Navarro-Frenk-White profiles may be excluded, depending on the degree to which the relevant properties of the different stellar populations are correlated. Finally, we derive a new optimal mass estimator that removes the residual biases and achieves a statistically significant reduction in the scatter to 20 percent overall for dispersion-dominated galaxies, allowing more precise and accurate mass estimates., Comment: Accepted for publication in MNRAS. 26 pages, 21 figures, 3 tables

Based on data aquired in 13 orbits of HST time, we present a detailed evolutionary history of the M31 dSph satellite Andromeda XVI, including its life-time star formation history, the spatial distribution of its stellar populations, and the properties of its variable stars. And XVI is characterized by prolonged star formation activity from the oldest epochs until star formation was quenched ~6 Gyr ago, and, notably, only half of the mass in stars of And XVI was in place 10 Gyr ago. And XVI appears to be a low mass galaxy for which the early quenching by either reionization or starburst feedback seems highly unlikely, and thus, is most likely due to an environmental effect (e.g., an interaction), possibly connected to a late infall in the densest regions of the Local Group. Studying the star formation history as a function of galactocentric radius, we detect a mild gradient in the star formation history: the star formation activity between 6 and 8 Gyr ago is significantly stronger in the central regions than in the external regions, although the quenching age appears to be the same, within 1 Gyr. We also report the discovery of 9 RR Lyrae stars, 8 of which belong to And XVI. The RR Lyrae stars allow a new estimate of the distance, (m-M)0= 23.72+/-0.09 mag, which is marginally larger than previous estimates based on the tip of the red giant branch., Comment: Accepted for publication on ApJ

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation., Comment: 10 pages, 7 figures, submitted to MNRAS

The Next Generation Virgo Cluster Survey (NGVS) is a deep $u^*giz$ survey targeting the Virgo cluster of galaxies at 16.5~Mpc. This survey provides high-quality photometry over an $\sim$ 100 deg$^2$ region straddling the constellations of Virgo and Coma Berenices. This sightline through the Milky Way is noteworthy in that it intersects two of the most prominent substructures in the Galactic halo: the Virgo Over-Density (VOD) and Sagittarius stellar stream (close to its bifurcation point). In this paper, we use deep $u^*gi$ imaging from the NGVS to perform tomography of the VOD and Sagittarius stream using main-sequence turnoff (MSTO) stars as a halo tracer population. The VOD, whose centroid is known to lie at somewhat lower declinations ($\alpha \sim 190^\circ$, $\delta \sim -5^\circ$) than is covered by the NGVS, is nevertheless clearly detected in the NGVS footprint at distances between $\sim$ 8 and 25~kpc. By contrast, the Sagittarius stream is found to slice directly across the NGVS field at distances between 25 and 40~kpc, with a density maximum at $\simeq$~35~kpc. No evidence is found for new substructures beyond the Sagittarius stream, at least out to a distance of $\sim$~90 kpc --- the largest distance to which we can reliably trace the halo using MSTO stars. We find clear evidence for a distance gradient in the Sagittarius stream across the $\sim 30$~deg of sky covered by the NGVS and its flanking fields. We compare our distance measurements along the stream to those predicted by leading stream models., Comment: 14 pages, 10 figures, accepted for publication in ApJ

We study the implications of galaxy formation on dark matter direct detection using high resolution hydrodynamic simulations of Milky Way-like galaxies simulated within the EAGLE and APOSTLE projects. We identify Milky Way analogues that satisfy observational constraints on the Milky Way rotation curve and total stellar mass. We then extract the dark matter density and velocity distribution in the Solar neighbourhood for this set of Milky Way analogues, and use them to analyse the results of current direct detection experiments. For most Milky Way analogues, the event rates in direct detection experiments obtained from the best fit Maxwellian distribution (with peak speed of 223 - 289 km/s) are similar to those obtained directly from the simulations. As a consequence, the allowed regions and exclusion limits set by direct detection experiments in the dark matter mass and spin-independent cross section plane shift by a few GeV compared to the Standard Halo Model, at low dark matter masses. For each dark matter mass, the halo-to-halo variation of the local dark matter density results in an overall shift of the allowed regions and exclusion limits for the cross section. However, the compatibility of the possible hints for a dark matter signal from DAMA and CDMS-Si and null results from LUX and SuperCDMS is not improved., Comment: 43 pages, 14 figures, 7 tables, accepted for publication in JCAP

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $\Lambda$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations., Comment: 21 pages, 16 figures

We use cosmological hydrodynamical simulations of the APOSTLE project along with high-quality rotation curve observations to examine the fraction of baryons in {\Lambda}CDM haloes that collect into galaxies. This 'galaxy formation efficiency' correlates strongly and with little scatter with halo mass, dropping steadily towards dwarf galaxies. The baryonic mass of a galaxy may thus be used to place a lower limit on total halo mass and, consequently, on its asymptotic maximum circular velocity. A number of observed dwarfs seem to violate this constraint, having baryonic masses up to ten times higher than expected from their rotation speeds, or, alternatively, rotating at only half the speed expected for their mass. Taking the data at face value, either these systems have formed galaxies with extraordinary efficiency - highly unlikely given their shallow potential wells - or their dark matter content is much lower than expected from {\Lambda}CDM haloes. This 'missing dark matter' is reminiscent of the inner mass deficit of galaxies with slowly-rising rotation curves, but cannot be explained away by star formation-induced 'cores' in the dark mass profile, since the anomalous deficit applies to regions larger than the luminous galaxies themselves. We argue that explaining the structure of these galaxies would require either substantial modification of the standard Lambda cold dark matter paradigm or else significant revision to the uncertainties in their inferred mass profiles, which should be much larger than reported. Systematic errors in inclination may provide a simple resolution to what would otherwise be a rather intractable problem for the current paradigm., Comment: 15 pages, 5 figures, 2 tables, MNRAS accepted version

We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few examples where stars form in two populations distinct in age in order to elucidate the origin of these gradients. The spatial distributions of the two components reflect their diverse origin; the old stellar component is assembled through mergers, but the young population forms largely in situ. The older component results from a first episode of star formation that begins early but is quickly shut off by the combined effects of stellar feedback and reionization. The younger component forms when a late accretion event adds gas and reignites star formation. The effect of mergers is to disperse the old stellar population, increasing their radius and decreasing their central density relative to the young population. We argue that dwarf-dwarf mergers offer a plausible scenario for the formation of systems with multiple distinct populations and, more generally, for the origin of age and metallicity gradients in dwarf spheroidals., Comment: 10 pages, 8 figures, Accepted for publication in MNRAS

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" ({\Lambda}CDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in {\Lambda}CDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals., Comment: 15 pages, submitted to MNRAS. arXiv admin note: text overlap with arXiv:1412.2748

We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies (APOSTLE project) carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed in the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White (NFW) profile between $r\approx1.5$ kpc and $r\approx10$ kpc. At smaller radii, $r\lesssim1.5$ kpc, the halos develop a flatter than NFW slope. This unexpected feature may be specific to our particular choice of subgrid physics model but nevertheless the dark matter density profiles agree within 30% as the mass resolution is increased by a factor 150. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.97$ within $r=1$ kpc) and there is no offset larger than 45 pc between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $\gamma$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three., Comment: 11 pages, 5 figures, accepted for publication in MNRAS, corrected typos

We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group of Galaxies (LG). We find, in agreement with the "timing argument" and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of $\sim 5\times 10^{12} \,M_{\odot}$. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of $\sim 2$ times smaller than suggested by the timing argument, but with large dispersion. Guided by these results, we have selected $12$ halo pairs with total mass in the range $1.6$-$3.6 \times 10^{12}\,M_{\odot}$ for the APOSTLE project (A Project Of Simulating The Local Environment), a suite of hydrodynamical resimulations at various numerical resolution levels (reaching up to $\sim10^{4}\,M_{\odot}$ per gas particle) that use the subgrid physics developed for the EAGLE project. These simulations reproduce, by construction, the main kinematics of the MW-M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. The APOSTLE candidate systems thus provide an excellent testbed to confront directly many of the predictions of the $\Lambda$CDM cosmology with observations of our local Universe., Comment: replaced with the accepted version

We examine the circular velocity profiles of galaxies in {\Lambda}CDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density "core". The "cusp vs core" issue is thus better characterized as an "inner mass deficit" problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of "inner mass deficit" galaxies inferred from kinematic data are incorrect., Comment: 17 pages, 6 figures, 2 tables, MNRAS accepted

The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the "missing satellites" problem), galaxies seem to avoid the largest substructures (the "too big to fail" problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three "problems" are resolved within the LCDM paradigm., Comment: 14 pages, 4 figures

We use cosmological N-body simulations from the Aquarius Project to study the tidal effects of a dark matter halo on the shape and orientation of its substructure. Although tides are often assumed to enhance asphericity and to stretch subhaloes tangentially, these effects are short lived: as in earlier work, we find that subhaloes affected by tides become substantially more spherical and show a strong radial alignment toward the centre of the host halo. These results, combined with a semi-analytic model of galaxy formation, may be used to assess the effect of Galactic tides on the observed population of dwarf spheroidal (dSph) satellites of the Milky Way and Andromeda galaxies. If, as the model suggests, the relatively low dark matter content of luminous dSphs such as Fornax and Leo I is due to tidal stripping, then their gravitational potential must be substantially more spherical than that of more heavily dark matter-dominated systems such as Draco or Carina. The model also predicts a tidally-induced statistical excess of satellites whose major axis aligns with the direction to the central galaxy. We find tantalizing evidence of this in the M31 satellite population, which suggests that tides may have played an important role in its evolution., Comment: 16 pages, 13 figures, submitted to MNRAS

We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope., Comment: 20 pages, 13 figures, minor changes to match version in press in MNRAS

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation., Comment: Accepted for publication in MNRAS. V2: Minor changes. For images and videos, see http://eagle.strw.leidenuniv.nl/ and http://icc.dur.ac.uk/Eagle/

Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating factors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution., Comment: 12 pages, 12 figures, submitted to MNRAS, comments welcome

We explore the impact of cosmic reionization on nearby isolated dwarf galaxies using a compilation of SFHs estimated from deep HST data and a cosmological hydrodynamical simulation of the Local Group. The nearby dwarfs show a wide diversity of star formation histories; from ancient systems that have largely completed their star formation $\sim 10$ Gyr ago to young dwarfs that have formed the majority of their stars in the past $\sim 5$ Gyr to two-component systems characterized by the overlap of comparable numbers of old and young stars. Taken as an ensemble, star formation in nearby dwarfs dips to lower-than-average rates at intermediate times ($4

We present the first comparison between the lifetime star formation histories (SFHs) of M31 and Milky Way (MW) satellites. Using the Advanced Camera for Surveys aboard the Hubble Space Telescope, we obtained deep optical imaging of Andromeda II (M$_{V} = -$12.0; log(M$_{\star}$/M$_{\odot}$) $\sim$ 6.7) and Andromeda XVI (M$_{V} = -$7.5; log(M$_{\star}$/M$_{\odot}$) $\sim$ 4.9) yielding color-magnitude diagrams that extend at least 1 magnitude below the oldest main sequence turnoff, and are similar in quality to those available for the MW companions. And II and And XVI show strikingly similar SFHs: both formed 50-70% of their total stellar mass between 12.5 and 5 Gyr ago (z$\sim$5-0.5) and both were abruptly quenched $\sim$ 5 Gyr ago (z$\sim$0.5). The predominance of intermediate age populations in And XVI makes it qualitatively different from faint companions of the MW and clearly not a pre-reionization fossil. Neither And II nor And XVI appears to have a clear analog among MW companions, and the degree of similarity in the SFHs of And II and And XVI is not seen among comparably faint-luminous pairs of MW satellites. These findings provide hints that satellite galaxy evolution may vary substantially among hosts of similar stellar mass. Although comparably deep observations of more M31 satellites are needed to further explore this hypothesis, our results underline the need for caution when interpreting satellite galaxies of an individual system in a broader cosmological context., Comment: 7 pages, 4 Figures. Accepted for publication in ApJ

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation., Comment: 8 pages, submitted to MNRAS, comments welcome

We reveal the highly structured nature of the Milky Way stellar halo within the footprint of the PAndAS photometric survey from blue main sequence and main sequence turn-off stars. We map no fewer than five stellar structures within a heliocentric range of ~5 to 30 kpc. Some of these are known (the Monoceros Ring, the Pisces/Triangulum globular cluster stream), but we also uncover three well-defined stellar structures that could be, at least partly, responsible for the so-called Triangulum/Andromeda and Triangulum/Andromeda 2 features. In particular, we trace a new faint stellar stream located at a heliocentric distance of ~17 kpc. With a surface brightness of \Sigma_V ~ 32-32.5 mag/arcsec^2, it follows an orbit that is almost parallel to the Galactic plane north of M31 and has so far eluded surveys of the Milky Way halo as these tend to steer away from regions dominated by the Galactic disk. Investigating our follow-up spectroscopic observations of PAndAS, we serendipitously uncover a radial velocity signature from stars that have colors and magnitudes compatible with the stream. From the velocity of eight likely member stars, we show that this stellar structure is dynamically cold, with an unresolved velocity dispersion that is lower than 7.1 km/s at the 90-percent confidence level. Along with the width of the stream (300-650 pc), its dynamics points to a dwarf-galaxy-accretion origin. The numerous stellar structures we can map in the Milky Way stellar halo between 5 and 30 kpc and their varying morphology is a testament to the complex nature of the stellar halo at these intermediate distances., Comment: 11 pages, 8 figures, accepted for publication in the ApJ, Figure 3 is the money plot

We present an analysis of the star formation history (SFH) of a field near the half light radius in the Local Group dwarf irregular galaxy IC 1613 based on deep Hubble Space Telescope Advanced Camera for Surveys imaging. Our observations reach the oldest main sequence turn-off, allowing a time resolution at the oldest ages of ~1 Gyr. Our analysis shows that the SFH of the observed field in IC 1613 is consistent with being constant over the entire lifetime of the galaxy. These observations rule out an early dominant episode of star formation in IC 1613. We compare the SFH of IC 1613 with expectations from cosmological models. Since most of the mass is in place at early times for low mass halos, a naive expectation is that most of the star formation should have taken place at early times. Models in which star formation follows mass accretion result in too many stars formed early and gas mass fractions which are too low today (the "over-cooling problem"). The depth of the present photometry of IC 1613 shows that, at a resolution of ~1 Gyr, the star formation rate is consistent with being constant, at even the earliest times, which is difficult to achieve in models where star formation follows mass assembly., Comment: 13 pages, 12 figures, accepted for publication in the ApJ

We analyze the Aquarius simulations to characterize the shape of dark matter halos with peak circular velocity in the range 8~0.75 and ~0.60 at r~1 kpc. Their velocity ellipsoid become strongly tangentially biased in the outskirts as a consequence of tidal stripping., Comment: 12 pages, 11 figures. MNRAS in press

We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift'' of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height'' mass parameter, $\nu(M,z)=\delta_{\rm crit}(z)/\sigma(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos., Comment: 13 pages, 10 figures, submitted to MNRAS

Counterrotating stars in disk galaxies are a puzzling dynamical feature whose origin has been ascribed to either satellite accretion events or to disk instabilities triggered by deviations from axisymmetry. We use a cosmological simulation of the formation of a disk galaxy to show that counterrotating stellar disk components may arise naturally in hierarchically-clustering scenarios even in the absence of merging. The simulated disk galaxy consists of two coplanar, overlapping stellar components with opposite spins: an inner counterrotating bar-like structure made up mostly of old stars surrounded by an extended, rotationally-supported disk of younger stars. The opposite-spin components originate from material accreted from two distinct filamentary structures which at turn around, when their net spin is acquired, intersect delineating a "V"-like structure. Each filament torques the other in opposite directions; the filament that first drains into the galaxy forms the inner counterrotating bar, while material accreted from the other filament forms the outer disk. Mergers do not play a substantial role and most stars in the galaxy are formed in situ; only 9% of all stars are contributed by accretion events. The formation scenario we describe here implies a significant age difference between the co- and counterrotating components, which may be used to discriminate between competing scenarios for the origin of counterrotating stars in disk galaxies., Comment: 7 pages, 7 figures. Accepted for publication in MNRAS

We use the very high resolution, fully cosmological simulations from the Aquarius project, coupled to a semi-analytical model of galaxy formation, to study the phase-space distribution of halo stars in "solar neighbourhood"-like volumes. We find that this distribution is very rich in substructure in the form of stellar streams for all five stellar haloes we have analysed. These streams can be easily identified in velocity space, as well as in spaces of pseudo-conserved quantities such as E vs. Lz. In our best-resolved local volumes, the number of identified streams ranges from ~ 300 to 600, in very good agreement with previous analytical predictions, even in the presence of chaotic mixing. The fraction of particles linked to (massive) stellar streams in these volumes can be as large as 84%. The number of identified streams is found to decrease as a power-law with galactocentric radius. We show that the strongest limitation to the quantification of substructure in our poorest-resolved local volumes is particle resolution rather than strong diffusion due to chaotic mixing., Comment: 13 pages, 10 figures. MNRAS accepted. Revised to reflect final version

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(\rho_{crit}(z)), resembles that of the enclosed density profile, M(<\rho >), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M(<\rho >) and M(\rho_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here., Comment: 12 pages, 8 figures, MNRAS 432 1103L (2013)

Current models of galaxy formation predict that galaxy pairs of comparable magnitudes should become increasingly rare with decreasing luminosity. This seems at odds with the relatively high frequency of pairings among dwarf galaxies in the Local Group. We use literature data to show that ~30% of all satellites of the Milky Way and Andromeda galaxies brighter than M_V=-8 are found in likely physical pairs of comparable luminosity. Besides the previously recognised pairings of the Magellanic Clouds and of NGC 147/NGC 185, other candidate pairs include the Ursa Minor and Draco dwarf spheroidals, as well as the And I/And III satellites of M31. These pairs are much closer than expected by chance if the radial and angular distributions of satellites were uncorrelated; in addition, they have very similar line-of-sight velocities and luminosities that differ by less than three magnitudes. In contrast, the same criteria pair fewer than 4% of satellites in N-body/semi-analytic models that match the radial distribution and luminosity function of Local Group satellites. If confirmed in studies of larger samples, the high frequency of dwarf galaxy pairings may provide interesting clues to the formation of faint galaxies in the current cosmological paradigm., Comment: 6 pages, 3 figures; updated version accepted for publication in MNRAS. minor changes, conclusions unchanged

We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the process. The dark matter content of these halos is unaffected but their gas can be efficiently removed by ram-pressure. The loss of gas is especially pronounced in low-mass halos due to their lower binding energy and has a dramatic effect on the star formation history of affected systems. This "cosmic web stripping" may help to explain the scarcity of dwarf galaxies compared with the numerous low-mass halos expected in \Lambda CDM and the large diversity of star formation histories and morphologies characteristic of faint galaxies. Although our results are based on a single high-resolution simulation, it is likely that the hydrodynamical interaction of dwarf galaxies with the cosmic web is a crucial ingredient so far missing from galaxy formation models., Comment: Submitted to ApJL. 6 pages, 4 figures. A set of movies showing the interaction between dwarf galaxies and the Cosmic Web can be found at mirror 1 http://www.astro.uvic.ca/~mario/dwarf-web/ or at mirror 2 http://www.iate.oac.uncor.edu/~alejandro/dwarf-web/ . Comments are welcomed

We study the abundance of satellite galaxies as a function of primary stellar mass using the SDSS/DR7 spectroscopic catalogue. In contrast with previous studies, which focussed mainly on bright primaries, our central galaxies span a wide range of stellar mass, 10^7.5 < M_*^pri/M_sun < 10^11, from dwarfs to central cluster galaxies. Our analysis confirms that the average number of satellites around bright primaries, when expressed in terms of satellite-to-primary stellar mass ratio (m_*^sat/M_*^pri), is a strong function of M_*^pri. On the other hand, satellite abundance is largely independent of primary mass for dwarf primaries (M_*^pri<10^10 M_sun). These results are consistent with galaxy formation models in the LCDM scenario. We find excellent agreement between SDSS data and semi-analytic mock galaxy catalogues constructed from the Millennium-II Simulation. Satellite galaxies trace dark matter substructure in LCDM, so satellite abundance reflects the dependence on halo mass, M_200, of both substructure and galaxy stellar mass (M_*). Since dark matter substructure is almost scale-free, the dependence of satellite abundance on primary mass results solely from the well-defined characteristic mass in the galaxy mass-halo mass relation. On dwarf galaxy scales, where models predict a power-law scaling, M_* \propto M_200^2.5, similarity is preserved and satellite abundance is independent of primary mass. For primaries brighter than the characteristic mass of the M_*-M_200 relation, satellite abundance increases strongly with primary mass. Our results provide strong support for the steep, approximately power-law dependence of dwarf galaxy mass on halo mass envisioned in LCDM galaxy formation models., Comment: 7 pages, 3 figures, updated version accepted for MNRAS main journal. Added extended discussion with respect to first submission. Conclusions unchanged

[Abbreviated] We have investigated the color-magnitude diagram of Omega Centauri and find that the blue main sequence (bMS) can be reproduced only by models that have a of helium abundance in the range Y=0.35-$0.40. To explain the faint subgiant branch of the reddest stars ("MS-a/RG-a" sequence), isochrones for the observed metallicity ([Fe/H]\approx0.7) appear to require both a high age (~13Gyr) and enhanced CNO abundances ([CNO/Fe]\approx0.9$). Y~0.35 must also be assumed in order to counteract the effects of high CNO on turnoff colors, and thereby to obtain a good fit to the relatively blue turnoff of this stellar population. This suggest a short chemical evolution period of time (<1Gyr) for Omega Cen. Our intermediate-mass (super-)AGB models are able to reproduce the high helium abundances, along with [N/Fe]~2 and substantial O depletions if uncertainties in the treatment of convection are fully taken into account. These abundance features distinguish the bMS stars from the dominant [Fe/H] $\approx1.7$ population. The most massive super-AGB stellar models (M_zams>=6.8M_sun, M_He,core>=1.245M_sun) predict too large N-enhancements, which limits their role in contributing to the extreme populations. We show quantitatively that highly He- and N-enriched AGB ejecta have particularly efficient cooling properties. Based on these results and on the reconstruction of the orbit of Omega Cen with respect to the Milky Way we propose the galactic plane passage gas purging scenario for the chemical evolution of this cluster. Our model addresses the formation and properties of the bMS population (including their central location in the cluster). We follow our model descriptively through four passage events, which could explain not only some key properties of the bMS, but also of the MS-a/RGB-a and the s-enriched stars., Comment: Accepted for publication in ApJ. 10 figures, 5 tables, 21 pages

We use the Millennium Simulation series to study how the dynamical state of dark matter halos affects the relation between mass and concentration. We find that a large fraction of massive systems are identified when they are substantially out of equilibrium and in a particular phase of their dynamical evolution: the more massive the halo, the more likely it is found at a transient stage of high concentration. This state reflects the recent assembly of massive halos and corresponds to the first pericentric passage of recently-accreted material when, before virialization, the kinetic and potential energies reach maximum and minimum values, respectively. This result explains the puzzling upturn in the mass-concentration relation reported in recent work for massive halos; indeed, the upturn disappears when only dynamically-relaxed systems are considered in the analysis. Our results warn against applying simple equilibrium models to describe the structure of rare, massive galaxy clusters and urges caution when extrapolating scaling laws calibrated on lower-mass systems, where such deviations from equilibrium are less common. The evolving dynamical state of galaxy clusters ought to be carefully taken into account if cluster studies are to provide precise cosmological constraints., Comment: 8 Pages. Minor changes to match published version

We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations., Comment: 20 pages, 16 figures. MNRAS in press. Minor changes in this version to match the version in press

Recent studies suggest that only three of the twelve brightest satellites of the Milky Way (MW) inhabit dark matter halos with maximum circular velocity, V_max, exceeding 30km/s. This is in apparent contradiction with the LCDM simulations of the Aquarius Project, which suggest that MW-sized halos should have at least 8 subhalos with V_max>30km/s. The absence of luminous satellites in such massive subhalos is thus puzzling and may present a challenge to the LCDM paradigm. We note, however, that the number of massive subhalos depends sensitively on the (poorly-known) virial mass of the Milky Way, and that their scarcity makes estimates of their abundance from a small simulation set like Aquarius uncertain. We use the Millennium Simulation series and the invariance of the scaled subhalo velocity function (i.e., the number of subhalos as a function of \nu, the ratio of subhalo V_max to host halo virial velocity, V_200) to secure improved estimates of the abundance of rare massive subsystems. In the range 0.1<\nu<0.5, N_sub(>\nu) is approximately Poisson-distributed about an average given by =10.2x(\nu/0.15)^(-3.11). This is slightly lower than in Aquarius halos, but consistent with recent results from the Phoenix Project. The probability that a LCDM halo has 3 or fewer subhalos with V_max above some threshold value, V_th, is then straightforward to compute. It decreases steeply both with decreasing V_th and with increasing halo mass. For V_th=30km/s, ~40% of M_halo=10^12 M_sun halos pass the test; fewer than 5% do so for M_halo>= 2x10^12 M_sun; and the probability effectively vanishes for M_halo>= 3x 10^12 M_sun. Rather than a failure of LCDM, the absence of massive subhalos might simply indicate that the Milky Way is less massive than is commonly thought., Comment: 5 figures, 8 pages. Accepted for publication in MNRAS. The number of massive satellites within a bigger radius r100 and in a WMAP7 cosmology model is considered now

In the simplest scenario, disk galaxies form predominantly in halos with high angular momentum and quiet recent assembly history, whereas spheroids are the slowly-rotating remnants of repeated merging events. We explore these assumptions using one hundred systems with halo masses similar to that of the Milky Way, identified in a series of cosmological gasdynamical simulations GIMIC. At z=0, the simulated galaxies exhibit a wide variety of morphologies, from dispersion-dominated spheroids to pure disk galaxies. Surprisingly, these morphological features are very poorly correlated with their halo properties: disks form in halos with high and low net spin, and mergers play a negligible role in the formation of spheroid stars, most of which form in-situ. More important to morphology is the coherent alignment of the angular momentum of baryons that accrete over time to form a galaxy. Spheroids tend to form when the spin of newly-accreted gas is misaligned with that of the extant galaxy, leading to the episodic formation of stars with different kinematics that cancel out the net rotation of the system. Disks, on the other hand, form out of gas that flows in with similar angular momentum to that of earlier-accreted material. Gas accretion from a hot corona thus favours disk formation, whereas gas that flows "cold", often along separate, misaligned filaments, favours the formation of spheroids. In this scenario, most spheroids consist of superpositions of stellar components with distinct kinematics, age, and metallicity, an arrangement that might survive to the present day given the paucity of major mergers. Since angular momentum is acquired largely at turnaround, morphology is imprinted early by the interplay of the tidal field and the shape of the material destined to form the galaxy., Comment: 13 pages, 10 figures. Extended discussion, conclusions unchanged. Accepted for publication in MNRAS

The cold dark matter halo mass function is much steeper than the galaxy stellar mass function on galactic and subgalactic scales. This difference is usually reconciled by assuming that the galaxy formation efficiency drops sharply with decreasing halo mass, so that virtually no dwarf galaxies form in halos less massive than ~ 10^10 M_sun. In turn, this implies that, at any given radius, the dark mass enclosed by a galaxy must exceed a certain minimum. We use rotation curves of dwarf galaxies compiled from the literature to explore whether their enclosed mass is consistent with these constraints. We find that almost one half of the dwarfs in our sample with stellar mass between 10^6-10^7 M_sun are at odds with this restriction: either they live in halos with masses substantially below 10^10 M_sun or there is a mechanism capable of reducing the dark mass enclosed by some of the faintest dwarfs. Neither possibility is easily accommodated within the standard LCDM scenario. Extending galaxy formation to halos well below 10^10 M_sun would lead to large numbers of dwarf galaxies in excess of current estimates; at the same time, the extremely low stellar mass of the systems involved makes it unlikely that baryonic effects can reduce their dark matter content. Resolving this challenge seems to require new insights into dwarf galaxy formation, or perhaps a radical revision of the prevailing paradigm., Comment: 8 pages, 4 Figures, MNRAS accepted

We report the discovery of eight lithium-rich field giants found in a high resolution spectroscopic sample of over 700 metal-poor stars ([Fe/H]<-0.5) selected from the RAVE survey. The majority of the Li-rich giants in our sample are very metal-poor ([Fe/H]<-1.9), and have a Li abundance (in the form of 7Li), A(Li)=log(n(Li)/n(H))+12, between 2.30 and 3.63, well above the typical upper red giant branch limit, A(Li)<0.5, while two stars, with A(Li)~1.7-1.8, show similar lithium abundances to normal giants at the same gravity. We further included two metal-poor, Li-rich globular cluster giants in our sample, namely the previously discovered M3-IV101 and newly discovered (in this work) M68-A96. This comprises the largest sample of metal-poor Li-rich giants to date. We performed a detailed abundance analysis of all stars, finding that the majority our sample stars have elemental abundances similar to that of Li-normal halo giants. Although the evolutionary phase of each Li-rich giant cannot be definitively determined, the Li-rich phase is likely connected to extra mixing at the red giant branch bump or early asymptotic giant branch that triggers cool bottom processing in which the bottom of the outer convective envelope is connected to the H-burning shell in the star. The surface of a star becomes Li-enhanced as 7Be (which burns to 7Li) is transported to the stellar surface via the Cameron-Fowler mechanism. We discuss and discriminate among several models for the extra mixing that can cause Li-production, given the detailed abundances of the Li-rich giants in our sample., Comment: 17 pages, 8 figures, accepted for publication in ApJ

We have undertaken the study of the elemental abundances and kinematic properties of a metal-poor sample of candidate thick-disk stars selected from the RAVE spectroscopic survey of bright stars to differentiate among the present scenarios of the formation of the thick disk. In this paper, we report on a sample of 214 red giant branch, 31 red clump/horizontal branch, and 74 main-sequence/sub-giant branch metal-poor stars, which serves to augment our previous sample of only giant stars. We find that the thick disk [alpha/Fe] ratios are enhanced, and have little variation (<0.1 dex), in agreement with our previous study. The augmented sample further allows, for the first time, investigation of the gradients in the metal-poor thick disk. For stars with [Fe/H] < -1.2, the thick disk shows very small gradients, <0.03 +/- 0.02 dex/kpc, in alpha-enhancement, while we find a +0.01 +/- 0.04 dex/kpc radial gradient and a -0.09 +/- 0.05 dex/kpc vertical gradient in iron abundance. In addition, we show that the peak of the distribution of orbital eccentricities for our sample agrees better with models in which the stars that comprise the thick disk were formed primarily in the Galaxy, with direct accretion of stars contributing little. Our results thus disfavor direct accretion of stars from dwarf galaxies into the thick disk as a major contributor to the thick disk population, but cannot discriminate between alternative models for the thick disk, such as those that invoke high-redshift (gas-rich) mergers, heating of a pre-existing thin stellar disk by a minor merger, or efficient radial migration of stars., Comment: 21 pages, 31 figures, and 5 tables. Accepted for publication in ApJ

We use cosmological simulations from the Aquarius Project to study the orbital history of the Large Magellanic Cloud (LMC) and its potential association with other satellites of the Milky Way (MW). We search for dynamical analogs to the LMC and find a subhalo that matches the LMC position and velocity at either of its two most recent pericentric passages. This suggests that the LMC is not necessarily on its first approach to the MW, provided that the virial mass of the Milky Way is M_200 ~1.8e12 M_sun. The simulation results yield specific predictions for the position and velocity of systems associated with the LMC prior to infall. If on first approach, most should lie close to the LMC because the Galactic tidal field has not yet had enough time to disperse them. If on second approach, the list of potential associates increases substantially. Interestingly, our analysis rules out an LMC association for Draco and Ursa Minor, two of the dwarf spheroidals suggested by Lynden-Bell & Lynden-Bell to form part of the "Magellanic Ghostly Stream". Our results also indicate that the direction of the orbital angular momentum is a powerful test of LMC association. This test, however, requires precise proper motions, which are unavailable for most MW satellites. Of the 4 satellites with published proper motions, only the Small Magellanic Cloud is clearly associated with the LMC. Taken at face value, the proper motions of Carina, Fornax and Sculptor rule them out as potential associates, but this conclusion should be revisited when better data become available. The dearth of satellites clearly associated with the Clouds might be solved by wide-field imaging surveys that target its surroundings, a region that may prove a fertile hunting ground for faint, previously unnoticed MW satellites., Comment: 13 pages, 8 figures, MNRAS in press. Minor corrections added to match published version

We use the high resolution cosmological N-body simulations from the Aquarius project to investigate in detail the mechanisms that determine the shape of Milky Way-type dark matter haloes. We find that, when measured at the instantaneous virial radius, the shape of individual haloes changes with time, evolving from a typically prolate configuration at early stages to a more triaxial/oblate geometry at the present day. This evolution in halo shape correlates well with the distribution of the infalling material: prolate configurations arise when haloes are fed through narrow filaments, which characterizes the early epochs of halo assembly, whereas triaxial/oblate configurations result as the accretion turns more isotropic at later times. Interestingly, at redshift z=0, clear imprints of the past history of each halo are recorded in their shapes at different radii, which also exhibit a variation from prolate in the inner regions to triaxial/oblate in the outskirts. Provided that the Aquarius haloes are fair representatives of Milky Way-like 10^12 Msun objects, we conclude that the shape of such dark matter haloes is a complex, time-dependent property, with each radial shell retaining memory of the conditions at the time of collapse., Comment: 16 pages, 14 figures. Accepted for publication in MNRAS. Minor changes to match published version

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations., Comment: The paper contains 23 pages, 13 figures. Accepted to MNRAS (2011, June 28)

The Monoceros ring is a collection of stars in nearly-circular orbits at roughly 18 kpc from the Galactic center. It may have originated (i) as the response of the disc to perturbations excited by satellite companions or (ii) from the tidal debris of a disrupted dwarf galaxy. The metallicity of Monoceros stars differs from that of disc stars at comparable Galactocentric distances, an observation that disfavours the first scenario. On the other hand, circular orbits are difficult to accommodate in the tidal-disruption scenario, since it requires a satellite which at the time of disruption was itself in a nearly circular orbit. Such satellite could not have formed at the location of the ring and, given its low mass, dynamical friction is unlikely to have played a major role in its orbital evolution. We search cosmological simulations for low-mass satellites in nearly-circular orbits and find that they result, almost invariably, from orbital changes induced by collisions with more massive satellites: the radius of the circular orbit thus traces the galactocentric distance of the collision. Interestingly, the Sagittarius dwarf, one of the most luminous satellites of the Milky Way, is in a polar orbit that crosses the Galactic plane at roughly the same Galactocentric distance as Monoceros. We use idealized simulations to demonstrate that an encounter with Sagittarius might well have led to the circularization and subsequent tidal demise of the progenitor of the Monoceros ring., Comment: 6 pages, 4 figures, to match version published in MNRAS Letters (http://onlinelibrary.wiley.com/doi/10.1111/j.1745-3933.2011.01035.x/abstract)

The Pan Andromeda Archeological Survey (PAndAS) CFHT Megaprime survey of the M31-M33 system has found a star stream which extends about 120 kpc NW from the center of M31. The great length of the stream, and the likelihood that it does not significantly intersect the disk of M31, means that it is unusually well suited for a measurement of stream gaps and clumps along its length as a test for the predicted thousands of dark matter sub-halos. The main result of this paper is that the density of the stream varies between zero and about three times the mean along its length on scales of 2 to 20 kpc. The probability that the variations are random fluctuations in the star density is less than 10^-5. As a control sample we search for density variations at precisely the same location in stars with metallicity higher than the stream, [Fe/H]=[0, -0.5] and find no variations above the expected shot noise. The lumpiness of the stream is not compatible with a low mass star stream in a smooth galactic potential, nor is it readily compatible with the disturbance caused by the visible M31 satellite galaxies. The stream's density variations appear to be consistent with the effects of a large population of steep mass function dark matter sub-halos, such as found in LCDM simulations, acting on an approximately 10Gyr old star stream. The effects of a single set of halo substructure realizations are shown for illustration, reserving a statistical comparison for another study., Comment: ApJ revised version submitted

Cosmological N-body simulations indicate that the spherically-averaged density profiles of cold dark matter halos are accurately described by Einasto profiles, where the logarithmic slope is a power-law of adjustable exponent, \gamma =dln\rho /dlnr ~ r^\alpha $. The pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable regularity, and are well approximated by simple power laws, Q(r)=\rho /\sigma ^3 ~ r^-\chi . We show that this is expected from dynamical equilibrium considerations, since Jeans' equations predict that the pseudo-phase-space density profiles of Einasto halos should resemble power laws over a wide range of radii. For the values of \alpha typical of CDM halos, the inner Q profiles of equilibrium halos deviate significantly from a power law only very close to the center, and simulations of extremely high-resolution would be needed to detect such deviations unambiguously. We use an ensemble of halos drawn from the Millennium-II simulation to study which of these two alternatives describe best the mass profile of CDM halos. Our analysis indicates that at the resolution of the best available simulations, both Einasto and power-law PPSD profiles (with adjustable exponents \alpha and \chi, respectively) provide equally acceptable fits to the simulations. A full account of the structure of CDM halos requires understanding how the shape parameters that characterize departures from self-similarity, like \alpha or \chi, are determined by evolutionary history, environment or initial conditions., Comment: 9 pages, 7 figures, submitted to MNRAS