We use the Cluster-EAGLE (C-EAGLE) hydrodynamical simulations to investigate the effects of self-interacting dark matter (SIDM) on galaxies as they fall into clusters. We find that SIDM galaxies follow similar orbits to their Cold Dark Matter (CDM) c
ounterparts, but end up with ${sim}$25 per cent less mass by the present day. One in three SIDM galaxies are entirely disrupted, compared to one in five CDM galaxies. However, the excess stripping will be harder to observe than suggested by previous DM-only simulations because the most stripped galaxies form cores and also lose stars: the most discriminating objects become unobservable. The best test will be to measure the stellar-to-halo mass relation (SHMR) for galaxies with stellar mass $10^{10-11},mathrm{M}_{odot}$. This is 8 times higher in a cluster than in the field for a CDM universe, but 13 times higher for an SIDM universe. Given intrinsic scatter in the SHMR, these models could be distinguished with noise-free galaxy-galaxy strong lensing of ${sim}32$ cluster galaxies.
We make the most precise determination to date of the number density of extragalactic 21-cm radio sources as a function of their spectral line widths - the HI velocity width function (HIWF) - based on 22832 sources from the final 7000 deg$^2$ data re
lease of the Arecibo Legacy Fast ALFA (ALFALFA) survey. The number density of sources as a function of their neutral hydrogen masses - the HI mass function (HIMF) - has previously been reported to have a significantly different low-mass slope and knee mass in the two sky regions surveyed during ALFALFA. In contrast with this, we find that the shape of the HIWF in the same two sky regions is remarkably similar, consistent with being identical within the confidence intervals implied by the data (but the overall normalisation differs). The spatial uniformity of the HIWF implies that it is likely a stable tracer of the mass function of dark matter haloes, in spite of the environmental processes to which the measured variation in the HIMF are attributed, at least for galaxies containing enough neutral hydrogen to be detected. This insensitivity of the HIWF to galaxy formation and evolution can be exploited to turn it into a powerful constraint on cosmological models as future surveys yield increasingly precise measurements. We also report on the possible influence of a previously overlooked systematic error affecting the HIWF, which may plausibly see its low-velocity slope steepen by $sim$40 per cent in analyses of future, deeper surveys.
We present measurements of the radial gravitational acceleration around isolated galaxies, comparing the expected gravitational acceleration given the baryonic matter with the observed gravitational acceleration, using weak lensing measurements from
the fourth data release of the Kilo-Degree Survey. These measurements extend the radial acceleration relation (RAR) by 2 decades into the low-acceleration regime beyond the outskirts of the observable galaxy. We compare our RAR measurements to the predictions of two modified gravity (MG) theories: MOND and Verlindes emergent gravity. We find that the measured RAR agrees well with the MG predictions. In addition, we find a difference of at least $6sigma$ between the RARs of early- and late-type galaxies (split by S{e}rsic index and $u-r$ colour) with the same stellar mass. Current MG theories involve a gravity modification that is independent of other galaxy properties, which would be unable to explain this behaviour. The difference might be explained if only the early-type galaxies have significant ($M_{gas} approx M_*$) circumgalactic gaseous haloes. The observed behaviour is also expected in $Lambda$CDM models where the galaxy-to-halo mass relation depends on the galaxy formation history. We find that MICE, a $Lambda$CDM simulation with hybrid halo occupation distribution modelling and abundance matching, reproduces the observed RAR but significantly differs from BAHAMAS, a hydrodynamical cosmological galaxy formation simulation. Our results are sensitive to the amount of circumgalactic gas; current observational constraints indicate that the resulting corrections are likely moderate. Measurements of the lensing RAR with future cosmological surveys will be able to further distinguish between MG and $Lambda$CDM models if systematic uncertainties in the baryonic mass distribution around galaxies are reduced.
We compare the star-forming properties of satellites around Milky Way (MW) analogs from the Stage~II release of the Satellites Around Galactic Analogs Survey (SAGA-II) to those from the APOSTLE and Auriga cosmological zoom-in simulation suites. We us
e archival GALEX UV imaging as a star-formation indicator for the SAGA-II sample and derive star-formation rates (SFRs) to compare with those from APOSTLE and Auriga. We compare our detection rates from the NUV and FUV bands to the SAGA-II H$alpha$ detections and find that they are broadly consistent with over $85%$ of observed satellites detected in all three tracers. We apply the same spatial selection criteria used around SAGA-II hosts to select satellites around the MW-like hosts in APOSTLE and Auriga. We find very good overall agreement in the derived SFRs for the star-forming satellites as well as the number of star-forming satellites per host in observed and simulated samples. However, the number and fraction of quenched satellites in the SAGA-II sample are significantly lower than those in APOSTLE and Auriga below a stellar mass of $M_*sim10^{8},M_{odot}$, even when the SAGA-II incompleteness and interloper corrections are included. This discrepancy is robust with respect to the resolution of the simulations and persists when alternative star-formation tracers are employed. We posit that this disagreement is not readily explained by vagaries in the observed or simulated samples considered here, suggesting a genuine discrepancy that may inform the physics of satellite populations around MW analogs.
We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is loc
ated at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the `stellar splashback feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (mu >> 32 mag/arcsec^2) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid.
We combine orbital information from N-body simulations with an analytic model for star formation quenching and SDSS observations to infer the differential effect of the group/cluster environment on star formation in satellite galaxies. We also consid
er a model for gas stripping, using the same input supplemented with HI fluxes from the ALFALFA survey. The models are motivated by and tested on the Hydrangea cosmological hydrodynamical simulation suite. We recover the characteristic times when satellite galaxies are stripped and quenched. Stripping in massive ($M_mathrm{ vir}sim 10^{14.5},mathrm{M}_odot$) clusters typically occurs at or just before the first pericentric passage. Lower mass ($sim10^{13.5},mathrm{M}_odot$) groups strip their satellites on a significantly longer (by $sim3,mathrm{Gyr}$) timescale. Quenching occurs later: Balmer emission lines typically fade $sim3.5,mathrm{Gyr}$ ($5.5,mathrm{Gyr}$) after first pericentre in clusters (groups), followed a few hundred $mathrm{Myr}$ later by reddenning in $(g-r)$ colour. These `delay timescales are remarkably constant across the entire satellite stellar mass range probed ($sim10^{9.5}-10^{11},mathrm{M}_odot$), a feature closely tied to our treatment of `group pre-processing. The lowest mass groups in our sample ($sim10^{12.5},mathrm{M}_odot$) strip and quench their satellites extremely inefficiently: typical timescales may approach the age of the Universe. Our measurements are qualitatively consistent with the `delayed-then-rapid quenching scenario advocated for by several other studies, but we find significantly longer delay times. Our combination of a homogeneous analysis and input catalogues yields new insight into the sequence of events leading to quenching across wide intervals in host and satellite mass.
The radial acceleration relation (RAR) locally relates the `observed acceleration inferred from the dynamics of a system to the acceleration implied by its baryonic matter distribution. The relation as traced by galaxy rotation curves is one-to-one w
ith remarkably little scatter, implying that the dynamics of a system can be predicted simply by measuring its density profile as traced by e.g. stellar light or gas emission lines. Extending the relation to accelerations below those usually probed by practically observable kinematic tracers is challenging, especially once accounting for faintly emitting baryons, such as the putative warm-hot intergalactic medium, becomes important. We show that in the low-acceleration regime, the (inverted) RAR predicts an unphysical, declining enclosed baryonic mass profile for systems with `observed acceleration profiles steeper than $g_{rm obs}propto r^{-1}$ (corresponding to density profiles steeper than isothermal - $rho(r)propto r^{-2}$). If the RAR is tantamount to a natural law, such acceleration profiles cannot exist. We apply this argument to test the compatibility of an extrapolation of the rotation curve-derived RAR to low accelerations with data from galaxy-galaxy weak lensing, dwarf spheroidal galaxy stellar kinematic, and outer Milky~Way dynamical measurements, fully independent of the uncertainties inherent in direct measurements of the baryonic matter distribution. In all cases we find that the data weakly favour a break to a steeper low-acceleration slope. Improvements in measurements and modelling of the outer Milky~Way, and weak lensing, seem like the most promising path toward stronger constraints on the low-acceleration behaviour of the RAR.
We study the gas kinematics of a sample of six isolated gas-rich low surface brightness galaxies, of the class called ultra-diffuse galaxies (UDGs). These galaxies have recently been shown to be outliers from the baryonic Tully-Fisher relation (BTFR)
, as they rotate much slower than expected given their baryonic mass, and to have baryon fractions similar to the cosmological mean. By means of a 3D kinematic modelling fitting technique, we show that the HI in our UDGs is distributed in thin regularly rotating discs and we determine their rotation velocity and gas velocity dispersion. We revisit the BTFR adding galaxies from other studies. We find a previously unknown trend between the deviation from the BTFR and the disc scale length valid for dwarf galaxies with circular speeds < 45 km/s, with our UDGs being at the extreme end. Based on our findings, we suggest that the high baryon fractions of our UDGs may originate due to the fact that they have experienced weak stellar feedback, likely due to their low star formation rate surface densities, and as a result they did not eject significant amounts of gas out of their discs. At the same time, we find indications that our UDGs may have higher-than-average stellar specific angular momentum, which can explain their large optical scale lengths.
The APOSTLE cosmological hydrodynamical simulation suite is a collection of twelve regions $sim 5$ Mpc in diameter, selected to resemble the Local Group of galaxies in terms of kinematics and environment, and re-simulated at high resolution (minimum
gas particle mass of $10^4,{rm M}_odot$) using the galaxy formation model and calibration developed for the EAGLE project. I select a sample of dwarf galaxies ($60 < V_{rm max}/{rm km},{rm s}^{-1} < 120$) from these simulations and construct synthetic spatially- and spectrally-resolved observations of their 21-cm emission. Using the $^{3{rm D}}$BAROLO tilted-ring modelling tool, I extract rotation curves from the synthetic data cubes. In many cases, non-circular motions present in the gas disc hinder the recovery of a rotation curve which accurately traces the underlying mass distribution; a large central deficit of dark matter, relative to the predictions of cold dark matter N-body simulations, may then be erroneously inferred.
The stellar halos of galaxies encode their accretion histories. In particular, the median metallicity of a halo is determined primarily by the mass of the most massive accreted object. We use hydrodynamical cosmological simulations from the APOSTLE p
roject to study the connection between the stellar mass, the metallicity distribution, and the stellar age distribution of a halo and the identity of its most massive progenitor. We find that the stellar populations in an accreted halo typically resemble the old stellar populations in a present-day dwarf galaxy with a stellar mass $sim 0.2-0.5$ dex greater than that of the stellar halo. This suggest that had they not been accreted, the primary progenitors of stellar halos would have evolved to resemble typical nearby dwarf irregulars.
