Source author record

Andrew Benson

Andrew Benson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

26 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's Euclid, and NASA'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.

preprint2016arXiv

The Astropy Problem

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

preprint2016arXiv

The connection between the host halo and the satellite galaxies of the Milky Way

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

preprint2015arXiv

An analytical model for galaxy metallicity: What do metallicity relations tell us about star formation and outflow?

We develop a simple analytical model that tracks galactic metallicities governed by star formation and feedback to gain insight from the observed galaxy stellar mass-metallicity relations over a large range of stellar masses and redshifts. The model reveals the following implications of star formation and feedback processes in galaxy formation. First, the observed metallicity relations provide a stringent upper limit for the averaged outflow mass-loading factors of local galaxies, which is ~20 for M_*~10^9Msun galaxies and monotonically decreases to ~1 for M_*~10^{11}Msun galaxies. Second, the inferred upper-limit for the outflow mass-loading factor sensitively depends on whether the outflow is metal-enriched with respect to the ISM metallicity. If half of the metals ejected from SNe leave the galaxy in metal-enriched winds, the outflow mass-loading factor for galaxies at any mass can barely be higher than ~10, which puts strong constraints on galaxy formation models. Third, the relatively lower stellar-phase to gas-phase metallicity ratio for lower-mass galaxies indicate that low-mass galaxies are still rapidly enriching their metallicities in recent times, while high-mass galaxies are more settled, which seems to show a downsizing effect in the metallicity evolution of galaxies. The analysis presented in the paper demonstrates the importance of accurate measurements of galaxy metallicities and the cold gas fraction of galaxies at different redshifts for constraining star formation and feedback processes, and demonstrates the power of these relations in constraining the physics of galaxy formation.

preprint2015arXiv

nIFTy Cosmology: Comparison of Galaxy Formation Models

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

preprint2014arXiv

A semi-analytic model comparison: testing cooling models against hydrodynamical simulations

We compare predictions of cooled masses and cooling rates from three stripped-down Semi-Analytic Models (SAMs) of galaxy formation with the results of N-body+SPH simulations with gas particle mass of 3.9x10^6 Msun, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ~2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z=0 to within ~20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved halos is subtracted out, SAMs tend to under-predict by ~30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses SAMs tend to over-predict cooling rates, though the numerical result may be affected by the use of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z~2, cold gas accumulates very quickly in the simulated halos. This accumulation is reproduced by SAMs with varying degrees of accuracy.

preprint2014arXiv

Metal distribution in the ICM - a comprehensive numerical study of twelve galaxy clusters

We present a simulation setup for studying the dynamical and chemical evolution of the intracluster medium (ICM) and analyze a sample of 12 galaxy clusters that are diverse both kinetically (pre-merger, merging, virialized) and in total mass (M vir = 1.17 x 10^14 - 1.06 x 10^15 M). We analyzed the metal mass fraction in the ICM as a function of redshift and discuss radial trends as well as projected 2D metallicity maps. The setup combines high mass resolution N-body simulations with the semi-analytical galaxy formation model Galacticus for consistent treatment of the subgrid physics (such as galactic winds and ram-pressure stripping) in the cosmological hydrodynamical simulations. The interface between Galacticus and the hydro simulation of the ICM with FLASH is discussed with respect to observations of star formation rate histories, radial star formation trends in galaxy clusters, and the metallicity at different redshifts. As a test for the robustness of the wind model, we compare three prescriptions from different approaches. For the wind model directly taken from Galacticus, we find mean ICM metallicities between 0.2 - 0.8Z within the inner 1Mpc at z = 0. The main contribution to the metal mass fraction comes from galactic winds. The outflows are efficiently mixed in the ICM, leading to a steady homogenization of metallicities until ram-pressure stripping becomes effective at low redshifts. We find a very peculiar and yet common drop in metal mass fractions within the inner ~200kpc of the cool cores, which is due to a combination of wind suppression by outer pressure within our model and a lack of mixing after the formation of these dense regions.

preprint2014arXiv

Oscillations and stability of polytropic filaments

We study the oscillations and stability of self-gravitating cylindrically symmetric fluid systems and collisionless systems. This is done by studying small perturbations to the equilibrium system and finding the normal modes, using methods similar to those used in astroseismology. We find that there is a single sequence of purely radial modes that become unstable if the adiabatic exponent is less than 1. Nonradial modes can be divided into p modes, which are stable and pressure-driven, and g modes, which are are gravity driven. The g modes become unstable if the adiabatic exponent is greater than the polytrope index. These modes are analogous to the modes of a spherical star, but their behavior is somewhat different because a cylindrical geometry has less symmetry than a spherical geometry. This implies that perturbations are classified by a radial quantum number, an azimuthal quantum number, and wavelength in the z direction, which can become arbitrarily large. We find that decreasing this wavelength increases the frequency of stable modes and increases the growth rate of unstable modes. We use use variational arguments to demonstrate that filaments of collisionless matter with ergodic distribution functions are stable to purely radial perturbations, and that filaments with ergodic power-law distribution functions are stable to all perturbations.

preprint2013arXiv

Constraining the WDM Particle Mass with Milky Way Satellites

Well-motivated particle physics theories predict the existence of particles (such as sterile neutrinos) which acquire non-negligible thermal velocities in the early universe. These particles could behave as warm dark matter (WDM) and generate a small-scale cutoff in the linear density power spectrum which scales approximately inversely with the particle mass. If this mass is of order a keV, the cutoff occurs on the scale of dwarf galaxies. Thus, in WDM models the abundance of small galaxies, such as the satellites that orbit in the halo of the Milky Way, depends on the mass of the warm particle. The abundance also scales with the mass of the host galactic halo. We use the \galform semi-analytic model of galaxy formation to calculate the properties of galaxies in universes in which the dark matter is warm. Using this method, we can compare the predicted satellite luminosity functions to the observed data for the Milky Way dwarf spheroidals, and determine a lower bound on the thermally produced WDM particle mass. This depends strongly on the value of the Milky Way halo mass and, to some extent, on the baryonic physics assumed; we examine both of these dependencies. For our fiducial model we find that for a particle mass of 3.3 keV (the 2$σ$ lower limit found by Viel et al. from a recent analysis of the Lyman-$α$ forest) the Milky Way halo mass is required to be $> 1.4 \times 10^{12}$ \msun. For this same fiducial model, we also find that all WDM particle masses are ruled out (at 95% confidence) if the halo of the Milky Way has a mass smaller than $1.1 \times 10^{12}$ \msun, while if the mass of the Galactic halo is greater than 1.8 $\times 10^{12}$ \msun, only WDM particle masses larger than 2 keV are allowed.

preprint2013arXiv

The Dwarfs Beyond: The Stellar-to-Halo Mass Relation for a New Sample of Intermediate Redshift Low Mass Galaxies

A number of recent challenges to the standard Lambda-CDM paradigm relate to discrepancies that arise in comparing the abundance and kinematics of local dwarf galaxies with the predictions of numerical simulations. Such arguments rely heavily on the assumption that the local dwarf and satellite galaxies form a representative distribution in terms of their stellar-to-halo mass ratios. To address this question, we present new, deep spectroscopy using DEIMOS on Keck for 82 low mass (10^7-10^9 solar masses) star-forming galaxies at intermediate redshift (z=0.2-1). For 50 percent of these we are able to determine resolved rotation curves using nebular emission lines and thereby construct the stellar mass Tully-Fisher relation to masses as low as 10^7 solar masses. Using scaling relations determined from weak lensing data, we convert this to a stellar-to-halo mass (SHM) relation for comparison with abundance matching predictions. We find a discrepancy between the propagated predictions from simulations compared to our observations, and suggest possible reasons for this as well as future tests that will be more effective.

preprint2013arXiv

The Escape Fraction of Ionizing Radiation from Galaxies

The escape of ionizing radiation from galaxies plays a critical role in the evolution of gas in galaxies, and the heating and ionization history of the intergalactic medium. We present semi-analytic calculations of the escape fraction of ionizing radiation for both hydrogen and helium from galaxies ranging from primordial systems to disk-type galaxies that are not heavily dust-obscured. We consider variations in the galaxy density profile, source type, location, and spectrum, and gas overdensity/distribution factors. For sufficiently hard first-light sources, the helium ionization fronts closely track or advance beyond that of hydrogen. Key new results in this work include calculations of the escape fractions for He I and He II ionizing radiation, and the impact of partial ionization from X-rays from early AGN or stellar clusters on the escape fractions from galaxy halos. When factoring in frequency-dependent effects, we find that X-rays play an important role in boosting the escape fractions for both hydrogen and helium, but especially for He II. We briefly discuss the implications of these results for recent observations of the He II reionization epoch at low redshifts, as well as the UV data and emission-line signatures from early galaxies anticipated from future satellite missions.

preprint2012arXiv

The SLUGGS Survey: Kinematics for over 2500 Globular Clusters in Twelve Early-type Galaxies

We present a spectro-photometric survey of 2522 extragalactic globular clusters (GCs) around twelve early-type galaxies, nine of which have not been published previously. Combining space-based and multi-colour wide field ground-based imaging, with spectra from the Keck DEIMOS instrument, we obtain an average of 160 GC radial velocities per galaxy, with a high velocity precision of 15 km/s per GC. After studying the photometric properties of the GC systems, such as their spatial and colour distributions, we focus on the kinematics of metal-poor (blue) and metal-rich (red) GC subpopulations to an average distance of ~8 effective radii from the galaxy centre. Our results show that for some systems the bimodality in GC colour is also present in GC kinematics. The kinematics of the red GC subpopulations are strongly coupled with the host galaxy stellar kinematics. The blue GC subpopulations are more dominated by random motions, especially in the outer regions, and decoupled from the red GCs. Peculiar GC kinematic profiles are seen in some galaxies: the blue GCs in NGC 821 rotate along the galaxy minor axis, whereas the GC system of the lenticular galaxy NGC 7457 appears to be strongly rotation supported in the outer region. We supplement our galaxy sample with data from the literature and carry out a number of tests to study the kinematic differences between the two GC subpopulations. We confirm that the GC kinematics are coupled with the host galaxy properties and find that the velocity kurtosis and the slope of their velocity dispersion profiles is different between the two GC subpopulations in more massive galaxies.

preprint2011arXiv

Bridging the gap between low and high mass dwarf galaxies

While the dark matter content within the most massive giant and smallest dwarf galaxies has been probed -- spanning a range of over one million in mass -- an important observational gap remains for galaxies of intermediate mass. This gap covers K band magnitudes of approximately -16 > M_K > -18 (for which dwarf galaxies have B--K ~ 2). On the high mass side of the gap are dwarf elliptical (dE) galaxies, that are dominated by stars in their inner regions. While the low mass side includes dwarf spheroidal (dSph) galaxies that are dark matter-dominated and ultra compact dwarf (UCD) objects that are star-dominated. Evolutionary pathways across the gap have been suggested but remain largely untested because the `gap' galaxies are faint, making dynamical measurements very challenging. With long exposures on the Keck telescope using the ESI instrument we have succeeded in bridging this gap by measuring the dynamical mass for five dwarf galaxies with M_K ~ -17.5 (M_B ~ --15.5). With the exception of our brightest dwarf galaxy, they possess relatively flat velocity dispersion profiles of around 20 km/s. By examining their 2D scaling relations and 3D fundamental manifold, we found that the sizes and velocity dispersions of these gap galaxies reveal continuous trends from dE to dSph galaxies. We conclude that low-luminosity dwarf elliptical galaxies are dominated by stars, not by dark matter, within their half light radii. This finding can be understood if internal feedback processes are operating most efficiently in gap galaxies, gravitationally heating the centrally-located dark matter to larger radii. Whereas external environmental processes, which can strip away stars, have a greater influence on dSph galaxies resulting in their higher dark matter fractions. Abridged.

preprint2011arXiv

New constraints on the evolution of the stellar-to-dark matter connection: a combined analysis of galaxy-galaxy lensing, clustering, and stellar mass functions from z=0.2 to z=1

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback.

preprint2009arXiv

Hot gas halos around disk galaxies: Confronting cosmological simulations with observations

Models of disk galaxy formation commonly predict the existence of an extended reservoir of accreted hot gas surrounding massive spirals at low redshift. As a test of these models, we use X-ray and H-alpha data of the two massive, quiescent edge-on spirals NGC 5746 and NGC 5170 to investigate the amount and origin of any hot gas in their halos. Contrary to our earlier claim, the Chandra analysis of NGC 5746, employing more recent calibration data, does not reveal any significant evidence for diffuse X-ray emission outside the optical disk, with a 3-sigma upper limit to the halo X-ray luminosity of 4e39 erg/s. An identical study of the less massive NGC 5170 also fails to detect any extraplanar X-ray emission. By extracting hot halo properties of disk galaxies formed in cosmological hydrodynamical simulations, we compare these results to expectations for cosmological accretion of hot gas by spirals. For Milky Way-sized galaxies, these high-resolution simulations predict hot halo X-ray luminosities which are lower by a factor of ~2 compared to our earlier results reported by Toft et al. (2002). We find the new simulation predictions to be consistent with our observational constraints for both NGC 5746 and NGC 5170, while also confirming that the hot gas detected so far around more actively star-forming spirals is in general probably associated with stellar activity in the disk. Observational results on quiescent disk galaxies at the high-mass end are nevertheless providing powerful constraints on theoretical predictions, and hence on the assumed input physics in numerical studies of disk galaxy formation and evolution.

preprint2009arXiv

The Evolutionary History of Lyman Break Galaxies Between Redshift 4 and 6: Observing Successive Generations of Massive Galaxies in Formation

We present new measurements of the evolution in the Lyman break galaxy (LBG) population between z~4 and z~6. By utilizing the extensive multiwavelength datasets available in the GOODS fields, we identify 2443 B, 506 V, and 137 i'-band dropout galaxies likely to be at z~4, 5, and 6. With the goal of understanding the duration of typical star formation episodes in galaxies at z>4, we examine the distribution of stellar masses and ages as a function of cosmic time. We find that at a fixed rest-UV luminosity, the average stellar masses and ages of galaxies do not increase significantly between z~6 and 4. In order to maintain this near equilibrium in the average properties of high redshift LBGs, we argue that there must be a steady flux of young, newly-luminous objects at each successive redshift. When considered along with the short duty cycles inferred from clustering measurements, these results may suggest that galaxies are undergoing star formation episodes lasting only several hundred million years. In contrast to the unchanging relationship between the average stellar mass and rest-UV luminosity, we find that the number density of massive galaxies increases considerably with time over 4<z<6. Given this rapid increase of UV luminous massive galaxies, we explore the possibility that a significant fraction of massive (M*>1e11 Msun) z~2-3 distant red galaxies (DRGs) were in part assembled in an LBG phase at earlier times. Integrating the growth in the stellar mass function of actively forming LBGs over 4<z<6 down to z~2, we find that z>3 LBGs could have contributed significantly to the quiescent DRG population, indicating that the intense star-forming systems probed by current sub-millimeter observations are not the only route toward the assembly of DRGs at z~2.