Researcher profile

Andrew Robertson

Andrew Robertson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2026arXiv

Accelerated calibration of semi-analytic galaxy formation models

We present an accelerated calibration framework for semi-analytic galaxy formation models, demonstrated with Galacticus. Rather than fitting directly to properties such as the low-redshift stellar mass function (SMF) - which requires evolving thousands of halos per likelihood evaluation - we construct a fast likelihood from the stellar-to-halo mass relation (SHMR; mean and scatter) evaluated at a small set of target halo masses, reducing each evaluation to simulating only tens of galaxies. We sample the posterior over Galacticus parameters with Markov Chain Monte Carlo and show that the resulting calibration reproduces the low-redshift SMF. We then extend the method to additional datasets, using a higher-redshift SHMR and the low-redshift stellar mass-size relation as examples, and assess performance for large scale structure survey-relevant properties: stellar masses, sizes, and emission-line strengths. The SMF matches data well at low redshift, but toward higher redshift the model yields too few massive galaxies and too many low-mass galaxies. Size evolution with redshift is approximately correct, but the mass-size relation is too flat, producing massive galaxies that are too small. The H$α$ luminosity function is well reproduced at z~2, but by z~0.4 the model overproduces highly star-forming, H$α$-bright systems. These discrepancies suggest the model lacks sufficient flexibility (e.g. in gas cooling/recycling or feedback) to reconcile all datasets simultaneously. Our strategy complements emulator-based methods for calibrating semi-analytic models by enabling rapid, low-cost scans of model choices and parameterisations - a capability we envision leveraging to supply calibrated starting points for more detailed follow-up inference.

preprint2023arXiv

Automated galaxy-galaxy strong lens modelling: no lens left behind

The distribution of dark and luminous matter can be mapped around galaxies that gravitationally lens background objects into arcs or Einstein rings. New surveys will soon observe hundreds of thousands of galaxy lenses, and current, labour-intensive analysis methods will not scale up to this challenge. We instead develop a fully automatic, Bayesian method which we use to fit a sample of 59 lenses imaged by the Hubble Space Telescope in uniform conditions. We set out to \textit{leave no lens behind} and focus on ways in which automated fits fail in a small handful of lenses, describing adjustments to the pipeline that allows us to infer accurate lens models. Our pipeline ultimately fits {\em all} 59 lenses in our sample, with a high success rate key because catastrophic outliers would bias large samples with small statistical errors. Machine Learning techniques might further improve the two most difficult steps: subtracting foreground lens light and initialising a first, approximate lens model. After that, increasing model complexity is straightforward. We find a mean $\sim1\%$ measurement precision on the measurement of the Einstein radius across the lens sample which {\em does not degrade with redshift} up to at least $z=0.7$ -- in stark contrast to other techniques used to study galaxy evolution, like stellar dynamics. Our \texttt{PyAutoLens} software is open source, and is also installed in the Science Data Centres of the ESA Euclid mission.

preprint2022arXiv

A forward-modelling method to infer the dark matter particle mass from strong gravitational lenses

A fundamental prediction of the cold dark matter (CDM) model of structure formation is the existence of a vast population of dark matter haloes extending to subsolar masses. By contrast, other dark matter models, such as a warm thermal relic (WDM), predict a cutoff in the mass function at a mass which, for popular models, lies approximately between $10^7$ and $10^{10}~{\rm M}_\odot$. We use mock observations to demonstrate the viability of a forward modelling approach to extract information about low-mass dark haloes lying along the line-of-sight to galaxy-galaxy strong lenses. This can be used to constrain the mass of a thermal relic dark matter particle, $m_\mathrm{DM}$. With 50 strong lenses at Hubble Space Telescope resolution and a maximum pixel signal-to-noise ratio of $\sim50$, the expected median 2$σ$ constraint for a CDM-like model (with a halo mass cutoff at $10^{7}~{\rm M}_\odot$) is $m_\mathrm{DM} > 4.10 \, \mathrm{keV}$ (50% chance of constraining $m_{\rm DM}$ to be better than 4.10 keV). If, however, the dark matter is a warm particle of $m_\mathrm{DM}=2.2 \, \mathrm{keV}$, our 'Approximate Bayesian Computation' method would result in a median estimate of $m_\mathrm{DM}$ between 1.43 and 3.21 keV. Our method can be extended to the large samples of strong lenses that will be observed by future telescopes, and could potentially rule out the standard CDM model of cosmogony. To aid future survey design, we quantify how these constraints will depend on data quality (spatial resolution and integration time) as well as on the lensing geometry (source and lens redshifts).

preprint2022arXiv

Cosmological simulations with rare and frequent dark matter self-interactions

Dark matter (DM) with self-interactions is a promising solution for the small-scale problems of the standard cosmological model. Here we perform the first cosmological simulation of frequent DM self-interactions, corresponding to small-angle DM scatterings. The focus of our analysis lies in finding and understanding differences to the traditionally assumed rare DM (large-angle) self scatterings. For this purpose, we compute the distribution of DM densities, the matter power spectrum, the two-point correlation function and the halo and subhalo mass functions. Furthermore, we investigate the density profiles of the DM haloes and their shapes. We find that overall large-angle and small-angle scatterings behave fairly similarly with a few exceptions. In particular, the number of satellites is considerably suppressed for frequent compared to rare self-interactions with the same cross-section. Overall we observe that while differences between the two cases may be difficult to establish using a single measure, the degeneracy may be broken through a combination of multiple ones. For instance, the combination of satellite counts with halo density or shape profiles could allow discriminating between rare and frequent self-interactions. As a by-product of our analysis, we provide - for the first time - upper limits on the cross-section for frequent self-interactions.

preprint2022arXiv

Unequal-mass mergers of dark matter haloes with rare and frequent self-interactions

Dark matter (DM) self-interactions have been proposed to solve problems on small length scales within the standard cold DM cosmology. Here, we investigate the effects of DM self-interactions in merging systems of galaxies and galaxy clusters with equal and unequal mass ratios. We perform N-body DM-only simulations of idealized setups to study the effects of DM self-interactions that are elastic and velocity-independent. We go beyond the commonly adopted assumption of large-angle (rare) DM scatterings, paying attention to the impact of small-angle (frequent) scatterings on astrophysical observables and related quantities. Specifically, we focus on DM-galaxy offsets, galaxy--galaxy distances, halo shapes, morphology, and the phase--space distribution. Moreover, we compare two methods to identify peaks: one based on the gravitational potential and one based on isodensity contours. We find that the results are sensitive to the peak finding method, which poses a challenge for the analysis of merging systems in simulations and observations, especially for minor mergers. Large DM-galaxy offsets can occur in minor mergers, especially with frequent self-interactions. The subhalo tends to dissolve quickly for these cases. While clusters in late merger phases lead to potentially large differences between rare and frequent scatterings, we believe that these differences are non-trivial to extract from observations. We therefore study the galaxy/star populations which remain distinct even after the DM haloes have coalesced. We find that these collisionless tracers behave differently for rare and frequent scatterings, potentially giving a handle to learn about the micro-physics of DM.

preprint2021arXiv

A high-resolution cosmological simulation of a strong gravitational lens

We present a cosmological hydrodynamical simulation of a 10^13 Msun galaxy group and its environment (out to 10 times the virial radius) carried out using the EAGLE model of galaxy formation. Exploiting a novel technique to increase the resolution of the dark matter calculation independently of that of the gas, the simulation resolves dark matter haloes and subhaloes of mass 5x10^6 Msun . It is therefore useful for studying the abundance and properties of the haloes and subhaloes targeted in strong lensing tests of the cold dark matter model. We estimate the halo and subhalo mass functions and discuss how they are affected both by the inclusion of baryons in the simulation and by the environment. We find that the halo and subhalo mass functions have lower amplitude in the hydrodynamical simulation than in its dark matter only counterpart. This reflects the reduced growth of haloes in the hydrodynamical simulation due to the early loss of gas by reionisation and galactic winds and, additionally, in the case of subhaloes, disruption by enhanced tidal effects within the host halo due to the presence of a massive central galaxy. The distribution of haloes is highly anisotropic reflecting the filamentary character of mass accretion onto the cluster. As a result, there is significant variation in the number of structures with viewing direction. The median number of structures near the centre of the halo, when viewed in projection, is reduced by a factor of two when baryons are included.

preprint2021arXiv

From dwarf galaxies to galaxy clusters: Self-Interacting Dark Matter over 7 orders of magnitude in halo mass

In this paper we study the density profiles of self-interacting dark matter (SIDM) haloes spanning the full observable mass range, from dwarf galaxies to galaxy clusters. Using realistic simulations that model the baryonic physics relevant for galaxy formation, we compare the density profiles of haloes simulated with either SIDM or cold and collisionless dark matter (CDM) to those inferred from observations of stellar velocity dispersion, gas rotation curves, weak and strong gravitational lensing, and/or X-ray maps. We make our comparison in terms of the maximal surface density of haloes, circumventing the need for semi-analytic or parametric models for dark matter density profiles. We find that the maximal surface density as a function of halo mass is well reproduced by CDM simulations that include baryons, while for SIDM with a velocity-independent cross-section of 1 cm$^2/$g, the simulated galaxy clusters have mean maximal surface densities that are below those of observed systems by an amount greater than the standard deviation of the observed maximal surface density at fixed mass. For less massive systems both CDM and SIDM agree with the observation equally well.

preprint2021arXiv

Systematic errors induced by the elliptical power-law model in galaxy-galaxy strong lens modeling

The elliptical power-law (EPL) model of the mass in a galaxy is widely used in strong gravitational lensing analyses. However, the distribution of mass in real galaxies is more complex. We quantify the biases due to this model mismatch by simulating and then analysing mock {\it Hubble Space Telescope} imaging of lenses with mass distributions inferred from SDSS-MaNGA stellar dynamics data. We find accurate recovery of source galaxy morphology, except for a slight tendency to infer sources to be more compact than their true size. The Einstein radius of the lens is also robustly recovered with 0.1% accuracy, as is the global density slope, with 2.5% relative systematic error, compared to the 3.4% intrinsic dispersion. However, asymmetry in real lenses also leads to a spurious fitted `external shear' with typical strength, $γ_{\rm ext}=0.015$. Furthermore, time delays inferred from lens modelling without measurements of stellar dynamics are typically underestimated by $\sim$5%. Using such measurements from a sub-sample of 37 lenses would bias measurements of the Hubble constant $H_0$ by $\sim$9%. Although this work is based on a particular set of MaNGA galaxies, and the specific value of the detected biases may change for another set of strong lenses, our results strongly suggest the next generation cosmography needs to use more complex lens mass models.

preprint2021arXiv

The effects of self-interacting dark matter on the stripping of galaxies that fall into clusters

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 ${\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.

preprint2021arXiv

The surprising accuracy of isothermal Jeans modelling of self-interacting dark matter density profiles

Recent claims of observational evidence for self-interacting dark matter (SIDM) have relied on a semi-analytic method for predicting the density profiles of galaxies and galaxy clusters containing SIDM. We present a thorough description of this method, known as isothermal Jeans modelling, and then test it with a large ensemble of haloes taken from cosmological simulations. Our simulations were run with cold and collisionless dark matter (CDM) as well as two different SIDM models, all with dark matter only variants as well as versions including baryons and relevant galaxy formation physics. Using a mix of different box sizes and resolutions, we study haloes with masses ranging from 3e10 to 3e15 Msun. Overall, we find that the isothermal Jeans model provides as accurate a description of simulated SIDM density profiles as the Navarro-Frenk-White profile does of CDM halos. We can use the model predictions, compared with the simulated density profiles, to determine the input DM-DM scattering cross-sections used to run the simulations. This works especially well for large cross-sections, while with CDM our results tend to favour non-zero (albeit fairly small) cross-sections, driven by a bias against small cross-sections inherent to our adopted method of sampling the model parameter space. The model works across the whole halo mass range we study, although including baryons leads to DM profiles of intermediate-mass (10^12 - 10^13 Msun) haloes that do not depend strongly on the SIDM cross-section. The tightest constraints will therefore come from lower and higher mass haloes: dwarf galaxies and galaxy clusters.

preprint2020arXiv

Baryonic clues to the puzzling diversity of dwarf galaxy rotation curves

We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue.

preprint2020arXiv

Exploring extensions to the standard cosmological model and the impact of baryons on small scales

It has been claimed that the standard model of cosmology (LCDM) cannot easily account for a number of observations on relatively small scales, motivating extensions to the standard model. Here we introduce a new suite of cosmological simulations that systematically explores three plausible extensions: warm dark matter, self-interacting dark matter, and a running of the scalar spectral index of density fluctuations. Current observational constraints are used to specify the additional parameters that come with these extensions. We examine a large range of observable metrics on small scales, including the halo mass function, density and circular velocity profiles, the abundance of satellite subhaloes, and halo concentrations. For any given metric, significant degeneracies can be present between the extensions. In detail, however, the different extensions have quantitatively distinct mass and radial dependencies, suggesting that a multi-probe approach over a range of scales can be used to break the degeneracies. We also demonstrate that the relative effects on the radial density profiles in the different extensions (compared to the standard model) are converged down to significantly smaller radii than are the absolute profiles. We compare the derived cosmological trends with the impact of baryonic physics using the EAGLE and ARTEMIS simulations. Significant degeneracies are also present between baryonic physics and cosmological variations (with both having similar magnitude effects on some observables). Given the inherent uncertainties both in the modelling of galaxy formation physics and extensions to LCDM, a systematic and simultaneous exploration of both is strongly warranted.

preprint2020arXiv

Mapping dark matter and finding filaments: calibration of lensing analysis techniques on simulated data

We quantify the performance of mass mapping techniques on mock imaging and gravitational lensing data of galaxy clusters. The optimum method depends upon the scientific goal. We assess measurements of clusters' radial density profiles, departures from sphericity, and their filamentary attachment to the cosmic web. We find that mass maps produced by direct (KS93) inversion of shear measurements are unbiased, and that their noise can be suppressed via filtering with MRLens. Forward-fitting techniques, such as Lenstool, suppress noise further, but at a cost of biased ellipticity in the cluster core and over-estimation of mass at large radii. Interestingly, current searches for filaments are noise-limited by the intrinsic shapes of weakly lensed galaxies, rather than by the projection of line-of-sight structures. Therefore, space-based or balloon-based imaging surveys that resolve a high density of lensed galaxies, could soon detect one or two filaments around most clusters.

preprint2020arXiv

Observable tests of self-interacting dark matter in galaxy clusters: cosmological simulations with SIDM and baryons

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

preprint2020arXiv

On building a cluster watch-list for identifying strongly lensed supernovae, gravitational waves and kilonovae

Motivated by discovering strongly-lensed supernovae, gravitational waves, and kilonovae in the 2020s, we investigate whether to build a watch-list of clusters based on observed cluster properties (i.e. lens-plane selection) or on the detectability of strongly-lensed background galaxies (i.e. source-plane selection). First, we estimate the fraction of high-redshift transient progenitors that reside in galaxies that are themselves too faint to be detected as being strongly-lensed. We find $\sim15-50$ per cent of transient progenitors reside in $z = 1-2$ galaxies too faint to be detected in surveys that reach ${\rm AB}\simeq23$, such as the Dark Energy Survey. This falls to $\ls10$ per cent at depths that will be probed by early data releases of LSST (${\rm AB}\simeq25$). Second, we estimate a conservative lower limit on the fraction of strong lensing clusters that will be missed by magnitude limited searches for multiply-imaged galaxies and giant arcs due to the faintness of such images. We find that DES-like surveys will miss $\sim75$ per cent of $10^{15}$M$_\odot$ strong lensing clusters, rising to $\sim100$ per cent of $10^{14}$M$_\odot$ clusters. Deeper surveys, such as LSST, will miss $\sim40$ per cent at $10^{15}$M$_\odot$, and $\sim95$ per cent at $10^{14}$M$_\odot$. Our results motivate building a cluster watch-list for strongly-lensed transients that includes those found by lens-plane selection.

preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

preprint2020arXiv

The distribution of dark matter and gas spanning six megaparsecs around the post-merger galaxy cluster MS0451-03

Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions.

preprint2020arXiv

To beta or not to beta: can higher-order Jeans analysis break the mass-anisotropy degeneracy in simulated dwarfs?

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 $χ^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.

preprint2020arXiv

Understanding the large inferred Einstein radii of observed low-mass galaxy clusters

We assess a claim that observed galaxy clusters with mass ~10^14 Msun are more centrally concentrated than predicted in LCDM. We generate mock strong gravitational lensing observations, taking the lenses from a cosmological hydrodynamical simulation, and analyse them in the same way as the real Universe. The observed and simulated lensing arcs are consistent with one another, with three main effects responsible for the previously claimed inconsistency. First, galaxy clusters containing baryonic matter have higher central densities than their counterparts simulated with only dark matter. Second, a sample of clusters selected because of the presence of pronounced gravitational lensing arcs preferentially finds centrally concentrated clusters with large Einstein radii. Third, lensed arcs are usually straighter than critical curves, and the chosen image analysis method (fitting circles through the arcs) overestimates the Einstein radii. After accounting for these three effects, LCDM predicts that galaxy clusters should produce giant lensing arcs that match those in the observed Universe.

preprint2020arXiv

Velocity-dependent Self-interacting Dark Matter from Groups and Clusters of Galaxies

We probe the self-interactions of dark matter using observational data of relaxed galaxy groups and clusters. Our analysis uses the Jeans formalism and considers a wider range of systematic effects than in previous work, including adiabatic contraction and stellar anisotropy, to robustly constrain the self-interaction cross section. For both groups and clusters, our results show a mild preference for a nonzero cross section compared with cold collisionless dark matter. Our groups result, $σ/m=0.5\pm0.2~\mathrm{cm}^2/\mathrm{g}$, places the first constraint on self-interacting dark matter (SIDM) at an intermediate scale between galaxies and massive clusters. Our clusters result is $σ/m=0.19\pm0.09~\mathrm{cm}^2/\mathrm{g}$, with an upper limit of $σ/ m < 0.35~\mathrm{cm}^2/\mathrm{g}$ (95% CL). Thus, our results disfavor a velocity-independent cross section of order $1~\mathrm{cm}^2/\mathrm{g}$ or larger needed to address small scale structure problems in galaxies, but are consistent with a velocity-dependent cross section that decreases with increasing scattering velocity. Comparing the cross sections with and without the effect of adiabatic contraction, we find that adiabatic contraction produces slightly larger values for our data sample, but they are consistent at the $1σ$ level. Finally, to validate our approach, we apply our Jeans analysis to a sample of mock data generated from SIDM-plus-baryons simulations with $σ/m = 1~\mathrm{cm}^2/\mathrm{g}$. This is the first test of the Jeans model at the level of stellar and lensing observables directly measured from simulations. We find our analysis gives a robust determination of the cross section, as well as consistently inferring the true baryon and dark matter density profiles.

preprint2020arXiv

What does strong gravitational lensing? The mass and redshift distribution of high-magnification lenses

Many distant objects can only be detected, or become more scientifically valuable, if they have been highly magnified by strong gravitational lensing. We use EAGLE and BAHAMAS, two recent cosmological hydrodynamical simulations, to predict the probability distribution for both the lens mass and lens redshift when point sources are highly magnified by gravitational lensing. For sources at a redshift of two, we find the distribution of lens redshifts to be broad, peaking at z=0.6. The contribution of different lens masses is also fairly broad, with most high-magnification lensing due to lenses with halo masses between 10^12 and 10^14 solar masses. Lower mass haloes are inefficient lenses, while more massive haloes are rare. We find that a simple model in which all haloes have singular isothermal sphere density profiles can approximately reproduce the simulation predictions, although such a model over-predicts the importance of haloes with mass <10^12 solar masses for lensing. We also calculate the probability that point sources at different redshifts are strongly lensed. At low redshift, high magnifications are extremely unlikely. Each z=0.5 source produces, on average, 5x10^-7 images with magnification greater than ten; for z =2 this increases to about 2x10^-5. Our results imply that searches for strongly lensed optical transients, including the optical counterparts to strongly lensed gravitational waves, can be optimized by monitoring massive galaxies, groups and clusters rather than concentrating on an individual population of lenses.

preprint2019arXiv

Subhalo destruction in the Apostle and Auriga simulations

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.