Search

We forward-model mass-weighted stellar ages (MWAs) and quiescent fractions in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in $\log(M_{\mathrm{vir}}/\mathrm{M}_{\odot})>14$ galaxy clusters at $z\sim 0$. We find the average deviation in MWA from the MWA-$M_\star$ relation depends on position in PPS, with a maximum difference between the inner cluster and infalling interloper galaxies of $\sim 1$ Gyr. Our model employs infall information from N-body simulations and stochastic star-formation histories from the UniverseMachine model. We find total quenching times of $t_\mathrm{Q}=3.7\pm 0.4$ Gyr and $t_\mathrm{Q}=4.0\pm 0.2$ Gyr after first pericentre, for $9, 21 pages, 13 figures, submitted to MNRAS

It is routinely assumed that galaxy rotation curves are equal to their circular velocity curves (modulo some corrections) such that they are good dynamical mass tracers. We take a visualisation-driven approach to exploring the limits of the validity of this assumption for a sample of $33$ low-mass galaxies ($60, Comment: MNRAS accepted. Fig. 2 is a video that will play on compatible software (e.g. Okular and Adobe Acrobat, but not Preview or browser viewers). The document may need to be 'trusted' before the video plays. The same video is also available as an ancillary file. The complete collection of supplementary materials are available from the publisher at https://doi.org/10.1093/mnras/stad868

We investigate the effect of observational constraints such as signal-to-noise, resolution and column density level on the HI morphological asymmetry ($\mathrm{A}_\mathrm{mod}$) and the effect of noise on the HI global profile ($\mathrm{A}_\mathrm{flux}$) asymmetry indices. Using mock galaxies from the EAGLE simulations we find an optimal combination of the observational constraints that are required for robust measurement of the $\mathrm{A}_\mathrm{mod}$ value of a galaxy: a column density threshold of $5\times10^{19}cm^{-2}$ or lower at a minimal signal-to-noise of 3 and a galaxy resolved with at least 11 beams. We also use mock galaxies to investigate the effect of noise on the $\mathrm{A}_\mathrm{flux}$ values and conclude that a global profile with signal-to-noise ratio greater than 5.5 is required to achieve a robust measurement of asymmetry. We investigate the relation between $\mathrm{A}_\mathrm{mod}$ and $\mathrm{A}_\mathrm{flux}$ indices and find them to be uncorrelated which implies that $\mathrm{A}_\mathrm{flux}$ values cannot be used to predict morphological asymmetries in galaxies., Comment: 19 pages, 12 figures and 2 appendices. Accepted for publication in the Monthly Notices of the Royal Astronomical Soceity

The number density of extragalactic 21-cm radio sources as a function of their spectral line-widths -- the HI width function (HIWF) -- is a sensitive tracer of the dark matter halo mass function (HMF). The $\Lambda$ cold dark matter model predicts that the HMF should be identical everywhere provided it is sampled in sufficiently large volumes, implying that the same should be true of the HIWF. The ALFALFA 21-cm survey measured the HIWF in northern and southern Galactic fields and found a systematically higher number density in the north. At face value, this is in tension with theoretical predictions. We use the Sibelius-DARK N-body simulation and the semi-analytical galaxy formation model GALFORM to create a mock ALFALFA survey. We find that the offset in number density has two origins: the sensitivity of the survey is different in the two fields, which has not been correctly accounted for in previous measurements; and the $1/V_{\mathrm{eff}}$ algorithm used for completeness corrections does not fully account for biases arising from spatial clustering in the galaxy distribution. The latter is primarily driven by a foreground overdensity in the northern field within $30\,\mathrm{Mpc}$, but more distant structure also plays a role. We provide updated measurements of the ALFALFA HIWF (and HIMF) correcting for the variations in survey sensitivity. Only when systematic effects such as these are understood and corrected for can cosmological models be tested against the HIWF., Comment: MNRAS accepted

We use $\Lambda$CDM cosmological hydrodynamical simulations to explore the kinematics of gaseous discs in late-type dwarf galaxies. We create high-resolution 21-cm 'observations' of simulated dwarfs produced in two variations of the EAGLE galaxy formation model: one where supernova-driven gas flows redistribute dark matter and form constant-density central 'cores', and another where the central 'cusps' survive intact. We 'observe' each galaxy along multiple sight lines and derive a rotation curve for each observation using a conventional tilted-ring approach to model the gas kinematics. We find that the modelling process introduces systematic discrepancies between the recovered rotation curve and the actual circular velocity curve driven primarily by (i) non-circular gas orbits within the discs; (ii) the finite thickness of gaseous discs, which leads to overlap of different radii in projection; and (iii) departures from dynamical equilibrium. Dwarfs with dark matter cusps often appear to have a core, whilst the inverse error is less common. These effects naturally reproduce an observed trend which other models struggle to explain: late-type dwarfs with more steeply-rising rotation curves appear to be dark matter-dominated in the inner regions, whereas the opposite seems to hold in galaxies with core-like rotation curves. We conclude that if similar effects affect the rotation curves of observed dwarfs, a late-type dwarf population in which all galaxies have sizeable dark matter cores is most likely incompatible with current measurements., Comment: 21 pages, 13 figures, 1 table. Accepted for publication in MNRAS. Appendices included as ancillary file

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) counterparts, but end up with ∼25 per cent less mass by the present day. One in three SIDM galaxies is 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}\!-\!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 ∼32 cluster galaxies.

Measurements of the rotation curves of dwarf galaxies are often interpreted as requiring a constant density core at the centre, at odds with the "cuspy" inner profiles predicted by $N$-body simulations of cold dark matter (CDM) haloes. It has been suggested that this conflict could be resolved by fluctuations in the inner gravitational potential caused by the periodic removal of gas following bursts of star formation. Earlier work has suggested that core formation requires a bursty and extended star formation history (SFH). Here we investigate the structure of CDM haloes of dwarf galaxies ($M_{{\rm DM}} \sim 10^9-5\times10^{10}\,{\rm M}_\odot$) formed in the APOSTLE ('A Project of Simulating the Local Environment') and AURIGA cosmological hydrodynamic simulations. Our simulations have comparable or better resolution than others that make cores ($M_{{\rm gas}} \sim 10^4\,{\rm M}_\odot$, gravitational softening $\sim 150$ pc). Yet, we do not find evidence of core formation at {\it any} mass or any correlation between the inner slope of the DM density profile and temporal variations in the SFH. APOSTLE and AURIGA dwarfs display a similar diversity in their cumulative SFHs to available data for Local Group dwarfs. Dwarfs in both simulations are DM-dominated on all resolved scales at all times, likely limiting the ability of gas outflows to alter significantly the central density profiles of their haloes. We conclude that recurrent bursts of star formation are not sufficient to cause the formation of cores, and that other conditions must also be met for baryons to be able to modify the central DM cusp., Comment: 15 pages, 9 figures, 2 tables. Updated version accepted for publication in MNRAS

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 Verlinde's emergent gravity. We find that the measured RAR agrees well with the MG predictions. In addition, we find a difference of at least $6\sigma$ 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., Comment: 31 pages, 15 figures, published in A&A

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 use 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., Comment: Accepted for publication in ApJ Letters after minor changes to text

We determine the dark matter pair-wise relative velocity distribution in a set of Milky Way-like halos in the Auriga and APOSTLE simulations. Focusing on the smooth halo component, the relative velocity distribution is well-described by a Maxwell-Boltzmann distribution over nearly all radii in the halo. We explore the implications for velocity-dependent dark matter annihilation, focusing on four models which scale as different powers of the relative velocity: Sommerfeld, s-wave, p-wave, and d-wave models. We show that the J-factors scale as the moments of the relative velocity distribution, and that the halo-to-halo scatter is largest for d-wave, and smallest for Sommerfeld models. The J-factor is strongly correlated with the dark matter density in the halo, and is very weakly correlated with the velocity dispersion. This implies that if the dark matter density in the Milky Way can be robustly determined, one can accurately predict the dark matter annihilation signal, without the need to identify the dark matter velocity distribution in the Galaxy., 25 pages, 11 figures, 2 tables; added figure 11; JCAP accepted version

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 located 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., 12 pages, 9 figures. Accepted for publication in MNRAS

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 consider 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., MNRAS accepted version. 25 pages, 11 figures. Appendices included in the source tarball (Download: Other formats: source)

We test a non-parametric higher-order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from the APOSTLE suite of cosmological hydrodynamics simulations with Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We find that, for CDM dwarfs, the recovered enclosed mass profiles have a bias of no more than 10 per cent, with a 50 per cent scatter in the inner regions and a 20 per cent scatter near the half-light radius, consistent with standard mass estimators. The density profiles are also recovered with a bias of no more than 10 per cent and a scatter of 30 per cent in the inner regions. For SIDM dwarfs, the mass and density profiles are recovered within our 95 per cent confidence intervals, but are biased towards cuspy dark matter distributions. This is mainly due to a lack of sufficient constraints from the data. We explore the sources of scatter in the accuracy of the recovered profiles and suggest a $\chi^2$ statistic to separate successful models from biased ones. Finally, we show that the uncertainties on the mass profiles obtained with GravSphere are smaller than those for comparable Jeans methods, and that they can be further improved if stronger priors, motivated by cosmological simulations, are placed on the velocity anisotropy. We conclude that GravSphere is a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies., Comment: 22 pages, 13 figures. Accepted in MNRAS

We use cosmological simulations of isolated Milky Way-mass galaxies, as well as Local Group analogues, to define the "edge" -- a caustic manifested in a drop in density or radial velocity -- of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ~1.4r_200m corresponding to the "splashback" radius, and a second caustic located at ~0.6r_200m which likely corresponds to the edge of the virialized material which has completed at least two pericentric passages. The splashback radius is ill defined in Local Group type environments where the halos of the two galaxies overlap. However, the second caustic is less affected by the presence of a companion, and is a more useful definition for the boundary of the Milky Way halo. Curiously, the stellar distribution also has a clearly defined caustic, which, in most cases, coincides with the second caustic of the dark matter. This can be identified in both radial density and radial velocity profiles, and should be measurable in future observational programmes. Finally, we show that the second caustic can also be identified in the phase-space distribution of dwarf galaxies in the Local Group. Using the current dwarf galaxy population, we predict the edge of the Milky Way halo to be 292 +/- 61 kpc., Comment: 14 pages, 11 figures. MNRAS in press

We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the Gaia DR2 Galactic rotation curve and other data. Using various hydrodynamical simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, r, Comment: 21 pages, 14 figures, main results: figs 3, 7, 10 and 11; MNRAS accepted version, code and resources available at https://github.com/MariusCautun/Milky_Way_mass_profile

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 a baryon fraction similar to the cosmological mean. By means of a 3D kinematic modelling fitting technique, we show that the H I 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 exponential disc scale length valid for dwarf galaxies with circular speeds ≲ 45 km s−1, 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., Monthly Notices of the Royal Astronomical Society, 495 (4), ISSN:0035-8711, ISSN:1365-2966, ISSN:1365-8711

N-body simulations make unambiguous predictions for the abundance of substructures within dark matter halos. However, the inclusion of baryons in the simulations changes the picture because processes associated with the presence of a large galaxy in the halo can destroy subhalos and substantially alter the mass function and velocity distribution of subhalos. We compare the effect of galaxy formation on subhalo populations in two state-of-the-art sets of hydrodynamical CDM simulations of Milky Way mass halos, APOSTLE and AURIGA. We introduce a new method for tracking the orbits of subhalos between simulation snapshots that gives accurate results down to a few kiloparsecs from the centre of the halo. Relative to a dark matter-only simulation, the abundance of subhalos in APOSTLE is reduced by 50% near the centre and by 10% within r200. In AURIGA the corresponding numbers are 80% and 40%. The velocity distributions of subhalos are also affected by the presence of the galaxy, much more so in AURIGA than in APOSTLE . The differences on subhalo properties in the two simulations can be traced back to the mass of the central galaxies, which in AURIGA are typically twice as massive as those in APOSTLE . We show that some of the results from previous studies are inaccurate due to systematic errors in the modelling of subhalo orbits near the centre of halos., Submitted to MNRAS

We investigated the disc-halo connection in massive (M⋆ > 5 × 1010 M⊙) disc galaxies from the cosmological hydrodynamical simulations EAGLE and IllustrisTNG, and compared this connection with that inferred from the study of H I rotation curves in nearby massive spirals from the Spitzer Photometry and Accurate Rotation Curves dataset. We find that discrepancies between the simulated and observed discs arise both on global and local scales. Globally, the simulated discs inhabit halos that are a factor ∼4 (in EAGLE) and ∼2 (in IllustrisTNG) more massive than those derived from the rotation curve analysis of the observed dataset. We also used synthetic rotation curves of the simulated discs to demonstrate that the recovery of the halo masses from rotation curves are not systematically biased. We find that the simulations predict systems dominated by dark matter with stellar-to-total enclosed mass ratios that are a factor of 1.5−2 smaller than real galaxies at all radii. This is an alternative manifestation of the ‘failed feedback problem,’ since it indicates that simulated halos hosting massive discs have been too inefficient at converting their baryons into stars, possibly owing to an overly efficient stellar and/or AGN feedback implementation.

We study the shape and kinematics of simulated dwarf galaxy discs in the APOSTLE suite of $\Lambda$CDM cosmological hydrodynamical simulations. We find that a large fraction of these gas-rich, star-forming discs show weak bars in their stellar component, despite being dark matter-dominated systems. The bar pattern shape and orientation reflect the ellipticity of the dark matter potential, and its rotation is locked to the slow figure rotation of the triaxial dark halo. The bar-like nature of the potential induces non-circular motions in the gas component, including strong bisymmetric flows that can be readily seen as m=3 harmonic perturbations in the HI line-of-sight velocity fields. Similar bisymmetric flows are seen in many galaxies of the THINGS and LITTLE THINGS surveys, although on average their amplitudes are a factor of ~2 weaker than in our simulated discs. Our results indicate that bar-like patterns may arise even when baryons are not dominant, and that they are common enough to warrant careful consideration when analyzing the gas kinematics of dwarf galaxy discs., Comment: 10 pages, 8 figures, submitted to MNRAS

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., To appear in the proceedings of IAUS 334: Rediscovering our Galaxy, July 10-14 2017, Telegrafenberg, Potsdam, Germany, Eds. C. Chiappini, I. Minchev, E. Starkenburg & M. Valentini

The predicted abundance and properties of the low-mass substructures embedded inside larger dark matter haloes differ sharply among alternative dark matter models. Too small to host galaxies themselves, these subhaloes may still be detected via gravitational lensing, or via perturbations of the Milky Way's globular cluster streams and its stellar disk. Here we use the Apostle cosmological simulations to predict the abundance and the spatial and velocity distributions of subhaloes in the range 10^6.5-10^8.5 solar masses inside haloes of mass ~ 10^12 solar masses in LCDM. Although these subhaloes are themselves devoid of baryons, we find that baryonic effects are important. Compared to corresponding dark matter only simulations, the loss of baryons from subhaloes and stronger tidal disruption due to the presence of baryons near the centre of the main halo, reduce the number of subhaloes by ~ 1/4 to 1/2, independently of subhalo mass, but increasingly towards the host halo centre. We also find that subhaloes have non-Maxwellian orbital velocity distributions, with centrally rising velocity anisotropy and positive velocity bias which reduces the number of low-velocity subhaloes, particularly near the halo centre. We parameterise the predicted population of subhaloes in terms of mass, galactocentric distance, and velocities. We discuss implications of our results for the prospects of detecting dark matter substructures and for possible inferences about the nature of dark matter., Comment: 19 pages, submitted to MNRAS, comments welcome

The Ophiuchus stream is a short arc-like stellar feature of uncertain origin located $\sim 5$ kpc North of the Galactic centre. New proper motions from the second $Gaia$ data release reconcile the direction of motion of stream members with the stream arc, resolving a puzzling mismatch reported in earlier work. We use N-body simulations to show that the stream is likely only on its second pericentric passage, and thus was formed recently. The simulations suggest that the entire disrupted progenitor is visible in the observed stream today, and that little further tidal debris lies beyond the ends of the stream. The luminosity, length, width, and velocity dispersion of the stream suggest a globular cluster (GC) progenitor substantially fainter and of lower surface brightness than estimated in previous work, and unlike any other known globulars in the Galaxy. This result suggests the existence of clusters that would extend the known GC population to fainter and more weakly bound systems than hitherto known. How such a weakly-bound cluster of old stars survived until it was disrupted so recently, however, remains a mystery. Integrating backwards in time, we find that the orbits of Sagittarius and Ophiuchus passed within $\sim 5$ kpc of each other about $\sim 100$ Myrs ago, an interaction that might help resolve this puzzle., 11 pages, 7 figures, submitted to MNRAS

We use integral field spectroscopy from the SAMI Galaxy Survey to identify galaxies that show evidence for recent quenching of star formation. The galaxies exhibit strong Balmer absorption in the absence of ongoing star formation in more than 10% of their spectra within the SAMI field of view. These $\rm{H}{\delta}$-strong galaxies (HDSGs) are rare, making up only $\sim 2$% (25/1220) of galaxies with stellar mass ${\rm log(}M_*/M_{\odot})>10$. The HDSGs make up a significant fraction of non-passive cluster galaxies (15%; 17/115) and a smaller fraction (2.0%; 8/387) of the non-passive population in low-density environments. The majority (9/17) of cluster HDSGs show evidence for star formation at their centers, with the HDS regions found in the outer parts of the galaxy. Conversely, the $\rm{H}{\delta}$-strong signal is more evenly spread across the galaxy for the majority (6/8) of HDSGs in low-density environments, and is often associated with emission lines that are not due to star formation. We investigate the location of the HDSGs in the clusters, finding that they are exclusively within 0.6$R_{200}$ of the cluster centre, and have a significantly higher velocity dispersion relative to the cluster population. Comparing their distribution in projected-phase-space to those derived from cosmological simulations indicates that the cluster HDSGs are consistent with an infalling population that have entered the central 0.5$r_{200, 3D}$ cluster region within the last $\sim 1\,$Gyr. In the 8/9 cluster HDSGs with central star formation, the extent of star formation is consistent with that expected of outside-in quenching by ram-pressure stripping. Our results indicate that the cluster HDSGs are currently being quenched by ram-pressure stripping on their first passage through the cluster., Comment: 44 pages (including 8 pages of appendices), 18 figures, 3 tables. Accepted for publicatioin in ApJ

Accurately predicting the shape of the HI velocity function of galaxies is regarded widely as a fundamental test of any viable dark matter model. Straightforward analyses of cosmological $N$-body simulations imply that the $\Lambda$CDM model predicts an overabundance of low circular velocity galaxies when compared to observed HI velocity functions. More nuanced analyses that account for the relationship between galaxies and their host haloes suggest that how we model the influence of baryonic processes has a significant impact on HI velocity function predictions. We explore this in detail by modelling HI emission lines of galaxies in the SHARK semi-analytic galaxy formation model, built on the SURFS suite of $\Lambda$CDM $N$-body simulations. We create a simulated ALFALFA survey, in which we apply the survey selection function and account for effects such as beam confusion, and compare simulated and observed HI velocity width distributions, finding differences of $\lesssim 50$%, orders of magnitude smaller than the discrepancies reported in the past. This is a direct consequence of our careful treatment of survey selection effects and, importantly, how we model the relationship between galaxy and halo circular velocity - the HI mass-maximum circular velocity relation of galaxies is characterised by a large scatter. These biases are complex enough that building a velocity function from the observed HI line widths cannot be done reliably., Comment: accepted to MNRAS 19 pages, 14 figures

We use mock interferometric HI measurements and a conventional tilted-ring modelling procedure to estimate circular velocity curves of dwarf galaxy discs from the APOSTLE suite of {\Lambda}CDM cosmological hydrodynamical simulations. The modelling yields a large diversity of rotation curves for an individual galaxy at fixed inclination, depending on the line-of-sight orientation. The diversity is driven by non-circular motions in the gas; in particular, by strong bisymmetric fluctuations in the azimuthal velocities that the tilted-ring model is ill-suited to account for and that are difficult to detect in model residuals. Large misestimates of the circular velocity arise when the kinematic major axis coincides with the extrema of the fluctuation pattern, in some cases mimicking the presence of kiloparsec-scale density 'cores', when none are actually present. The thickness of APOSTLE discs compounds this effect: more slowly-rotating extra-planar gas systematically reduces the average line-of-sight speeds. The recovered rotation curves thus tend to underestimate the true circular velocity of APOSTLE galaxies in the inner regions. Non-circular motions provide an appealing explanation for the large apparent cores observed in galaxies such as DDO 47 and DDO 87, where the model residuals suggest that such motions might have affected estimates of the inner circular velocities. Although residuals from tilted ring models in the simulations appear larger than in observed galaxies, our results suggest that non-circular motions should be carefully taken into account when considering the evidence for dark matter cores in individual galaxies., Comment: 29 pages, 19 figures, 4 tables, supplementary materials available in arXiv tarball. MNRAS accepted version

We study the gas kinematics traced by the 21-cm emission of a sample of six HI$-$rich low surface brightness galaxies classified as ultra-diffuse galaxies (UDGs). Using the 3D kinematic modelling code $\mathrm{^{3D}}$Barolo we derive robust circular velocities, revealing a startling feature: HI$-$rich UDGs are clear outliers from the baryonic Tully-Fisher relation, with circular velocities much lower than galaxies with similar baryonic mass. Notably, the baryon fraction of our UDG sample is consistent with the cosmological value: these UDGs are compatible with having no "missing baryons" within their virial radii. Moreover, the gravitational potential provided by the baryons is sufficient to account for the amplitude of the rotation curve out to the outermost measured point, contrary to other galaxies with similar circular velocities. We speculate that any formation scenario for these objects will require very inefficient feedback and a broad diversity in their inner dark matter content., Accepted for publication by The Astrophysical Journal Letters (ApJL). V2: Acknowledgments have been updated, and a typo has been corrected

A number of Local Group dwarf galaxies are known to have two spatially segregated stellar metallicity populations, a centrally concentrated metal-rich population and a more extended metal-poor population. In this work we discuss mechanisms that lead to the formation of two spatially segregated metallicity populations. Using a set of high-resolution hydrodynamical simulations of Local Group-like environments, we select a sample of satellite and field galaxies, spanning the stellar mass range $10^6-10^9$M$_{\odot}$, that exhibit bimodality in their metallicity distributions. Among those, we identify a subsample with a strong spatial segregation in the two populations. We find three distinct mechanisms for their formation. In field dwarfs and in a small fraction of satellites, a merger causes the metal-poor stars to migrate to larger radii and encourages the available gas to sink to the centre of the dwarf. Most of the gas is subsequently blown out of the halo through star formation feedback, but the remaining gas is consumed in the formation of a metal-rich population. In the exclusive case of satellites that have retained some of their gas at infall, it is the compression of this gas by ram pressure near pericentre that triggers the formation of metal-rich stars, whilst simultaneously preventing star formation at larger radii through stripping. Additionally, in a small number of field and satellite dwarfs, interactions with gaseous filaments and other galaxies can result in the formation of a metal-rich population. Regardless of the formation mechanism, a history of mergers typically enhances the spatial segregation., Comment: 20 pages, 19 figures, accepted in MNRAS

We examine the origin of radial and vertical gradients in the age/metallicity of the stellar component of a galaxy disc formed in the APOSTLE cosmological hydrody- namical simulations. Some of these gradients resemble those in the Milky Way, where they have sometimes been interpreted as due to internal evolution, such as scattering off giant molecular clouds, radial migration driven by spiral patterns, or orbital reso- nances with a bar. Secular processes play a minor role in the simulated galaxy, which lacks strong spiral or bar patterns, and where such gradients arise as a result of the gradual enrichment of a gaseous disc that is born thick but thins as it turns into stars and settles into centrifugal equilibrium. The settling is controlled by the feedback of young stars; which links the star formation, enrichment, and equilibration timescales, inducing radial and vertical gradients in the gaseous disc and its descendent stars. The kinematics of coeval stars evolve little after birth and provide a faithful snapshot of the gaseous disc structure at the time of their formation. In this interpretation, the age-velocity dispersion relation would reflect the gradual thinning of the disc rather than the importance of secular orbit scattering; the outward flaring of stars would result from the gas disc flare rather than from radial migration; and vertical gradients would arise because the gas disc gradually thinned as it enriched. Such radial and vertical trends might just reflect the evolving properties of the parent gaseous disc, and are not necessarily the result of secular evolutionary processes., Resubmitted to MNRAS

The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Penarrubia have used this approach to argue that data for these galaxies are consistent with constant-density `cores' in their inner regions and rule out `cuspy' Navarro-Frenk-White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE (A Project Of Simulating The Local Environment) Lambda cold dark matter cosmological hydrodynamics simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Penarrubia to a sample of 53 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed., 15 pages, 13 figures, accepted for publication in MNRAS

The shallow faint-end slope of the galaxy mass function is usually reproduced in $\Lambda$CDM galaxy formation models by assuming that the fraction of baryons that turns into stars drops steeply with decreasing halo mass and essentially vanishes in haloes with maximum circular velocities $V_{\rm max}, Comment: replaced to match the accepted version. 2 new figures added compared to the previous version

We use the APOSTLE (A Project Of Simulating The Local Environment) cosmological hydrodynamic simulations to examine the effects of tidal stripping on cold dark matter subhaloes that host three of the most luminous Milky Way dwarf satellite galaxies: Fornax, Sculptor and Leo I. We identify simulated satellites that match the observed spatial and kinematic distributions of stars in these galaxies, and track their evolution after infall. We find ∼30 per cent of subhaloes hosting satellites with present-day stellar mass 106–108 M⊙ experience >20 per cent stellar mass-loss after infall. Fornax analogues have earlier infall times compared to Sculptor and Leo I analogues. Star formation in Fornax analogues continues for ∼3–6 Gyr after infall, whereas Sculptor and Leo I analogues stop forming stars

We examine the baryon content of low-mass A cold dark matter (ACDM) haloes (10(8)

We examine the spatial distribution of the oldest and most metal poor stellar populations of Milky Way-sized galaxies using the APOSTLE cosmological hydrodynamical simulations of the Local Group. In agreement with earlier work, we find strong radial gradients in the fraction of the oldest (tform < 0.8 Gyr) and most metal poor ([Fe/H]< -2.5) stars, both of which increase outwards. The most metal poor stars form over an extended period of time; half of them form after z = 5.3, and the last 10% after z = 2.8. The age of the metal poor stellar population also shows significant variation with environment; a high fraction of them are old in the galaxy's central regions and an even higher fraction in some individual dwarf galaxies, with substantial scatter from dwarf to dwarf. Overall, over half of the stars that belong to both the oldest and most metal-poor population are found outside the solar circle. Somewhat counter-intuitively, we find that dwarf galaxies with a large fraction of metal poor stars that are very old are systems where metal poor stars are relatively rare, but where a substantial old population is present. Our results provide guidance for interpreting the results of surveys designed to hunt for the earliest and most pristine stellar component of our Milky Way., 14 pages, 9 figures. Submitted to MNRAS, revised after referee comments

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

We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum aceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in $\Lambda$CDM. This is because: (i) disk galaxies in $\Lambda$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $\Lambda$CDM galaxy formation model. Explaining the MDAR in $\Lambda$CDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process., Comment: 8 pages, 3 figures, MNRAS submitted

In order to understand the processes that quench star formation within rich clusters, we construct a library of subhalo orbits drawn from $\Lambda$CDM cosmological N-body simulations of four rich clusters. The orbits are combined with models of star formation followed by quenching in the cluster environment. These are compared with observed bulge and disc colours and stellar absorption linestrength indices of satellite galaxies. Models in which the bulge stellar populations depend only on the galaxy subhalo mass while the disc quenching depends on the cluster environment are acceptable fits to the data. An exponential disc quenching timescale of 3 - 3.5 Gyr is preferred. Models with short ($\lesssim 1$ Gyr) quenching timescales yield cluster-centric gradients in disc colours and Balmer line indices that are too steep compared to observations. We also examine models in which there is quenching in lower mass groups prior to cluster infall ("pre-processing"), finding that such models are a better fit to the data than models without pre-processing and require similar quenching times. The data slightly prefer models where quenching occurs only for galaxies falling within about 0.5 $r_{200}$. Finally, we have examined models with short quenching timescales of 1 Gyr, but a long delay time of 3 Gyr prior to quenching. All models with short quenching timescales, even such "delayed-then-rapid" quenching models, produce excessively red galaxies near the cluster core and are strongly disfavoured by the data. These results imply that the environments of rich clusters must impact star formation rates of infalling galaxies on relatively long timescales -- several times longer than a typical halo spends within the virial radius of a cluster. This scenario favours gentler quenching mechanisms such as slow "strangulation" over more rapid ram-pressure stripping., Comment: MNRAS submitted. 18 pages, 9 figures, 1 table

Physical processes regulating star formation in satellite galaxies represent an area of ongoing research, but the projected nature of observed coordinates makes separating different populations of satellites (with different processes at work) difficult. The orbital history of a satellite galaxy leads to its present-day phase space coordinates; we can also work backwards and use these coordinates to statistically infer information about the orbital history. We use merger trees from the MultiDark Run 1 N-body simulation to compile a catalog of the orbits of satellite haloes in cluster environments. We parameterize the orbital history by the time since crossing within 2.5 rvir of the cluster centre and use our catalog to estimate the probability density over a range of this parameter given a set of present-day projected (i.e. observable) phase space coordinates. We show that different populations of satellite haloes, e.g. infalling, backsplash and virialized, occupy distinct regions of phase space, and semi-distinct regions of projected phase space. This will allow us to probabilistically determine the time since infall of a large sample of observed satellite galaxies, and ultimately to study the effect of orbital history on star formation history (the topic of a future paper). We test the accuracy of our method and find that we can reliably recover this time within +/-2.58 Gyr in 68 per cent of cases by using all available phase space coordinate information, compared to +/-2.64 Gyr using only position coordinates and +/-3.10 Gyr guessing 'blindly', i.e. using no coordinate information, but with knowledge of the overall distribution of infall times. In some regions of phase space, the accuracy of the infall time estimate improves to +/-1.85 Gyr. Although we focus on time since infall, our method is easily generalizable to other orbital parameters (e.g. pericentric distance and time)., Comment: 10 pages, 9 figures, submitted to MNRAS

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

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

We measure the star formation quenching efficiency and timescale in cluster environments. Our method uses N-body simulations to estimate the probability distribution of possible orbits for a sample of observed SDSS galaxies in and around clusters based on their position and velocity offsets from their host cluster. We study the relationship between their star formation rates and their likely orbital histories via a simple model in which star formation is quenched once a delay time after infall has elapsed. Our orbit library method is designed to isolate the environmental effect on the star formation rate due to a galaxy's present-day host cluster from `pre-processing' in previous group hosts. We find that quenching of satellite galaxies of all stellar masses in our sample ($10^{9}-10^{11.5}\,{\rm M}_\odot$) by massive ($> 10^{13}\,{\rm M}_\odot$) clusters is essentially $100$ per cent efficient. Our fits show that all galaxies quench on their first infall, approximately at or within a Gyr of their first pericentric passage. There is little variation in the onset of quenching from galaxy-to-galaxy: the spread in this time is at most $\sim 2$ Gyr at fixed $M_*$. Higher mass satellites quench earlier, with very little dependence on host cluster mass in the range probed by our sample., 14 pages, 10 figures, 1 table. MNRAS submitted. Comments welcome

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

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