Bright Satellites Research Articles

Halo densities and pericenter distances of the bright Milky Way satellites as a test of dark matter physics

ABSTRACT We provide new constraints on the dark matter halo density profile of Milky Way (MW) dwarf spheroidal galaxies (dSphs) using the phase-space distribution function (DF) method. After assessing the systematics of the approach against mock data from the Gaia Challenge project, we apply the DF analysis to the entire kinematic sample of well-measured MW dwarf satellites for the first time. Contrary to previous findings for some of these objects, we find that the DF analysis yields results consistent with the standard Jeans analysis. In particular, in this study we rediscover (i) a large diversity in the inner halo densities of dSphs (bracketed by Draco and Fornax), and (ii) an anticorrelation between inner halo density and pericenter distance of the bright MW satellites. Regardless of the strength of the anticorrelation, we find that the distribution of these satellites in density versus pericenter space is inconsistent with the results of the high-resolution N-body simulations that include a disc potential. Our analysis motivates further studies on the role of internal feedback and dark matter microphysics in these dSphs.

LIGHTS. Survey Overview and a Search for Low Surface Brightness Satellite Galaxies

We present an overview of the LBT Imaging of Galactic Halos and Tidal Structures survey, which currently includes 25 nearby galaxies that are on average ∼1 mag fainter than the Milky Way, and a catalog of 54 low central surface brightness (24 < μ 0,g /mag arcsec−2 < 28) satellite galaxy candidates, most of which were previously uncatalogued. The depth of the imaging exceeds the full 10 yr depth of the Rubin Observatory’s Legacy Survey of Space and Time. We find, after applying completeness corrections, rising numbers of candidate satellites as we approach the limiting luminosity (M r ∼ −8 mag) and central surface brightness (μ 0,g ∼ 28 mag arcsec−2). Over the parameter range we explore, each host galaxy (excluding those that are in overdense regions, apparently groups) has nearly four such candidate satellites to a projected radius of ∼100 kpc. These objects are mostly just at or beyond the reach of spectroscopy unless they are H i rich or have ongoing star formation. We identify three, possibly four, ultra-diffuse satellite galaxies (effective radius >1.5 kpc). This incidence rate falls within expectations of the extrapolation of the published relationship between the number of ultra-diffuse satellite galaxies and host halo mass. Last, we visually identify 12 candidate satellites that host a nuclear star cluster (NSC). The NSC occupation fraction for the sample (12/54) matches that published for satellites of early-type galaxies, suggesting that the parent’s morphological type plays at most a limited role in determining the NSC occupation fraction.

Forecasts for Galaxy Formation and Dark Matter Constraints from Dwarf Galaxy Surveys

The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)–mass cosmological zoom-in simulations and an empirical galaxy–halo connection model, and assess sensitivity to galaxy formation and DM signals when marginalizing over galaxy–halo connection uncertainties. We find that a survey of all satellites around one MW-mass host can constrain a galaxy formation cutoff at peak virial masses of M50=108M⊙ at the 1σ level; however, a tail toward low M50 prevents a 2σ measurement. In this scenario, combining hosts with differing bright satellite abundances significantly reduces uncertainties on M50 at the 1σ level, but the 2σ tail toward low M50 persists. We project that observations of one (two) complete satellite populations can constrain warm DM models with m WDM ≈ 10 keV (20 keV). Subhalo mass function (SHMF) suppression can be constrained to ≈70%, 60%, and 50% that in cold dark matter (CDM) at peak virial masses of 108, 109, and 1010 M ⊙, respectively; SHMF enhancement constraints are weaker (≈20, 4, and 2 times that in CDM, respectively) due to galaxy–halo connection degeneracies. These results motivate searches for faint dwarf galaxies beyond the MW and indicate that ongoing missions like Euclid and upcoming facilities including the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope will probe new galaxy formation and DM physics.

Satellite Visibility During the 2024 April Total Eclipse

On 2024 April 8, tens of millions of people across North America will be able to view a total solar eclipse. Such astronomical events have been important throughout history, but with nearly 10,000 satellites in orbit, we question whether total eclipses will now reveal a sky full of satellites, fundamentally changing this experience for humanity. Using the current population of Starlink satellites, we find that the brightest satellites would be naked-eye visible in dark skies, but the high sky brightness during totality will make them undetectable to the unaided eye. Our model does not take into account the effects of chance reflections from large, flat surfaces like solar panels, which we expect will cause glints and flares that could be visible from large satellites and abandoned rocket bodies. Time-lapse all-sky imaging might reveal satellites during the eclipse.

Near-Earth Object Observations using Synthetic Tracking

Synthetic tracking (ST) has emerged as a potent technique for observing fast-moving near-Earth objects (NEOs), offering enhanced detection sensitivity and astrometric accuracy by avoiding trailing loss. This approach also empowers small telescopes to use prolonged integration times to achieve high sensitivity for NEO surveys and follow-up observations. In this study, we present the outcomes of ST observations conducted with Pomona College’s 1 m telescope at the Table Mountain Facility and JPL’s robotic telescopes at the Sierra Remote Observatory. The results showcase astrometric accuracy statistics comparable to stellar astrometry, irrespective of an object’s rate of motion, and the capability to detect faint asteroids beyond 20.5th magnitude using 11 inch telescopes. Furthermore, we detail the technical aspects of data processing, including the correction of differential chromatic refraction in the atmosphere and accurate timing for image stacking, which contribute to achieving precise astrometry. We also provide compelling examples that showcase the robustness of ST even when asteroids closely approach stars or bright satellites cause disturbances. Moreover, we illustrate the proficiency of ST in recovering NEO candidates with highly uncertain ephemerides. As a glimpse of the potential of NEO surveys utilizing small robotic telescopes with ST, we present significant statistics from our NEO survey conducted for testing purposes. These findings underscore the promise and effectiveness of ST as a powerful tool for observing fast-moving NEOs, offering valuable insights into their trajectories and characteristics. Overall, the adoption of ST stands to revolutionize fast-moving NEO observations for planetary defense and studying these celestial bodies.

New dwarf galaxy candidates in the sphere of influence of the Sombrero galaxy

ABSTRACT We report the discovery of 40 new satellite dwarf galaxy candidates in the sphere of influence of the Sombrero Galaxy (M104), the most luminous galaxy in the Local Volume. Using the Subaru Hyper Suprime-Cam, we surveyed 14.4 deg2 of its surroundings, extending to the virial radius. Visual inspection of the deep images and galfit modelling yielded a galaxy sample highly complete down to Mg ∼ −9 ($L_{g}\sim 3\times 10^{5}\ \mathrm{ L}_\odot$) and spanning magnitudes −16.4 &amp;lt; Mg &amp;lt; −8 and half-light radii 50 pc &amp;lt; re &amp;lt; 1600 pc assuming the distance of M104. These 40 new candidates, out of which 27 are group members with high confidence, double the number of potential satellites of M104 within the virial radius, placing it among the richest hosts in the Local Volume. Using a principal component analysis, we find that the entire sample of candidates is consistent with an almost circular on-sky distribution, more circular than any comparable environment found in the Illustris TNG100-1 (The Next Generation) simulation. However, the distribution of the high-probability sample is more oblate and consistent with the simulation. The cumulative satellite luminosity function is broadly consistent with analogues from the simulation, albeit it contains no bright satellite with Mg &amp;lt; −16.4 ($L_{g}\sim 3 \times 10^{8}\ \mathrm{ L}_\odot$), a $2.3\, \sigma$ occurrence. Follow-up spectroscopy to confirm group membership will begin to demonstrate how these systems can act as probes of the structure and formation history of the halo of M104.

Galaxy And Mass Assembly: the xSAGA galaxy complement in nearby galaxy groups

ABSTRACT Groups of galaxies are the intermediate density environment in which much of the evolution of galaxies is thought to take place. In spectroscopic redshift surveys, one can identify these as close spatial-redshift associations. However, spectroscopic surveys will always be more limited in luminosity and completeness than imaging ones. Here, we combine the Galaxy And Mass Assembly (GAMA) group catalogue with the extended Satellites Around Galactic Analogues (xSAGA) catalogue of machine learning identified low-redshift satellite galaxies. We find 1825 xSAGA galaxies within the bounds of the GAMA equatorial fields (mr &amp;lt; 21), 1562 of which could have a counterpart in the GAMA spectroscopic catalogue (mr &amp;lt; 19.8). Of these, 1326 do have a GAMA counterpart with 974 below z = 0.03 (true positives) and 352 above z = 0.03 (false positives). By cross-correlating the GAMA group catalogue with the xSAGA catalogue, we can extend and characterize the satellite content of GAMA galaxy groups. We find that most groups have &amp;lt;5 xSAGA galaxies associated with them, but richer groups may have more. Each additional xSAGA galaxy contributes only a small fraction of the group’s total stellar mass (${\ll}10$ per cent). Selecting GAMA groups that resemble the Milky Way halo, with a few (&amp;lt;4) bright galaxies, we find that xSAGA can add a magnitude fainter sources to a group and that the Local Group does not stand out in the number of bright satellites. We explore the quiescent fraction of xSAGA galaxies in GAMA groups and find a good agreement with the literature.

The high optical brightness of the BlueWalker 3 satellite

Large constellations of bright artificial satellites in low Earth orbit pose significant challenges to ground-based astronomy1. Current orbiting constellation satellites have brightnesses between apparent magnitudes 4 and 6, whereas in the near-infrared Ks band, they can reach magnitude 2 (ref. 2). Satellite operators, astronomers and other users of the night sky are working on brightness mitigation strategies3,4. Radio emissions induce further potential risk to ground-based radio telescopes that also need to be evaluated. Here we report the outcome of an international optical observation campaign of a prototype constellation satellite, AST SpaceMobile’s BlueWalker 3. BlueWalker 3 features a 64.3 m2 phased-array antenna as well as a launch vehicle adaptor (LVA)5. The peak brightness of the satellite reached an apparent magnitude of 0.4. This made the new satellite one of the brightest objects in the night sky. Additionally, the LVA reached an apparent V-band magnitude of 5.5, four times brighter than the current International Astronomical Union recommendation of magnitude 7 (refs. 3,6); it jettisoned on 10 November 2022 (Universal Time), and its orbital ephemeris was not publicly released until 4 days later. The expected build-out of constellations with hundreds of thousands of new bright objects1 will make active satellite tracking and avoidance strategies a necessity for ground-based telescopes.

Dark Skies and Bright Satellites

Dark Skies and Bright Satellites

Galactic satellite systems in CDM, WDM and SIDM

ABSTRACT We investigate the population of bright satellites ($M_{*} \ge 10^{5} \, \mathrm{M}_{\odot }$) of haloes of mass comparable to that of the Milky Way in cosmological simulations in which the dark matter (DM) is either cold, warm, or self-interacting (CDM, WDM, and SIDM, respectively). The nature of the DM gives rise to differences in the abundance and structural properties of field haloes. In WDM, the main feature is a reduction in the total number of galaxies that form, reflecting a suppression of low-mass DM haloes and lower galaxy formation efficiency compared to CDM. For SIDM, the changes are structural, restricted to the central regions of haloes and dependent on the assumed self-interaction cross-section. We also consider different baryonic subgrid physics models for galaxy formation, in which supernova gas blowouts can or cannot induce the formation of a core in dwarf galaxies. Overall, the inclusion of baryons lessen the differences in the halo properties in the different DM models compared to DM-only simulations. This affects the satellite properties at infall and therefore their subsequent tidal stripping and survival rates. None the less, we find slightly less concentrated satellite radial distributions as the SIDM cross-section increases. Unfortunately, we also find that the satellite populations in simulations with baryon-induced cores in CDM and WDM can mimic the results found in SIDM, making the satellite stellar mass and maximum circular velocity functions heavily degenerate on the assumed nature of the DM and the adopted subgrid modelling. These degeneracies preclude using the brightest satellites of the Milky Way to constrain the nature of DM.

The Lensed Lyman-Alpha MUSE Arcs Sample (LLAMAS)

Aims. We present the Lensed Lyman-Alpha MUSE Arcs Sample (LLAMAS) selected from MUSE and HST observations of 17 lensing clusters. The sample consists of 603 continuum-faint (−23 &lt; MUV &lt; −14) lensed Lyman-α emitters (producing 959 images) with secure spectroscopic redshifts between 2.9 and 6.7. Combining the power of cluster magnification with 3D spectroscopic observations, we were able to reveal the resolved morphological properties of 268 Lyman-α emitters. Methods. We used a forward-modeling approach to model both Lyman-α and rest-frame UV continuum emission profiles in the source plane and measure spatial extent, ellipticity, and spatial offsets between UV and Lyman-α emission. Results. We find a significant correlation between UV continuum and Lyman-α spatial extent. Our characterization of the Lyman-α halos indicates that the halo size is linked to the physical properties of the host galaxy (SFR, Lyman-α equivalent width, Lyman-α line FWHM). We find that 48% of Lyman-α halos are best fit by an elliptical emission distribution with a median axis ratio of q = 0.48. We observe that 60% of galaxies detected both in UV and Lyman-α emission show a significant spatial offset (ΔLyα − UV). We measure a median offset of ΔLyα − UV = 0.58 ± 0.14 kpc for the entire sample. By comparing the spatial offset values with the size of the UV component, we show that 40% of the offsets could be due to star-forming sub-structures in the UV component, while the larger offsets (60%) are more likely due to greater-distance processes such as scattering effects inside the circumgalactic medium or emission from faint satellites or merging galaxies. Comparisons with a zoom-in radiative hydrodynamics simulation of a typical Lyman-α emitting galaxy show a very good agreement with LLAMAS galaxies and indicate that bright star-formation clumps and satellite galaxies could produce a similar spatial offset distribution.

Plaskett 1.8 m Observations of Starlink Satellites

Abstract We present observations of 23 Starlink satellites in the g′ bandpass, obtained from the Dominion Astrophysical Observatory’s Plaskett 1.8 m telescope. The targets include a mixture of satellites with and without brightness mitigation measures (i.e., visors). At the time of the observations (2021 July 16), Starlink satellites were sunlit throughout the night, and, even with strict elevation and azimuth limits, there were over 800 candidate Starlink arcs. The satellites altogether have a median range corrected brightness (550 km) of H ¯ g 550 = 5.3 mag. Dividing the targets into those without and with visors, their median range corrected magnitudes are H ¯ g 550 ( no visor ) = 5.1 and H ¯ g 550 ( visor ) = 5.7 mag, respectively. While the visor sample is dimmer in aggregate, the range corrected brightness distribution ranged from H g 550 = 4.3 mag to 9.4 mag. The two brightest satellites have similar H g 550 to within the uncertainty, but one is visored and the other nonvisored. The dimmest satellite we observed is one without mitigations. The intrinsic brightness dispersion among the full sample is σ g = 0.5 mag.

Extending the SAGA Survey (xSAGA). I. Satellite Radial Profiles as a Function of Host-galaxy Properties

Abstract We present “Extending the Satellites Around Galactic Analogs Survey” (xSAGA), a method for identifying low-z galaxies on the basis of optical imaging and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify z &lt; 0.03 galaxies from more-distant objects using image cutouts from the DESI Legacy Imaging Surveys. From the sample of &gt;100,000 CNN-selected low-z galaxies, we identify &gt;20,000 probable satellites located between 36–300 projected kpc from NASA-Sloan Atlas central galaxies in the stellar-mass range 9.5 &lt; log ( M ⋆ / M ⊙ ) &lt; 11 . We characterize the incompleteness and contamination for CNN-selected samples and apply corrections in order to estimate the true number of satellites as a function of projected radial distance from their hosts. Satellite richness depends strongly on host stellar mass, such that more-massive host galaxies have more satellites, and on host morphology, such that elliptical hosts have more satellites than disky hosts with comparable stellar masses. We also find a strong inverse correlation between satellite richness and the magnitude gap between a host and its brightest satellite. The normalized satellite radial distribution between 36–300 kpc does not depend on host stellar mass, morphology, or magnitude gap. The satellite abundances and radial distributions we measure are in reasonable agreement with predictions from hydrodynamic simulations. Our results deliver unprecedented statistical power for studying satellite galaxy populations and highlight the promise of using machine-learning for extending galaxy samples of wide-area surveys.

The formation of M101-alike galaxies in the cold dark matter model

ABSTRACT The population of satellite galaxies in a host galaxy is a combination of the cumulative accretion of subhaloes and their associated star formation efficiencies; therefore, the luminosity distribution of satellites provides valuable information of both dark matter properties and star formation physics. Recently, the luminosity function of satellites in nearby Milky Way-mass galaxies has been well measured to satellites as faint as Leo I with MV ∼ −8. In addition to the finding of the diversity in the satellite luminosity functions, it has been noticed that there is a big gap among the magnitude of satellites in some host galaxies, such as M101, where the gap is around 5 in magnitude, noticeably larger than the prediction from the halo abundance matching method. The reason of this gap is still unknown. In this paper, we use a semi-analytical model of galaxy formation, combined with high-resolution N-body simulation, to investigate the probability and origin of such big gap in M101-alike galaxies. We found that, although M101 analogues are very rare with probability of $\sim 0.1\ \mathrm{ to}\ 0.2\,{{\rm per\, cent}}$ in the local universe, their formation is a natural outcome of the cold dark matter model. The gap in magnitude is mainly due to the mass of the accreted subhaloes, not from the stochastic star formation in them. We also found that the gap is correlated with the total satellite mass and host halo mass. By tracing the formation history of M101-type galaxies, we find that they were likely formed after z ∼ 1 due to the newly accreted bright satellites. The gap is not in a stable state, and it will disappear in 7 Gyr due to mergers of bright satellites with the central galaxy.

Luminosity Functions and Host-to-host Scatter of Dwarf Satellite Systems in the Local Volume

Abstract Low-mass satellites around Milky Way (MW)-like galaxies are important probes of small-scale structure and galaxy formation. However, confirmation of satellite candidates with distance measurements remains a key barrier to fast progress in the Local Volume (LV). We measure the surface brightness fluctuation distances to recently cataloged candidate dwarf satellites around 10 massive hosts within D &lt; 12 Mpc to confirm association. The satellite systems of these hosts are complete and mostly cleaned of contaminants down to M g ∼ −9 to −10, within the area of the search footprints. Joining this sample with hosts surveyed to comparable or better completeness in the literature, we explore how well cosmological simulations combined with common stellar to halo mass relations (SHMR) match observed satellite luminosity functions in the classical satellite luminosity regime. Adopting an SHMR that matches hydrodynamic simulations, we find that the predicted overall satellite abundance agrees well with the observations. The MW is remarkably typical in its luminosity function among LV hosts. We find that the host-to-host scatter predicted by the model is in close agreement with the scatter between the observed systems, once the different masses of the observed systems are taken into account. However, we find significant evidence that the observed systems have more bright and fewer faint satellites than the SHMR model predicts, possibly necessitating a higher normalization of the SHMR around halo masses of 1011 M ⊙ or significantly greater scatter than present in common SHMRs. These results demonstrate the utility of nearby satellite systems in inferring the galaxy–subhalo connection in the low-mass regime.

The orbit of asteroid (317) Roxane’s satellite Olympias from Gemini, Keck, VLT and the SOR, and (22) Kalliope’s Linus from the SOR

A definitive orbit is derived for asteroid (317) Roxane’s satellite Olympias [S/2009 (317)1] by combining the 2009 discovery images from Gemini North (Merline et al. 2009) with images from Keck and the VLT obtained in 2012, as well as images from its 2016-2017 apparition from the Starfire Optical Range. The orbit is retrograde with respect to the ecliptic but in the same sense as Roxane’s spin. Olympias has a period of P=11.9440±0.0005 days, a semi-major axis of a=245±3 km, and an orbital pole at RA=97∘, Dec=−71∘, or ecliptic coordinates λ=245∘, β=−85∘, close to the south ecliptic pole. This satellite orbital pole is only 3∘ from Roxane’s orbital pole (but in a retrograde sense) and restricts all observations of Olympias from Earth to within 4∘ of the satellite’s orbital plane. By fitting the brightness ratios between Roxane (rotational period of 8.16961±0.00005 h) and Olympias as a Fourier series, we find a rotational period for Olympias of 8.2587±0.0001 h, making this an asynchronous wide binary. From the brightness ratios, and with the average infrared modeling diameter found in the literature of 19.16±0.39 km (error of the mean), we estimate triaxial ellipsoid radii of 14.5×8.5×7.2 km for Roxane and 3.6×2.5×2.0 km for Olympias. We can then apportion the mass between the two objects and find a density for both (assumed to be the same) of 2.16±0.18 g/cm3. There are only a few E-type binaries known and this is the first direct determination of E-type density from a binary. We suggest that the system was formed by the Escaping Ejecta Binary (EEB) mechanism of Durda et al. (2004a), probably forming closer together, and then undergoing the complex evolution steps described by Jacobson et al. (2014) involving synchronization, BYORP orbit expansion, loss of tidal lock, and then YORP spinup. Roxane and Olympias may be the only known EEB system to date.From the same 2016-2017 apparition the orbit of Linus around asteroid (22) Kalliope is derived from the SOR. This well-observed bright satellite is found to have a circular orbit with a period of P=3.5956±0.0004 days, in good agreement with the latest elements of Vachier et al. (2012) of P=3.5957±0.0001 days, and a semi-major axis of a=1099±6 km, somewhat greater than their a=1082±11 km for a slightly eccentric orbit (e=0.007±0.010). With a diameter for Kalliope of 161±6 km (Hanuš et al. 2017), we derive a density for Kalliope of 3.72±0.25 g/cm3 from our one apparition study, the same as Hanuš et al. (2017) but greater than the 3.24±0.16 of Vachier et al. (2012).

The twisted dark matter halo of the Milky Way

ABSTRACT We analyse systems analogous to the Milky Way (MW) in the eagle cosmological hydrodynamics simulation in order to deduce the likely structure of the MW’s dark matter (DM) halo. We identify MW mass haloes in the simulation whose satellite galaxies have similar kinematics and spatial distribution to those of the bright satellites of the MW, specifically systems in which the majority of the satellites (8 out of 11) have nearly coplanar orbits that are also perpendicular to the central stellar disc. We find that the normal to the common orbital plane of the coplanar satellites is well aligned with the minor axis of the host DM halo, with a median misalignment angle of only 17.3°. Based on this result, we infer that the minor axis of the Galactic DM halo points towards (l, b) = (182°, −2°), with an angular uncertainty at the 68 and 95 percentile confidence levels of 22° and 43°, respectively. Thus, the inferred minor axis of the MW halo lies in the plane of the stellar disc. The halo, however, is not homologous and its flattening and orientation vary with radius. The inner parts of the halo are rounder than the outer parts and well aligned with the stellar disc (that is the minor axis of the halo is perpendicular to the disc). Further out, the halo twists and the minor axis changes direction by 90°. This twist occurs over a very narrow radial range and reflects variations in the filamentary network along which mass was accreted into the MW.

New dwarfs around the curly spiral galaxy M 63

We present a deep (50‐hr exposed) image of the nearby spiral galaxy M 63 (NGC 5055), taken with a 0.14‐m aperture telescope. The galaxy halo exhibits the known, very faint system of stellar streams extending across 110 kpc. We found five very low‐surface‐brightness dwarf galaxies around M 63. Assuming they are satellites of M 63, their median parameters are: absolute B‐magnitude −8.8 mag, linear diameter 1.3 kpc, surface brightness ∼27.8 mag/sq. arcsec, and linear projected separation 93 kpc. Based on four brighter satellites with measured radial velocities, we found the low orbital mass estimate of M 63 to be (5.1 ± 1.8)1011M⊙ on a scale of ∼216 kpc. The specific property of M 63 is its declining rotation curve. Taking into account the declining rotation curves of the M 63 and three nearby massive galaxies, NGC 2683, NGC 2903, and NGC 3521, we recognize their low mean orbital mass‐to‐K‐band luminosity ratio, (4.8 ± 1.1) M⊙/L⊙, which is only ∼1/6 of the corresponding ratio for the Milky Way and M 31.

Diverse dark matter density at sub-kiloparsec scales in Milky Way satellites: Implications for the nature of dark matter

Milky Way (MW) satellites reside within dark matter (DM) subhalos with a broad distribution of circular velocity profiles. This diversity is enhanced with the inclusion of ultra-faint satellites, which seemingly have very high DM densities, albeit with large systematic uncertainties. We argue that if confirmed, this large diversity in the MW satellite population poses a serious test for the structure formation theory with possible implications for the DM nature. For the Cold Dark Matter model, the diversity might be a signature of the combined effects of subhalo tidal disruption by the MW disk and strong supernova feedback. For models with a dwarf-scale cutoff in the power spectrum, the diversity is a consequence of the lower abundance of dwarf-scale halos. This diversity is most challenging for Self-Interacting Dark Matter (SIDM) models with cross sections $\sigma/m_\chi\gtrsim1~$cm$^2$g$^{-1}$ where subhalos have too low densities to explain the ultra-faint galaxies. We propose a novel solution to explain the diversity of MW satellites based on the gravothermal collapse of SIDM haloes. This solution requires a velocity-dependent cross section that predicts a bimodal distribution of cuspy dense (collapsed) subhaloes consistent with the ultra-faint satellites, and cored lower density subhaloes consistent with the brighter satellites.

The Milky Way’s halo and subhaloes in self-interacting dark matter

ABSTRACT We perform high-resolution simulations of an MW-like galaxy in a self-interacting cold dark matter model with elastic cross-section over mass of $1~\rm cm^2\, g^{-1}$ (SIDM) and compare to a model without self-interactions (CDM). We run our simulations with and without a time-dependent embedded potential to capture effects of the baryonic disc and bulge contributions. The CDM and SIDM simulations with the embedded baryonic potential exhibit remarkably similar host halo profiles, subhalo abundances, and radial distributions within the virial radius. The SIDM host halo is denser in the centre than the CDM host and has no discernible core, in sharp contrast to the case without the baryonic potential (core size ${\sim}7 \, \rm kpc$). The most massive subhaloes (with $V_{\mathrm{peak}}\gt 20 \, \rm km\, s^{-1}$) in our SIDM simulations, expected to host the classical satellite galaxies, have density profiles that are less dense than their CDM analogues at radii less than 500 pc but the deviation diminishes for less massive subhaloes. With the baryonic potential included in the CDM and SIDM simulations, the most massive subhaloes do not display the too-big-to-fail problem. However, the least dense among the massive subhaloes in both these simulations tend to have the smallest pericenter values, a trend that is not apparent among the bright MW satellite galaxies.