Showing 1–50 of 134 results for author: Benson, A J

JWST Lensed quasar dark matter survey II: Strongest gravitational lensing limit on the dark matter free streaming length to date

Authors: Ryan E. Keeley, Anna M. Nierenberg, Daniel Gilman, Charles Gannon, Simon Birrer, Tommaso Treu, Andrew J. Benson, Xiaolong Du, K. N. Abazajian, T. Anguita, V. N. Bennert, S. G. Djorgovski, K. K. Gupta, S. F. Hoenig, A. Kusenko, C. Lemon, M. Malkan, V. Motta, L. A. Moustakas, M. S. H. Oh, D. Sluse, D. Stern, R. H. Wechsler

Abstract: This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to pertu… ▽ More This is the second in a series of papers in which we use JWST MIRI multiband imaging to measure the warm dust emission in a sample of 31 multiply imaged quasars, to be used as a probe of the particle nature of dark matter. We present measurements of the relative magnifications of the strongly lensed warm dust emission in a sample of 9 systems. The warm dust region is compact and sensitive to perturbations by populations of halos down to masses $\sim 10^6$ M$_{\odot}$. Using these warm dust flux-ratio measurements in combination with 5 previous narrow-line flux-ratio measurements, we constrain the halo mass function. In our model, we allow for complex deflector macromodels with flexible third and fourth-order multipole deviations from ellipticity, and we introduce an improved model of the tidal evolution of subhalos. We constrain a WDM model and find an upper limit on the half-mode mass of $10^{7.6} M_\odot$ at posterior odds of 10:1. This corresponds to a lower limit on a thermally produced dark matter particle mass of 6.1 keV. This is the strongest gravitational lensing constraint to date, and comparable to those from independent probes such as the Ly$α$ forest and Milky Way satellite galaxies. △ Less

Submitted 2 May, 2024; originally announced May 2024.

LATIS: Constraints on the Galaxy-halo Connection at $z \sim 2.5$ from Galaxy-galaxy and Galaxy-Ly$α$ Clustering

Authors: Andrew B. Newman, Mahdi Qezlou, Nima Chartab, Gwen C. Rudie, Guillermo A. Blanc, Simeon Bird, Andrew J. Benson, Daniel D. Kelson, Brian C. Lemaux

Abstract: The connection between galaxies and dark matter halos is often quantified using the stellar mass-halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at $z \gtrsim 2$ can also constrain the SMHM relation via the g… ▽ More The connection between galaxies and dark matter halos is often quantified using the stellar mass-halo mass (SMHM) relation. Optical and near-infrared imaging surveys have led to a broadly consistent picture of the evolving SMHM relation based on measurements of galaxy abundances and angular correlation functions. Spectroscopic surveys at $z \gtrsim 2$ can also constrain the SMHM relation via the galaxy autocorrelation function and through the cross-correlation between galaxies and Ly$α$ absorption measured in transverse sightlines; however, such studies are very few and have produced some unexpected or inconclusive results. We use $\sim$3000 spectra of $z\sim2.5$ galaxies from the Lyman-alpha Tomography IMACS Survey (LATIS) to measure the galaxy-galaxy and galaxy-Ly$α$ correlation functions in four bins of stellar mass spanning $10^{9.2} \lesssim M_* / M_{\odot} \lesssim 10^{10.5}$. Parallel analyses of the MultiDark N-body and ASTRID hydrodynamic cosmological simulations allow us to model the correlation functions, estimate covariance matrices, and infer halo masses. We find that results of the two methods are mutually consistent and are broadly in accord with standard SMHM relations. This consistency demonstrates that we are able to accurately measure and model Ly$α$ transmission fluctuations $δ_F$ in LATIS. We also show that the galaxy-Ly$α$ cross-correlation, a free byproduct of optical spectroscopic galaxy surveys at these redshifts, can constrain halo masses with similar precision to galaxy-galaxy clustering. △ Less

Submitted 31 October, 2023; originally announced November 2023.

Comments: Accepted to ApJ

JWST lensed quasar dark matter survey I: Description and First Results

Authors: A. M. Nierenberg, R. E. Keeley, D. Sluse, D. Gilman, S. Birrer, T. Treu, K. N. Abazajian, T. Anguita, A. J. Benson, V. N. Bennert, S. G. Djorgovski, X. Du, C. D. Fassnacht, S. F. Hoenig, A. Kusenko, C. Lemon, M. Malkan, V. Motta, L. A. Moustakas, D. Stern, R. H. Wechsler

Abstract: The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light… ▽ More The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light is emitted from a physical region of the quasar with an angular scale of milli-arcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10$μ$m. The JWST Mid-Infrared Instrument (MIRI) is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3% precision. Our simulations suggest that this precision makes it feasible to detect the presence of 10$^7$ M$_\odot$ dark matter halos at cosmological distances. Such halos are expected to be completely dark in Cold Dark Matter models. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Devouring The Centaurus A Satellites: Modeling Dwarf Galaxies with Galacticus

Authors: Sachi Weerasooriya, Mia Sauda Bovill, Matthew A. Taylor, Andrew J. Benson, Cameron Leahy

Abstract: For the first time, systematic studies of dwarf galaxies are being conducted throughout the Local Volume, including Centaurus A (NGC 5128), which is the nearest elliptical galaxy. Given Centaurus As mass (roughly ten times that of the Milky Way), AGN activity, and recent major mergers, investigating these dwarfs and their star formation physics is imperative. However, simulating the faintest dwarf… ▽ More For the first time, systematic studies of dwarf galaxies are being conducted throughout the Local Volume, including Centaurus A (NGC 5128), which is the nearest elliptical galaxy. Given Centaurus As mass (roughly ten times that of the Milky Way), AGN activity, and recent major mergers, investigating these dwarfs and their star formation physics is imperative. However, simulating the faintest dwarfs in a massive galaxy like Centaurus A with sufficient resolution in a hydrodynamic simulation is computationally expensive and currently unfeasible. In this study, we seek to reproduce Centaurus A dwarfs using the same star formation physics as the Milky Way. We employ the semi-analytic model Galacticus to model dwarfs within a 600 kpc region. Utilizing astrophysical prescriptions and parameters matching the Milky Way satellites, we explore predictions for various properties and star formation histories (SFHs) to investigate environmental effects. We also reproduce cumulative luminosity and luminosity metallicity relations consistent with observations for the overall Centaurus A satellite population, while predicting half-light radii, velocity dispersion, and SFHs for the dwarf galaxies in Centaurus A. The agreement between our predicted SFHs for Centaurus A dwarfs and those of the Milky Way implies the presence of universal processes governing star formation in these galaxies. Overall, our findings shed light on the star formation physics of dwarf galaxies in the Centaurus A system, revealing insights into their properties and dependence on the host environment. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 12 pages, 8 figures

Do assumptions about the central density of subhaloes affect dark matter annihilation and lensing calculations?

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: A growing body of evidence suggests that the central density of cuspy dark matter subhaloes is conserved in minor mergers. However, empirical models of subhalo evolution, calibrated from simulations, often assume a drop in the central density. Since empirical models of subhaloes are used in galaxy-galaxy lensing studies and dark matter annihilation calculations, we explore the consequences of assu… ▽ More A growing body of evidence suggests that the central density of cuspy dark matter subhaloes is conserved in minor mergers. However, empirical models of subhalo evolution, calibrated from simulations, often assume a drop in the central density. Since empirical models of subhaloes are used in galaxy-galaxy lensing studies and dark matter annihilation calculations, we explore the consequences of assuming different subhalo models. We find that dark matter annihilation calculations are very sensitive to the assumed subhalo mass profile, and different models can give more than a magnitude difference in the J-factor and boost factor in individual haloes. On the other hand, the shear and convergence profiles used in galaxy-galaxy lensing are sensitive to the initial profile assumed (e.g., NFW versus Einato) but are otherwise well-approximated by a simple model in which the original profile is sharply truncated. We conclude that since the innermost parts of haloes are difficult to resolve in simulations, it is important to have a theoretical understanding of how subhaloes evolve to make accurate predictions of the dark matter annihilation signal. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Comments: 15 pages, 13 figures. Submitted to MNRAS

Characterizing the ELG luminosity functions in the nearby Universe

Authors: G. Favole, V. Gonzalez-Perez, Y. Ascasibar, P. Corcho-Caballero, A. D. Montero-Dorta, A. J. Benson, J. Comparat, S. A. Cora, D. Croton, H. Guo, D. Izquierdo-Villalba, A. Knebe, Á. Orsi, D. Stoppacher, C. A. Vega-Martínez

Abstract: Nebular emission lines are powerful diagnostics for the physical processes at play in galaxy formation and evolution. Moreover, emission-line galaxies (ELGs) are one of the main targets of current and forthcoming spectroscopic cosmological surveys. We investigate the contributions to the line luminosity functions (LFs) of different galaxy populations in the local Universe, providing a benchmark fo… ▽ More Nebular emission lines are powerful diagnostics for the physical processes at play in galaxy formation and evolution. Moreover, emission-line galaxies (ELGs) are one of the main targets of current and forthcoming spectroscopic cosmological surveys. We investigate the contributions to the line luminosity functions (LFs) of different galaxy populations in the local Universe, providing a benchmark for future surveys of earlier cosmic epochs. The large statistics of the observations from the SDSS DR7 Main galaxy sample and the MPA-JHU spectral catalogue enabled us to precisely measure the H$α$, H$β$, [OII], [OIII], and, for the first time, the [NII], and [SII] emission-line LFs over ~2.4 Gyrs in the low-z Universe, 0.02<z<0.22. We present a generalized 1/Vmax LF estimator capable of simultaneously correcting for spectroscopic, r-band magnitude, and emission-line incompleteness. We studied the contribution to the LF of different types of ELGs classified using two methods: (i) the value of the specific star formation rate (sSFR), and (ii) the line ratios on the Baldwin-Phillips-Terlevich (BPT) and the WHAN (i.e., H$α$ equivalent width versus the [NII]/H$α$ line ratio) diagrams. The ELGs in our sample are mostly star forming, with 84 per cent having sSFR>10$^{-11}$/yr. When classifying ELGs using the BPT+WHAN diagrams, we find that 63.3 per cent are star forming, only 0.03 are passively evolving, and 1.3 have nuclear activity (Seyfert). The rest are low-ionization narrow emission-line regions (LINERs) and composite ELGs. We found that a Saunders function is the most appropriate to describe all of the emission-line LFs, both observed and dust-extinction-corrected. They are dominated by star-forming regions, except for the bright end of the [OIII] and [NII] LFs (i.e. L[NII]>10$^{42}$ erg/s, L[OIII]>10$^{43}$ erg/s), where the contribution of Seyfert galaxies is not negligible. △ Less

Submitted 14 December, 2023; v1 submitted 20 March, 2023; originally announced March 2023.

Evaluating the sustainability of a de facto harvest strategy for British Columbia's Spot Prawn (Pandalus platyceros) fishery in the presence of environmental drivers of recruitment and hyperstable catch rates

Authors: Steven P. Rossi, Sean P. Cox, Samuel D. N. Johnson, Ashleen J. Benson

Abstract: The Spot Prawn trap fishery off the west coast of British Columbia (BC) is managed using a fixed escapement strategy that aims to prevent recruitment overfishing while maximizing expected long-term yield by closing the fishery when the catch rate of spawners, projected to the following spring, drops below 1.7 spawners per trap (the de jure rule). We develop a management strategy evaluation framewo… ▽ More The Spot Prawn trap fishery off the west coast of British Columbia (BC) is managed using a fixed escapement strategy that aims to prevent recruitment overfishing while maximizing expected long-term yield by closing the fishery when the catch rate of spawners, projected to the following spring, drops below 1.7 spawners per trap (the de jure rule). We develop a management strategy evaluation framework for BC's Spot Prawn fishery that examines the expected performance of the management procedure implemented in practice (the de facto rule), which was significantly more conservative than the de jure rule, usually closing the fishery when spawner catch rates were at least twice as high as specified by the de jure rule. Simulations indicate that the de facto spawner index rule using average empirical March 31st targets from 2000 to 2019 maintains most stocks near or above 0.8 BMSY with or without accounting for environmental effects and/or increasing future SST on recruitment. Abundance indices were found to be strongly hyperstable, with fishing efficiency 1.5 to 3.0 times higher under low biomass than high biomass. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 86 pages, 12 figures, 33 appendix figures

A universal model for the evolution of tidally stripped systems

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: Accurate models of the structural evolution of dark matter subhaloes, as they orbit within larger systems, are fundamental to understanding the detailed distribution of dark matter at the present day. Numerical simulations of subhalo evolution support the idea that the mass loss associated with tidal stripping is most naturally understood in energy space, with the particles that are the least boun… ▽ More Accurate models of the structural evolution of dark matter subhaloes, as they orbit within larger systems, are fundamental to understanding the detailed distribution of dark matter at the present day. Numerical simulations of subhalo evolution support the idea that the mass loss associated with tidal stripping is most naturally understood in energy space, with the particles that are the least bound being removed first. Starting from this premise, we recently proposed a zero-parameter "energy-truncation model" for subhalo evolution. We tested this model with simulations of tidal stripping of satellites with initial NFW profiles, and showed that the energy-truncation model accurately predicts both the mass loss and density profiles. In this work, we apply the model to a variety of Hernquist, Einasto and King profiles. We show that it matches the simulation results quite closely in all cases, indicating that it may serve as a universal model to describe tidally stripped collisionless systems. A key prediction of the energy-truncation model is that the central density of dark matter subhaloes is conserved as they lose mass; this has important implications for dark matter annihilation calculations, and for other observational tests of dark matter. △ Less

Submitted 29 July, 2022; originally announced July 2022.

Comments: 19 pages, 14 figures. Accepted for publication in MNRAS

Tidal Tracks and Artificial Disruption of Cold Dark Matter Halos

Authors: Andrew J. Benson, Xiaolong Du

Abstract: We describe a simple extension to existing models for the tidal heating of dark matter subhalos which takes into account second order terms in the impulse approximation for tidal heating. We show that this revised model can accurately match the tidal tracks along which subhalos evolve as measured in high-resolution N-body simulations. We further demonstrate that, when a constant density core is in… ▽ More We describe a simple extension to existing models for the tidal heating of dark matter subhalos which takes into account second order terms in the impulse approximation for tidal heating. We show that this revised model can accurately match the tidal tracks along which subhalos evolve as measured in high-resolution N-body simulations. We further demonstrate that, when a constant density core is introduced into a subhalo, this model is able to quantitatively reproduce the evolution and artificial disruption of N-body subhalos arising from finite resolution effects. Combining these results we confirm prior work indicating that artificial disruption in N-body simulations can result in a factor two underestimate of the subhalo mass function in the inner regions of host halos, and a 10--20% reduction over the entire virial volume. △ Less

Submitted 22 September, 2022; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 9 pages, accepted for publication in MNRAS

Constraints on sterile neutrino models from strong gravitational lensing, Milky Way satellites, and Lyman-$α$ forest

Authors: Ioana A. Zelko, Tommaso Treu, Kevork N. Abazajian, Daniel Gilman, Andrew J. Benson, Simon Birrer, Anna M. Nierenberg, Alexander Kusenko

Abstract: The nature of dark matter is one of the most important unsolved questions in science. Some dark matter candidates do not have sufficient nongravitational interactions to be probed in laboratory or accelerator experiments. It is thus important to develop astrophysical probes which can constrain or lead to a discovery of such candidates. We illustrate this using state-of-the-art measurements of stro… ▽ More The nature of dark matter is one of the most important unsolved questions in science. Some dark matter candidates do not have sufficient nongravitational interactions to be probed in laboratory or accelerator experiments. It is thus important to develop astrophysical probes which can constrain or lead to a discovery of such candidates. We illustrate this using state-of-the-art measurements of strong gravitationally-lensed quasars to constrain four of the most popular sterile neutrino models, and also report the constraints for other independent methods that are comparable in procedure. First, we derive effective relations to describe the correspondence between the mass of a thermal relic warm dark matter particle and the mass of sterile neutrinos produced via Higgs decay and GUT-scale scenarios, in terms of large-scale structure and galaxy formation astrophysical effects. Second, we show that sterile neutrinos produced through the Higgs decay mechanism are allowed only for mass $>26$ keV, and GUT-scale scenario $>5.3$ keV. Third, we show that the single sterile neutrino model produced through active neutrino oscillations is allowed for mass $>92$ keV, and the 3 sterile neutrino minimal standard model ($ν$MSM) for mass $>16$ keV. These are the most stringent experimental limits on these models. △ Less

Submitted 18 July, 2023; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: Published in Physical Review Letters (PRL); 10 pages, 3 figures

Journal ref: Phys. Rev. Lett. 129, 191301. Published 4 November 2022

A Random Walk Model for Halo Triaxiality

Authors: Paul Menker, Andrew J. Benson

Abstract: We describe a semi-analytic model to predict the triaxial shapes of dark matter halos utilizing the sequences of random merging events captured in merger trees to follow the evolution of each halo's energy tensor. When coupled with a simple model for relaxation toward a spherical shape, we find that this model predicts distributions of halo axis length ratios which approximately agree with those m… ▽ More We describe a semi-analytic model to predict the triaxial shapes of dark matter halos utilizing the sequences of random merging events captured in merger trees to follow the evolution of each halo's energy tensor. When coupled with a simple model for relaxation toward a spherical shape, we find that this model predicts distributions of halo axis length ratios which approximately agree with those measured from cosmological N-body simulations once constrained to match the median axis ratio at a single halo mass. We demonstrate the predictive and explanatory power of this model by considering conditioned distributions of axis length ratios, and the mass-dependence of halo shapes, finding these to be in good agreement with N-body results. This model provides both insight into the physics driving the evolution of halo triaxial shapes, and rapid quantitative predictions for the statistics of triaxiality connected directly to the formation history of the halo. △ Less

Submitted 29 August, 2022; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: 9 pages, 5 figures

Probing Hot Gas Components of Circumgalactic Medium in Cosmological Simulations with the Thermal Sunyaev-Zel'dovich Effect

Authors: Junhan Kim, Sunil Golwala, James G. Bartlett, Stefania Amodeo, Nicholas Battaglia, Andrew J. Benson, J. Colin Hill, Philip F. Hopkins, Cameron B. Hummels, Emily Moser, Matthew E. Orr

Abstract: The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful tool with the potential for constraining directly the properties of the hot gas that dominates dark matter halos because it measures pressure and thus thermal energy density. Studying this hot component of the circumgalactic medium (CGM) is important because it is strongly impacted by star-formation and active galactic nuclei (AGN) activity… ▽ More The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful tool with the potential for constraining directly the properties of the hot gas that dominates dark matter halos because it measures pressure and thus thermal energy density. Studying this hot component of the circumgalactic medium (CGM) is important because it is strongly impacted by star-formation and active galactic nuclei (AGN) activity in galaxies, participating in the feedback loop that regulates star and black hole mass growth in galaxies. We study the tSZ effect across a wide halo mass range using three cosmological hydrodynamical simulations: Illustris-TNG, EAGLE, and FIRE-2. Specifically, we present the scaling relation between tSZ signal and halo mass and radial profiles of gas density, temperature, and pressure for all three simulations. The analysis includes comparisons to Planck tSZ observations and to the thermal pressure profile inferred from the Atacama Cosmology Telescope (ACT) measurements. We compare these tSZ data to the simulations to interpret the measurements in terms of feedback and accretion processes in the CGM. We also identify as-yet unobserved potential signatures of these processes that may be visible in future measurements, which will have the capability of measuring tSZ signals to even lower masses. We also perform internal comparisons between runs with different physical assumptions. We conclude: (1) there is strong evidence for the impact of feedback at $R_{500}$ but that this impact decreases by $5R_{500}$, and (2) the thermodynamic profiles of the CGM are highly dependent on the implemented model, such as cosmic-ray or AGN feedback prescriptions. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: 21 pages, 10 figures, submitted to ApJ, comments welcome

Validating Synthetic Galaxy Catalogs for Dark Energy Science in the LSST Era

Authors: Eve Kovacs, Yao-Yuan Mao, Michel Aguena, Anita Bahmanyar, Adam Broussard, James Butler, Duncan Campbell, Chihway Chang, Shenming Fu, Katrin Heitmann, Danila Korytov, François Lanusse, Patricia Larsen, Rachel Mandelbaum, Christopher B. Morrison, Constantin Payerne, Marina Ricci, Eli Rykoff, F. Javier Sánchez, Ignacio Sevilla-Noarbe, Melanie Simet, Chun-Hao To, Vinu Vikraman, Rongpu Zhou, Camille Avestruz , et al. (14 additional authors not shown)

Abstract: Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation proc… ▽ More Large simulation efforts are required to provide synthetic galaxy catalogs for ongoing and upcoming cosmology surveys. These extragalactic catalogs are being used for many diverse purposes covering a wide range of scientific topics. In order to be useful, they must offer realistically complex information about the galaxies they contain. Hence, it is critical to implement a rigorous validation procedure that ensures that the simulated galaxy properties faithfully capture observations and delivers an assessment of the level of realism attained by the catalog. We present here a suite of validation tests that have been developed by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). We discuss how the inclusion of each test is driven by the scientific targets for static ground-based dark energy science and by the availability of suitable validation data. The validation criteria that are used to assess the performance of a catalog are flexible and depend on the science goals. We illustrate the utility of this suite by showing examples for the validation of cosmoDC2, the extragalactic catalog recently released for the LSST DESC second Data Challenge. △ Less

Submitted 13 January, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 46 pages, 33 figures

The Galaxy Evolution Probe

Authors: Jason Glenn, Charles M. Bradford, Erik Rosolowsky, Rashied Amini, Katherine Alatalo, Lee Armus, Andrew J. Benson, Tzu-Ching Chang, Jeremy Darling, Peter K. Day, Jeanette Domber, Duncan Farrah, Brandon Hensley, Sarah Lipscy, Bradley Moore, Seb Oliver, Joanna Perido, David Redding, Michael Rodgers, Raphael Shirley, Howard A. Smith, John B. Steeves, Carole Tucker, Jonas Zmuidzinas

Abstract: The Galaxy Evolution Probe (GEP) is a concept for a mid- and far-infrared space observatory to measure key properties of large samples of galaxies with large and unbiased surveys. GEP will attempt to achieve zodiacal light and Galactic dust emission photon background-limited observations by utilizing a 6 Kelvin, 2.0 meter primary mirror and sensitive arrays of kinetic inductance detectors. It will… ▽ More The Galaxy Evolution Probe (GEP) is a concept for a mid- and far-infrared space observatory to measure key properties of large samples of galaxies with large and unbiased surveys. GEP will attempt to achieve zodiacal light and Galactic dust emission photon background-limited observations by utilizing a 6 Kelvin, 2.0 meter primary mirror and sensitive arrays of kinetic inductance detectors. It will have two instrument modules: a 10 - 400 micron hyperspectral imager with spectral resolution R = 8 (GEP-I) and a 24 - 193 micron, R = 200 grating spectrometer (GEP-S). GEP-I surveys will identify star-forming galaxies via their thermal dust emission and simultaneously measure redshifts using polycyclic aromatic hydrocarbon emission lines. Galaxy luminosities derived from star formation and nuclear supermassive black hole accretion will be measured for each source, enabling the cosmic star formation history to be measured to much greater precision than previously possible. Using optically thin far-infrared fine-structure lines, surveys with GEP-S will measure the growth of metallicity in the hearts of galaxies over cosmic time and extraplanar gas will be mapped in spiral galaxies in the local universe to investigate feedback processes. The science case and mission architecture designed to meet the science requirements are described, and the kinetic inductance detector and readout electronics state of the art and needed developments are described. This paper supersedes the GEP concept study report cited in it by providing new content, including: a summary of recent mid-infrared KID development, a discussion of microlens array fabrication for mid-infrared KIDs, and additional context for galaxy surveys. The reader interested in more technical details may want to consult the concept study report. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 84 pages, 24 figures, published in JATIS

Journal ref: J. Astron. Telesc. Instrum. Syst. 7(3), 034004 (2021)

The Origin of the Orbital Parameter Distribution of Merging Halos

Authors: A. J. Benson

Abstract: We describe a simple model which explains the qualitative and (approximate) quantitative features of the distribution of orbital velocities of merging pairs of dark matter halos. Our model considers a primary dark matter halo as a perturber in a background of secondary halos with velocities described by linear theory. By evaluating the ensemble of secondary halos on orbits within the perturbing ha… ▽ More We describe a simple model which explains the qualitative and (approximate) quantitative features of the distribution of orbital velocities of merging pairs of dark matter halos. Our model considers a primary dark matter halo as a perturber in a background of secondary halos with velocities described by linear theory. By evaluating the ensemble of secondary halos on orbits within the perturbing halo's "loss cone" we derive the distribution of orbital parameters of these captured halos. This model is able provide qualitative explanations for the features of this distribution as measured from N-body simulations, and is in approximate quantitative agreement with those measurements. As the velocity dispersion of the background halos is larger on smaller scales our model predicts an overall increase in the characteristic velocities of merging halos, relative to the virial velocities of those halos, in lower mass systems. Our model also provides a simple explanation for the measured independence of the orbital velocity distribution function on redshift when considered at fixed peak height. By connecting the orbital parameter distribution to the underlying power spectrum our model also allows for estimates to be made of the effect of modifying that power spectrum, for example by including a truncation at large wavenumber. For plausible warm dark matter models we find that this truncation has only a small effect on the predicted distributions. △ Less

Submitted 13 May, 2021; v1 submitted 5 April, 2021; originally announced April 2021.

Comments: 11 pages; accepted in MNRAS

Constraints on the properties of warm dark matter using the satellite galaxies of the Milky Way

Authors: Oliver Newton, Matteo Leo, Marius Cautun, Adrian Jenkins, Carlos S. Frenk, Mark R. Lovell, John C. Helly, Andrew J. Benson

Abstract: The satellite galaxies of the Milky Way (MW) are effective probes of the underlying dark matter (DM) substructure, which is sensitive to the nature of the DM particle. In particular, a class of DM models have a power spectrum cut-off on the mass scale of dwarf galaxies and thus predict only small numbers of substructures below the cut-off mass. This makes the MW satellite system appealing to const… ▽ More The satellite galaxies of the Milky Way (MW) are effective probes of the underlying dark matter (DM) substructure, which is sensitive to the nature of the DM particle. In particular, a class of DM models have a power spectrum cut-off on the mass scale of dwarf galaxies and thus predict only small numbers of substructures below the cut-off mass. This makes the MW satellite system appealing to constrain the DM properties: feasible models must produce enough substructure to host the number of observed Galactic satellites. Here, we compare theoretical predictions of the abundance of DM substructure in thermal relic warm DM (WDM) models with estimates of the total satellite population of the MW. This produces conservative robust lower limits on the allowed mass, $m_\mathrm{th}$, of the thermal relic WDM particle. As the abundance of satellite galaxies depends on the MW halo mass, we marginalize over the corresponding uncertainties and rule out $m_\mathrm{th} \leq 2.02\, \mathrm{keV}$ at 95 per cent confidence independently of assumptions about galaxy formation processes. Modelling some of these - in particular, the effect of reionization, which suppresses the formation of dwarf galaxies - strengthens our constraints on the DM properties and excludes models with $m_\mathrm{th} \leq 3.99\, \mathrm{keV}$ in our fiducial model. We also find that thermal relic models cannot produce enough satellites if the MW halo mass is $M_{200}\leq 0.6\times 10^{12}\, \mathrm{M_\odot}$, which imposes a lower limit on the MW halo mass in CDM. We address several observational and theoretical uncertainties and discuss how improvements in these will strengthen the DM mass constraints. △ Less

Submitted 1 September, 2021; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: 29 pages (JCAP format), 3 appendices, 10 figures, 1 table. Version accepted by JCAP and published by IOP Publishing under Gold Open Access - Transformative Agreement with a CC-BY 4.0 licence

A Random Walk Model for Dark Matter Halo Concentrations

Authors: Turner Johnson, Andrew J. Benson, Daniel Grin

Abstract: For idealized (spherical, smooth) dark matter halos described by single-parameter density profiles (such as the NFW profile) there exists a one-to-one mapping between the energy of the halo and the scale radius of its density profile. The energy therefore uniquely determines the concentration parameter of such halos. We exploit this fact to predict the concentrations of dark matter halos via a ran… ▽ More For idealized (spherical, smooth) dark matter halos described by single-parameter density profiles (such as the NFW profile) there exists a one-to-one mapping between the energy of the halo and the scale radius of its density profile. The energy therefore uniquely determines the concentration parameter of such halos. We exploit this fact to predict the concentrations of dark matter halos via a random walk in halo energy space. Given a full merger tree for a halo, the total internal energy of each halo in that tree is determined by summing the internal and orbital energies of progenitor halos. We show that, when calibrated, this model can accurately reproduce the mean of the concentration--mass relation measured in N-body simulations, and reproduces more of the scatter in that relation than previous models. We further test this model by examining both the autocorrelation of scale radii across time, and the correlations between halo concentration and spin, and comparing to results measured from cosmological N-body simulations. In both cases we find that our model closely matches the N-body results. Our model is implemented within the open source Galacticus toolkit. △ Less

Submitted 13 January, 2021; v1 submitted 26 June, 2020; originally announced June 2020.

Comments: 22 pages, 12 figures, ApJ in press

A new calibration method of sub-halo orbital evolution for semi-analytic models

Authors: Shengqi Yang, Xiaolong Du, Andrew J. Benson, Anthony R. Pullen, Annika H. G. Peter

Abstract: Understanding the non-linear dynamics of satellite halos (a.k.a. "sub-halos") is important for predicting the abundance and distribution of dark matter substructures and satellite galaxies, and for distinguishing among microphysical dark matter models using observations. Typically, modeling these dynamics requires large N-body simulations with high resolution. Semi-analytic models can provide a mo… ▽ More Understanding the non-linear dynamics of satellite halos (a.k.a. "sub-halos") is important for predicting the abundance and distribution of dark matter substructures and satellite galaxies, and for distinguishing among microphysical dark matter models using observations. Typically, modeling these dynamics requires large N-body simulations with high resolution. Semi-analytic models can provide a more efficient way to describe the key physical processes such as dynamical friction, tidal mass loss, and tidal heating, with only a few free parameters. In this work, we present a fast Monte Carlo Markov Chain fitting approach to explore the parameter space of such a sub-halo non-linear evolution model. We use the dynamical models described in an earlier work and calibrate the models to two sets of high-resolution cold dark matter N-body simulations, ELVIS and Caterpillar. Compared to previous calibrations that used manual parameter tuning, our approach provides a more robust way to determine the best-fit parameters and their posterior probabilities. We find that jointly fitting for the sub-halo mass and maximum velocity functions can break the degeneracy between tidal stripping and tidal heating parameters, as well as providing better constraints on the strength of dynamical friction. We show that our semi-analytic simulation can accurately reproduce N-body simulations statistics, and that the calibration results for the two sets of N-body simulations agree at 95% confidence level. Dynamical models calibrated in this work will be important for future dark matter substructure studies. △ Less

Submitted 23 September, 2020; v1 submitted 23 March, 2020; originally announced March 2020.

Comments: 13 pages, 8 figures

Mass loss in tidally stripped systems; the energy-based truncation method

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does… ▽ More The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does not provide a physical understanding of the processes involved. Recently, we demonstrated that truncating NFW distribution functions sharply in energy results in density profiles that resemble those of tidally stripped systems, offering a path to constructing physically motivated models of tidal mass loss. In this work, we review calculations of mass loss based on energy truncation alone, and then consider what secondary effects may modulate mass loss beyond this. We find that a combination of dependence on additional orbital parameters and variations in individual particle energies over an orbit results in a less abrupt truncation in energy space as a subhalo loses mass. Combining the energy truncation approach with a simple prediction for the mass-loss rate, we construct a full model of mass loss that can accurately predict the evolution of a subhalo in terms of a single parameter $η_{\rm eff}$. This parameter can be fully determined from the initial orbital and halo properties, and does not require calibration with numerical simulations. △ Less

Submitted 20 March, 2020; originally announced March 2020.

Comments: 20 pages, 14 figures, accepted to MNRAS

LATIS: The Ly$α$ Tomography IMACS Survey

Authors: Andrew B. Newman, Gwen C. Rudie, Guillermo A. Blanc, Daniel D. Kelson, Sunny Rhoades, Tyson Hare, Victoria Pérez, Andrew J. Benson, Alan Dressler, Valentino Gonzalez, Juna A. Kollmeier, Nicholas P. Konidaris, John S. Mulchaey, Michael Rauch, Olivier Le Fèvre, Brian C. Lemaux, Olga Cucciati, Simon J. Lilly

Abstract: We introduce LATIS, the Ly$α$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$α$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightl… ▽ More We introduce LATIS, the Ly$α$ Tomography IMACS Survey, a spectroscopic survey at Magellan designed to map the z=2.2-2.8 intergalactic medium (IGM) in three dimensions by observing the Ly$α$ forest in the spectra of galaxies and QSOs. Within an area of 1.7 deg${}^2$, we will observe approximately half of $\gtrsim L^*$ galaxies at z=2.2-3.2 for typically 12 hours, providing a dense network of sightlines piercing the IGM with an average transverse separation of 2.5 $h^{-1}$ comoving Mpc (1 physical Mpc). At these scales, the opacity of the IGM is expected to be closely related to the dark matter density, and LATIS will therefore map the density field in the $z \sim 2.5$ universe at $\sim$Mpc resolution over the largest volume to date. Ultimately LATIS will produce approximately 3800 spectra of z=2.2-3.2 galaxies that probe the IGM within a volume of $4 \times 10^6 h^{-3}$ Mpc${}^3$, large enough to contain a representative sample of structures from protoclusters to large voids. Observations are already complete over one-third of the survey area. In this paper, we describe the survey design and execution. We present the largest IGM tomographic maps at comparable resolution yet made. We show that the recovered matter overdensities are broadly consistent with cosmological expectations based on realistic mock surveys, that they correspond to galaxy overdensities, and that we can recover structures identified using other tracers. LATIS is conducted in Canada-France-Hawaii Telescope Legacy Survey fields, including COSMOS. Coupling the LATIS tomographic maps with the rich data sets collected in these fields will enable novel studies of environment-dependent galaxy evolution and the galaxy-IGM connection at cosmic noon. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: accepted to ApJ; Fig. 26 contains animated rendering of the IGM maps

The Normalization and Slope of the Dark Matter (Sub-)Halo Mass Function on Sub-Galactic Scales

Authors: Andrew J. Benson

Abstract: Simulations of cold dark matter make robust predictions about the slope and normalization of the dark matter halo and subhalo mass functions on small scales. Recent observational advances utilizing strong gravitational lensing have demonstrated the ability of this technique to place constraints on these quantities on subgalactic scales corresponding to dark matter halo masses of $10^6$--… ▽ More Simulations of cold dark matter make robust predictions about the slope and normalization of the dark matter halo and subhalo mass functions on small scales. Recent observational advances utilizing strong gravitational lensing have demonstrated the ability of this technique to place constraints on these quantities on subgalactic scales corresponding to dark matter halo masses of $10^6$--$10^9\mathrm{M}_\odot$. On these scales the physics of baryons, which make up around 17% of the matter content of the Universe but which are not included in pure dark matter N-body simulations, are expected to affect the growth of structure and the collapse of dark matter halos. In this work we develop a semi-analytic model to predict the amplitude and slope of the dark matter halo and subhalo mass functions on subgalactic scales in the presence of baryons. We find that the halo mass function is suppressed by up to 25%, and the slope is modified, ranging from $-1.916$ to $-1.868$ in this mass range. These results are consistent with current measurements, but differ sufficiently from the expectations for a dark matter only universe that it may be testable in the near future. △ Less

Submitted 6 February, 2020; v1 submitted 11 November, 2019; originally announced November 2019.

Comments: 9 pages, MNRAS accepted

Gravity and the Nonlinear Growth of Structure in the Carnegie-Spitzer-IMACS Redshift Survey

Authors: Daniel D. Kelson, Louis E. Abramson, Andrew J. Benson, Shannon G. Patel, Stephen A. Shectman, Alan Dressler, Patrick J. McCarthy, John S. Mulchaey, Rik J. Williams

Abstract: A key obstacle to developing a satisfying theory of galaxy evolution is the difficulty in extending analytic descriptions of early structure formation into full nonlinearity, the regime in which galaxy growth occurs. Extant techniques, though powerful, are based on approximate numerical methods whose Monte Carlo-like nature hinders intuition building. Here, we develop a new solution to this proble… ▽ More A key obstacle to developing a satisfying theory of galaxy evolution is the difficulty in extending analytic descriptions of early structure formation into full nonlinearity, the regime in which galaxy growth occurs. Extant techniques, though powerful, are based on approximate numerical methods whose Monte Carlo-like nature hinders intuition building. Here, we develop a new solution to this problem and its empirical validation. We first derive closed-form analytic expectations for the evolution of fixed percentiles in the real-space cosmic density distribution, {\it averaged over representative volumes observers can track cross-sectionally\}. Using the Lagrangian forms of the fluid equations, we show that percentiles in $δ$---the density relative to the median---should grow as $δ(t)\proptoδ_{0}^α\,t^β$, where $α\equiv2$ and $β\equiv2$ for Newtonian gravity at epochs after the overdensities transitioned to nonlinear growth. We then use 9.5 sq. deg. of Carnegie-Spitzer-IMACS Redshift Survey data to map {\it galaxy\} environmental densities over $0.2<z<1.5$ ($\sim$7 Gyr) and infer $α=1.98\pm0.04$ and $β=2.01\pm0.11$---consistent with our analytic prediction. These findings---enabled by swapping the Eulerian domain of most work on density growth for a Lagrangian approach to real-space volumetric averages---provide some of the strongest evidence that a lognormal distribution of early density fluctuations indeed decoupled from cosmic expansion to grow through gravitational accretion. They also comprise the first exact, analytic description of the nonlinear growth of structure extensible to (arbitrarily) low redshift. We hope these results open the door to new modeling of, and insight-building into, the galaxy growth and its diversity in cosmological contexts. △ Less

Submitted 21 January, 2020; v1 submitted 23 August, 2019; originally announced August 2019.

Comments: 13 pages, 12 figures, submitted to Monthly Notices of the Royal Astronomical Society, 23 August 2019, and accepted for publication, 9 January 2020. Some more typo(s) fixed

[OII] emitters in MultiDark-Galaxies and DEEP2

Authors: G. Favole, V. Gonzalez-Perez, A. Orsi, D. Stoppacher, J. Comparat, S. A. Cora, C. A. Vega-Martinez, A. R. H. Stevens, C. Maraston, D. Croton, A. Knebe, A. J. Benson, A. D. Montero-Dorta, N. Padilla, F. Prada, D. Thomas

Abstract: We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] lumi… ▽ More We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] luminosity, L[OII], using simple relations derived both from the models and observations and also using a public code. This code ideally uses as input instantaneous star formation rates (SFRs), which are only provided by one of the SAMs under consideration. We use this SAM to study the feasibility of inferring galaxies' L[OII] for models that only provide average SFRs. We find that the post-processing computation of L[OII] from average SFRs is accurate for model galaxies with dust attenuated L[OII]$\lesssim10^{42.2}$erg s$^{-1}$ ($<5\%$ discrepancy). We also explore how to derive the [OII] luminosity from simple relations using global properties usually output by SAMs. Besides the SFR, the model L[OII] is best correlated with the observed-frame $u$ and $g$ broad-band magnitudes. These correlations have coefficients (r-values) above 0.64 and a dispersion that varies with L[OII]. We use these correlations and an observational one based on SFR and metallicity to derive L[OII]. These relations result in [OII] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used. Nevertheless, for all the considered models, the amplitude of the clustering at scales above 1$h^{-1}$Mpc remains unchanged independently of the method used to derive L[OII]. △ Less

Submitted 12 February, 2021; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 23 pages, 19 figures, published in MNRAS. Data available at: http://popia.ft.uam.es/MultiDarkEmissionLines/. Dust attenuation code: https://github.com/gfavole/dust

CosmoDC2: A Synthetic Sky Catalog for Dark Energy Science with LSST

Authors: Danila Korytov, Andrew Hearin, Eve Kovacs, Patricia Larsen, Esteban Rangel, Joseph Hollowed, Andrew J. Benson, Katrin Heitmann, Yao-Yuan Mao, Anita Bahmanyar, Chihway Chang, Duncan Campbell, Joseph Derose, Hal Finkel, Nicholas Frontiere, Eric Gawiser, Salman Habib, Benjamin Joachimi, François Lanusse, Nan Li, Rachel Mandelbaum, Christopher Morrison, Jeffrey A. Newman, Adrian Pope, Eli Rykoff , et al. (5 additional authors not shown)

Abstract: This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `… ▽ More This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data. △ Less

Submitted 27 July, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: 27 pages, 17 figures, submitted to APJS

A semi-analytical perspective on massive galaxies at $z\sim0.55$

Authors: D. Stoppacher, F. Prada, A. D. Montero-Dorta, S. Rodríguez-Torres, A. Knebe, G. Favole, W. Cui, A. J. Benson, C. Behrens, A. A. Klypin

Abstract: The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMAS… ▽ More The most massive and luminous galaxies in the Universe serve as powerful probes to study the formation of structure, the assembly of mass, and cosmology. However, their detailed formation and evolution is still barely understood. Here we extract a sample of massive mock galaxies from the semi-analytical model of galaxy formation (SAM) GALACTICUS from the MultiDark-Galaxies, by replicating the CMASS photometric selection from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The comparison of the GALACTICUS CMASS-mock with BOSS-CMASS data allows us to explore different aspects of the massive galaxy population at $0.5<z<0.6$, including the galaxy-halo connection and the galaxy clustering. We find good agreement between our modelled galaxies and observations regarding the galaxy-halo connection, but our CMASS-mock over-estimates the clustering amplitude of the 2-point correlation function, due to a smaller number density compared to BOSS, a lack of blue objects, and a small intrinsic scatter in stellar mass at fixed halo mass of $<0.1$ dex. To alleviate this problem, we construct an alternative mock catalogue mimicking the CMASS colour-magnitude distribution by randomly down-sampling the SAM catalogue. This CMASS-mock reproduces the clustering of CMASS galaxies within 1$σ$ and shows some environmental dependency of star formation properties that could be connected to the quenching of star formation and the assembly bias. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Comments: 15 pages, 10 figures, 2 tables, submitted to MNRAS

Halo concentrations from extended Press-Schechter merger histories

Authors: Andrew J. Benson, Aaron Ludlow, Shaun Cole

Abstract: We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the… ▽ More We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the correct scatter, we postulate that the missing scatter is due to the lack of any environmental dependence in merger trees derived from excursion set theory. We show that a simple modification to the merger tree construction algorithm, which makes merger rates dependent on environment, can increase the scatter by the required amount, and simultaneously provide a qualitatively correct correlation between environment and formation epoch in the excursion set merger trees. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 10 pages, 7 figures, submitted to MNRAS

The effect of dark matter-dark radiation interactions on halo abundance -- a Press-Schechter approach

Authors: Omid Sameie, Andrew J. Benson, Laura V. Sales, Hai-Bo Yu, Leonidas A. Moustakas, Peter Creasey

Abstract: We study halo mass functions with the Press-Schechter formalism for interacting dark matter models, where matter power spectra are damped due to dark acoustic oscillations in the early universe. After adopting a smooth window function, we calibrate the analytical model with numerical simulations from the "effective theory of structure formation" (ETHOS) project and fix the model parameters in the… ▽ More We study halo mass functions with the Press-Schechter formalism for interacting dark matter models, where matter power spectra are damped due to dark acoustic oscillations in the early universe. After adopting a smooth window function, we calibrate the analytical model with numerical simulations from the "effective theory of structure formation" (ETHOS) project and fix the model parameters in the high mass regime, $M_{\rm h}\gtrsim3\times10^{10}\;{\rm M}_{\odot}$. We also perform high-resolution cosmological simulations with halo masses down to $M_{\rm h}\sim10^8\;{\rm M}_{\odot}$ to cover a wide mass range for comparison. Although the model is calibrated with ETHOS1 and CDM simulations for high halo masses at redshift $z=0$, it successfully reproduces simulations for two other ETHOS models in the low mass regime at low and high redshifts. As an application, we compare the cumulative number density of haloes to that of observed galaxies at $z=6$, and find the interacting dark matter models with a kinetic decoupling temperature below $0.5\ \rm{keV}$ is disfavored. We also perform the abundance-matching analysis and derive the stellar-halo mass relation for these models at $z=4$. Suppression in halo abundance leads to less massive haloes that host observed galaxies in the stellar mass range $M_*\simeq 10^5-10^7\ {\rm M}_{\odot}$. △ Less

Submitted 16 April, 2019; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: 9 pages, 7 figures, published version

A compendium of extinction curves for simple galactic geometries

Authors: Andrew J. Benson

Abstract: We calculate the net extinction of galactic light as a function of wavelength, inclination, central optical depth, and morphology for simple galactic geometries using the Hyperion radiative transfer code. Compared to previous, similar works we tabulate extinction over a much broader range of galactic properties, and using a much finer grid in the model parameters. We expect these results to be use… ▽ More We calculate the net extinction of galactic light as a function of wavelength, inclination, central optical depth, and morphology for simple galactic geometries using the Hyperion radiative transfer code. Compared to previous, similar works we tabulate extinction over a much broader range of galactic properties, and using a much finer grid in the model parameters. We expect these results to be useful for constructing dust-extinguished spectra and luminosities of model galaxies and, therefore, for synthetic survey building. Results are made available as an HDF5 file at https://doi.org/10.5281/zenodo.1442826 △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: Accepted by RNAAS

Quantifying the power spectrum of small-scale structure in semi-analytic galaxies

Authors: Sean Brennan, Andrew J. Benson, Francis-Yan Cyr-Racine, Charles R. Keeton, Leonidas A. Moustakas, Anthony R. Pullen

Abstract: In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested th… ▽ More In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested that the full population of subhalos could be probed using a power spectrum analysis. In this paper we quantify the power spectrum of small-scale structure in simulated galaxies, with the goal of understanding theoretical predictions and setting the stage for using measurements of the power spectrum to test dark matter models. We use a sample of simulated galaxies generated from the \texttt{Galacticus} semi-analytic model to determine the power spectrum distribution first in the CDM paradigm and then in a warm dark matter scenario. We find that a measurement of the slope and amplitude of the power spectrum on galaxy strong lensing scales ($k\sim 1$ kpc$^{-1}$) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhalos can be directly detected via gravitational imaging. △ Less

Submitted 29 August, 2019; v1 submitted 10 August, 2018; originally announced August 2018.

Comments: 8 pages, 10 figures

Covariances of Galaxy Stellar Mass Functions and Correlation Functions

Authors: Andrew J. Benson

Abstract: We compute covariance matrices for many observed estimates of the stellar mass function of galaxies from $z=0$ to $z\approx 4$, and for one estimate of the projected correlation function of galaxies split by stellar mass at $z\lesssim 0.5$. All covariance matrices include contributions due to large scale structure, the preference for galaxies to be found in groups and clusters, and for shot noise.… ▽ More We compute covariance matrices for many observed estimates of the stellar mass function of galaxies from $z=0$ to $z\approx 4$, and for one estimate of the projected correlation function of galaxies split by stellar mass at $z\lesssim 0.5$. All covariance matrices include contributions due to large scale structure, the preference for galaxies to be found in groups and clusters, and for shot noise. These covariance matrices are made available for use in constraining models of galaxy formation and the galaxy-halo connection. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: Submitted to MNRAS. Associated data files will be made available when (and if) the paper is accepted for publication

DESCQA: An Automated Validation Framework for Synthetic Sky Catalogs

Authors: Yao-Yuan Mao, Eve Kovacs, Katrin Heitmann, Thomas D. Uram, Andrew J. Benson, Duncan Campbell, Sofía A. Cora, Joseph DeRose, Tiziana Di Matteo, Salman Habib, Andrew P. Hearin, J. Bryce Kalmbach, K. Simon Krughoff, François Lanusse, Zarija Lukić, Rachel Mandelbaum, Jeffrey A. Newman, Nelson Padilla, Enrique Paillas, Adrian Pope, Paul M. Ricker, Andrés N. Ruiz, Ananth Tenneti, Cristian Vega-Martínez, Risa H. Wechsler , et al. (2 additional authors not shown)

Abstract: The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline developm… ▽ More The use of high-quality simulated sky catalogs is essential for the success of cosmological surveys. The catalogs have diverse applications, such as investigating signatures of fundamental physics in cosmological observables, understanding the effect of systematic uncertainties on measured signals and testing mitigation strategies for reducing these uncertainties, aiding analysis pipeline development and testing, and survey strategy optimization. The list of applications is growing with improvements in the quality of the catalogs and the details that they can provide. Given the importance of simulated catalogs, it is critical to provide rigorous validation protocols that enable both catalog providers and users to assess the quality of the catalogs in a straightforward and comprehensive way. For this purpose, we have developed the DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science Collaboration as well as for the broader community. The goal of DESCQA is to enable the inspection, validation, and comparison of an inhomogeneous set of synthetic catalogs via the provision of a common interface within an automated framework. In this paper, we present the design concept and first implementation of DESCQA. In order to establish and demonstrate its full functionality we use a set of interim catalogs and validation tests. We highlight several important aspects, both technical and scientific, that require thoughtful consideration when designing a validation framework, including validation metrics and how these metrics impose requirements on the synthetic sky catalogs. △ Less

Submitted 8 February, 2018; v1 submitted 27 September, 2017; originally announced September 2017.

Comments: 16+4 pages, 9 figures, 2+1 tables. Published in the ApJ Supplement. Code repository at https://github.com/LSSTDESC/descqa

Journal ref: ApJS (2018), 234, 36

Constraining the Noise-Free Distribution of Halo Spin Parameters

Authors: Andrew J. Benson

Abstract: Any measurement made using an N-body simulation is subject to noise due to the finite number of particles used to sample the dark matter distribution function, and the lack of structure below the simulation resolution. This noise can be particularly significant when attempting to measure intrinsically small quantities, such as halo spin. In this work we develop a model to describe the effects of p… ▽ More Any measurement made using an N-body simulation is subject to noise due to the finite number of particles used to sample the dark matter distribution function, and the lack of structure below the simulation resolution. This noise can be particularly significant when attempting to measure intrinsically small quantities, such as halo spin. In this work we develop a model to describe the effects of particle noise on halo spin parameters. This model is calibrated using N-body simulations in which the particle noise can be treated as a Poisson process on the underlying dark matter distribution function, and we demonstrate that this calibrated model reproduces measurements of halo spin parameter error distributions previously measured in N-body convergence studies. Utilizing this model, along with previous measurements of the distribution of halo spin parameters in N-body simulations, we place constraints on the noise-free distribution of halo spins. We find that the noise-free median spin is 3% lower than that measured directly from the N-body simulation, corresponding to a shift of approximately 40 times the statistical uncertainty in this measurement arising purely from halo counting statistics. We also show that measurement of the spin of an individual halo to 10% precision requires at least $4\times 10^4$ particles in the halo - for halos containing 200 particles the fractional error on spins measured for individual halos is of order unity. △ Less

Submitted 4 May, 2017; originally announced May 2017.

Comments: 12 pages, 7 figures, submitted to MNRAS, comments welcome

The phase-space structure of tidally stripped halos

Authors: Nicole E. Drakos, James E. Taylor, Andrew J. Benson

Abstract: We propose a new method for generating equilibrium models of spherical systems of collisionless particles that are finite in extent, but whose central regions resemble dark matter halos from cosmological simulations. This method involves iteratively removing unbound particles from a Navarro-Frenk-White profile truncated sharply at some radius. The resulting models are extremely stable, and thus pr… ▽ More We propose a new method for generating equilibrium models of spherical systems of collisionless particles that are finite in extent, but whose central regions resemble dark matter halos from cosmological simulations. This method involves iteratively removing unbound particles from a Navarro-Frenk-White profile truncated sharply at some radius. The resulting models are extremely stable, and thus provide a good starting point for N-body simulations of isolated halos. We provide a code to generate such models for NFW and a variety of other common density profiles. We then develop an analytic approximation to this truncated distribution function. Our method proceeds by analogy with the King model, truncating and shifting the original distribution function of an infinitely extended Navarro-Frenk-White profile in energy space. We show that the density profiles of our models closely resemble the tidally truncated density profiles seen previously in studies of satellite evolution. Pursuing this analogy further with a series of simulations of tidal mass loss, we find that our models provide a good approximation to the full distribution function of tidally stripped systems, thus allowing theoretically motivated phase-space calculations for such systems. △ Less

Submitted 22 March, 2017; originally announced March 2017.

On the Origin and Evolution of the Galaxy Stellar Mass Function

Authors: Daniel D. Kelson, Andrew J. Benson, Louis E. Abramson

Abstract: Here we explore the evolution of galaxy ensembles at early times by writing the in situ stellar mass growth of galaxies purely as a stationary stochastic (e.g., quasi-steady state) process. By combining the mathematics of such processes with Newtonian gravity and a mean local star formation efficiency, we show that the stellar mass evolution of galaxy ensembles is directly related to the average a… ▽ More Here we explore the evolution of galaxy ensembles at early times by writing the in situ stellar mass growth of galaxies purely as a stationary stochastic (e.g., quasi-steady state) process. By combining the mathematics of such processes with Newtonian gravity and a mean local star formation efficiency, we show that the stellar mass evolution of galaxy ensembles is directly related to the average acceleration of baryons onto dark matter halos at the onset of star formation, with explicit dependencies on initial local matter densities and halo mass. The density term specifically implies more rapid average rates of growth in higher density regions of the universe compared to low density regions, i.e., assembly bias. With this framework, using standard cosmological parameters, a mean star formation efficiency derived by other authors, and knowledge of the shape of the cosmological matter power spectrum at small scales, we analytically derive (1) the characteristic stellar masses of galaxies (M*), (2) the power-law low-mass slope (alpha) and normalization (phi*) of the stellar mass function, and (3) the evolution of the stellar mass function in time over 12.5 > z > 2. Correspondingly, the rise in the cosmic star formation rate density over these epochs, while the universe can sustain unabated fueling of star formation, also emerges naturally. All of our findings are consistent with the deepest available data, including the expectation of alpha~-7/5; i.e., a stellar mass function low-mass slope that is notably shallower than that of the halo mass function, and with no systematic deviations from a mean star formation efficiency with density or mass, nor any explicit, additional feedback mechanisms. These derivations yield a compelling richness and complexity but also show that very few astrophysical details are required to understand the evolution of cosmic ensemble of galaxies at early times. △ Less

Submitted 20 October, 2016; originally announced October 2016.

Comments: 19 pages, 3 figures, submitted to the Astrophysical Journal

The Mass Function of Unprocessed Dark Matter Halos and Merger Tree Branching Rates

Authors: Andrew J. Benson

Abstract: A common approach in semi-analytic modeling of galaxy formation is to construct Monte Carlo realizations of merger histories of dark matter halos whose masses are sampled from a halo mass function. Both the mass function itself, and the merger rates used to construct merging histories are calibrated to N-body simulations. Typically, "backsplash" halos (those which were once subhalos within a large… ▽ More A common approach in semi-analytic modeling of galaxy formation is to construct Monte Carlo realizations of merger histories of dark matter halos whose masses are sampled from a halo mass function. Both the mass function itself, and the merger rates used to construct merging histories are calibrated to N-body simulations. Typically, "backsplash" halos (those which were once subhalos within a larger halo, but which have since moved outside of the halo) are counted in both the halo mass function, and in the merger rates (or, equivalently, progenitor mass functions). This leads to a double-counting of mass in Monte Carlo merger histories which will bias results relative to N-body results. We measure halo mass functions and merger rates with this double-counting removed in a large, cosmological N-body simulation with cosmological parameters consistent with current constraints. Furthermore, we account for the inherently noisy nature of N-body halo mass estimates when fitting functions to N-body data, and show that ignoring these errors leads to a significant systematic bias given the precision statistics available from state-of-the-art N-body cosmological simulations. △ Less

Submitted 4 October, 2016; originally announced October 2016.

Comments: 13 pages, submitted to MNRAS, comments welcome

Achieving Convergence in Galaxy Formation Models by Augmenting N-body Merger Trees

Authors: Andrew J. Benson, Chris Cannella, Shaun Cole

Abstract: Accurate modeling of galaxy formation in a hierarchical, cold dark matter universe requires the use of sufficiently high-resolution merger trees to obtain convergence in the predicted properties of galaxies. When semi-analytic galaxy formation models are applied to cosmological N-body simulation merger trees, it is often the case that those trees have insufficient resolution to give converged gala… ▽ More Accurate modeling of galaxy formation in a hierarchical, cold dark matter universe requires the use of sufficiently high-resolution merger trees to obtain convergence in the predicted properties of galaxies. When semi-analytic galaxy formation models are applied to cosmological N-body simulation merger trees, it is often the case that those trees have insufficient resolution to give converged galaxy properties. We demonstrate a method to augment the resolution of N-body merger trees by grafting in branches of Monte Carlo merger trees with higher resolution, but which are consistent with the pre-existing branches in the N-body tree. We show that this approach leads to converged galaxy properties. △ Less

Submitted 16 August, 2016; v1 submitted 7 April, 2016; originally announced April 2016.

Comments: Accepted for publication in Computational Astrophysics and Cosmology. Version with full resolution figures available at https://users.obs.carnegiescience.edu/abenson/galacticus/augment.pdf

A unified multi-wavelength model of galaxy formation

Authors: Cedric G. Lacey, Carlton M. Baugh, Carlos S. Frenk, Andrew J. Benson, Richard G. Bower, Shaun Cole, Violeta Gonzalez-Perez, John C. Helly, Claudia D. P. Lagos, Peter D. Mitchell

Abstract: We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy di… ▽ More We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 mu. This was not previously achieved by a single physical model within the LambdaCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z=0 to z=6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z=0, as well as far-IR number counts, and far-UV luminosity functions at z ~ 3-6. [abridged] △ Less

Submitted 6 August, 2016; v1 submitted 28 September, 2015; originally announced September 2015.

Comments: Accepted for publication in MNRAS. Revised in response to referee; minor changes only (including some figures), no change in results or conclusions. 62 pages (44 without appendices); 57 figures (32 without appendices)

Sussing Merger Trees: A proposed Merger Tree data format

Authors: Peter A. Thomas, Julian Onions, Dylan Tweed, Andrew J. Benson, Darren Croton, Pascal Elahi, Bruno Henriques, Ilian T. Iliev, Alexander Knebe, Hanni Lux, Yao-Yuan Mao, Mark Neyrinck, Frazer R. Pearce, Vicente Rodriguez-Gomez, Aurel Schneider, Chaichalit Srisawat

Abstract: We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON. We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON. △ Less

Submitted 21 August, 2015; originally announced August 2015.

Comments: Technical report: a proposed merger tree data format. Not submitted to any journal. Source code available at https://bitbucket.org/ProfPAThomas/mergertree

The scaling relations of early-type dwarf galaxies across a range of environments

Authors: Samantha J. Penny, Joachim Janz, Duncan A. Forbes, Andrew J. Benson, Jeremy Mould

Abstract: We present the results of a Keck-ESI study of dwarf galaxies across a range of environment: the Perseus Cluster, the Virgo Cluster, the NGC 1407 group, and the NGC 1023 group. Eighteen dEs are targeted for spectroscopy, three for the first time. We confirm cluster membership for one Virgo dE, and group membership for one dE in the NGC 1023 group, and one dE in the NGC 1407 group for the first time… ▽ More We present the results of a Keck-ESI study of dwarf galaxies across a range of environment: the Perseus Cluster, the Virgo Cluster, the NGC 1407 group, and the NGC 1023 group. Eighteen dEs are targeted for spectroscopy, three for the first time. We confirm cluster membership for one Virgo dE, and group membership for one dE in the NGC 1023 group, and one dE in the NGC 1407 group for the first time. Regardless of environment, the dEs follow the same size-magnitude and $σ$-luminosity relation. Two of the Virgo dwarfs, VCC 1199 and VCC 1627, have among the highest central velocity dispersions ($σ_{0}$ = 58.4 km s$^{-1}$ and 49.2 km s$^{-1}$) measured for dwarfs of their luminosity ($M_{R}\approx -17$). Given their small sizes ($R_{e} < 300$ pc) and large central velocity dispersions, we classify these two dwarfs as compact ellipticals rather than dEs. Group dEs typically have higher mean dynamical-to-stellar mass ratios than the cluster dEs, with $M_{dyn}/M_{\star} = 5.1\pm0.6$ for the group dwarfs, vs. $M_{dyn}/M_{\star} = 2.2\pm0.5$ for the cluster sample, which includes two cEs. We also search for trends in $M_{dyn}/M_{\star}$ vs. distance from M87 for the Virgo Cluster population, and find no preference for dwarfs with high values of $M_{dyn}/M_{\star}$ to reside in the cluster outskirts vs. centre. △ Less

Submitted 18 August, 2015; originally announced August 2015.

Comments: 15 pages, 8 figures, 5 tables. Accepted for publication in MNRAS

The Dark Matter Halos of Moderate Luminosity X-ray AGN as Determined from Weak Gravitational Lensing and Host Stellar Masses

Authors: A. Leauthaud, A. J. Benson, F. Civano, A. L. Coil, K. Bundy, R. Massey, M. Schramm, A. Schulze, P. Capak, M. Elvis, A. Kulier, J. Rhodes

Abstract: Understanding the relationship between galaxies hosting active galactic nuclei (AGN) and the dark matter halos in which they reside is key to constraining how black-hole fueling is triggered and regulated. Previous efforts have relied on simple halo mass estimates inferred from clustering, weak gravitational lensing, or halo occupation distribution modeling. In practice, these approaches remain un… ▽ More Understanding the relationship between galaxies hosting active galactic nuclei (AGN) and the dark matter halos in which they reside is key to constraining how black-hole fueling is triggered and regulated. Previous efforts have relied on simple halo mass estimates inferred from clustering, weak gravitational lensing, or halo occupation distribution modeling. In practice, these approaches remain uncertain because AGN, no matter how they are identified, potentially live a wide range of halo masses with an occupation function whose general shape and normalization are poorly known. In this work, we show that better constraints can be achieved through a rigorous comparison of the clustering, lensing, and cross-correlation signals of AGN hosts to a fiducial stellar-to-halo mass relation (SHMR) derived for all galaxies. Our technique exploits the fact that the global SHMR can be measured with much higher accuracy than any statistic derived from AGN samples alone. Using 382 moderate luminosity X-ray AGN at z<1 from the COSMOS field, we report the first measurements of weak gravitational lensing from an X-ray selected sample. Comparing this signal to predictions from the global SHMR, we find that, contrary to previous results, most X-ray AGN do not live in medium size groups ---nearly half reside in relatively low mass halos with Mh~10^12.5 Msun. The AGN occupation function is well described by the same form derived for all galaxies but with a lower normalization---the fraction of halos with AGN in our sample is a few percent. By highlighting the relatively "normal" way in which moderate luminosity X-ray AGN hosts occupy halos, our results suggest that the environmental signature of distinct fueling modes for luminous QSOs compared to moderate luminosity X-ray AGN is less obvious than previously claimed. △ Less

Submitted 21 October, 2014; originally announced October 2014.

Comments: 15 pages, accepted to MNRAS

Trends in Dwarf Early-Type Kinematics with Cluster-centric Radius Driven By Tidal Stirring

Authors: A. J. Benson, E. Toloba, L. Mayer, J. D. Simon, P. Guhathakurta

Abstract: We model the dynamics of dwarf early-type galaxies in the Virgo cluster when subject to a variety of environmental processes. We focus on how these processes imprint trends in dynamical state (rotational vs. pressure support as measured by the $λ^*_{\rm Re/2}$ statistic) with projected distance from the cluster center, and compare these results to observational estimates. We find a large scatter i… ▽ More We model the dynamics of dwarf early-type galaxies in the Virgo cluster when subject to a variety of environmental processes. We focus on how these processes imprint trends in dynamical state (rotational vs. pressure support as measured by the $λ^*_{\rm Re/2}$ statistic) with projected distance from the cluster center, and compare these results to observational estimates. We find a large scatter in the gradient of $λ^*_{\rm Re/2}$ with projected radius. A statistical analysis shows that models with no environmental effects produce gradients as steep as those observed in none of the 100 cluster realizations we consider, while in a model incorporating tidal stirring by the cluster potential 34% of realizations produce gradients as steep as that observed. Our results suggest that tidal stirring may be the cause of the observed radial dependence of dwarf early-type dynamics in galaxy clusters. △ Less

Submitted 20 October, 2014; originally announced October 2014.

Comments: 5 pages, accepted for publication in ApJ

Nonlinear evolution of dark matter subhalos and applications to warm dark matter

Authors: Anthony R. Pullen, Andrew J. Benson, Leonidas A. Moustakas

Abstract: We describe the methodology to include nonlinear evolution, including tidal effects, in the computation of subhalo distribution properties in both cold (CDM) and warm (WDM) dark matter universes. Using semi-analytic modeling, we include effects from dynamical friction, tidal stripping, and tidal heating, allowing us to dynamically evolve the subhalo distribution. We calibrate our nonlinear evoluti… ▽ More We describe the methodology to include nonlinear evolution, including tidal effects, in the computation of subhalo distribution properties in both cold (CDM) and warm (WDM) dark matter universes. Using semi-analytic modeling, we include effects from dynamical friction, tidal stripping, and tidal heating, allowing us to dynamically evolve the subhalo distribution. We calibrate our nonlinear evolution scheme to the CDM subhalo mass function in the Aquarius N-body simulation, producing a subhalo mass function within the range of simulations. We find tidal effects to be the dominant mechanism of nonlinear evolution in the subhalo population. Finally, we compute the subhalo mass function for $m_χ=1.5$ keV WDM including the effects of nonlinear evolution, and compare radial number densities and mass density profiles of subhalos in CDM and WDM models. We show that all three signatures differ between the two dark matter models, suggesting that probes of substructure may be able to differentiate between them. △ Less

Submitted 14 August, 2014; v1 submitted 30 July, 2014; originally announced July 2014.

Comments: Added a couple of references

Galaxy Formation as a Cosmological Tool. I: The Galaxy Merger History as a Measure of Cosmological Parameters

Authors: Christopher J. Conselice, Asa F. L. Bluck, Alice Mortlock, David Palamara, Andrew J. Benson

Abstract: As galaxy formation and evolution over long cosmic time-scales depends to a large degree on the structure of the universe, the assembly history of galaxies is potentially a powerful approach for learning about the universe itself. In this paper we examine the merger history of dark matter halos based on the Extended Press-Schechter formalism as a function of cosmological parameters, redshift and h… ▽ More As galaxy formation and evolution over long cosmic time-scales depends to a large degree on the structure of the universe, the assembly history of galaxies is potentially a powerful approach for learning about the universe itself. In this paper we examine the merger history of dark matter halos based on the Extended Press-Schechter formalism as a function of cosmological parameters, redshift and halo mass. We calculate how major halo mergers are influenced by changes in the cosmological values of $Ω_{\rm m}$, $Ω_Λ$, $σ_{8}$, the dark matter particle temperature (warm vs. cold dark matter), and the value of a constant and evolving equation of state parameter $w(z)$. We find that the merger fraction at a given halo mass varies by up to a factor of three for halos forming under the assumption of Cold Dark Matter, within different underling cosmological parameters. We find that the current measurements of the merger history, as measured through observed galaxy pairs as well as through structure, are in agreement with the concordance cosmology with the current best fit giving $1 - Ω_{\rm m} = Ω_{\rm Λ} = 0.84^{+0.16}_{-0.17}$. To obtain a more accurate constraint competitive with recently measured cosmological parameters from Planck and WMAP requires a measured merger accuracy of $δf_{\rm m} \sim 0.01$, implying surveys with an accurately measured merger history over 2 - 20 deg$^{2}$, which will be feasible with the next generation of imaging and spectroscopic surveys such as Euclid and LSST. △ Less

Submitted 14 July, 2014; originally announced July 2014.

Comments: MNRAS, accepted for publication, 22 pages

Building a Predictive Model of Galaxy Formation - I: Phenomenological Model Constrained to the $z=0$ Stellar Mass Function

Authors: A. J. Benson

Abstract: We constrain a highly simplified semi-analytic model of galaxy formation using the $z\approx 0$ stellar mass function of galaxies. Particular attention is paid to assessing the role of random and systematic errors in the determination of stellar masses, to systematic uncertainties in the model, and to correlations between bins in the measured and modeled stellar mass functions, in order to constru… ▽ More We constrain a highly simplified semi-analytic model of galaxy formation using the $z\approx 0$ stellar mass function of galaxies. Particular attention is paid to assessing the role of random and systematic errors in the determination of stellar masses, to systematic uncertainties in the model, and to correlations between bins in the measured and modeled stellar mass functions, in order to construct a realistic likelihood function. We derive constraints on model parameters and explore which aspects of the observational data constrain particular parameter combinations. We find that our model, once constrained, provides a remarkable match to the measured evolution of the stellar mass function to $z=1$, although fails dramatically to match the local galaxy HI mass function. Several "nuisance parameters" contribute significantly to uncertainties in model predictions. In particular, systematic errors in stellar mass estimate are the dominant source of uncertainty in model predictions at $z\approx 1$, with additional, non-negligble contributions arising from systematic uncertainties in halo mass functions and the residual uncertainties in cosmological parameters. Ignoring any of these sources of uncertainties could lead to viable models being erroneously ruled out. Additionally, we demonstrate that ignoring the significant covariance between bins the observed stellar mass function leads to significant biases in the constraints derived on model parameters. Careful treatment of systematic and random errors in the constraining data, and in the model being constrained, are crucial if this methodology is to be used to test hypotheses relating to the physics of galaxy formation. △ Less

Submitted 21 May, 2014; originally announced May 2014.

Comments: 40 pages, submitted to MNRAS, comments welcome, version with full resolution figures available at http://users.obs.carnegiescience.edu/abenson/galacticus/galacticusConstraints1.pdf

The Circular Velocity Function of Group Galaxies

Authors: Louis E. Abramson, Rik J. Williams, Andrew J. Benson, Juna A. Kollmeier, John S. Mulchaey

Abstract: A robust prediction of $Λ{\rm CDM}$ cosmology is the halo circular velocity function (CVF), a dynamical cousin of the halo mass function. The correspondence between theoretical and observed CVFs is uncertain, however: cluster galaxies are reported to exhibit a power-law CVF consistent with $N$-body simulations, but that of the field is distinctly Schechter-like, flattened relative to $Λ{\rm CDM}$… ▽ More A robust prediction of $Λ{\rm CDM}$ cosmology is the halo circular velocity function (CVF), a dynamical cousin of the halo mass function. The correspondence between theoretical and observed CVFs is uncertain, however: cluster galaxies are reported to exhibit a power-law CVF consistent with $N$-body simulations, but that of the field is distinctly Schechter-like, flattened relative to $Λ{\rm CDM}$ expectations at circular velocities $v_c \lesssim 200\, {\rm km\, s^{-1}}$. Groups offer a powerful probe of the role environment plays in this discrepancy as they bridge the field and clusters. Here, we construct the CVF for a large, mass- and multiplicity-complete sample of group galaxies from the Sloan Digital Sky Survey. Using independent photometric $v_c$ estimators, we find no transition from a field- to $Λ{\rm CDM}$-shaped CVF above $v_c = 50\, {\rm km\, s^{-1}}$ as a function of group halo mass. All groups with $12.4 \lesssim \log M_{\rm halo} / M_{\odot} \lesssim 15.1$ (Local Group analogs to rich clusters) display similar Schechter-like CVFs marginally suppressed at low-$v_c$ compared to that of the field. Conversely, some agreement with $N$-body results emerges for samples saturated with late-type galaxies, with isolated late-types displaying a CVF similar in shape to $Λ{\rm CDM}$ predictions. We conclude that the flattening of the low-$v_c$ slope in groups is due to their depressed late-type fractions -- environment affecting the CVF only to the extent that it correlates with this quantity -- and that previous cluster analyses may suffer from interloper contamination. These results serve as useful benchmarks for cosmological simulations of galaxy formation. △ Less

Submitted 5 August, 2014; v1 submitted 24 July, 2013; originally announced July 2013.

Comments: 16 pages, 16 figures; Accepted to ApJ August 2014. ArXiv version updated to reflect changes made during refereeing; results unchanged

Constraints on black hole fuelling modes from the clustering of X-ray AGN

Authors: N. Fanidakis, A. Georgakakis, G. Mountrichas, M. Krumpe, C. M. Baugh, C. G. Lacey, C. S. Frenk, T. Miyaji, A. J. Benson

Abstract: We present a clustering analysis of X-ray selected AGN by compiling X-ray samples from the literature and re-estimating the dark matter (DM) halo masses of AGN in a uniform manner. We find that moderate luminosity AGN (Lx(2-10 keV)=10^42-10^44 erg/sec) in the z=0-1.3 Universe are typically found in DM haloes with masses of ~10^13 Msun. We then compare our findings to the theoretical predictions of… ▽ More We present a clustering analysis of X-ray selected AGN by compiling X-ray samples from the literature and re-estimating the dark matter (DM) halo masses of AGN in a uniform manner. We find that moderate luminosity AGN (Lx(2-10 keV)=10^42-10^44 erg/sec) in the z=0-1.3 Universe are typically found in DM haloes with masses of ~10^13 Msun. We then compare our findings to the theoretical predictions of the coupled galaxy and black hole formation model GALFORM. We find good agreement when our calculation includes the hot-halo mode of accretion onto the central black hole. This type of accretion, which is additional to the common cold accretion during disk instabilities and galaxy mergers, is tightly coupled to the AGN feedback in the model. The hot-halo mode becomes prominent in DM haloes with masses greater than ~10^12.5 Msun, where AGN feedback typically operates, giving rise to a distinct class of moderate luminosity AGN that inhabit rich clusters and superclusters. Cold gas fuelling of the black hole cannot produce the observationally inferred DM halo masses of X-ray AGN. Switching off AGN feedback in the model results in a large population of luminous quasars (Lx(2-10 keV) > 10^44 erg/sec) in DM haloes with masses up to ~10^14 Msun, which is inconsistent with the observed clustering of quasars. The abundance of hot-halo AGN decreases significantly in the z~3-4 universe. At such high redshifts, the cold accretion mode is solely responsible for shaping the environment of moderate luminosity AGN. Our analysis supports two accretion modes (cold and hot) for the fuelling of supermassive black holes and strongly underlines the importance of AGN feedback in cosmological models both of galaxy formation and black hole growth. △ Less

Submitted 9 May, 2013; originally announced May 2013.

Comments: 10 pages, 4 figures, submitted to MNRAS (comments are welcome)

Excursion Set Theory for Correlated Random Walks

Authors: Arya Farahi, Andrew J. Benson

Abstract: We present a new method to compute the first crossing distribution in excursion set theory for the case of correlated random walks. We use a combination of the path integral formalism of Maggiore & Riotto, and the integral equation solution of Zhang & Hui, and Benson et al. to find a numerically robust and convenient algorithm to derive the first crossing distribution in terms of a perturbative ex… ▽ More We present a new method to compute the first crossing distribution in excursion set theory for the case of correlated random walks. We use a combination of the path integral formalism of Maggiore & Riotto, and the integral equation solution of Zhang & Hui, and Benson et al. to find a numerically robust and convenient algorithm to derive the first crossing distribution in terms of a perturbative expansion around the limit of an uncorrelated random walk. We apply this methodology to the specific case of a Gaussian random density field filtered with a Gaussian smoothing function. By comparing our solutions to results from Monte Carlo calculations of the first crossing distribution we demonstrate that our method accurate for power spectra $P(k)\propto k^n$ for $n=1$, becoming less accurate for smaller values of $n$. It is therefore complementary to the method of Musso & Sheth, which will therefore be more useful for standard $Λ$CDM power spectra. Our approach is quite general, and can be adapted to other smoothing functions, and also to non-Gaussian density fields. △ Less

Submitted 8 June, 2013; v1 submitted 1 March, 2013; originally announced March 2013.

Comments: 15 pages, 7 figures, Accepted for publication in MNRAS. Some corrections following comments from referee

MSC Class: 85-08 ACM Class: J.2

Journal ref: MNRAS; 2013; 433 (4): 3428-3439

A Research Note on the Implementation of Star Formation and Stellar Feedback in Semi-Analytic Models

Authors: Fabio Fontanot, Gabriella De Lucia, Andrew J. Benson, Pierluigi Monaco, Michael Boylan-Kolchin

Abstract: We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, disc instabilities or AGN feedback. We run these versions… ▽ More We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, disc instabilities or AGN feedback. We run these versions on identical samples of dark matter (DM) haloes extracted from high-resolution N-body simulations in order to perform both statistical analysis and object-by-object comparisons. We compare our results with previous work based on stripped-down versions of the same SAMs including only gas cooling, and show that all feedback models provide coherent modifications in the distribution of baryons between the various gas phases. In particular, we find that the predicted hot gas fractions are considerably increased by up to a factor of three, while the corresponding cold gas fractions are correspondingly decreased, and a significant amount of mass is ejected from the DM halo. Nonetheless, we also find relevant differences in the predicted properties of model galaxies among the three SAMs: these deviations are more relevant at mass scales comparable to that of our own Galaxy, and are reduced at larger masses, confirming the varying impact of stellar feedback at different mass scales. We also check the effect of enhanced star formation events (i.e. starbursts modes), defined in connection with galaxy mergers. We find that, in general, these episodes have a limited impact in the overall star formation histories of model galaxies, even in massive DM halos where merger-driven star formation has often been considered very important. △ Less

Submitted 17 January, 2013; originally announced January 2013.

Comments: 19 pages, 7 figures, Research Note, Comments welcome

Dark Matter Halo Merger Histories Beyond Cold Dark Matter: I - Methods and Application to Warm Dark Matter

Authors: Andrew J. Benson, Arya Farahi, Shaun Cole, Leonidas A. Moustakas, Adrian Jenkins, Mark Lovell, Rachel Kennedy, John Helly, Carlos Frenk

Abstract: We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle propertie… ▽ More We describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos. △ Less

Submitted 7 October, 2012; v1 submitted 13 September, 2012; originally announced September 2012.

Comments: 19 pages, 14 figures, MNRAS in press. Comments welcome

Threshing in Action - The tidal disruption of a dwarf galaxy by the Hydra I Cluster

Authors: Andreas Koch, Andreas Burkert, R. Michael Rich, Michelle L. M. Collins, Christine S. Black, Michael Hilker, Andrew J. Benson

Abstract: We report on the discovery of strong tidal features around a dwarf spheroidal galaxy in the Hydra I galaxy cluster, indicating its ongoing tidal disruption. This very low surface brightness object, HCC-087, was originally classified as an early-type dwarf in the Hydra Cluster Catalogue (HCC), but our re-analysis of the ESO-VLT/FORS images of the HCC unearthed a clear indication of an S-shaped morp… ▽ More We report on the discovery of strong tidal features around a dwarf spheroidal galaxy in the Hydra I galaxy cluster, indicating its ongoing tidal disruption. This very low surface brightness object, HCC-087, was originally classified as an early-type dwarf in the Hydra Cluster Catalogue (HCC), but our re-analysis of the ESO-VLT/FORS images of the HCC unearthed a clear indication of an S-shaped morphology and a large spatial extent. Its shape, luminosity (M_V=-11.6 mag), and physical size (at a half-light radius of 3.1 kpc and a full length of ~5.9 kpc) are comparable to the recently discovered NGC 4449B and the Sagittarius dwarf spheroidal, all of which are undergoing clear tidal disruption. Aided by N-body simulations we argue that HCC-087 is currently at its first apocenter, at 150 kpc, around the cluster center and that it is being tidally disrupted by the galaxy cluster's potential itself. An interaction with the near-by (50 kpc) S0 cluster galaxy HCC-005, at M* ~ 3 x 10^10 M_sun is rather unlikely, as this constellation requires a significant amount of dynamical friction and thus low relative velocities. The S-shaped morphology and large spatial extent of the satellite would, however, also appear if HCC-087 would orbit the cluster center. These features appear to be characteristic properties of satellites that are seen in the process of being tidally disrupted, independent of the environment of the destruction. An important finding of our simulations is an orientation of the tidal tails perpendicular to the orbit. △ Less

Submitted 11 July, 2012; originally announced July 2012.

Comments: 5 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters. Some figure sizes reduced