Luminous Satellites Research Articles

Exploring the diversity and similarity of radially anisotropic Milky Way-like stellar haloes: implications for disrupted dwarf galaxy searches

ABSTRACT We investigate the properties of mergers comparable to the Gaia–Sausage–Enceladus (GSE) using cosmological hydrodynamical simulations of Milky Way-like galaxies. The merger progenitors span an order of magnitude in their peak stellar mass ($3\times 10^8\lt M_{\star }/\rm {M}_{\odot }\lt 4\times 10^9$) and include both rotation and pressure-supported galaxies (0.10 < D/T < 0.77). In a minority of cases, the GSE-like debris is comprised of stars from more than one merger progenitor. However, there is a close similarity in their chemodynamical properties and the triaxial shapes of their debris, and so it is not always possible to distinguish them. The merger progenitors host a variety of luminous satellites (0 and 8 with $M_{\star }\gt 10^6\, \rm {M}_{\odot }$), but most of these do not follow the merger to low orbital energies. Between 0 and 1 of these satellites may survive to z = 0, but with no clear signatures of their past association. We show that the fraction of stars originating from GSE-like mergers is reduced for lower metallicities (reaching a minimum around [Fe/H] = −2), and also within 5 kpc of the Galactic Centre. Whilst these central regions are dominated by in-situ stars, the ex-situ fraction trends towards a 100 per cent asymptote when considering the most metal-poor stars ([Fe/H] ≪ −2.5). Considering this, its near proximity, and its small volume on the sky, the Galactic Centre lends itself as a prime environment in the search for the stars from the earliest galaxies, whilst avoiding contamination from GSE stars.

The VMC survey

Context. The Large Magellanic Cloud (LMC) is the most luminous satellite galaxy of the Milky Way and, owing to its companion, the Small Magellanic Cloud (SMC), represents an excellent laboratory to study the interaction of dwarf galaxies. Aims. The aim of this study is to investigate the kinematics of the outer regions of the LMC by using stellar proper motions to understand the impact of interactions, for example with the SMC about 250 Myr ago. Methods. We calculate proper motions using multi-epoch Ks-band images from the VISTA survey of the Magellanic Cloud system (VMC). Observations span a time baseline of 2−5 yr. We combine the VMC data with data from the Gaia Early Data Release 3 and introduce a new method to distinguish between Magellanic and Milky Way stars based on a machine learning algorithm. This new technique enables a larger and cleaner sample selection of fainter sources as it reaches below the red clump of the LMC. Results. We investigate the impact of the SMC on the rotational field of the LMC and find hints of stripped SMC debris. The south-eastern region of the LMC shows a slow rotational speed compared to the overall rotation. N-body simulations suggest that this could be caused by a fraction of stripped SMC stars located in that particular region that move opposite to the expected rotation.

Tracing stars in Milky Way satellites with a-sloth

ABSTRACT We study the stellar mass-to-halo mass (SMHM) relation at z = 0 in 30 Milky Way (MW)-like systems down to the ultrafaint ($M_{\rm *}\lt 10^5\, {\rm M}_\odot$) regime using the semi-analytic model a-sloth. A new model allows us to follow star formation (SF) and the stochastic stellar feedback from individually sampled Pop II stars. Our fiducial model produces consistent results with the SMHM relation derived from abundance matching and the observed cumulative stellar mass function above the observational completeness. We find a plateau in the SMHM relation in the ultrafaint regime. The stellar mass of this plateau tells us how many stars formed before supernovae occur and regulate further SF, which is determined by the Pop II SF efficiency. We also find that the number of luminous satellites increases rapidly as M* decreases until $M_{\rm *}\approx 10^4\, {\rm M}_\odot$. Finally, we find that the relative streaming velocity between baryons and dark matter at a high redshift is important in determining the number of ultrafaint dwarf galaxies at z = 0. The new model in a-sloth provides a framework to study the stellar properties and the formation history of metal-poor stars in MW and its satellites.

The spatial distribution of Milky Way satellites, gaps in streams, and the nature of dark matter

ABSTRACT The spatial distribution of Milky Way (MW) subhaloes provides an important set of observables for testing cosmological models. These include the radial distribution of luminous satellites, planar configurations, and the abundance of dark subhaloes whose existence or absence is key to distinguishing among dark matter models. We use the cocoN-body simulations of cold dark matter (CDM) and 3.3 keV thermal relic warm dark matter (WDM) to predict the satellite spatial distribution in the limit that the impact of baryonic physics is minimal. We demonstrate that the radial distributions of CDM and 3.3 keV-WDM luminous satellites are identical if the minimum pre-infall halo mass to form a galaxy is >108.5 ${\, \rm M_\odot }$. The distribution of dark subhaloes is significantly more concentrated in WDM due to the absence of low mass, recently accreted substructures that typically inhabit the outer parts of a MW halo in CDM. We show that subhaloes of mass [107, 108] ${\, \rm M_\odot }$ and within 30 kpc of the centre are the stripped remnants of larger haloes in both models. Therefore, their abundance in WDM is 3× higher than one would anticipate from the overall WDM subhalo population. We estimate that differences between CDM and WDM concentration–mass relations can be probed for subhalo–stream impact parameters <2 kpc. Finally, we find that the impact of WDM on planes of satellites is likely negligible. Comprehensive comparisons with observations will require further work with high resolution, self-consistent hydrodynamical simulations.

Joint constraints on thermal relic dark matter from strong gravitational lensing, the Ly α forest, and Milky Way satellites

ABSTRACT We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Ly α forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of λhm = 0.089 Mpc h−1 at the 95 per cent confidence level, which we show to be stable for a broad range of prior choices. Assuming a Planck cosmology and that WDM particles are thermal relics, this corresponds to an upper limit on the half-mode mass of Mhm < 3 × 107 M⊙ h−1, and a lower limit on the particle mass of mth > 6.048 keV, both at the 95 per cent confidence level. We find that models with λhm > 0.223 Mpc h−1 (corresponding to mth > 2.552 keV and Mhm < 4.8 × 108 M⊙ h−1) are ruled out with respect to the maximum likelihood model by a factor ≤1/20. For lepton asymmetries L6 > 10, we rule out the 7.1 keV sterile neutrino dark matter model, which presents a possible explanation to the unidentified 3.55 keV line in the Milky Way and clusters of galaxies. The inferred 95 percentiles suggest that we further rule out the ETHOS-4 model of self-interacting DM. Our results highlight the importance of extending the current constraints to lower half-mode scales. We address important sources of systematic errors and provide prospects for how the constraints of these probes can be improved upon in the future.

Magellanic satellites in ΛCDM cosmological hydrodynamical simulations of the Local Group

ABSTRACT We use the APOSTLE ΛCDM cosmological hydrodynamical simulations of the Local Group to study the recent accretion of massive satellites into the halo of Milky Way (MW)-sized galaxies. These systems are selected to be close analogues to the Large Magellanic Cloud (LMC), the most massive satellite of the MW. The simulations allow us to address, in a cosmological context, the impact of the Clouds on the MW, including the contribution of Magellanic satellites to the MW satellite population, and the constraints placed on the Galactic potential by the motion of the LMC. We show that LMC-like satellites are twice more common around Local Group-like primaries than around isolated haloes of similar mass; these satellites come from large turnaround radii and are on highly eccentric orbits whose velocities at first pericentre are comparable with the primary’s escape velocity. This implies $V_{\rm esc}^{\rm MW} (50$ kpc) ∼ 365 km s−1, a strong constraint on Galactic potential models. LMC analogues contribute about two satellites with $M_*\gt 10^5\, \mathrm{ M}_\odot$, having thus only a mild impact on the luminous satellite population of their hosts. At first pericentre, LMC-associated satellites are close to the LMC in position and velocity, and are distributed along the LMC’s orbital plane. Their orbital angular momenta roughly align with the LMC’s, but, interestingly, they may appear to ‘counter-rotate’ the MW in some cases. These criteria refine earlier estimates of the LMC association of MW satellites: only the SMC, Hydrus1, Car3, Hor1, Tuc4, Ret2, and Phoenix2 are compatible with all criteria. Carina, Grus2, Hor2, and Fornax are less probable associates given their large LMC relative velocity.

Constraining the primordial magnetic field with dwarf galaxy simulations

Using a set of cosmological hydro-dynamical simulations, we constrained the properties of primordial magnetic fields by studying their impact on the formation and evolution of dwarf galaxies. We performed a large set of simulations (8 dark matter only and 72 chemo-hydrodynamical) including primordial magnetic fields through the extra density fluctuations they induce at small length scales (k ≥ 10 h Mpc−1) in the matter power spectrum. Our sample of dwarfs includes nine systems selected out of the initial (3.4 Mpc h−1)3 parent box, resimulated from z = 200 to z = 0 using a zoom-in technique and including the physics of baryons. We explored a wide variety of primordial magnetic fields with strength Bλ ranging from 0.05 to 0.50 nG and magnetic energy spectrum slopes nB from −2.9 to −2.1. Strong magnetic fields characterized by a high amplitude (Bλ = 0.50, 0.20 nG with nB = −2.9) or by a steep initial power spectrum slope (nB = −2.1, −2.4, with Bλ = 0.05 nG) induce perturbations on mass scales from 107 to 109 M⊙. In this context emerging galaxies see their star formation rates strongly boosted. They become more luminous and metal rich than their counterparts without primordial magnetic fields. Such strong fields are ruled out by their inability to reproduce the observed scaling relations of dwarf galaxies. They predict that dwarf galaxies are at the origin of an unrealistically early reionization of the Universe and that they also overproduce luminous satellites in the Local Group. Weaker magnetic fields impacting the primordial density field at corresponding masses ≲106 M⊙, produce a large number of mini dark matter halos orbiting the dwarfs, however out of reach for current lensing observations. This study allows us, for the first time, to constrain the properties of primordial magnetic fields based on realistic cosmological simulations of dwarf galaxies.

SHARP – VII. New constraints on the dark matter free-streaming properties and substructure abundance from gravitationally lensed quasars

ABSTRACT We present an analysis of seven strongly gravitationally lensed quasars and the corresponding constraints on the properties of dark matter. Our results are derived by modelling the lensed image positions and flux-ratios using a combination of smooth macro-models and a population of low-mass haloes within the mass range of 106–109 M⊙. Our lens models explicitly include higher order complexity in the form of stellar discs and luminous satellites, as well as low-mass haloes located along the observed lines of sight for the first time. Assuming a cold dark matter (CDM) cosmology, we infer an average total mass fraction in substructure of $f_{\rm sub} = 0.012^{+0.007}_{-0.004}$ (68 per cent confidence limits), which is in agreement with the predictions from CDM hydrodynamical simulations to within 1σ. This result is closer to the predictions than those from previous studies that did not include line-of-sight haloes. Under the assumption of a thermal relic dark matter model, we derive a lower limit on the particle relic mass of mth > 5.58 keV (95 per cent confidence limits), which is consistent with a value of mth > 5.3 keV from the recent analysis of the Ly α forest. We also identify two main sources of possible systematic errors and conclude that deeper investigations in the complex structure of lens galaxies as well as the size of the background sources should be a priority for this field.

Under the FIRElight: Stellar Tracers of the Local Dark Matter Velocity Distribution in the Milky Way

The Gaia era opens new possibilities for discovering the remnants of disrupted satellite galaxies in the Solar neighborhood. If the population of local accreted stars is correlated with the dark matter sourced by the same mergers, one can then map the dark matter distribution directly. Using two cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies from the Latte suite of Fire-2 simulations, we find a strong correlation between the velocity distribution of stars and dark matter at the solar circle that were accreted from luminous satellites. This correspondence holds for dark matter that is either relaxed or in kinematic substructure called debris flow, and is consistent between two simulated hosts with different merger histories. The correspondence is more problematic for streams because of possible spatial offsets between the dark matter and stars. We demonstrate how to reconstruct the dark matter velocity distribution from the observed properties of the accreted stellar population by properly accounting for the ratio of stars to dark matter contributed by individual mergers. After demonstrating this method using the Fire-2 simulations, we apply it to the Milky Way and use it to recover the dark matter velocity distribution associated with the recently discovered stellar debris field in the Solar neighborhood. Based on results from Gaia, we estimate that $42 ^{+26}_{-22}\%$ of the local dark matter that is accreted from luminous mergers is in debris flow.

Dark and luminous satellites of LMC-mass galaxies in the FIRE simulations

ABSTRACT Within lambda cold dark matter ($\Lambda$CDM), dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhaloes, 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$$\gt 10^5$ M$_{\odot }$; Carina and Fornax) as well as several ultrafaint 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$$\sim$6$\times 10^ {3}$ M$_{\odot }$) substructure population of isolated LMC-mass hosts ($M_{\text{200m}}$ = 1–3$\times 10^ {11}$ M$_{\odot }$) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhaloes in simulations with and without baryons, we find that, within 0.2 $r_{\text{200m}}$, LMC-mass hosts deplete $\sim$30 per cent of their substructure, significantly lower than the $\sim$70 per cent of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$ = 880 M$_{\odot }$), $\sim 5\!-\!10$ subhaloes form galaxies with $M_\ast$$\ge 10^{4}$ M$_{\odot }$ , in agreement with the seven 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 ultrafaint 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.

Proper motion of the Magellanic Bridge: Removal of foreground stars

AbstractThe Large and Small Magellanic Cloud (LMC and SMC) are the most luminous dwarf galaxy satellites of the Milky Way. Thanks to their close proximity (50-60 kpc), they provide one of the best opportunities to study in detail the kinematics of resolved stellar populations in an interacting pair of galaxies. Large photometric surveys like the ongoing Gaia mission and the near-infrared VISTA survey of the Magellanic Cloud system (VMC) will have a significant impact on our insight into the Magellanic system. We have combined the individual strengths of VMC and Gaia DR2 data to improve our understanding of the internal kinematics of the galaxies. In this study, we present results from our ongoing project dedicated to measure and analyse the proper motions of large samples of stars across the Magellanic Clouds, efficiently removing Milk Way foreground stars utilising distances derived with the StarHorse code.

Inferred Evidence for Dark Matter Kinematic Substructure with SDSS–Gaia

Abstract We use the distribution of accreted stars in Sloan Digital Sky Survey–Gaia DR2 to demonstrate that a nontrivial fraction of the dark matter halo within galactocentric radii of 7.5–10 kpc and is in substructure and thus may not be in equilibrium. Using a mixture likelihood analysis, we separate the contributions of an old, isotropic stellar halo and a younger anisotropic population. The latter dominates and is uniform within the region studied. It can be explained as the tidal debris of a disrupted massive satellite on a highly radial orbit and is consistent with mounting evidence from recent studies. Simulations that track the tidal debris from such mergers find that the dark matter traces the kinematics of its stellar counterpart. If so, our results indicate that a component of the nearby dark matter halo that is sourced by luminous satellites is in kinematic substructure referred to as debris flow. These results challenge the Standard Halo Model, which is discrepant with the distribution recovered from the stellar data, and have important ramifications for the interpretation of direct detection experiments.

Missing Satellites Problem: Completeness Corrections to the Number of Satellite Galaxies in the MilkyWay are Consistent with Cold Dark Matter Predictions.

A critical challenge to the cold dark matter (CDM) paradigm is that there are fewer satellites observed around the MilkyWay than found in simulations of dark matter substructure. We show that there is a match between the observed satellite counts corrected by the detection efficiency of the Sloan Digital Sky Survey (for luminosities L≳340 L_{⊙}) and the number of luminous satellites predicted by CDM, assuming an empirical relation between stellar mass and halo mass. The "missing satellites problem," cast in terms of number counts, is thus solved. We also show that warm dark matter models with a thermal relic mass smaller than 4keV are in tension with satellite counts, putting pressure on the sterile neutrino interpretation of recent x-ray observations. Importantly, the total number of MilkyWay satellites depends sensitively on the spatial distribution of satellites, possibly leading to a "too many satellites" problem. Measurements of completely dark halos below 10^{8} M_{⊙}, achievable with substructure lensing and stellar stream perturbations, are the next frontier for tests of CDM.

No assembly required: mergers are mostly irrelevant for the growth of low-mass dwarf galaxies

We investigate the merger histories of isolated dwarf galaxies based on a suite of 15 high-resolution cosmological zoom-in simulations, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_{\odot}$ (and M$_\star\sim10^5-10^7\,{\rm M}_{\odot}$) at $z=0$, from the Feedback in Realistic Environments (FIRE) project. The stellar populations of these dwarf galaxies at $z=0$ are formed essentially entirely "in situ": over 90$\%$ of the stellar mass is formed in the main progenitor in all but two cases, and all 15 of the galaxies have >70$\%$ of their stellar mass formed in situ. Virtually all galaxy mergers occur prior to $z\sim3$, meaning that accreted stellar populations are ancient. On average, our simulated dwarfs undergo 5 galaxy mergers in their lifetimes, with typical pre-merger galaxy mass ratios that are less than 1:10. This merger frequency is generally comparable to what has been found in dissipationless simulations when coupled with abundance matching. Two of the simulated dwarfs have a luminous satellite companion at $z=0$. These ultra-faint dwarfs lie at or below current detectability thresholds but are intriguing targets for next-generation facilities. The small contribution of accreted stars make it extremely difficult to discern the effects of mergers in the vast majority of dwarfs either photometrically or using resolved-star color-magnitude diagrams (CMDs). The important implication for near-field cosmology is that star formation histories of comparably massive galaxies derived from resolved CMDs should trace the build-up of stellar mass in one main system across cosmic time as opposed to reflecting the contributions of many individual star formation histories of merged dwarfs.

The Strong Gravitationally Lensed Herschel Galaxy HLock01: Optical Spectroscopy Reveals a Close Galaxy Merger with Evidence of Inflowing Gas

Abstract The submillimeter galaxy (SMG) HERMES J105751.1+573027 (hereafter HLock01) at z = 2.9574 ± 0.0001 is one of the brightest gravitationally lensed sources discovered in the Herschel Multi-tiered Extragalactic Survey. Apart from the high flux densities in the far-infrared, it is also extremely bright in the rest-frame ultraviolet (UV), with a total apparent magnitude m UV ≃ 19.7 mag. We report here deep spectroscopic observations with the Gran Telescopio Canarias of the optically bright lensed images of HLock01. Our results suggest that HLock01 is a merger system composed of the Herschel-selected SMG and an optically bright Lyman break-like galaxy (LBG), separated by only 3.3 kpc in projection. While the SMG appears very massive (M * ≃ 5 × 1011 M ⊙), with a highly extinguished stellar component (A V ≃ 4.3 ), the LBG is a young, lower-mass (M * ≃ 1 × 1010 M ⊙), but still luminous ( ) satellite galaxy. Detailed analysis of the high signal-to-noise ratio (S/N) rest-frame UV spectrum of the LBG shows complex kinematics of the gas, exhibiting both blueshifted and redshifted absorption components. While the blueshifted component is associated with strong galactic outflows from the massive stars in the LBG, as is common in most star-forming galaxies, the redshifted component may be associated with gas inflow seen along a favorable sightline to the LBG. We also find evidence of an extended gas reservoir around HLock01 at an impact parameter of 110 kpc, through the detection of C ii λλ1334 absorption in the red wing of a bright Lyα emitter at z ≃ 3.327. The data presented here highlight the power of gravitational lensing in high S/N studies to probe deeply into the physics of high-z star-forming galaxies.

The predicted luminous satellite populations around SMC- and LMC-mass galaxies – a missing satellite problem around the LMC?

Recent discovery of many dwarf satellite galaxies in the direction of the Small and Large Magellanic Clouds (SMC and LMC) provokes questions of their origins, and what they can reveal about galaxy evolution theory. Here, we predict the satellite stellar mass function of Magellanic Cloud-mass host galaxies using abundance matching and reionization models applied to the \textit{Caterpillar} simulations. Specifically focusing on the volume within $50$~kpc of the LMC, we predict a mean of 4-8 satellites with stellar mass $M_* > 10^4\, \mathrm{M_\odot}$, and 3-4 satellites with $80 < M_* \leq 3000\, \mathrm{M_\odot}$. Surprisingly, all $12$ currently known satellite candidates have stellar masses of $80 < M_* \leq 3000\, \mathrm{M_\odot}$. Reconciling the dearth of large satellites and profusion of small satellites is challenging and may require a combination of a major modification of the \mstarmhalo relationship (steep, but with an abrupt flattening at $10^3\, \mathrm{M_\odot}$), late reionization for the Local Group ($z_{\rm{reion}} \lesssim 9$ preferred), and/or strong tidal stripping. We can more robustly predict that $\sim 53\%$ of satellites within this volume were accreted together with the LMC and SMC, and $\sim 47\%$ were only ever Milky Way satellites. Observing satellites of isolated LMC-sized field galaxies is essential to placing the LMC in context, and to better constrain the \mstarmhalo relationship. Modeling known LMC-sized galaxies within $8$ Mpc, we predict 1-6 (2-12) satellites with $M_* > 10^5\, \mathrm{M_\odot}$ ($M_* > 10^4\, \mathrm{M_\odot}$) within the virial volume of each, and 1-3 (1-7) within a single $1.5^{\circ}$ diameter field of view, making their discovery likely.

An observer's guide to the (Local Group) dwarf galaxies: predictions for their own dwarf satellite populations

A recent surge in the discovery of new ultrafaint dwarf satellites of the Milky Way has inspired the idea of searching for faint satellites, $10^3\, \mathrm{M_{\odot}} < M_* < 10^6\, \mathrm{M_{\odot}}$, around less massive field galaxies in the Local Group. Such satellites would be subject to weaker environmental influences than Milky Way satellites, and could lead to new insights on low mass galaxy formation. In this paper, we predict the number of luminous satellites expected around field dwarf galaxies by applying several abundance matching models and a reionization model to the dark-matter only \textit{Caterpillar} simulation suite. For three of the four abundance matching models used, we find a $>99\%$ chance that at least one satellite with stellar mass $M_*> 10^5\, \mathrm{M_{\odot}}$ exists around the combined five Local Group field dwarf galaxies with the largest stellar mass. When considering satellites with $M_*> 10^4\, \mathrm{M_{\odot}}$, we predict a combined $5-25$ satellites for the five largest field dwarfs, and $10-50$ for the whole Local Group field dwarf population. Because of the relatively small number of predicted dwarfs, and their extended spatial distribution, a large fraction each Local Group dwarf's virial volume will need to be surveyed to guarantee discoveries. We compute the predicted number of satellites in a given field of view of specific Local Group galaxies, as a function of minimum satellite luminosity, and explicitly obtain such values for the Solitary Local dwarfs survey. Uncertainties in abundance matching and reionization models are large, implying that comprehensive searches could lead to refinements of both models.

Dark census: Statistically detecting the satellite populations of distant galaxies

In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Galaxy-scale strong gravitational lenses have been successfully used to study mass substructures lying close to lensed images of bright background sources. However, in typical galaxy-scale lenses, the strong lensing region only covers a small projected area of the lens's dark matter halo, implying that the vast majority of subhalos cannot be directly detected in lensing observations. In this paper, we point out that this large population of dark satellites can collectively affect gravitational lensing observables, hence possibly allowing their statistical detection. Focusing on the region of the galactic halo outside the strong lensing area, we compute from first principles the statistical properties of perturbations to the gravitational time delay and position of lensed images in the presence of a mass substructure population. We find that in the standard cosmological scenario, the statistics of these lensing observables are well approximated by Gaussian distributions. The formalism developed as part of this calculation is very general and can be applied to any halo geometry and choice of subhalo mass function. Our results significantly reduce the computational cost of including a large substructure population in lens models and enable the use of Bayesian inference techniques to detect and characterize the distributed satellite population of distant lens galaxies.

Satellite galaxies in semi-analytic models of galaxy formation with sterile neutrino dark matter

The sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the lepton asymmetry, $L_6$, that is present in the primordial plasma and leads to resonant sterile neutrino production. We compare these power spectra with those of thermally produced dark matter particles and show that resonantly produced sterile neutrinos are much colder than their thermal relic counterparts. We also demonstrate that the shape of these power spectra is not determined by the free-streaming scale alone. We then use the power spectra as an input for semi-analytic models of galaxy formation in order to predict the number of luminous satellite galaxies in a Milky Way-like halo. By assuming that the mass of the Milky Way halo must be no more than $2\times10^{12}M_{\odot}$ (the adopted upper bound based on current astronomical observations) we are able to constrain the value of $L_6$ for $M_s\le 8$~keV. We also show that the range of $L_6$ that is in best agreement with the 3.5~keV line (if produced by decays of 7~keV sterile neutrino) requires that the Milky Way halo has a mass no smaller than $1.5\times10^{12}M_{\odot}$. Finally, we compare the power spectra obtained by direct integration of the Boltzmann equations for a non-resonantly produced sterile neutrino with the fitting formula of Viel~et~al. and find that the latter significantly underestimates the power amplitude on scales relevant to satellite galaxies.

Are rotating planes of satellite galaxies ubiquitous?

We compare the dynamics of satellite galaxies in the Sloan Digital Sky Survey to simple models in order to test the hypothesis that a large fraction of satellites co-rotate in coherent planes. We confirm the previously-reported excess of co-rotating satellite pairs located near diametric opposition with respect to the host, but show that this signal is unlikely to be due to rotating discs (or planes) of satellites. In particular, no overabundance of co-rotating satellites pairs is observed within $\sim 20^{\circ}-50^{\circ}$ of direct opposition, as would be expected for planar distributions inclined relative to the line-of-sight. Instead, the excess co-rotation for satellite pairs within $\sim 10^{\circ}$ of opposition is consistent with random noise associated with undersampling of an underlying isotropic velocity distribution. We conclude that at most $10\%$ of the hosts in our sample harbor co-rotating satellite planes (as traced by the luminous satellite population).