Researcher profile

Andrew Benson

Andrew Benson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
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

11 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.

preprint2024arXiv

Constrain the Dark Matter Distribution of Ultra-diffuse Galaxies with Globular-Cluster Mass Segregation: A Case Study with NGC5846-UDG1

The properties of globular clusters (GCs) contain valuable information of their host galaxies and dark-matter halos. In the remarkable example of ultra-diffuse galaxy, NGC5846-UDG1, the GC population exhibits strong radial mass segregation, indicative of dynamical-friction-driven orbital decay, which opens the possibility of using imaging data alone to constrain the dark-matter content of the galaxy. To explore this possibility, we develop a semi-analytical model of GC evolution, which starts from the initial mass function, the initial structure-mass relation, and the initial spatial distribution of the GC progenitors, and follows the effects of dynamical friction, tidal evolution, and two-body relaxation. Using Markov Chain Monte Carlo, we forward-model the GCs in a NGC5846-UDG1-like potential to match the observed GC mass, size, and spatial distributions, and to constrain the profile of the host halo and the origin of the GCs. We find that, with the assumptions of zero mass segregation when the star clusters were born, NGC5846-UDG1 is relatively dark-matter poor compared to what is expected from stellar-to-halo-mass relations, and its halo concentration is lower than the cosmological average, irrespective of having a cuspy or a cored profile. Its GC population has an initial spatial distribution more extended than the smooth stellar distribution. We discuss the results in the context of scaling laws of galaxy-halo connections, and warn against naively using the GC-abundance-halo-mass relation to infer the halo mass of UDGs. Our model is generally applicable to GC-rich dwarf galaxies, and is publicly available at https://github.com/JiangFangzhou/GCevo.

preprint2022arXiv

Core-collapse, evaporation and tidal effects: the life story of a self-interacting dark matter subhalo

Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass halos, accelerated if halos are subhalos of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of $\lesssim 10^8 M_\odot$ subhalos in lens-mass systems ($\sim 10^{13}M_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhalos with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhalos ($\lesssim 1/1000$ host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhalos to core-collapse, finding that it is nearly impossible to drive core-collapse in subhalos in SIDM models with constant cross sections. Any discovery of ultra-compact dark substructures with future substructure lensing observations favors additional degrees of freedom, such as velocity-dependence, in the cross section.

preprint2022arXiv

The primordial matter power spectrum on sub-galactic scales

The primordial matter power spectrum quantifies fluctuations in the distribution of dark matter immediately following inflation. Over cosmic time, over-dense regions of the primordial density field grow and collapse into dark matter halos, whose abundance and density profiles retain memory of the initial conditions. By analyzing the image magnifications in eleven strongly-lensed and quadruply-imaged quasars, we infer the abundance and concentrations of low-mass halos, and cast the measurement in terms of the amplitude of the primordial matter power spectrum. We anchor the power spectrum on large scales, isolating the effect of small-scale deviations from the $Λ$CDM prediction. Assuming an analytic model for the power spectrum and accounting for several sources of potential systematic uncertainty, including three different models for the halo mass function, we obtain correlated inferences of $\log_{10}\left(P / P_{Λ\rm{CDM}}\right)$, the power spectrum amplitude relative to the predictions of the concordance cosmological model, of $0.0_{-0.4}^{+0.5}$, $0.1_{-0.6}^{+0.7}$, and $0.2_{-0.9}^{+1.0}$ at k = 10, 25 and 50 $\rm{Mpc^{-1}}$ at $68 \%$ confidence, consistent with cold dark matter and single-field slow-roll inflation.

preprint2021arXiv

The High Latitude Spectroscopic Survey on the Nancy Grace Roman Space Telescope

The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0-1.93 $μ$m, $R=435-865$) and the 0.28 deg$^2$ wide field camera. We present a reference HLSS which maps 2000 deg$^2$ and achieves an emission line flux limit of 10$^{-16}$ erg/s/cm$^2$ at 6.5$σ$, requiring $\sim$0.6 yrs of observing time. We summarize the flowdown of the Roman science objectives to the science and technical requirements of the HLSS. We construct a mock redshift survey over the full HLSS volume by applying a semi-analytic galaxy formation model to a cosmological N-body simulation, and use this mock survey to create pixel-level simulations of 4 deg$^2$ of HLSS grism spectroscopy. We find that the reference HLSS would measure $\sim$ 10 million H$α$ galaxy redshifts that densely map large scale structure at $z=1-2$ and 2 million [OIII] galaxy redshifts that sparsely map structures at $z=2-3$. We forecast the performance of this survey for measurements of the cosmic expansion history with baryon acoustic oscillations and the growth of large scale structure with redshift space distortions. We also study possible deviations from the reference design, and find that a deep HLSS at $f_{\rm line}>7\times10^{-17}$erg/s/cm$^2$ over 4000 deg$^2$ (requiring $\sim$1.5 yrs of observing time) provides the most compelling stand-alone constraints on dark energy from Roman alone. This provides a useful reference for future optimizations. The reference survey, simulated data sets, and forecasts presented here will inform community decisions on the final scope and design of the Roman HLSS.

preprint2020arXiv

A random-walk model for dark matter halo spins

We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the angular momenta of infalling subhalos. We further show that this model is in approximate agreement with the correlation of the spin magnitude over time as determined from N-body simulations, while it slightly underpredicts the correlation in the direction of the spin vector measured from the same simulations. This model is useful for predicting spins from merger histories derived from non-N-body sources, thereby circumventing the need for very high resolution simulations to permit accurate measurements of spins. It may be particularly relevant to modeling systems which accumulate angular momentum from halos over time (such as galactic discs) - we show that this model makes small but significant changes in the distribution of galactic disc sizes computed using the Galacticus semi-analytic galaxy formation model.

preprint2020arXiv

Clustering in the Simulated H$α$ Galaxy Redshift Survey from Nancy Grace Roman Space Telescope

We present a realistic 2000 deg$^{2}$ H$α$ galaxy mock catalog with $1<z<2$ for the Nancy Grace Roman Space Telescope galaxy redshift survey, the High Latitude Spectroscopic Survey (HLSS), created using Galacticus, a semi-analytical galaxy formation model, and high resolution cosmological N-body simulations. Galaxy clustering can probe dark energy and test gravity via baryon acoustic oscillation (BAO) and redshift space distortion (RSD) measurements. Using our realistic mock as the simulated Roman HLSS data, and a covariance matrix computed using a large set of approximate mocks created using EZmock, we investigate the expected precision and accuracy of the BAO and RSD measurements using the same analysis techniques used in analyzing real data. We find that the Roman H$α$ galaxy survey alone can measure the angular diameter distance with 2\% uncertainty, the Hubble parameter with 3-6\% uncertainty, and the linear growth parameter with 7\% uncertainty, in each of four redshift bins. Our realistic forecast illustrates the power of the Roman galaxy survey in probing the nature of dark energy and testing gravity.

preprint2020arXiv

Report from the Tri-Agency Cosmological Simulation Task Force

The Tri-Agency Cosmological Simulations (TACS) Task Force was formed when Program Managers from the Department of Energy (DOE), the National Aeronautics and Space Administration (NASA), and the National Science Foundation (NSF) expressed an interest in receiving input into the cosmological simulations landscape related to the upcoming DOE/NSF Vera Rubin Observatory (Rubin), NASA/ESA&#39;s Euclid, and NASA&#39;s Wide Field Infrared Survey Telescope (WFIRST). The Co-Chairs of TACS, Katrin Heitmann and Alina Kiessling, invited community scientists from the USA and Europe who are each subject matter experts and are also members of one or more of the surveys to contribute. The following report represents the input from TACS that was delivered to the Agencies in December 2018.

preprint2019arXiv

Constraints on the mass-concentration relation of cold dark matter halos with 11 strong gravitational lenses

The mass-concentration relation of dark matter halos reflects the assembly history of objects in hierarchical structure formation scenarios, and depends on fundamental quantities in cosmology such as the slope of the primordial matter power-spectrum. This relation is unconstrained by observations on sub-galactic scales. We derive the first measurement of the mass-concentration relation using the image positions and flux ratios from eleven quadruple-image strong gravitational lenses (quads) in the mass range $10^{6} - 10^{10} M_{\odot}$, assuming cold dark matter. Our analysis framework includes both subhalos and line of sight halos, marginalizes over nuisance parameters describing the lens macromodel, accounts for finite source effects on lensing observables, and simultaneously constrains the normalization and logarithmic slope of the mass-concentration relation, and the normalization of the subhalo mass function. At $z=0$, we constrain the concentration of $10^{8} M_{\odot}$ halos $c=12_{-5}^{+6}$ at $68 \%$ CI, and $c=12_{-9}^{+15}$ at $95 \%$ CI. For a $10^{7} M_{\odot}$ halo, we obtain $68 \%$ ($95 \%$) constraints $c=15_{-8}^{+9}$ ($c=15_{-11}^{+18}$), while for $10^{9} M_{\odot}$ halos $c=10_{-4}^{+7}$ ($c=10_{-7}^{+14}$). These results are consistent with the theoretical predictions from mass-concentration relations in the literature, and establish strong lensing by galaxies as a powerful probe of halo concentrations on sub-galactic scales across cosmological distance.

preprint2019arXiv

Generating Synthetic Cosmological Data with GalSampler

As part of the effort to meet the needs of the Large Synoptic Survey Telescope Dark Energy Science Collaboration (LSST DESC) for accurate, realistically complex mock galaxy catalogs, we have developed GalSampler, an open-source python package that assists in generating large volumes of synthetic cosmological data. The key idea behind GalSampler is to recast hydrodynamical simulations and semi-analytic models as physically-motivated galaxy libraries. GalSampler populates a new, larger-volume halo catalog with galaxies drawn from the baseline library; by using weighted sampling guided by empirical modeling techniques, GalSampler inherits statistical accuracy from the empirical model and physically-motivated complexity from the baseline library. We have recently used GalSampler to produce the cosmoDC2 extragalactic catalog made for the LSST DESC Data Challenge 2. Using cosmoDC2 as a guiding example, we outline how GalSampler can continue to support ongoing and near-future galaxy surveys such as the Dark Energy Survey (DES), the Dark Energy Spectroscopic Instrument (DESI), WFIRST, and Euclid.

preprint2019arXiv

Warm dark matter chills out: constraints on the halo mass function and the free-streaming length of dark matter with 8 quadruple-image strong gravitational lenses

The free-streaming length of dark matter depends on fundamental dark matter physics, and determines the abundance and concentration of dark matter halos on sub-galactic scales. Using the image positions and flux ratios from eight quadruply-imaged quasars, we constrain the free-streaming length of dark matter and the amplitude of the subhalo mass function (SHMF). We model both main deflector subhalos and halos along the line of sight, and account for warm dark matter (WDM) free-streaming effects on the mass function and mass-concentration relation. By calibrating the scaling of the SHMF with host halo mass and redshift using a suite of simulated halos, we infer a global normalization for the SHMF. We account for finite-size background sources, and marginalize over the mass profile of the main deflector. Parameterizing dark matter free-streaming through the half-mode mass $m_{\rm{hm}}$, we constrain the thermal relic particle mass $m_{\rm{DM}}$ corresponding to $m_{\rm{hm}}$. At $95 \%$ CI: $m_{\rm{hm}} < 10^{7.8} M_{\odot}$ ($m_{\rm{DM}} > 5.2 \ \rm{keV}$). We disfavor $m_{\rm{DM}} = 4.0 \rm{keV}$ and $ m_{\rm{DM}} = 3.0 \rm{keV}$ with likelihood ratios of 7:1 and 30:1, respectively, relative to the peak of the posterior distribution. Assuming cold dark matter, we constrain the projected mass in substructure between $10^6 - 10^{9} M_{\odot}$ near lensed images. At $68 \%$ CI, we infer $2.0 - 6.1 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to mean projected mass fraction $\bar{f}_{\rm{sub}} = 0.035_{-0.017}^{+0.021}$. At $95 \%$ CI, we obtain a lower bound on the projected mass of $0.6 \times 10^{7} M_{\odot} \rm{kpc^{-2}}$, corresponding to $\bar{f}_{\rm{sub}} > 0.005$. These results agree with the predictions of cold dark matter.