Search

We study gas inflows onto supermassive black holes using hydrodynamics simulations of isolated galaxies and idealized galaxy mergers with an explicit, multiphase interstellar medium (ISM). Our simulations use the recently developed ISM and stellar evolution model called Stars and MUltiphase Gas in GaLaxiEs (SMUGGLE). We implement a novel super-Lagrangian refinement scheme that increases the gas mass resolution in the immediate neighborhood of the black holes (BHs) to accurately resolve gas accretion. We do not include black hole feedback in our simulations. We find that the complex and turbulent nature of the SMUGGLE ISM leads to highly variable BH accretion. BH growth in SMUGGLE converges at gas mass resolutions $\lesssim3\times10^3{\rm M_\odot}$. We show that the low resolution simulations combined with the super-Lagrangian refinement scheme are able to produce central gas dynamics and BH accretion rates very similar to that of the uniform high resolution simulations. We further explore BH fueling by simulating galaxy mergers. The interaction between the galaxies causes an inflow of gas towards the galactic centres and results in elevated and bursty star formation. The peak gas densities near the BHs increase by orders of magnitude resulting in enhanced accretion. Our results support the idea that galaxy mergers can trigger AGN activity, although the instantaneous accretion rate depends strongly on the local ISM. We also show that the level of merger-induced enhancement of BH fueling predicted by the SMUGGLE model is much smaller compared to the predictions by simulations using an effective equation of state model of the ISM., 16 pages, 10 figures

Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-alpha (Ly{\alpha}), and Balmer line (H{\alpha} and H\b{eta}) radiative transfer study of a high-resolution isolated Milky-Way simulation using the Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry, and dust grain evolution in the interstellar medium (ISM). We extend our Cosmic Ly{\alpha} Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust (27.5%), helium ionization (8.7%), and anisotropic escape fractions (7.9%), as these reduce the available budget for hydrogen line emission (55.9%). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas phase structure, and spatial, spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local H{\alpha} surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Ly{\alpha} emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve, Comment: 28 pages, 32 figures, MNRAS, in press

The arrangement of dwarf galaxies in a thin plane surrounding the Milky Way has been thought to contradict the prevailing cosmological model of cold dark matter in the Universe. New work suggests that this arrangement may just be a temporary alignment, bringing our galaxy back into agreement with theoretical expectations once the radial distribution of satellites is taken into account., Comment: Invited 'News and Views' contribution for Nature Astronomy on paper by Sawala et al. 2022 "The Milky Way's plane of satellites is consistent with {\Lambda}CDM"

Ultra-diffuse galaxies (UDGs) are the lowest-surface brightness galaxies known, with typical stellar masses of dwarf galaxies but sizes similar to larger galaxies like the Milky Way. The reason for their extended sizes is debated, with suggested internal processes like angular momentum, feedback or mergers versus external mechanisms or a combination of both. Observationally, we know that UDGs are red and quiescent in groups and clusters while their counterparts in the field are blue and star-forming. This dichotomy suggests environmental effects as main culprit. However, this scenario is challenged by recent observations of isolated quiescent UDGs in the field. Here we use $\Lambda$CDM cosmological hydrodynamical simulation to show that isolated quenched UDGs are formed as backsplash galaxies that were once satellites of another galactic, group or cluster halo but are today a few Mpc away from them. These interactions, albeit brief, remove the gas and tidally strip the outskirts of the dark matter haloes of the now quenched seemingly-isolated UDGs, which are born as star-forming field UDGs occupying dwarf-mass dark matter haloes. Quiescent UDGs may therefore be found in non-negligible numbers in filaments and voids, bearing the mark of past interactions as stripped outer haloes devoid of dark matter and gas compared to dwarfs with similar stellar content., Comment: Author's version for the main article and supplementary information (13 pages, 10 figures). Accepted for publication in Nature Astronomy

We use the IllustrisTNG cosmological hydrodynamical simulation to study the formation of ultra-diffuse galaxies (UDGs) in galaxy clusters. We supplement the simulations with a realistic mass-size relation for galaxies at the time of infall into the cluster, as well as an analytical model to describe the tidally-induced evolution of their stellar mass, velocity dispersion and size. The model assumes "cuspy" NFW halos and, contrary to recent claims, has no difficulty reproducing the observed number of UDGs in clusters. Our results further suggest that the UDG population consists of a mixture of "normal" low surface brightness galaxies such as those found in the field ("born" UDGs, or B-UDGs), as well as a distinct population that owe their large size and low surface brightness to the effects of cluster tides ("tidal", or T-UDGs). The simulations indicate that T-UDGs entered the cluster earlier and should be more prevalent than B-UDGs near the cluster centres. T-UDGs should also have, at given stellar mass, lower velocity dispersion, higher metallicities, and lower dark matter content than B-UDGs. Our results suggest that systems like DF-44 are consistent with having been born as UDGs, while others such as DF2, DF4 and VLSB-D are possibly extreme T-UDG examples., Comment: 12 pages, 9 figures. Accepted version for publication in MNRAS, March 23rd 2020

Galaxies and their dark-matter halos have posed several challenges to the Dark Energy plus Cold Dark Matter (LCDM) cosmological model. These discrepancies between observations and theory intensify for the lowest-mass (`dwarf') galaxies. LCDM predictions for the number, spatial distribution, and internal structure of low-mass dark-matter halos have historically been at odds with observed dwarf galaxies, but this is partially expected, because many predictions modeled only the dark-matter component. Any robust LCDM prediction must include, hand-in-hand, a model for galaxy formation to understand how baryonic matter populates and affects dark-matter halos. In this article, we review the most notable challenges to LCDM regarding dwarf galaxies, and we discuss how recent cosmological numerical simulations have pinpointed baryonic solutions to these challenges. We identify remaining tensions, including the diversity of the inner dark-matter content, planes of satellites, stellar morphologies, and star-formation quenching. Their resolution, or validation as actual problems to LCDM, will likely require both refining galaxy formation models and improving numerical accuracy in simulations., Comment: Invited Review for Nature Astronomy, 25 pages (12 main article), 5 figures. Accepted version

Impulsive supernova feedback and non-standard dark matter models, such as self-interacting dark matter (SIDM), are the two main contenders for the role of the dominant core formation mechanism at the dwarf galaxy scale. Here we show that the impulsive supernova cycles that follow episodes of bursty star formation leave distinct features in the distribution function of stars: groups of stars with similar ages and metallicities develop overdense shells in phase space. If cores are formed through supernova feedback, we predict the presence of such features in star-forming dwarf galaxies with cored host halos. Their systematic absence would favor alternative dark matter models, such as SIDM, as the dominant core formation mechanism., Comment: 16 pages, 13 figures, arxiv version contains both the letter and the supplemental material \c{opyright} 2022 American Physical Society

We present a post-processing catalog of globular clusters (GCs) for the $39$ most massive groups and clusters in the TNG50 simulation of the IlllustrisTNG project (virial masses $M_{200} =[5\times 10^{12} \rm - 2 \times 10^{14}$] M$_{\odot}$). We tag GC particles to all galaxies with stellar mass $M_* \geq 5\times10^6$ M$_{\odot}$, and we calibrate their masses to reproduce the observed power-law relation between GC mass and halo mass for galaxies with $M_{200} \geq 10^{11}$ M$_{\odot}$ (corresponding to $M_* \sim 10^9$ $M_{\odot}$). Here we explore whether an extrapolation of this $M_{\rm GC}$-$M_{200}$ relation to lower-mass dwarfs is consistent with current observations. We find a good agreement between our predicted number and specific frequency of GCs in dwarfs with $\rm M_*=[5 \times 10^6 \rm - 10^9]$ M$_{\odot}$ and observations. Moreover, we predict a steep decline in the GC occupation fraction for dwarfs with $M_*, 16 pages, 8 figures (main body of paper). Resubmitted to MNRAS after positive referee report. The globular cluster catalogs for the 39 most massive groups and clusters in TNG50 are available for download from the IllustrisTNG public data release at www.tng-project.org/doppel22

Within $��$CDM, dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhalos, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC-system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$ > $10^5$ M$_{\odot}$; Carina and Fornax) as well as several ultra-faint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, $M_\ast$ ~ $6\times10^3$ M$_{\odot}$) substructure population of isolated LMC-mass hosts ($M_\text{200}$ = $1-3\times10^{11}$ M$_{\odot}$) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhalos in simulations with and without baryons, we find that, within 0.2 $r_\text{200}$, LMC-mass hosts deplete ~30% of their substructure, significantly lower than the ~70% of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$ = 880 M$_{\odot}$), ~5-10 subhalos form galaxies with $M_\ast$ > $10^4$ M$_{\odot}$, in agreement with the 7 observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faint end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultra faint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association., 19 pages, 9 figures, 4 tables. Accepted for publication in MNRAS

The abundance of the faintest galaxies provides insight into the nature of dark matter and the process of dwarf galaxy formation. In the LCDM scenario, low mass halos are so numerous that the efficiency of dwarf formation must decline sharply with decreasing halo mass in order to accommodate the relative scarcity of observed dwarfs and satellites in the Local Group. The nature of this decline contains important clues to the mechanisms regulating the onset of galaxy formation in the faintest systems. We explore here two possible models for the stellar mass ($M_*$)-halo mass ($M_{200}$) relation at the faint end, motivated by some of the latest LCDM cosmological hydrodynamical simulations. One model includes a sharp mass threshold below which no luminous galaxies form, as expected if galaxy formation proceeds only in systems above the Hydrogen-cooling limit. In the second model, $M_*$ scales as a steep power-law of $M_{200}$ with no explicit cutoff, as suggested by recent semianalytic work. Although both models predict satellite numbers around Milky Way-like galaxies consistent with current observations, they predict vastly different numbers of ultra-faint dwarfs and of satellites around isolated dwarf galaxies. Our results illustrate how the satellite mass function around dwarfs may be used to probe the $M_*$-$M_{200}$ relation at the faint end and to elucidate the mechanisms that determine which low-mass halos "light up" or remain dark in the LCDM scenario., Comment: 14 pages, 11 figures. Accepted version for publication in MNRAS

We present a suite of 16 high-resolution hydrodynamic simulations of an isolated dwarf galaxy (gaseous and stellar disk plus a stellar bulge) within an initially cuspy dark matter (DM) halo, including self-interactions between the DM particles (SIDM); as well as stochastic star formation and subsequent supernova feedback (SNF), implemented using the stellar feedback model SMUGGLE. The simulations start from identical initial conditions and we regulate the strength of SIDM and SNF by systematically varying the SIDM momentum transfer cross section and the gas density threshold for star formation. The DM halo forms a constant density core of similar size and shape for several combinations of those two parameters. Haloes with cores that are formed due to SIDM (adiabatic cusp-core transformation) have velocity dispersion profiles which are closer to isothermal than those of haloes with cores that are formed due to SNF in simulations with bursty star formation (impulsive cusp-core transformation). Impulsive SNF can generate positive stellar age gradients and increase random motion in the gas at the centre of the galaxy. Simulated galaxies in haloes with cores that were formed adiabatically are spatially more extended, with stellar metallicity gradients that are shallower (at late times) than those of galaxies in other simulations. Such observable properties of the gas and the stars, which indicate either an adiabatic or an impulsive evolution of the gravitational potential, may be used to determine whether observed cores in DM haloes are formed through self-interactions between the DM particles or in response to impulsive SNF., Comment: accepted for publication in MNRAS, contains 3 appendices; main paper: 18 pages, 11 figures; appendices: 8 pages, 7 figures; for questions/comments contact jburger@ucr.edu (Jan Burger)

While it is generally believed that supermassive black holes (SMBH) lie in most galaxies with bulges, few SMBHs have been confirmed in bulgeless galaxies. Identifying such a population could provide important insights to the BH seed population and secular BH growth. To this end, we obtained near-infrared spectroscopic observations of a sample of low-redshift bulgeless galaxies with mid-infrared colors suggestive of AGN. We find additional evidence of AGN activity (such as coronal lines and broad permitted lines) in 69$\%$ (9/13) of the sample, demonstrating that mid-infrared selection is a powerful tool to detect AGN. More than half of the galaxies with confirmed AGN activity show fast outflows in [O III] in the optical and/or [Si VI] in the NIR, with the latter generally having much faster velocities that are also correlated to their spatial extent. We are also able to obtain virial BH masses for some targets and find they fall within the scatter of other late-type galaxies in the $M_{\rm{BH}}$-$M_{\rm{stellar}}$ relation. The fact that they lack a significant bulge component indicates that secular processes, likely independent of major mergers, grew these BHs to supermassive sizes. Finally, we analyze the rotational gas kinematics and find two notable exceptions: two AGN hosts with outflows that appear to be rotating faster than expected. There is an indication that these two galaxies have stellar masses significantly lower than expected from their dark matter halo masses. This, combined with the observed AGN activity and strong gas outflows may be evidence of the effects of AGN feedback., Accepted by ApJ. 25 pages, 8 figures, and 6 tables

We study evolution of self-interacting dark matter (SIDM) subhalos in the Milky Way (MW) tidal field. The interaction between the subhalos and the MW's tides lead to more diverse dark matter distribution in the inner region, compared to their cold dark matter counterparts. We test this scenario with two MW satellite galaxies, Draco and Fornax, opposite extremes in the inner dark matter content, and find that they can be accommodated within the SIDM model proposed to explain the diverse rotation curves of spiral galaxies in the field., Comment: 6 pages, 3 figures. Updated figures and text. Accepted for publication in PRL

Globular clusters (GCs) are often used to estimate the dark matter content of galaxies, especially dwarfs, where other kinematic tracers are lacking. These estimates typically assume spherical symmetry and dynamical equilibrium, assumptions that may not hold for the sparse GC population of dwarfs in galaxy clusters. We use a catalog of GCs tagged onto the Illustris simulation to study the accuracy of GC-based mass estimates. We focus on galaxies in the stellar mass range 10$^{8} - 10^{11.8}$ M$_{\odot}$ identified in $9$ simulated Virgo-like clusters. Our results indicate that mass estimates are, on average, quite accurate in systems with GC numbers $N_{\rm GC} \geq 10$ and where the uncertainty of individual GC line-of-sight velocities is smaller than the inferred velocity dispersion, $\sigma_{\rm GC}$. In cases where $N_{\rm GC} \leq 10$, however, biases may result depending on how $\sigma_{\rm GC}$ is computed. We provide calibrations that may help alleviate these biases in methods widely used in the literature. As an application, we find a number of dwarfs with $M_{*} \sim 10^{8.5}\, M_{\odot}$ (comparable to the ultradiffuse galaxy DF2, notable for the low $\sigma_{GC}$ of its $10$ GCs) with $\sigma_{\rm GC} \sim 7$ - $15\; \rm km \rm s^{-1}$. These DF2 analogs correspond to relatively massive systems at their infall time ($M_{200} \sim 1$ - $3 \times 10^{11}$ $M_{\odot}$) which have retained only $3$-$17$ GCs and have been stripped of more than 95$\%$ of their dark matter. Our results suggest that extreme tidal mass loss in otherwise normal dwarf galaxies may be a possible formation channel for ultradiffuse objects like DF2., Comment: 19 pages, 15 figures. Accepted to MNRAS Dec. 11 2020

Recent studies using Gaia DR2 have identified a massive merger in the early history of the Milky Way (MW) whose debris is dominated by radial and counterrotating orbits. This event, dubbed the Gaia-Enceladus/Gaia-Sausage (GE/GS), is also hypothesized to have built the majority of the inner stellar halo. We use the cosmological hydrodynamic simulation Illustris to place this merger in the context of galaxy assembly within lambda cold dark matter. From ∼150 MW analogues, $\sim \!80 {{ \rm {per\ cent}}}$ have experienced at least one merger of similar mass and infall time as the GE/GS. Within this sample, 37 have debris as radial as the GE/GS, which we dub the ancient radial mergers (ARMs). Counterrotation is not rare among ARMs, with $43 {{ \rm {per\ cent}}}$ having $\gt 40 {{ \rm {per\ cent}}}$ of their debris in counterrotating orbits. However, the compactness inferred for the GE/GS debris given its large radial orbital anisotropy, β, and its substantial contribution to the stellar halo are difficult to reproduce. The median radius of ARM debris is r*,deb ≃ 45 kpc, while GE/GS is thought to be mostly contained within r ∼ 30 kpc. For most MW analogues, few mergers are required to build the inner stellar halo, and ARM debris only accounts for (median) $\sim \!12 {{ \rm {per\ cent}}}$ of inner accreted stars. Encouragingly, we find one ARM that is both compact and dominates the inner halo of its central, making it our best GE/GS analogue. Interestingly, this merger deposits a significant number of stars (M* ≃ 1.5 × 109 M⊙) in the outer halo, suggesting that an undiscovered section of GE/GS may await detection.

We use deep Hubble Space Telescope imaging to derive a distance to the Virgo Cluster ultradiffuse galaxy (UDG) VCC 615 using the tip of the red giant branch (TRGB) distance estimator. We detect 5,023 stars within the galaxy, down to a 50% completeness limit of F814W = 28.0, using counts in the surrounding field to correct for contamination due to background sources and Virgo intracluster stars. We derive an extinction-corrected F814W tip magnitude of $m_{\rm tip,0} = 27.19^{+0.07}_{-0.05}$, yielding a distance of $d=17.7^{+0.6}_{-0.4}$ Mpc. This places VCC 615 on the far side of the Virgo Cluster ($d_{\rm Virgo} = 16.5 Mpc$), at a Virgocentric distance of 1.3 Mpc and near the virial radius of the main body of Virgo. Coupling this distance with the galaxy's observed radial velocity, we find that VCC 615 is on an outbound trajectory, having survived a recent passage through the inner parts of the cluster. Indeed, our orbit modeling gives a 50% chance the galaxy passed inside the Virgo core (r, 18 pages, 10 figures, accepted for publication in The Astrophysical Journal

We present a summary of the predictions from numerical simulations to our understanding of dwarf galaxies. It centers the discussion around the Λ Cold Dark Matter scenario (ΛCDM) but discusses also implications for alternative dark matter models. Four key predictions are identified: the abundance of dwarf galaxies, their dark matter content, their relation with environment and the existence of dwarf satellites orbiting dwarf field galaxies. We discuss tensions with observations and identify the most exciting predictions expected from simulations in the future, including i) the existence of “dark galaxies” (dark matter halos without stars), ii) the ability to resolve the structure (size, morphology, dark matter distribution) in dwarfs and iii) the number of ultra-faint satellites around dwarf galaxies. All of these predictions shall inform future observations, not only the faintest galaxies to be discovered within the Local Volume but also distant dwarfs driving galaxy formation in the early universe.

We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift = 1.3, on timescales that depend sen- sitively on the strength of the pattern. Strong bars develop relatively quickly (in a few Gyr, = 10 disc rotation periods) in systems that are disc dominated, gas poor, and have declining rotation curves. Weak bars develop more slowly in systems where the disc is less gravitation- ally important, and are still growing at z = 0. Unbarred galaxies are comparatively gas-rich discs whose rotation speeds do not exceed the maximum circular velocity of the halos they inhabit. Bar lengths compare favourably with observations, ranging from 0.2 to 0.8 times the radius containing 90% of the stars. Bars slow down remarkably quickly as they grow, causing the inner regions of the surrounding dark halo to expand. At z = 0 strong bars have corota- tion radii roughly ten times the bar length. Such slow bars are inconsistent with the few cases where pattern speeds have been measured or inferred observationally, a discrepancy that, if confirmed, might prove a challenge for disc galaxy formation in LCDM., 11 pages, 10 figures

Galactic rotation curves are a fundamental constraint for any cosmological model. We use controlled N-body simulations of galaxies to study the gravitational effect of baryons in a scenario with collisionless cold dark matter (CDM) versus one with a self-interacting dark matter (SIDM) component. In particular, we examine the inner profiles of the rotation curves in the velocity range Vmax = 30-250 km/s, whose diversity has been found to be greater than predicted by the Lambda-CDM scenario. We find that the scatter in the observed rotation curves exceeds that predicted by dark matter only mass-concentration relations in either the CDM nor SIDM models. Allowing for realistic baryonic content and spatial distributions, however, helps create a large variety of rotation curve shapes; which is in better agreement with observations in the case of self-interactions due to the characteristic cored profiles being more accommodating to the slowly rising rotation curves than CDM. We find individual fits to model two of the most remarkable outliers of similar Vmax, UGC 5721 and IC 2574, the former a cusp-like rotation curve and the latter a seemingly 8 kpc cored profile. This diversity in SIDM arises as permutations of overly concentrated halos with compact baryonic distributions versus underdense halos with extended baryonic disks. The SIDM solution is promising and its feasibility ultimately depends on the sampling of the halo mass-concentration relation and its interplay with the baryonic profiles, emphasising the need for a better understanding of the frequency of extreme outliers present in current observational samples., Comment: 13 pages, 8 figures, accepted by MNRAS

We have used dedicated 0.7m telescopes in California and Israel to image the halos of ~200 galaxies in the Local Volume to 29 mag/sq arcsec, the sample mainly drawn from the 2MASS Large Galaxy Atlas (LGA). We supplement the LGA sample with dwarf galaxies and more distant giant ellipticals. Low surface brightness halos exceeding 50 kpc in diameter are found only in galaxies more luminous than L* and classic interaction signatures are relatively infrequent. Halo diameter is correlated with total galaxy luminosity. Extended low surface brightness halos are present even in galaxies as faint as M_V=-18. Edge-on galaxies with boxy bulges tend to lack extended spheroidal halos, while those with large classical bulges exhibit extended round halos, supporting the notions that boxy or barlike bulges originate from disks. Most face-on spiral galaxies present features that appear to be irregular extensions of spiral arms, although rare cases show smooth boundaries with no sign of star formation. Although we serendipitously discovered a dwarf galaxy undergoing tidal disruption in the halo of NGC 4449, we found no comparable examples in our general survey. A search for similar examples in the Local Volume identified hcc087, a tidally disrupting dwarf galaxy in the Hercules Cluster, but we do not confirm an anomalously large half-light radius reported for the dwarf VCC 1661., Comment: 4 pages, 4 figures

Observationally, the fraction of blue satellite galaxies decreases steeply with host halo mass, and their radial distribution around central galaxies is significantly shallower in massive (M_* >10e11M_sun) than in Milky Way like systems. Theoretical models, based primarily on semi-analytical techniques, have had a long-standing problem with reproducing these trends, instead predicting too few blue satellites in general but also estimating a radial distribution that is too shallow, regardless of primary mass. In this Letter, we use the Illustris cosmological simulation to study the properties of satellite galaxies around isolated primaries. For the first time, we find good agreement between theory and observations. We identify the main source of this success relative to earlier work to be a consequence of the large gas contents of satellites at infall, a factor ~5-10 times larger than in semi-analytical models. Because of their relatively large gas reservoirs, satellites can continue to form stars long after infall, with a typical timescale for star-formation to be quenched ~2 Gyr in groups but more than ~5 Gyr for satellites around Milky Way like primaries. The gas contents we infer are consistent with z=0 observations of HI gas in galaxies, although we find large discrepancies among reported values in the literature. A testable prediction of our model is that the gas-to-stellar mass ratio of satellite progenitors should vary only weakly with cosmic time., 5 pages, 3 figures. Accepted for publication in MNRAS Letters

We explore a scenario where metal poor globular clusters (GCs) form at the centres of their own dark matter halos in the early universe before reionization. This hypothesis leads to predictions about the abundance, distribution and kinematics of GCs today that we explore using cosmological N-body simulations and analytical modelling. We find that selecting the massive tail of collapsed objects at $z\gtrsim 9$ as GC formation sites leads to four main predictions: i) a highly clustered population of GCs around galaxies today, ii) a natural scaling between number of GCs and halo virial mass that follows roughly the observed trend, iii) a very low number of free floating GCs outside massive halos and iv) GCs should be embedded within massive and extended dark matter (sub)halos. We find that the strongest constraint to the model is given by the combination of (i) and (ii): a mass cut to tagged GCs halos which accounts for the number density of metal poor GCs today predicts a radial distribution that is too extended compared to recent observations. On the other hand, a mass cut sufficient to match the observed half number radius could only explain 60% of the metal poor population. In all cases, observations favour early redshifts for GC formation ($z\geq 15$)., Comment: Accepted in MNRAS. 13 pages, 9 figures

We employ isolated N-body simulations to study the response of self-interacting dark matter (SIDM) haloes in the presence of the baryonic potentials. Dark matter self-interactions lead to kinematic thermalization in the inner halo, resulting in a tight correlation between the dark matter and baryon distributions. A deep baryonic potential shortens the phase of SIDM core expansion and triggers core contraction. This effect can be further enhanced by a large selfscattering cross-section. We find the final SIDM density profile is sensitive to the baryonic concentration and the strength of dark matter self-interactions. Assuming a spherical initial halo, we also study evolution of the SIDM halo shape together with the density profile. The halo shape at later epochs deviates from spherical symmetry due to the influence of the non-spherical disc potential, and its significance depends on the baryonic contribution to the total gravitational potential, relative to the dark matter one. In addition, we construct a multicomponent model for the MilkyWay, including an SIDM halo, a stellar disc, and a bulge, and show it is consistent with observations from stellar kinematics and streams. ©2018 The Author(s). Published by Oxford University Press on behalf of The Royal Astronomical Society., NASA MUREP Institutional Research Opportunity (grant no. NNX15AP99A ), Hubble Space Telescope (grant no. HST-AR-14582), U. S. Department of Energy (grant no. de-sc0008541), Hellman Fellows Fund, MIT RSC award, NASA ATP (grant NNX17AG29G), Grant of Excellence Iceland Research Fund (grant no. 173929-051)

We study how optical galaxy morphology depends on mass and star formation rate (SFR) in the Illustris Simulation. To do so, we measure automated galaxy structures in 10808 simulated galaxies at z=0 with stellar masses 10^9.7 < M_*/M_sun < 10^12.3. We add observational realism to idealized synthetic images and measure non-parametric statistics in rest-frame optical and near-IR images from four directions. We find that Illustris creates a morphologically diverse galaxy population, occupying the observed bulge strength locus and reproducing median morphology trends versus stellar mass, SFR, and compactness. Morphology correlates realistically with rotation, following classification schemes put forth by kinematic surveys. Type fractions as a function of environment agree roughly with data. These results imply that connections among mass, star formation, and galaxy structure arise naturally from models matching global star formation and halo occupation functions when simulated with accurate methods. This raises a question of how to construct experiments on galaxy surveys to better distinguish between models. We predict that at fixed halo mass near 10^12 M_sun, disc-dominated galaxies have higher stellar mass than bulge-dominated ones, a possible consequence of the Illustris feedback model. While Illustris galaxies at M_* ~ 10^11 M_sun have a reasonable size distribution, those at M_* ~ 10^10 M_sun have half-light radii larger than observed by a factor of two. Furthermore, at M_* ~ 10^10.5-10^11 M_sun, a relevant fraction of Illustris galaxies have distinct "ring-like" features, such that the bright pixels have an unusually wide spatial extent., 26 pages, 15 figures, MNRAS accepted version

We use the SDSS/DR8 galaxy sample to study the radial distribution of satellite galaxies around isolated primaries, comparing to semi-analytic models of galaxy formation based on the Millennium and Millennium-II simulations. SDSS satellites behave differently around high- and low-mass primaries: those orbiting objects with $M_*>10^{11}M_\odot$ are mostly red and are less concentrated towards their host than the inferred dark matter halo, an effect that is very pronounced for the few blue satellites. On the other hand, less massive primaries have steeper satellite profiles that agree quite well with the expected dark matter distribution and are dominated by blue satellites, even in the inner regions where strong environmental effects are expected. In fact, such effects appear to be strong only for primaries with $M_* > 10^{11}M_\odot$. This behaviour is not reproduced by current semi-analytic simulations, where satellite profiles always parallel those of the dark matter and satellite populations are predominantly red for primaries of all masses. The disagreement with SDSS suggests that environmental effects are too efficient in the models. Modifying the treatment of environmental and star formation processes can substantially increase the fraction of blue satellites, but their radial distribution remains significantly shallower than observed. It seems that most satellites of low-mass primaries can continue to form stars even after orbiting within their joint halo for 5 Gyr or more., 19 pages, 11 figures, accepted by MNRAS

The relative impact of radiation pressure and photoionization feedback from young stars on surrounding gas is studied with hydrodynamic radiative transfer (RT) simulations. The calculations focus on the single-scattering (direct radiation pressure) and optically thick regime, and adopt a moment-based RT-method implemented in the moving-mesh code AREPO. The source luminosity, gas density profile and initial temperature are varied. At typical temperatures and densities of molecular clouds, radiation pressure drives velocities of order ~20 km/s over 1-5 Myr; enough to unbind the smaller clouds. However, these estimates ignore the effects of photoionization that naturally occur concurrently. When radiation pressure and photoionization act together, the latter is substantially more efficient, inducing velocities comparable to the sound speed of the hot ionized medium (10-15 km/s) on timescales far shorter than required for accumulating similar momentum with radiation pressure. This mismatch allows photoionization to dominate the feedback as the heating and expansion of gas lowers the central densities, further diminishing the impact of radiation pressure. Our results indicate that a proper treatment of the impact of young stars on the interstellar medium needs to primarily account for their ionization power whereas direct radiation pressure appears to be a secondary effect. This conclusion may change if extreme boosts of the radiation pressure by photon trapping are assumed., Comment: 18 pages, 19 figures (main results presented in 13 pages, 10 figures; extended appendix for RT tests with extra 9 figures). Accepted for publication in MNRAS after tiny changes

In the current concordance cosmology small halos are expected to be completely dark and can significantly perturb low-mass galaxies during minor merger interactions. These interactions may well contribute to the diversity of the dwarf galaxy population. Dwarf galaxies in the field are often observed to have peculiarities in their structure, morphology, and kinematics, as well as strong bursts of star formation without apparent cause. We aim to characterize the signatures of minor mergers of dwarf galaxies with dark satellites to aid their observational identification. We explore and quantify a variety of structural, morphological, and kinematic indicators of merging dwarf galaxies and their remnants using a suite of hydrodynamical simulations. The most sensitive indicators of mergers with dark satellites are large asymmetries in the gaseous and stellar distributions, enhanced central surface brightness and starbursts, and velocity offsets and misalignments between the cold gas and stellar components. In general, merging systems span a wide range of values of the most commonly used indicators, while isolated objects tend to have more confined values. Interestingly, we find in our simulations that a significantly off-centered burst of star formation can pinpoint the location of the dark satellite. Observational systems with such characteristics are perhaps the most promising for unveiling the presence of the hitherto, missing satellites., Comment: 9 pages, 7 figures. Accepted in A&A

The hierarchical nature of LCDM suggests that the Magellanic Clouds must have been surrounded by a number of satellites before their infall into the Milky Way. Many of those satellites should still be in close proximity to the Clouds, but some could have dispersed ahead/behind the Clouds along their Galactic orbit. Either way, prior association with the Clouds results in strong restrictions on the present-day positions and velocities of candidate Magellanic satellites: they must lie close to the nearly-polar orbital plane of the Magellanic stream, and their distances and radial velocities must follow the latitude dependence expected for a tidal stream with the Clouds at pericenter. We use a cosmological numerical simulation of the disruption of a massive subhalo in a Milky Way-sized LCDM halo to test whether any of the 20 dwarfs recently-discovered in the DES, SMASH, Pan-STARRS, and ATLAS surveys are truly associated with the Clouds. Of the 6 systems with kinematic data, only Hydra II and Hor 1 have distances and radial velocities consistent with a Magellanic origin. Of the remaining dwarfs, six (Hor 2, Eri 3, Ret 3, Tuc 4, Tuc 5, and Phx 2) have positions and distances consistent with a Magellanic origin, but kinematic data are needed to substantiate that possibility. Conclusive evidence for association would require proper motions to constrain the orbital angular momentum direction, which, for true Magellanic satellites, must coincide with that of the Clouds. We use this result to predict radial velocities and proper motions for all new dwarfs. Our results are relatively insensitive to the assumption of first or second pericenter for the Clouds., Comment: 10 pages, 6 figures, 3 tables. Submitted to MNRAS

Mergers and the spin of the dark matter halo are factors traditionally believed to determine the morphology of galaxies within a $\Lambda$CDM cosmology. We study this hypothesis by considering approximately 18,000 central galaxies at $z=0$ with stellar masses $M_{\ast} = 10^{9}-10^{12} \, {\rm M}_{\odot}$ selected from the Illustris cosmological hydrodynamic simulation. The fraction of accreted stars -- which measures the importance of massive, recent and dry mergers -- increases steeply with galaxy stellar mass, from less than 5 per cent in dwarfs to 80 per cent in the most massive objects, and the impact of mergers on galaxy morphology increases accordingly. For galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, mergers have the expected effect: if gas-poor they promote the formation of spheroidal galaxies, whereas gas-rich mergers favour the formation and survivability of massive discs. This trend, however, breaks at lower masses. For objects with $M_{\ast} \lesssim 10^{11} \, {\rm M}_{\odot}$, mergers do not seem to play any significant role in determining the morphology, with accreted stellar fractions and mean merger gas fractions that are indistinguishable between spheroidal and disc-dominated galaxies. On the other hand, halo spin correlates with morphology primarily in the least massive objects in the sample ($M_{\ast} \lesssim 10^{10} \, {\rm M}_{\odot}$), but only weakly for galaxies above that mass. Our results support a scenario where (1) mergers play a dominant role in shaping the morphology of massive galaxies, (2) halo spin is important for the morphology of dwarfs, and (3) the morphology of medium-sized galaxies -- including the Milky Way -- shows little dependence on galaxy assembly history or halo spin, at least when these two factors are considered individually., Comment: 16 pages, 14 figures, 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 haloes less massive than ∼1010 M⊙. 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 in the range of 106 < Mgal/M⊙ < 107 are at odds with this restriction: either they live in haloes with masses substantially below 1010 M⊙ 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 Λ cold dark matter scenario. Extending galaxy formation to haloes well below 1010 M⊙ 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.

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.

Semi-analytic models (SAMs) are currently one of the primary tools with which to model statistically significant ensembles of galaxies. The underlying physical prescriptions inherent to each SAM are, in many cases, different from one another. Several SAMs have been applied to the dark matter merger trees extracted from the Millennium Run, including those associated with the well-known Munich and Durham lineages. We compare the predicted luminosity distributions of galaxy groups using four publicly available SAMs (De Lucia et al. 2006; Bower et al. 2006; Bertone et al. 2007; Font et al. 2008), in order to explore a galactic environment in which the models have not been explored to the same degree as they have in the field or in rich clusters. We identify a characteristic "wiggle" in the group galaxy luminosity function generated using the De Lucia et al. (2006) SAM, that is not present in the Durham-based models, consistent to some degree with observations. However, a comparison between conditional luminosity functions of groups between the models and observations of Yang et al. (2007) suggest that neither model is a particularly good match. The luminosity function wiggle is interpreted as the result of the two-mode AGN feedback implementation used in the Munich models, which itself results in flattened magnitude gap distribution. An associated analysis of the magnitude gap distribution between first- and second-ranked group galaxies shows that while the Durham models yield distributions with approximately equal luminosity first- and second-ranked galaxies, in agreement with observations, the De Lucia et al. models favours the scenario in which the second-ranked galaxy is approximately one magnitude fainter than the primary,especially when the dynamic range of the mock data is limited to 3 magnitudes.

We study the properties of simulated high-redshift galaxies using cosmological N-body/gasdynamical runs from the OverWhelmingly Large Simulations (OWLS) project. The runs contrast several feedback implementations of varying effectiveness: from no-feedback, to supernova-driven winds to powerful AGN-driven outflows. These different feedback models result in large variations in the abundance and structural properties of bright galaxies at z=2. We find that feedback affects the baryonic mass of a galaxy much more severely than its spin, which is on average roughly half that of its surrounding dark matter halo in our runs. Feedback induces strong correlations between angular momentum content and galaxy mass that leave their imprint on galaxy scaling relations and morphologies. Encouragingly, we find that galaxy disks are common in moderate-feedback runs, making up typically ~50% of all galaxies at the centers of haloes with virial mass exceeding 1e11 M_sun. The size, stellar masses, and circular speeds of simulated galaxies formed in such runs have properties that straddle those of large star-forming disks and of compact early-type galaxies at z=2. Once the detailed abundance and structural properties of these rare objects are well established it may be possible to use them to gauge the overall efficacy of feedback in the formation of high redshift galaxies.

Ultra-diffuse galaxies (UDGs) are unusual galaxies with low luminosities, similar to classical dwarf galaxies, but sizes up to $\sim\!5$ larger than expected for their mass. Some UDGs have large populations of globular clusters (GCs), something unexpected in galaxies with such low stellar density and mass. We have carried out a comprehensive study of GCs in both UDGs and classical dwarf galaxies at comparable stellar masses using HST observations of the Coma cluster. We present new imaging for 33 Dragonfly UDGs with the largest effective radii ($>2$ kpc), and additionally include 15 UDGs and 54 classical dwarf galaxies from the HST/ACS Coma Treasury Survey and the literature. Out of a total of 48 UDGs, 27 have statistically significant GC systems, and 11 have candidate nuclear star clusters. The GC specific frequency ($S_N$) varies dramatically, with the mean $S_N$ being higher for UDGs than for classical dwarfs. At constant stellar mass, galaxies with larger sizes (or lower surface brightnesses) have higher $S_N$, with the trend being stronger at higher stellar mass. At lower stellar masses, UDGs tend to have higher $S_N$ when closer to the center of the cluster, i.e., in denser environments. The fraction of UDGs with a nuclear star cluster also depends on environment, varying from $\sim\!40$\% in the cluster core, where it is slightly lower than the nucleation fraction of classical dwarfs, to $\lesssim20\%$ in the outskirts. Collectively, we observe an unmistakable diversity in the abundance of GCs, and this may point to multiple formation routes., 15 pages, 10 figures, 2 tables, accepted for publication in ApJ

We use N-body simulations to study the effects of tides on the kinematical structure of satellite galaxies orbiting a Milky Way-like potential. Our work is motivated by observations of dwarf spheroidal galaxies in the Local Group, for which often a distinction is possible between a cold centrally concentrated metal rich and a hot, extended metal poor population. We find that an important attenuation of the initial differences in the distribution of the two stellar components occurs for orbits with small pericentric radii (r_per < 20 kpc). This is mainly due to: i) the loss of the gravitational support provided by the dark matter component after tidal stripping takes place, which forces a re-configuration of the luminous components, and ii) tides preferentially affect the more extended stellar component, leading to a net decrease in its velocity dispersion as a response for the mass loss, which thus shrinks the kinematical gap. We apply these ideas to the Sculptor and Carina dwarf spheroidals. Differences in their orbits might help to explain, under the assumption of similar initial configurations, why in the former a clear kinematical separation between metal poor and metal rich stars is apparent, while in Carina this segregation is significantly more subtle., 17 pages, 7 figures, Accepted for publication in Advances in Astronomy, special issue on "Dwarf-Galaxy Cosmology"

We analyse the projected radial distribution of satellites around bright primary galaxies in the 2dFGRS. We have considered several primary-satellite subsamples to search for dependences of the satellite number density profile, \rho(r_p), on properties of satellites and primaries. We find significant differences of the behaviour of \rho(r_p) depending on primary characteristics. In star-forming primaries, the satellite number density profile is consistent with power laws within projected distance 20100 kpc), the density profiles of all primaries is well described by power laws, although we notice that for red, early spectral type primaries, the outer slope obtained is steeper than that corresponding to blue, late spectral type ones. We have tested our results by control samples of galaxies identical to the samples of satellites in apparent magnitude and projected distance to the primary, but with a large relative velocity. This sample of unphysical primary-galaxy pairs shows a flat radial density beyond r_p=20 kpc indicating that our results are not biased toward a decrease of the true number of objects due to catalogue selection effects. Our results can be understood in terms of dynamical friction and tidal stripping on satellites in the primary haloes. These processes can effectively transfer energy to the dark matter, flattening the central steep profiles of the satellite distribution in evolved systems., Comment: 8 pages, 6 figures. Accepted for publication in MNRAS

We have analysed a sample of satellite and primary galaxies in the 2dF galaxy redshift survey. In our study we find a strong statistical evidence of the Holmberg effect (that is a tendency for satellites to avoid regions along the line defined by the primary plane) within 500 kpc of projected distance to the primary. This effect is present only when we restrict to objects with radial velocity relative to the primary $|\Delta v

We investigate the role of the delineated cosmic web/filaments on the star formation activity by exploring a sample of 425 narrow-band selected H{\alpha} emitters, as well as 2846 color-color selected underlying star-forming galaxies for a large scale structure (LSS) at z=0.84 in the COSMOS field from the HiZELS survey. Using the scale-independent Multi-scale Morphology Filter (MMF) algorithm, we are able to quantitatively describe the density field and disentangle it into its major components: fields, filaments and clusters. We show that the observed median star formation rate (SFR), stellar mass, specific star formation rate (sSFR), the mean SFR-Mass relation and its scatter for both H{\alpha} emitters and underlying star-forming galaxies do not strongly depend on different classes of environment, in agreement with previous studies. However, the fraction of H{\alpha} emitters varies with environment and is enhanced in filamentary structures at z~1. We propose mild galaxy-galaxy interactions as the possible physical agent for the elevation of the fraction of H{\alpha} star-forming galaxies in filaments. Our results show that filaments are the likely physical environments which are often classed as the "intermediate" densities, and that the cosmic web likely plays a major role in galaxy formation and evolution which has so far been poorly investigated., Comment: 15 pages, 11 figures, Accepted for publication in the ApJ

We study how feedback influences baryon infall onto galaxies using cosmological, zoom-in simulations of haloes with present mass $M_{vir}=6.9\times10^{11} M_{\odot}$ to $1.7\times10^{12} M_{\odot}$. Starting at z=4 from identical initial conditions, implementations of weak and strong stellar feedback produce bulge- and disc-dominated galaxies, respectively. Strong feedback favours disc formation: (1) because conversion of gas into stars is suppressed at early times, as required by abundance matching arguments, resulting in flat star formation histories and higher gas fractions; (2) because 50% of the stars form in situ from recycled disc gas with angular momentum only weakly related to that of the z=0 dark halo; (3) because late-time gas accretion is typically an order of magnitude stronger and has higher specific angular momentum, with recycled gas dominating over primordial infall; (4) because 25-30% of the total accreted gas is ejected entirely before z~1, removing primarily low angular momentum material which enriches the nearby inter-galactic medium. Most recycled gas roughly conserves its angular momentum, but material ejected for long times and to large radii can gain significant angular momentum before re-accretion. These processes lower galaxy formation efficiency in addition to promoting disc formation., Comment: 23 pages, 29 figures, accepted for publication in 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 a semianalytic galaxy catalogue constructed from the Millennium Simulation to study the satellites of isolated galaxies in the LCDM cosmogony. This sample (~80,000$ bright primaries, surrounded by ~178,000 satellites) allows the characterization, with minimal statistical uncertainty, of the dynamical properties of satellite/primary galaxy systems in a LCDM universe. We find that, overall, the satellite population traces the dark matter rather well: its spatial distribution and kinematics may be approximated by an NFW profile with a mildly anisotropic velocity distribution. Their spatial distribution is also mildly anisotropic, with a well-defined ``anti-Holmberg'' effect that reflects the misalignment between the major axis and angular momentum of the host halo. The isolation criteria for our primaries picks not only galaxies in sparse environments, but also a number of primaries at the centre of ''fossil'' groups. We find that the abundance and luminosity function of these unusual systems are in reasonable agreement with the few available observational constraints. We recover the expected L_{host} \sigma_{sat}^3 relation for LCDM models for truly-isolated primaries. Less strict primary selection, however, leads to substantial modification of the scaling relation. Our analysis also highlights a number of difficulties afflicting studies that rely on blind stacking of satellite systems to constrain the mean halo mass of the primary galaxies., Comment: 18 pages, 14 figures, MNRAS in press. Accepted version with minor changes. Version with high resolution figures available at: http://www.astro.uvic.ca/~lsales/SatPapers/SatPapers.html

We study the population of satellite galaxies formed in a suite of N-body/gasdynamical simulations of galaxy formation in a LCDM universe. We find little spatial or kinematic bias between the dark matter and the satellite population. The velocity dispersion of the satellites is a good indicator of the virial velocity of the halo: \sigma_{sat}/V_{vir}=0.9 +/- 0.2. Applied to the Milky Way and M31 this gives V_{vir}^{MW}=109 +/- 22$ km/s and V_{vir}^{M31} = 138 +/- 35 km/s, respectively, substantially lower than the rotation speed of their disk components. The detailed kinematics of simulated satellites and dark matter are also in good agreement. By contrast, the stellar halo of the simulated galaxies is kinematically and spatially distinct from the population of surviving satellites. This is because the survival of a satellite depends on mass and on time of accretion; surviving satellites are biased toward low-mass systems that have been recently accreted by the galaxy. Our results support recent proposals for the origin of the systematic differences between stars in the Galactic halo and in Galactic satellites: the elusive ``building blocks'' of the Milky Way stellar halo were on average more massive, and were accreted (and disrupted) earlier than the population of dwarfs that has survived self-bound until the present., Comment: 13 pages, 11 figures, MNRAS in press. Accepted version with minor changes. Version with high resolution figures available at: http://www.astro.uvic.ca/~lsales/SatPapers/SatPapers.html