Showing 1–50 of 101 results for author: Zentner, A R

Subhalos are Distributed Anisotropically About Their Hosts

Authors: Lorena Mezini, Andrew R. Zentner, Kuan Wang, Catherine Fielder

Abstract: We investigate the distributions of subhalos about their hosts in two suites of zoom-in N-body simulations of halo growth -- one suite focused on Milky Way Mass halos ($\sim 10^{12} \mathrm{M}_{\odot}$) and another focused on cluster ($\sim 10^{15} \mathrm{M}_{\odot}$) halos in the Symphony simulation suite. We find, in agreement with previous work on this subject, that subhalos are distributed an… ▽ More We investigate the distributions of subhalos about their hosts in two suites of zoom-in N-body simulations of halo growth -- one suite focused on Milky Way Mass halos ($\sim 10^{12} \mathrm{M}_{\odot}$) and another focused on cluster ($\sim 10^{15} \mathrm{M}_{\odot}$) halos in the Symphony simulation suite. We find, in agreement with previous work on this subject, that subhalos are distributed anisotropically about their host halos. In particular, the positions of subhalos lie preferentially near the major axes of their host halos, possibly implying that satellite galaxies will exhibit a similar alignment. Furthermore, we show that in two-dimensional projection subhalos are more likely to be observed near the halo center (where the central galaxy presumably resides) when the host halo is projected nearly along its major axis. This projection effect is significant. Within a projected radius of $5\%$ of the virial radius of the host halo, the fraction of mass in subhalos is $\sim 44\%$ larger for Milky Way mass halos and as much as $\sim 145\%$ larger for cluster halos when projected along the major axis as compared to the average from a random projection. This result has consequences for many applications including the interpretation of gravitational lenses. Finally, we find that the orbital angular momentum vector of subhalos is aligned with the angular momentum vector of their host halo, indicating that a significant component of a halo's angular momentum may be carried in its subhalos. This has consequences for galaxy formation models which use host halo angular momentum as a proxy for galaxy momentum. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: submitted to MNRAS

The DESI One-Percent Survey: Evidence for Assembly Bias from Low-Redshift Counts-in-Cylinders Measurements

Authors: Alan N. Pearl, Andrew R. Zentner, Jeffrey A. Newman, Rachel Bezanson, Kuan Wang, John Moustakas, Jessica N. Aguilar, Steven Ahlen, David Brooks, Todd Claybaugh, Shaun Cole, Kyle Dawson, Axel de la Macorra, Peter Doel, Jamie E. Forero-Romero, Satya Gontcho A Gontcho, Klaus Honscheid, Martin Landriau, Marc Manera, Paul Martini Aaron Meisner, Ramon Miquel, Jundan Nie, Will Percival, Francisco Prada, Mehdi Rezaie , et al. (6 additional authors not shown)

Abstract: We explore the galaxy-halo connection information that is available in low-redshift samples from the early data release of the Dark Energy Spectroscopic Instrument (DESI). We model the halo occupation distribution (HOD) from z=0.1-0.3 using Survey Validation 3 (SV3; a.k.a., the One-Percent Survey) data of the DESI Bright Galaxy Survey (BGS). In addition to more commonly used metrics, we incorporat… ▽ More We explore the galaxy-halo connection information that is available in low-redshift samples from the early data release of the Dark Energy Spectroscopic Instrument (DESI). We model the halo occupation distribution (HOD) from z=0.1-0.3 using Survey Validation 3 (SV3; a.k.a., the One-Percent Survey) data of the DESI Bright Galaxy Survey (BGS). In addition to more commonly used metrics, we incorporate counts-in-cylinders (CiC) measurements, which drastically tighten HOD constraints. Our analysis is aided by the Python package, galtab, which enables the rapid, precise prediction of CiC for any HOD model available in halotools. This methodology allows our Markov chains to converge with much fewer trial points, and enables even more drastic speedups due to its GPU portability. Our HOD fits constrain characteristic halo masses tightly and provide statistical evidence for assembly bias, especially at lower luminosity thresholds: the HOD of central galaxies in $z\sim0.15$ samples with limiting absolute magnitude $M_r < -20.0$ and $M_r < -20.5$ samples is positively correlated with halo concentration with a significance of 99.9% and 99.5%, respectively. Our models also favor positive central assembly bias for the brighter $M_r < -21.0$ sample at $z\sim0.25$ (94.8% significance), but there is no significant evidence for assembly bias with the same luminosity threshold at $z\sim0.15$. We provide our constraints for each threshold sample's characteristic halo masses, assembly bias, and other HOD parameters. These constraints are expected to be significantly tightened with future DESI data, which will span an area 100 times larger than that of SV3. △ Less

Submitted 15 September, 2023; originally announced September 2023.

The Influence of Subhaloes on Host Halo Properties

Authors: Lorena Mezini, Catherine E. Fielder, Andrew R. Zentner, Yao-Yuan Mao, Kuan Wang, Hao-Yi Wu

Abstract: Within the $Λ$CDM cosmology, dark matter haloes are comprised of both a smooth component and a population of smaller, gravitationally bound subhaloes. These components are often treated as a single halo when halo properties, such as density profiles, are extracted from simulations. Recent work has shown that density profiles change substantially when subhalo mass is excluded. In this paper, we exp… ▽ More Within the $Λ$CDM cosmology, dark matter haloes are comprised of both a smooth component and a population of smaller, gravitationally bound subhaloes. These components are often treated as a single halo when halo properties, such as density profiles, are extracted from simulations. Recent work has shown that density profiles change substantially when subhalo mass is excluded. In this paper, we expand on this result by analysing the change in three specific host halo properties -- concentration ($c_{\rm{NFW}}$), spin ($λ_{\rm Bullock}$), and shape ($c/a$), -- when calculated only from the smooth component of the halo. This analysis is performed on both Milky Way-mass haloes and cluster-mass haloes in high-resolution, zoom-in, $N$-body simulations. We find that when subhaloes are excluded the median value of (1) $c_{\rm{NFW}}$ is enhanced by $\approx 38 \pm 12\%$ and $\approx 88 \pm 7.7\%$ for Milky Way mass ($10^{12.1}\,\text{M}_\odot$) and cluster mass ($10^{14.8}\,\text{M}_\odot$) haloes respectively, (2) $λ_{\rm Bullock}$ is reduced for Milky Way mass by $\approx 16 \pm 6.8\%$ and cluster mass haloes by $\approx 32 \pm 8.9\%$. Additionally, with the removal of subhaloes, cluster mass haloes tend to become more spherical as the ratio of minor-to-major axis, $c/a$, increases by $\approx 12 \pm 4\%$, whereas Milky Way mass haloes remain approximately the same shape with $c/a$ changed by $\approx 1.2 \pm 5.6\%$. The fractional change of each of these properties depends primarily on the amount of mass that is removed from the halo system and, to a lesser extent, mass accretion history. Our findings demonstrate that the properties of the smooth components of dark matter haloes are biased relative to the total mass of the halo including subhaloes. △ Less

Submitted 4 October, 2023; v1 submitted 26 April, 2023; originally announced April 2023.

Comments: Accepted to MNRAS 09/26/2023

Evidence of Galaxy Assembly Bias in SDSS DR7 Galaxy Samples from Count Statistics

Authors: Kuan Wang, Yao-Yuan Mao, Andrew R. Zentner, Hong Guo, Johannes U. Lange, Frank C. van den Bosch, Lorena Mezini

Abstract: We present observational constraints on the galaxy-halo connection, focusing particularly on galaxy assembly bias, from a novel combination of counts-in-cylinders statistics, $P(N_{\rm{CIC}})$, with the standard measurements of the projected two-point correlation function, $w_{\rm{p}}(r_{\rm{p}})$, and number density, $n_{\rm{gal}}$, of galaxies. We measure $n_{\rm{gal}}$,… ▽ More We present observational constraints on the galaxy-halo connection, focusing particularly on galaxy assembly bias, from a novel combination of counts-in-cylinders statistics, $P(N_{\rm{CIC}})$, with the standard measurements of the projected two-point correlation function, $w_{\rm{p}}(r_{\rm{p}})$, and number density, $n_{\rm{gal}}$, of galaxies. We measure $n_{\rm{gal}}$, $w_{\rm{p}}(r_{\rm{p}})$ and $P(N_{\rm{CIC}})$ for volume-limited, luminosity-threshold samples of galaxies selected from SDSS DR7, and use them to constrain halo occupation distribution (HOD) models, including a model in which galaxy occupation depends upon a secondary halo property, namely halo concentration. We detect significant positive central assembly bias for the $M_r<-20.0$ and $M_r<-19.5$ samples. Central galaxies preferentially reside within haloes of high concentration at fixed mass. Positive central assembly bias is also favoured in the $M_r<-20.5$ and $M_r<-19.0$ samples. We find no evidence of central assembly bias in the $M_r<-21.0$ sample. We observe only a marginal preference for negative satellite assembly bias in the $M_r<-20.0$ and $M_r<-19.0$ samples, and non-zero satellite assembly bias is not indicated in other samples. Our findings underscore the necessity of accounting for galaxy assembly bias when interpreting galaxy survey data, and demonstrate the potential of count statistics in extracting information from the spatial distribution of galaxies, which could be applied to both galaxy-halo connection studies and cosmological analyses. △ Less

Submitted 20 September, 2022; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: 16+5 pages, 9+3 figures, 5+2 tables, Fig. 6 shows the main result. To be submitted to MNRAS, comments welcome

CLIMBER: Galaxy-Halo Connection Constraints from Next-Generation Surveys

Authors: Alan N. Pearl, Rachel Bezanson, Andrew R. Zentner, Jeffrey A. Newman, Andy D. Goulding, Katherine E. Whitaker, Sean D. Johnson, Jenny E. Greene

Abstract: In the coming decade, a new generation of massively multiplexed spectroscopic surveys, such as PFS, WAVES, and MOONS, will probe galaxies in the distant universe in vastly greater numbers than was previously possible. In this work, we generate mock catalogs for each of these three planned surveys to help quantify and optimize their scientific output. To assign photometry into the UniverseMachine e… ▽ More In the coming decade, a new generation of massively multiplexed spectroscopic surveys, such as PFS, WAVES, and MOONS, will probe galaxies in the distant universe in vastly greater numbers than was previously possible. In this work, we generate mock catalogs for each of these three planned surveys to help quantify and optimize their scientific output. To assign photometry into the UniverseMachine empirical model, we develop the Calibrating Light: Illuminating Mocks By Empirical Relations (CLIMBER) procedure using UltraVISTA photometry. Using the published empirical selection functions for each aforementioned survey, we quantify the mass completeness of each survey. We compare different targeting strategies by varying the area and targeting completeness, and quantify how these survey parameters affect the uncertainty of the two-point correlation function. We demonstrate that the PFS and MOONS measurements will be primarily dominated by cosmic variance, not shot noise, motivating the need for increasingly large survey areas. On the other hand, the WAVES survey, which covers a much larger area, will strike a good balance between cosmic variance and shot noise. For a fixed number of targets, a 5% increased survey area (and $\sim$5% decreased completeness) would decrease the uncertainty of the correlation function at intermediate scales by 0.15%, 1.2%, and 1.1% for our WAVES, PFS, and MOONS samples, respectively. Meanwhile, for a fixed survey area, 5% increased targeting completeness improves the same constraints by 0.7%, 0.25%, and 0.1%. All of the utilities used to construct our mock catalogs and many of the catalogs themselves are publicly available. △ Less

Submitted 30 November, 2021; originally announced December 2021.

The Effects of Asymmetric Dark Matter on Stellar Evolution I: Spin-Dependent Scattering

Authors: Troy J. Raen, Héctor Martínez-Rodríguez, Travis J. Hurst, Andrew R. Zentner, Carles Badenes, Rachel Tao

Abstract: Most of the dark matter (DM) search over the last few decades has focused on WIMPs, but the viable parameter space is quickly shrinking. Asymmetric Dark Matter (ADM) is a WIMP-like DM candidate with slightly smaller masses and no present day annihilation, meaning that stars can capture and build up large quantities. The captured ADM can transport energy through a significant volume of the star. We… ▽ More Most of the dark matter (DM) search over the last few decades has focused on WIMPs, but the viable parameter space is quickly shrinking. Asymmetric Dark Matter (ADM) is a WIMP-like DM candidate with slightly smaller masses and no present day annihilation, meaning that stars can capture and build up large quantities. The captured ADM can transport energy through a significant volume of the star. We investigate the effects of spin-dependent ADM energy transport on stellar structure and evolution in stars with $0.9 \leq M_{\star}/\mathrm{M}_{\odot} \leq 5.0$ in varying DM environments. We wrote a MESA module that calculates the capture of DM and the subsequent energy transport within the star. We fix the DM mass to 5 GeV and the cross section to $10^{-37} \mathrm{cm{^2}}$, and study varying environments by scaling the DM capture rate. For stars with radiative cores ($M_{\star} \lesssim 1.3\ \mathrm{M}_{\odot}$), the presence of ADM flattens the temperature and burning profiles in the core and increases MS ($X_c > 10^{-3}$) lifetimes by up to $\sim 20\%$. We find that strict requirements on energy conservation are crucial to the simulation of ADM's effects on these stars. In higher-mass stars, ADM energy transport shuts off core convection, limiting available fuel and shortening MS lifetimes by up to $\sim 40\%$. This may translate to changes in the luminosity and effective temperature of the MS turnoff in population isochrones. The tip of the red giant branch may occur at lower luminosities. The effects are largest in DM environments with high densities and/or low velocity dispersions, making dwarf and early forming galaxies most likely to display the effects. △ Less

Submitted 16 March, 2021; v1 submitted 8 October, 2020; originally announced October 2020.

Comments: 14 pages, 8 figures, submitted to Monthly Notices of the Royal Astronomical Society

Illuminating Dark Matter Halo Density Profiles Without Subhaloes

Authors: Catherine E. Fielder, Yao-Yuan Mao, Andrew R. Zentner, Jeffrey A. Newman, Hao-Yi Wu, Risa Wechsler

Abstract: Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all halo mass, including mass in subhaloes, in studies of halo density profiles. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host.… ▽ More Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all halo mass, including mass in subhaloes, in studies of halo density profiles. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host. This difference makes comparisons between theoretical and observed results difficult. In this work we investigate density profiles of the smooth components of host haloes by excluding mass contained within subhaloes. We find that the density profiles of the smooth halo component (without subhaloes) differs substantially from the conventional halo density profile. Smooth profiles decline more rapidly at large radii and are not well characterised by the standard NFW profile. We also find that concentrations derived from smooth density profiles exhibit less scatter at fixed mass and a weaker mass dependence than standard concentrations. Both smooth and standard halo profiles can be described by a generalised Einasto profile, an Einasto profile with a modified central slope, with smaller residuals than either an NFW or Einasto profile. These results hold for both Milky Way-mass and cluster-mass haloes. This new characterisation of smooth halo profiles can be useful for many analyses, such as lensing and dark matter annihilation, in which the smooth and clumpy components of a halo should be accounted for separately. △ Less

Submitted 6 July, 2020; originally announced July 2020.

Comments: 19 pages, 9 figures, appendices, submitted to MNRAS

Concentrations of Dark Haloes Emerge from Their Merger Histories

Authors: Kuan Wang, Yao-Yuan Mao, Andrew R. Zentner, Johannes U. Lange, Frank C. van den Bosch, Risa H. Wechsler

Abstract: The concentration parameter is a key characteristic of a dark matter halo that conveniently connects the halo's present-day structure with its assembly history. Using 'Dark Sky', a suite of cosmological $N$-body simulations, we investigate how halo concentration evolves with time and emerges from the mass assembly history. We also explore the origin of the scatter in the relation between concentra… ▽ More The concentration parameter is a key characteristic of a dark matter halo that conveniently connects the halo's present-day structure with its assembly history. Using 'Dark Sky', a suite of cosmological $N$-body simulations, we investigate how halo concentration evolves with time and emerges from the mass assembly history. We also explore the origin of the scatter in the relation between concentration and assembly history. We show that the evolution of halo concentration has two primary modes: (1) smooth increase due to pseudo-evolution; and (2) intense responses to physical merger events. Merger events induce lasting and substantial changes in halo structures, and we observe a universal response in the concentration parameter. We argue that merger events are a major contributor to the uncertainty in halo concentration at fixed halo mass and formation time. In fact, even haloes that are typically classified as having quiescent formation histories experience multiple minor mergers. These minor mergers drive small deviations from pseudo-evolution, which cause fluctuations in the concentration parameters and result in effectively irreducible scatter in the relation between concentration and assembly history. Hence, caution should be taken when using present-day halo concentration parameter as a proxy for the halo assembly history, especially if the recent merger history is unknown. △ Less

Submitted 27 September, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: 15 pages, 7 + 1 figures, accepted for publication in MNRAS, comments welcome

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 498, Issue 3, November 2020, Pages 4450-4464

The Clustering of DESI-like Luminous Red Galaxies Using Photometric Redshifts

Authors: Rongpu Zhou, Jeffrey A. Newman, Yao-Yuan Mao, Aaron Meisner, John Moustakas, Adam D. Myers, Abhishek Prakash, Andrew R. Zentner, David Brooks, Yutong Duan, Martin Landriau, Michael E. Levi, Francisco Prada, Gregory Tarle

Abstract: We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging dataset, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The photometric LRG sample in this study contains 2.7 million objects over the redshift range of $0.4 < z < 0.9$ over 5655 deg$^2$. We have developed new p… ▽ More We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging dataset, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The photometric LRG sample in this study contains 2.7 million objects over the redshift range of $0.4 < z < 0.9$ over 5655 deg$^2$. We have developed new photometric redshift (photo-$z$) estimates using the Legacy Survey DECam and WISE photometry, with $σ_{\mathrm{NMAD}} = 0.02$ precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise ratio while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a five-parameter HOD model fit using the MultiDark halo catalog, we find that there is little evolution in HOD parameters except at the highest redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of Markov chain Monte Carlo fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogs, and present a new technique for finding best-fit parameters in this situation. Accompanying this paper we have released the Photometric Redshifts for the Legacy Surveys (PRLS) catalog of photo-$z$'s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 dataset. This catalog provides accurate photometric redshifts for objects with $z < 21$ over more than 16,000 deg$^2$ of sky. △ Less

Submitted 2 April, 2021; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: 24 pages, 25 figures, published in MNRAS

Journal ref: MNRAS 501, 3309-3331 (2021)

Non-adiabatic cosmological production of ultra-light Dark Matter

Authors: Nathan Herring, Daniel Boyanovsky, Andrew R. Zentner

Abstract: We study the non-adiabatic cosmological production of ultra light dark matter (ULDM) under a minimal set of assumptions: a free ultra light real scalar as a spectator field in its Bunch-Davies vacuum state during inflation and instantaneous reheating into a radiation dominated era. For (ULDM) fields minimally coupled to gravity, non-adiabatic particle production yields a \emph{distribution functio… ▽ More We study the non-adiabatic cosmological production of ultra light dark matter (ULDM) under a minimal set of assumptions: a free ultra light real scalar as a spectator field in its Bunch-Davies vacuum state during inflation and instantaneous reheating into a radiation dominated era. For (ULDM) fields minimally coupled to gravity, non-adiabatic particle production yields a \emph{distribution function} peaked at \emph{low} comoving momentum $\mathcal{N}_k \propto 1/k^3$. The infrared behavior is a remnant of the infrared enhancement of light minimally coupled fields during inflation. We obtain the full energy momentum tensor, show explicity its equivalence with the fluid-kinetic one in the adiabatic regime, and extract the abundance, equation of state and free streaming length (cutoff in the matter power spectrum). Taking the upper bound on the scale of inflation from Planck, the (UDLM) saturates the dark matter abundance for $m \simeq 1.5\,\times 10^{-5}\mathrm{eV}$ with an equation of state parameter $w \simeq 10^{-14}$ and a free streaming length $λ_{fs} \simeq 70\,\mathrm{pc}$. Thus this cosmologically produced (ULDM) yields a \emph{cold} dark matter particle. We argue that the abundance from non-adiabatic production yields a \emph{lower bound} on generic (ULDM) and axion-like particles that must be included in any assessment of (ULDM) as a dark matter candidate. △ Less

Submitted 30 March, 2020; v1 submitted 20 December, 2019; originally announced December 2019.

Comments: Discussion on isocurvature perturbations. To appear in PRD

Journal ref: Phys. Rev. D 101, 083516 (2020)

On the Ejection of Dark Matter from Globular Clusters

Authors: Travis J. Hurst, Andrew R. Zentner

Abstract: We investigate analytically whether in a close encounter with a star, a Dark Matter particle can be accelerated above the escape speed of a Globular Cluster and be ejected. We find that this mechanism is not sufficient to eject a massive, extended Dark Matter halo by the present time. Combined with observations of isolated Globular Clusters that may not have had their halos tidally stripped, these… ▽ More We investigate analytically whether in a close encounter with a star, a Dark Matter particle can be accelerated above the escape speed of a Globular Cluster and be ejected. We find that this mechanism is not sufficient to eject a massive, extended Dark Matter halo by the present time. Combined with observations of isolated Globular Clusters that may not have had their halos tidally stripped, these results cast doubt on the scenario in which Globular Clusters formed in Dark Matter halos. △ Less

Submitted 21 May, 2020; v1 submitted 1 October, 2019; originally announced October 2019.

Cosmological Evidence Modelling: a new simulation-based approach to constrain cosmology on non-linear scales

Authors: Johannes U. Lange, Frank C. van den Bosch, Andrew R. Zentner, Kuan Wang, Andrew P. Hearin, Hong Guo

Abstract: Extracting accurate cosmological information from galaxy-galaxy and galaxy-matter correlation functions on non-linear scales ($\lesssim 10 h^{-1} \mathrm{Mpc}$) requires cosmological simulations. Additionally, one has to marginalise over several nuisance parameters of the galaxy-halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic poste… ▽ More Extracting accurate cosmological information from galaxy-galaxy and galaxy-matter correlation functions on non-linear scales ($\lesssim 10 h^{-1} \mathrm{Mpc}$) requires cosmological simulations. Additionally, one has to marginalise over several nuisance parameters of the galaxy-halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $\mathcal{O}(10^6)$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realisations in parameter space. Afterwards, this emulator is used as a surrogate model in an MCMC analysis. Here, we demonstrate a different method called Cosmological Evidence Modelling (CEM). First, for each simulation, we calculate the Bayesian evidence marginalised over the galaxy-halo connection by repeatedly populating the simulation with galaxies. We show that this Bayesian evidence is directly related to the posterior probability of cosmological parameters. Finally, we build a physically motivated model for how the evidence depends on cosmological parameters as sampled by the simulations. We demonstrate the feasibility of CEM by using simulations from the Aemulus simulation suite and forecasting cosmological constraints from BOSS CMASS measurements of redshift-space distortions. Our analysis includes an exploration of how galaxy assembly bias affects cosmological inference. Overall, CEM has several potential advantages over the more common approach of emulating summary statistics, including the ability to easily marginalise over highly complex models of the galaxy-halo connection and greater accuracy, thereby reducing the number of simulations required. △ Less

Submitted 6 September, 2019; originally announced September 2019.

Comments: 10 pages, 2 figures, submitted to MNRAS, comments welcome

BASILISK: Bayesian Hierarchical Inference of the Galaxy-Halo Connection using Satellite Kinematics--I. Method and Validation

Authors: Frank C. van den Bosch, Johannes U. Lange, Andrew R. Zentner

Abstract: We present a Bayesian hierarchical inference formalism (Basilisk) to constrain the galaxy-halo connection using satellite kinematics. Unlike traditional methods, Basilisk does not resort to stacking the kinematics of satellite galaxies in bins of central luminosity, and does not make use of summary statistics, such as satellite velocity dispersion. Rather, Basilisk leaves the data in its raw form… ▽ More We present a Bayesian hierarchical inference formalism (Basilisk) to constrain the galaxy-halo connection using satellite kinematics. Unlike traditional methods, Basilisk does not resort to stacking the kinematics of satellite galaxies in bins of central luminosity, and does not make use of summary statistics, such as satellite velocity dispersion. Rather, Basilisk leaves the data in its raw form and computes the corresponding likelihood. In addition, Basilisk can be applied to flux-limited, rather than volume-limited samples, greatly enhancing the quantity and dynamic range of the data. And finally, Basilisk is the only available method that simultaneously solves for halo mass and orbital anisotropy of the satellite galaxies, while properly accounting for scatter in the galaxy-halo connection. Basilisk uses the conditional luminosity function to model halo occupation statistics, and assumes that satellite galaxies are a relaxed tracer population of the host halo's potential with kinematics that obey the spherical Jeans equation. We test and validate Basilisk using mocks of varying complexity, and demonstrate that it yields unbiased constraints on the galaxy-halo connection and at a precision that rivals galaxy-galaxy lensing. In particular, Basilisk accurately recovers the full PDF of the relation between halo mass and central galaxy luminosity, and simultaneously constrains the orbital anisotropy of the satellite galaxies. Basilisk's inference is not affected by potential velocity bias of the central galaxies, or by slight errors in the inferred, radial profile of satellite galaxies that arise as a consequence of interlopers and sample impurity. △ Less

Submitted 20 August, 2019; originally announced August 2019.

Comments: 30 pages, 14 figures, accepted for publication in MNRAS

How to Optimally Constrain Galaxy Assembly Bias: Supplement Projected Correlation Functions with Count-in-cells Statistics

Authors: Kuan Wang, Yao-Yuan Mao, Andrew R. Zentner, Frank C. van den Bosch, Johannes U. Lange, Chad M. Schafer, Antonia S. Villarreal, Andrew P. Hearin, Duncan Campbell

Abstract: Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real U… ▽ More Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $ΔΣ(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+ΔΣ(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $ΔΣ(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection. △ Less

Submitted 6 September, 2019; v1 submitted 22 March, 2019; originally announced March 2019.

Comments: Figures 3 and 4 show the main results. Published in Monthly Notices of the Royal Astronomical Society

Astro2020: Empirically Constraining Galaxy Evolution

Authors: Peter Behroozi, Matthew R. Becker, Frank C. van den Bosch, Jarle Brinchmann, Charlie Conroy, Mark Dickinson, Christopher M. Hirata, Andrew Hearin, Alexie Leauthaud, Chun Ly, Yao-Yuan Mao, Benjamin P. Moster, Christine O'Donnell, Casey Papovich, Aldo Rodríguez-Puebla, Rachel Somerville, Erik Tollerud, Jeremy Tinker, Yun Wang, Risa H. Wechsler, Charity Woodrum, Ann Zabludoff, Dennis Zaritsky, Andrew R. Zentner, Huanian Zhang

Abstract: Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive bl… ▽ More Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope

Authors: Alex Drlica-Wagner, Yao-Yuan Mao, Susmita Adhikari, Robert Armstrong, Arka Banerjee, Nilanjan Banik, Keith Bechtol, Simeon Bird, Kimberly K. Boddy, Ana Bonaca, Jo Bovy, Matthew R. Buckley, Esra Bulbul, Chihway Chang, George Chapline, Johann Cohen-Tanugi, Alessandro Cuoco, Francis-Yan Cyr-Racine, William A. Dawson, Ana Díaz Rivero, Cora Dvorkin, Denis Erkal, Christopher D. Fassnacht, Juan García-Bellido, Maurizio Giannotti , et al. (75 additional authors not shown)

Abstract: Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We d… ▽ More Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ . △ Less

Submitted 24 April, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 96 pages, 22 figures, 1 table

Report number: FERMILAB-PUB-19-048-A-AE

Updated Results on the Galaxy-Halo Connection from Satellite Kinematics in SDSS

Authors: Johannes U. Lange, Frank C. van den Bosch, Andrew R. Zentner, Kuan Wang, Antonia S. Villarreal

Abstract: We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter… ▽ More We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter haloes, as well as the radial profile of satellite galaxies. Confirming previous results, we find that red central galaxies live in more massive haloes than blue galaxies at fixed luminosity. Additionally, our results suggest that satellite galaxies have a radial profile less centrally concentrated than dark matter but not as cored as resolved subhaloes in dark matter-only simulations. Compared to previous works using satellite kinematics by More et al., we find much more competitive constraints on the galaxy-halo connection, on par with those derived from a combination of galaxy clustering and galaxy-galaxy lensing. We compare our results on the galaxy-halo connection to other studies using galaxy clustering and group catalogues, showing very good agreement between these different techniques. We discuss future applications of satellite kinematics in the context of constraining cosmology and the relationship between galaxies and dark matter haloes. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 18 pages, 10 figures, submitted to MNRAS, comments welcome

Maturing Satellite Kinematics into a Competitive Probe of the Galaxy-Halo Connection

Authors: Johannes U. Lange, Frank C. van den Bosch, Andrew R. Zentner, Kuan Wang, Antonia S. Villarreal

Abstract: The kinematics of satellite galaxies moving in a dark matter halo are a direct probe of the underlying gravitational potential. Thus, the phase-space distributions of satellites represent a powerful tool to determine the galaxy-halo connection from observations. By stacking the signal of a large number of satellite galaxies this potential can be unlocked even for haloes hosting a few satellites on… ▽ More The kinematics of satellite galaxies moving in a dark matter halo are a direct probe of the underlying gravitational potential. Thus, the phase-space distributions of satellites represent a powerful tool to determine the galaxy-halo connection from observations. By stacking the signal of a large number of satellite galaxies this potential can be unlocked even for haloes hosting a few satellites on average. In this work, we test the impact of various modelling assumptions on constraints derived from analysing satellite phase-space distributions in the non-linear, 1-halo regime. We discuss their potential to explain the discrepancy between average halo masses derived from satellite kinematics and gravitational lensing previously reported. Furthermore, we develop an updated, more robust analysis to extract constraints on the galaxy-halo relation from satellite properties in spectroscopic galaxy surveys such as the SDSS. We test the accuracy of this approach using a large number of realistic mock catalogues. Furthermore, we find that constraints derived from such an analysis are complementary and competitive with respect to the commonly used galaxy clustering and galaxy-galaxy lensing observables. △ Less

Submitted 24 October, 2018; originally announced October 2018.

Comments: 24 pages, 15 figures; resubmitted to MNRAS after first referee report

Particle decay in post inflationary cosmology

Authors: Nathan Herring, Brian Pardo, Daniel Boyanovsky, Andrew R. Zentner

Abstract: We study scalar particle decay during the radiation and matter dominated epochs of a standard cosmological model. An adiabatic approximation is introduced that is valid for degrees of freedom with typical wavelengths much smaller than the particle horizon ($\propto$~Hubble radius) at a given time. We implement a non-perturbative method that includes the cosmological expansion and obtain a cosmolog… ▽ More We study scalar particle decay during the radiation and matter dominated epochs of a standard cosmological model. An adiabatic approximation is introduced that is valid for degrees of freedom with typical wavelengths much smaller than the particle horizon ($\propto$~Hubble radius) at a given time. We implement a non-perturbative method that includes the cosmological expansion and obtain a cosmological Fermi's Golden Rule that enables one to compute the decay law of a parent particle of mass $m_1$, along with the build up of the population of daughter particles of mass $m_2$. The survival probability of the decaying particle is $P(t)=e^{-\widetildeΓ_k(t)\,t}$ with $\widetildeΓ_k(t)$ being an \emph{effective momentum and time dependent decay rate}. It features a transition time scale $t_{nr}$ between the relativistic and non-relativistic regimes and for $k \neq 0$ is always smaller than the analogous rate in Minkowski spacetime, as a consequence of (local) time dilation and the cosmological redshift. For $t \ll t_{nr}$ the decay law is a "stretched exponential" $P(t) = e^{-(t/t^*)^{3/2}}$, whereas for the non-relativistic stage with $t \gg t_{nr}$, we find $P(t) = e^{-Γ_0 t}\,(t/t_{nr})^{Γ_0\,t_{nr}/2}$. The Hubble time scale $\propto 1/H(t)$ introduces an energy uncertainty $ΔE \sim H(t)$ which relaxes the constraints of kinematic thresholds. This opens new decay channels into heavier particles for $2πE_k(t) H(t) \gg 4m^2_2-m^2_1$, with $E_k(t)$ the (local) comoving energy of the decaying particle. As the expansion proceeds this channel closes and the usual two particle thresholds restrict the decay kinematics. △ Less

Submitted 2 October, 2018; v1 submitted 7 August, 2018; originally announced August 2018.

Comments: more comments. published version

Journal ref: Phys. Rev. D 98, 083503 (2018)

Predictably Missing Satellites: Subhalo Abundance in Milky Way-like Halos

Authors: Catherine E. Fielder, Yao-Yuan Mao, Jeffrey A. Newman, Andrew R. Zentner, Timothy C. Licquia

Abstract: On small scales there have been a number of claims of discrepancies between the standard Cold Dark Matter (CDM) model and observations. The 'missing satellites problem' infamously describes the overabundance of subhalos from CDM simulations compared to the number of satellites observed in the Milky Way. A variety of solutions to this discrepancy have been proposed; however, the impact of the speci… ▽ More On small scales there have been a number of claims of discrepancies between the standard Cold Dark Matter (CDM) model and observations. The 'missing satellites problem' infamously describes the overabundance of subhalos from CDM simulations compared to the number of satellites observed in the Milky Way. A variety of solutions to this discrepancy have been proposed; however, the impact of the specific properties of the Milky Way halo relative to the typical halo of its mass have yet to be explored. Motivated by recent studies that identified ways in which the Milky Way is atypical (e.g., Licquia et al. 2015), we investigate how the properties of dark matter halos with mass comparable to our Galaxy's --- including concentration, spin, shape, and scale factor of the last major merger --- correlate with the subhalo abundance. Using zoom-in simulations of Milky Way-like halos, we build two models of subhalo abundance as functions of host halo properties and conclude that the Milky Way should be expected to have 22%-44% fewer subhalos with low maximum rotation velocities ($V_{\rm max}^{\rm sat} \sim 10$kms$^{-1}$) at the 95% confidence level and up to 72% fewer than average subhalos with high rotation velocities ($V_{\rm max}^{\rm sat} \gtrsim 30$kms$^{-1}$, comparable to the Magellanic Clouds) than would be expected for a typical halo of the Milky Way's mass. Concentration is the most informative single parameter for predicting subhalo abundance. Our results imply that models tuned to explain the missing satellites problem assuming typical subhalo abundances for our Galaxy will be over-correcting. △ Less

Submitted 17 June, 2019; v1 submitted 13 July, 2018; originally announced July 2018.

Comments: 21 pages, 11 figures, submitted to Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 486, Issue 4, July 2019, Pages 4545-4568

General Relativistic Corrections to the Weak Lensing Convergence Power Spectrum

Authors: John T. Giblin Jr, James B. Mertens, Glenn D. Starkman, Andrew R. Zentner

Abstract: We compute the weak lensing convergence power spectrum, $C^{κκ}(θ)$, in a dust-filled universe using fully non-linear general relativistic simulations. The spectrum is then compared to more standard, approximate calculations by computing the Bardeen (Newtonian) potentials in linearized gravity and utilizing the Born approximation. We find corrections to the angular power spectrum amplitude of orde… ▽ More We compute the weak lensing convergence power spectrum, $C^{κκ}(θ)$, in a dust-filled universe using fully non-linear general relativistic simulations. The spectrum is then compared to more standard, approximate calculations by computing the Bardeen (Newtonian) potentials in linearized gravity and utilizing the Born approximation. We find corrections to the angular power spectrum amplitude of order ten percent at very large angular scales, $\ell ~ 2-3$, and percent-level corrections at intermediate angular scales of $\ell ~ 20-30$. △ Less

Submitted 20 July, 2017; originally announced July 2017.

Comments: 9 pages, 3 figures

Journal ref: Phys. Rev. D 96, 103530 (2017)

Updated Constraints on Self-Interacting Dark Matter from Supernova 1987A

Authors: Cameron Mahoney, Adam K. Leibovich, Andrew R. Zentner

Abstract: We revisit SN1987A constraints on light, hidden sector gauge bosons ("dark photons") that are coupled to the standard model through kinetic mixing with the photon. These constraints are realized because excessive bremsstrahlung radiation of the dark photon can lead to rapid cooling of the SN1987A progenitor core, in contradiction to the observed neutrinos from that event. The models we consider ar… ▽ More We revisit SN1987A constraints on light, hidden sector gauge bosons ("dark photons") that are coupled to the standard model through kinetic mixing with the photon. These constraints are realized because excessive bremsstrahlung radiation of the dark photon can lead to rapid cooling of the SN1987A progenitor core, in contradiction to the observed neutrinos from that event. The models we consider are of interest as phenomenological models of strongly self-interacting dark matter. We clarify several possible ambiguities in the literature and identify errors in prior analyses. We find constraints on the dark photon mixing parameter that are in rough agreement with the early estimates of Dent et al., but only because significant errors in their analyses fortuitously canceled. Our constraints are in good agreement with subsequent analyses by Rrapaj & Reddy and Hardy & Lasenby. We estimate the dark photon bremsstrahlung rate using one-pion exchange (OPE), while Rrapaj & Reddy use a soft radiation approximation (SRA) to exploit measured nuclear scattering cross sections. We find that the differences between mixing parameter constraints obtained through the OPE approximation or the SRA approximation are roughly a factor of $\sim 2-3$. Hardy & Laseby include plasma effects in their calculations finding significantly weaker constraints on dark photon mixing for dark photon masses below $\sim 10\, \mathrm{MeV}$. We do not consider plasma effects. Lastly, we point out that the properties of the SN1987A progenitor core remain somewhat uncertain and that this uncertainty alone causes uncertainty of at least a factor of $\sim 2-3$ in the excluded values of the dark photon mixing parameter. Further refinement of these estimates is unwarranted until either the interior of the SN1987A progenitor is more well understood or additional, large, and heretofore neglected effects, are identified. △ Less

Submitted 5 August, 2017; v1 submitted 20 June, 2017; originally announced June 2017.

Comments: 11 pages, 4 figures; fixed calculation error in generating main figure

Journal ref: Phys. Rev. D 96, 043018 (2017)

The Galaxy Clustering Crisis in Abundance Matching

Authors: Duncan Campbell, Frank C. van den Bosch, Nikhil Padmanabhan, Yao-Yuan Mao, Andrew R. Zentner, Johannes U. Lange, Fangzhou Jiang, Antonia Villarreal

Abstract: Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxi… ▽ More Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxies produced in these models compared to those based on $V_{\rm peak}$. Despite success in predicting clustering, a simple $V_{\rm peak}$ based SHAM model results in predictions for galaxy growth that are at odds with observations. We evaluate three possible remedies that could "save" mass-based SHAM: (1) SHAM models require a significant population of "orphan" galaxies as a result of artificial disruption/merging of sub-haloes in modern high resolution dark matter simulations; (2) satellites must grow significantly after their accretion; and (3) stellar mass is significantly affected by halo assembly history. No solution is entirely satisfactory. However, regardless of the particulars, we show that popular SHAM models based on $M_{\rm peak}$ cannot be complete physical models as presented. Either $V_{\rm peak}$ truly is a better predictor of stellar mass at $z\sim 0$ and it remains to be seen how the correlation between stellar mass and $V_{\rm peak}$ comes about, or SHAM models are missing vital component(s) that significantly affect galaxy clustering. △ Less

Submitted 17 May, 2017; originally announced May 2017.

Comments: 25 pages, 22 figures, submitted to MNRAS, comments welcome

Brightest galaxies as halo centre tracers in SDSS DR7

Authors: Johannes U. Lange, Frank C. van den Bosch, Andrew Hearin, Duncan Campbell, Andrew R. Zentner, Antonia Villarreal, Yao-Yuan Mao

Abstract: Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight veloc… ▽ More Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight velocity and spatial offsets between brightest galaxies and their neighbours. We compare those to detailed mock catalogues, constructed from high-resolution, dark-matter-only $N$-body simulations, in which it is assumed that satellite galaxies trace dark matter subhaloes. This allows us to place constraints on the fraction $f_{\rm BNC}$ of haloes in which the brightest galaxy is not the central. Compared to previous studies we explicitly take into account the unrelaxed state of the host haloes, velocity offsets of halo cores and correlations between $f_{\rm BNC}$ and the satellite occupation. We find that $f_{\rm BNC}$ strongly decreases with the luminosity of the brightest galaxy and increases with the mass of the host halo. Overall, in the halo mass range $10^{13} - 10^{14.5} h^{-1} M_\odot$ we find $f_{\rm BNC} \sim 30\%$, in good agreement with a previous study by Skibba et al. We discuss the implications of these findings for studies inferring the galaxy--halo connection from satellite kinematics, models of the conditional luminosity function and galaxy formation in general. △ Less

Submitted 28 November, 2017; v1 submitted 14 May, 2017; originally announced May 2017.

Comments: 24 pages, 15 figures. Accepted for publication in MNRAS

The Immitigable Nature of Assembly Bias: The Impact of Halo Definition on Assembly Bias

Authors: Antonia S. Villarreal, Andrew R. Zentner, Yao-Yuan Mao, Chris W. Purcell, Frank C. van den Bosch, Benedikt Diemer, Johannes U. Lange, Kuan Wang, Duncan Campbell

Abstract: Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $Δ$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo m… ▽ More Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $Δ$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $Δ$. For conventional halo definitions ($Δ\sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $Δ\sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $Δ$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $Δ\gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass. △ Less

Submitted 16 September, 2017; v1 submitted 11 May, 2017; originally announced May 2017.

Comments: 19 pages, 13 figures. Updated to published version. Main result summarized in Figure 10

Journal ref: MNRAS, 472, 1088-1105 (2017)

Beyond Assembly Bias: Exploring Secondary Halo Biases for Cluster-size Haloes

Authors: Yao-Yuan Mao, Andrew R. Zentner, Risa H. Wechsler

Abstract: Secondary halo bias, commonly known as 'assembly bias,' is the dependence of halo clustering on a halo property other than mass. This prediction of the Lambda-Cold Dark Matter cosmology is essential to modelling the galaxy distribution to high precision and interpreting clustering measurements. As the name suggests, different manifestations of secondary halo bias have been thought to originate fro… ▽ More Secondary halo bias, commonly known as 'assembly bias,' is the dependence of halo clustering on a halo property other than mass. This prediction of the Lambda-Cold Dark Matter cosmology is essential to modelling the galaxy distribution to high precision and interpreting clustering measurements. As the name suggests, different manifestations of secondary halo bias have been thought to originate from halo assembly histories. We show conclusively that this is incorrect for cluster-size haloes. We present an up-to-date summary of secondary halo biases of high-mass haloes due to various halo properties including concentration, spin, several proxies of assembly history, and subhalo properties. While concentration, spin, and the abundance and radial distribution of subhaloes exhibit significant secondary biases, properties that directly quantify halo assembly history do not. In fact, the entire assembly histories of haloes in pairs are nearly identical to those of isolated haloes. In general, a global correlation between two halo properties does not predict whether or not these two properties exhibit similar secondary biases. For example, assembly history and concentration (or subhalo abundance) are correlated for both paired and isolated haloes, but follow slightly different conditional distributions in these two cases. This results in a secondary halo bias due to concentration (or subhalo abundance), despite the lack of assembly bias in the strict sense for cluster-size haloes. Due to this complexity, caution must be exercised in using any one halo property as a proxy to study the secondary bias due to another property. △ Less

Submitted 28 August, 2018; v1 submitted 10 May, 2017; originally announced May 2017.

Comments: 14+1 pages, 9+1 figures. Figures 2 and 8 highlight the main results. Accepted by MNRAS. Fig 9 updated; main result unchanged

Journal ref: MNRAS, 474, 5143 (2018)

The Radial Acceleration Relation in Disk Galaxies in the MassiveBlack-II Simulation

Authors: Ananth Tenneti, Yao-Yuan Mao, Rupert A. C. Croft, Tiziana Di Matteo, Arthur Kosowsky, Fernando Zago, Andrew R. Zentner

Abstract: A strong correlation has been measured between the observed centripetal accelerations in galaxies and the accelerations implied by the baryonic components of galaxies. This empirical radial acceleration relation must be accounted for in any viable model of galaxy formation. We measure and compare the radial accelerations contributed by baryons and by dark matter in disk galaxies in the MassiveBlac… ▽ More A strong correlation has been measured between the observed centripetal accelerations in galaxies and the accelerations implied by the baryonic components of galaxies. This empirical radial acceleration relation must be accounted for in any viable model of galaxy formation. We measure and compare the radial accelerations contributed by baryons and by dark matter in disk galaxies in the MassiveBlack-II hydrodynamic galaxy formation simulation. The sample of 1594 galaxies spans three orders of magnitude in luminosity and four in surface brightness, comparable to the observed sample from the Spitzer Photometry & Accurate Rotation Curves (SPARC) dataset used by McGaugh et al. (2016). We find that radial accelerations contributed by baryonic matter only and by total matter are highly correlated, with only small scatter around their mean or median relation, despite the wide ranges of galaxy luminosity and surface brightness. We further find that the radial acceleration relation in this simulation differs from that of the SPARC sample, and can be described by a simple power law in the acceleration range we are probing. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Comments: 7 pages, 4 figures

Journal ref: MNRAS, 474, 3125 (2018)

Constraints on Assembly Bias from Galaxy Clustering

Authors: Andrew R. Zentner, Andrew Hearin, Frank C. van den Bosch, Johannes U. Lange, Antonia Villarreal

Abstract: We constrain the newly-introduced decorated Halo Occupation Distribution (HOD) model using SDSS DR7 measurements of projected galaxy clustering or r-band luminosity threshold samples. The decorated HOD is a model for the galaxy-halo connection that augments the HOD by allowing for the possibility of galaxy assembly bias: galaxy luminosity may be correlated with dark matter halo properties besides… ▽ More We constrain the newly-introduced decorated Halo Occupation Distribution (HOD) model using SDSS DR7 measurements of projected galaxy clustering or r-band luminosity threshold samples. The decorated HOD is a model for the galaxy-halo connection that augments the HOD by allowing for the possibility of galaxy assembly bias: galaxy luminosity may be correlated with dark matter halo properties besides mass, Mvir. We demonstrate that it is not possible to rule out galaxy assembly bias using DR7 measurements of galaxy clustering alone. Moreover, galaxy samples with Mr < -20 and Mr < -20.5 favor strong central galaxy assembly bias. These samples prefer scenarios in which high-concentration are more likely to host a central galaxy relative to low-concentration halos of the same mass. We exclude zero assembly bias with high significance for these samples. Satellite galaxy assembly bias is significant for the faintest sample, Mr < -19. We find no evidence for assembly bias in the Mr < -21 sample. Assembly bias should be accounted for in galaxy clustering analyses or attempts to exploit galaxy clustering to constrain cosmology. In addition to presenting the first constraints on HOD models that accommodate assembly bias, our analysis includes several improvements over previous analyses of these data. Therefore, our inferences supersede previously-published results even in the case of a standard HOD analysis. △ Less

Submitted 24 June, 2016; originally announced June 2016.

Comments: 15 pages, 8 figures. To be submitted to MNRAS. Comments Welcome. Python scripts to perform this analysis and MCMC chains will all be made publicly available

Introducing Decorated HODs: modeling assembly bias in the galaxy-halo connection

Authors: Andrew P. Hearin, Andrew R. Zentner, Frank C. van den Bosch, Duncan Campbell, Erik Tollerud

Abstract: The connection between galaxies and dark matter halos is often inferred from data using probabilistic models, such as the Halo Occupation Distribution (HOD). Conventional HOD formulations assume that only halo mass governs the galaxy-halo connection. Violations of this assumption, known as galaxy assembly bias, threaten the HOD program. We introduce decorated HODs, a new, flexible class of models… ▽ More The connection between galaxies and dark matter halos is often inferred from data using probabilistic models, such as the Halo Occupation Distribution (HOD). Conventional HOD formulations assume that only halo mass governs the galaxy-halo connection. Violations of this assumption, known as galaxy assembly bias, threaten the HOD program. We introduce decorated HODs, a new, flexible class of models designed to account for assembly bias. Decorated HODs minimally expand the parameter space and maximize the independence between traditional and novel HOD parameters. We use decorated HODs to quantify the influence of assembly bias on clustering and lensing statistics. For SDSS-like samples, the impact of assembly bias on galaxy clustering can be as large as a factor of two on r ~ 200 kpc scales and ~15% in the linear regime. Assembly bias can either enhance or diminish clustering on large scales, but generally increases clustering on scales r <~ 1 Mpc. We performed our calculations with Halotools, an open-source, community-driven python package for studying the galaxy-halo connection (http://halotools.readthedocs.org). We conclude by describing the use of decorated HODs to treat assembly bias in otherwise conventional likelihood analyses. △ Less

Submitted 9 December, 2015; originally announced December 2015.

Comments: Figure 1 provides a cartoon visualization of the analytical formalism. Figures 4 & 5 show the characteristic signature that assembly bias imprints on the clustering and lensing of SDSS-like galaxy samples. 21 pages, 8 figures, submitted to MNRAS

Incorporating Astrophysical Systematics into a Generalized Likelihood for Cosmology with Type Ia Supernovae

Authors: Kara A. Ponder, W. Michael Wood-Vasey, Andrew R. Zentner

Abstract: Traditional cosmological inference using Type Ia supernovae (SNeIa) have used stretch- and color-corrected fits of SN Ia light curves and assumed a resulting fiducial mean and symmetric intrinsic dispersion for the resulting relative luminosity. As systematics become the main contributors to the error budget, it has become imperative to expand supernova cosmology analyses to include a more general… ▽ More Traditional cosmological inference using Type Ia supernovae (SNeIa) have used stretch- and color-corrected fits of SN Ia light curves and assumed a resulting fiducial mean and symmetric intrinsic dispersion for the resulting relative luminosity. As systematics become the main contributors to the error budget, it has become imperative to expand supernova cosmology analyses to include a more general likelihood to model systematics to remove biases with losses in precision. To illustrate an example likelihood analysis, we use a simple model of two populations with a relative luminosity shift, independent intrinsic dispersions, and linear redshift evolution of the relative fraction of each population. Treating observationally viable two-population mock data using a one-population model results in an inferred dark energy equation of state parameter $w$ that is biased by roughly 2 times its statistical error for a sample of N $ \gtrsim$ 2500 SNeIa. Modeling the two-population data with a two-population model removes this bias at a cost of an approximately $\sim20\%$ increase in the statistical constraint on $w$. These significant biases can be realized even if the support for two underlying SNeIa populations, in the form of model selection criteria, is inconclusive. With the current observationally-estimated difference in the two proposed populations, a sample of N $ \gtrsim$ 10,000 SNeIa is necessary to yield conclusive evidence of two populations. △ Less

Submitted 14 October, 2016; v1 submitted 14 November, 2015; originally announced November 2015.

Comments: 14 pages, 11 figures, 3 tables. Updated to include acknowledgements

Journal ref: The Astrophysical Journal, Volume 825, Issue 1, article id. 35, 13 pp. (2016)

Mass assembly history and infall time of the Fornax dwarf spheroidal galaxy

Authors: Mei-Yu Wang, Louis E. Strigari, Mark R. Lovell, Carlos S. Frenk, Andrew R. Zentner

Abstract: We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold dark matter (CDM), warm dark matter (WDM), and decaying dark matter (DDM) simulations for our models of structure formation. The subhalo candidates in CDM typically have smaller mass and higher concentrations at z =… ▽ More We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold dark matter (CDM), warm dark matter (WDM), and decaying dark matter (DDM) simulations for our models of structure formation. The subhalo candidates in CDM typically have smaller mass and higher concentrations at z = 0 than the corresponding candidates in WDM and DDM. We examine the formation histories of the ~ 100 Fornax candidate subhalos identified in CDM simulations and, using approximate luminosity-mass relationships for subhalos, we find two of these subhalos that are consistent with both the Fornax luminosity and kinematics. These two subhalos have a peak mass over ten times larger than their z = 0 mass. We suggest that in CDM the dark matter halo hosting Fornax must have been severely stripped of mass and that it had an infall time into the Milky Way of ~ 9 Gyr ago. In WDM, we find that candidate subhalos consistent with the properties of Fornax have a similar infall time and a similar degree of mass loss, while in DDM we find a later infall time of ~ 3 - 4 Gyr ago and significantly less mass loss. We discuss these results in the context of the Fornax star formation history, and show that these predicted subhalo infall times can be linked to different star formation quenching mechanisms. This emphasizes the links between the properties of the dark matter and the mechanisms that drive galaxy evolution. △ Less

Submitted 22 September, 2015; v1 submitted 14 September, 2015; originally announced September 2015.

Comments: 14 pages, 12 figures, some typos are corrected, submitted to MNRAS

Connecting massive galaxies to dark matter halos in BOSS - I. Is galaxy color a stochastic process in high-mass halos?

Authors: Shun Saito, Alexie Leauthaud, Andrew P. Hearin, Kevin Bundy, Andrew R. Zentner, Peter S. Behroozi, Beth A. Reid, Manodeep Sinha, Jean Coupon, Jeremy L. Tinker, Martin White, Donald P. Schneider

Abstract: We use subhalo abundance matching (SHAM) to model the stellar mass function (SMF) and clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) "CMASS" sample at $z\sim0.5$. We introduce a novel method which accounts for the stellar mass incompleteness of CMASS as a function of redshift, and produce CMASS mock catalogs which include selection effects, reproduce the overall SMF, the projecte… ▽ More We use subhalo abundance matching (SHAM) to model the stellar mass function (SMF) and clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) "CMASS" sample at $z\sim0.5$. We introduce a novel method which accounts for the stellar mass incompleteness of CMASS as a function of redshift, and produce CMASS mock catalogs which include selection effects, reproduce the overall SMF, the projected two-point correlation function $w_{\rm p}$, the CMASS $dn/dz$, and are made publicly available. We study the effects of assembly bias above collapse mass in the context of "age matching" and show that these effects are markedly different compared to the ones explored by Hearin et al. (2013) at lower stellar masses. We construct two models, one in which galaxy color is stochastic ("AbM" model) as well as a model which contains assembly bias effects ("AgM" model). By confronting the redshift dependent clustering of CMASS with the predictions from our model, we argue that that galaxy colors are not a stochastic process in high-mass halos. Our results suggest that the colors of galaxies in high-mass halos are determined by other halo properties besides halo peak velocity and that assembly bias effects play an important role in determining the clustering properties of this sample. △ Less

Submitted 24 May, 2016; v1 submitted 1 September, 2015; originally announced September 2015.

Comments: 22 pages. Appendix. B added. Matches the version accepted by MNRAS. Mock galaxy catalog and HOD table are available at http://www.massivegalaxies.com

Report number: IPMU15-0129

Indirect Probes of Dark Matter and Globular Cluster Properties From Dark Matter Annihilation within the Coolest White Dwarfs

Authors: Travis J. Hurst, Andrew R. Zentner, Aravind Natarajan, Carles Badenes

Abstract: White Dwarfs (WD) capture Dark Matter (DM) as they orbit within their host halos. These captured particles may subsequently annihilate, heating the stellar core and preventing the WD from cooling. The potential wells of WDs are considerably deeper and core temperatures significantly cooler than those of main sequence stars. Consequently, DM evaporation is less important in WDs and DM with masses… ▽ More White Dwarfs (WD) capture Dark Matter (DM) as they orbit within their host halos. These captured particles may subsequently annihilate, heating the stellar core and preventing the WD from cooling. The potential wells of WDs are considerably deeper and core temperatures significantly cooler than those of main sequence stars. Consequently, DM evaporation is less important in WDs and DM with masses $M_χ \gtrsim 100\, \kev$ and annihilation cross-sections orders of magnitude below the canonical thermal cross-section ($\sigmav \gtrsim 10^{-46}\, \cm^3$/s) can significantly alter WD cooling in particular astrophysical environments. We consider WDs in globular clusters (GCs) and dwarf galaxies. If the parameters of the DM particle are known, then the temperature of the coolest WD in a GC can be used to constrain the DM density of the cluster's halo (potentially even ruling out the presence of a halo if the inferred density is of order the ambient Galactic density). Recently several direct detection experiments have seen signals whose origins might be due to low mass DM. In this paper, we show that if these claims from CRESST, DAMA, CDMS-Si, and CoGeNT could be interpreted as DM, then observations of NGC 6397 limit the fraction of DM in that cluster to be $f_{\mathrm{DM}} \lesssim 10^{-3}$. This would be an improvement over existing constraints of 3 orders of magnitude and clearly rule out a significant DM halo for this cluster... △ Less

Submitted 19 May, 2015; v1 submitted 15 October, 2014; originally announced October 2014.

Comments: 13 pages, 7 figures

Journal ref: Phys. Rev. D 91, 103514 (2015)

Cosmological Simulations of Decaying Dark Matter: Implications for Small-scale Structure of Dark Matter Halos

Authors: Mei-Yu Wang, Annika H. G. Peter, Louis E. Strigari, Andrew R. Zentner, Bryan Arant, Shea Garrison-Kimmel, Miguel Rocha

Abstract: We present a set of N-body simulations of a class of models in which an unstable dark matter particle decays into a stable non-interacting dark matter particle, with decay lifetime comparable to the Hubble time. We study the effects of the kinematic recoil velocity received by the stable dark matter on the structures of dark matter halos ranging from galaxy-cluster to Milky Way mass scales. For Mi… ▽ More We present a set of N-body simulations of a class of models in which an unstable dark matter particle decays into a stable non-interacting dark matter particle, with decay lifetime comparable to the Hubble time. We study the effects of the kinematic recoil velocity received by the stable dark matter on the structures of dark matter halos ranging from galaxy-cluster to Milky Way mass scales. For Milky Way-mass halos, we use high-resolution, zoom-in simulations to explore the effects of decays on Galactic substructure. In general, halos with circular velocities comparable to the magnitude of kick velocity are most strongly affected by decays. We show that decaying dark matter models with lifetimes comparable to Hubble time and recoil speeds about 20-40 km/s can significantly reduce both the abundance of Galactic subhalos and the internal densities of the subhalos. We also compare subhalo circular velocity profiles with observational constraints on the Milky Way dwarf satellite galaxies. Interestingly, we find that decaying dark matter models that do not violate current astrophysical constraints, can significantly mitigate both the well-documented "missing satellites problem" and the more recent "too big to fail problem" associated with the abundances and densities of Local Group dwarf satellite galaxies. A relatively unique feature of late decaying dark matter models is that they predict significant evolution of halos as a function of time. This is an important consideration because at high redshifts, prior to decays, decaying models exhibit the same sequence of structure formation as cold dark matter. We conclude that models of decaying dark matter make predictions that are relevant for the interpretation of observations of small galaxies in the Local Group and can be tested or constrained by the kinematics of Local Group dwarf galaxies as well as by forthcoming large-scale surveys. △ Less

Submitted 12 June, 2014; v1 submitted 2 June, 2014; originally announced June 2014.

Comments: 17 pages, 14 figures, references added, labels added in figure 1 & 2, submitted to MNRAS

Systematic errors in the measurement of neutrino masses due to baryonic feedback processes: Prospects for stage IV lensing surveys

Authors: Aravind Natarajan, Andrew R. Zentner, Nicholas Battaglia, Hy Trac

Abstract: We examine the importance of baryonic feedback effects on the matter power spectrum on small scales, and the implications for the precise measurement of neutrino masses through gravitational weak lensing. Planned large galaxy surveys such as the Large Synoptic Sky Telescope (LSST) and Euclid are expected to measure the sum of neutrino masses to extremely high precision, sufficient to detect non-ze… ▽ More We examine the importance of baryonic feedback effects on the matter power spectrum on small scales, and the implications for the precise measurement of neutrino masses through gravitational weak lensing. Planned large galaxy surveys such as the Large Synoptic Sky Telescope (LSST) and Euclid are expected to measure the sum of neutrino masses to extremely high precision, sufficient to detect non-zero neutrino masses even in the minimal mass normal hierarchy. We show that weak lensing of galaxies while being a very good probe of neutrino masses, is extremely sensitive to baryonic feedback processes. We use publicly available results from the Overwhelmingly Large Simulations (OWLS) project to investigate the effects of active galactic nuclei feedback, the nature of the stellar initial mass function, and gas cooling rates, on the measured weak lensing shear power spectrum. Using the Fisher matrix formalism and priors from CMB+BAO data, we show that when one does not account for feedback, the measured neutrino mass may be substantially larger or smaller than the true mass, depending on the dominant feedback mechanism, with the mass error |Δm_nu| often exceeding the mass m_nu itself. We also consider gravitational lensing of the cosmic microwave background (CMB) and show that it is not sensitive to baryonic feedback on scales l < 2000, although CMB experiments that aim for sensitivities sigma(m_nu) < 0.02 eV will need to include baryonic effects in modeling the CMB lensing potential. A combination of CMB lensing and galaxy lensing can help break the degeneracy between neutrino masses and baryonic feedback processes. We conclude that future large galaxy lensing surveys such as LSST and Euclid can only measure neutrino masses accurately if the matter power spectrum can be measured to similar accuracy. △ Less

Submitted 3 September, 2014; v1 submitted 23 May, 2014; originally announced May 2014.

Comments: No major changes. References added, and text expanded. Published in Phys. Rev. D, and replaced to reflect the published version

Journal ref: Phys. Rev. D 90, 063516 (2014)

Predicting Galaxy Star Formation Rates via the Co-evolution of Galaxies and Halos

Authors: Douglas F. Watson, Andrew P. Hearin, Andreas A. Berlind, Matthew R. Becker, Peter S. Behroozi, Ramin A. Skibba, Reinabelle Reyes, Andrew R. Zentner, Frank C. van den Bosch

Abstract: In this paper, we test the age matching hypothesis that the star formation rate (SFR) of a galaxy of fixed stellar mass is determined by its dark matter halo formation history, and as such, that more quiescent galaxies reside in older halos. This simple model has been remarkably successful at predicting color-based galaxy statistics at low redshift as measured in the Sloan Digital Sky Survey (SDSS… ▽ More In this paper, we test the age matching hypothesis that the star formation rate (SFR) of a galaxy of fixed stellar mass is determined by its dark matter halo formation history, and as such, that more quiescent galaxies reside in older halos. This simple model has been remarkably successful at predicting color-based galaxy statistics at low redshift as measured in the Sloan Digital Sky Survey (SDSS). To further test this method with observations, we present new SDSS measurements of the galaxy two-point correlation function and galaxy-galaxy lensing as a function of stellar mass and SFR, separated into quenched and star-forming galaxy samples. We find that our age matching model is in excellent agreement with these new measurements. We also employ a galaxy group finder and show that our model is able to predict: (1) the relative SFRs of central and satellite galaxies, (2) the SFR-dependence of the radial distribution of satellite galaxy populations within galaxy groups, rich groups, and clusters and their surrounding larger scale environments, and (3) the interesting feature that the satellite quenched fraction as a function of projected radial distance from the central galaxy exhibits an ~ r^-.15 slope, independent of environment. The accurate prediction for the spatial distribution of satellites is intriguing given the fact that we do not explicitly model satellite-specific processes after infall, and that in our model the virial radius does not mark a special transition region in the evolution of a satellite, contrary to most galaxy evolution models. The success of the model suggests that present-day galaxy SFR is strongly correlated with halo mass assembly history. △ Less

Submitted 6 October, 2014; v1 submitted 6 March, 2014; originally announced March 2014.

Comments: 14 pages, 6 figures, and one Appendix of tabulated data of new SFR-dependent galaxy clustering and galaxy-galaxy lensing measurements from SDSS. Accepted to MNRAS

Galaxy Assembly Bias: A Significant Source of Systematic Error in the Galaxy-Halo Relationship

Authors: Andrew R. Zentner, Andrew P. Hearin, Frank C. van den Bosch

Abstract: It is common practice for methods that use galaxy clustering to constrain the galaxy-halo relationship, such as the halo occupation distribution (HOD) and/or conditional luminosity function (CLF), to assume that halo mass alone suffices to determine a halo's resident galaxy population. Yet, the clustering strength of cold dark matter halos depends upon halo properties in addition to mass, such as… ▽ More It is common practice for methods that use galaxy clustering to constrain the galaxy-halo relationship, such as the halo occupation distribution (HOD) and/or conditional luminosity function (CLF), to assume that halo mass alone suffices to determine a halo's resident galaxy population. Yet, the clustering strength of cold dark matter halos depends upon halo properties in addition to mass, such as formation time, an effect referred to as assembly bias. If galaxy characteristics are correlated with any of these auxiliary halo properties, the basic assumption of HOD/CLF methods is violated. We estimate the potential for assembly bias to induce systematic errors in inferred halo occupation statistics. We use halo abundance matching and age matching to construct fiducial mock galaxy catalogs that exhibit assembly bias as well as additional mock catalogs with identical HODs, but with assembly bias removed. We fit a parameterized HOD to the projected two-point clustering of mock galaxies in each catalog to assess the systematic errors induced by reasonable levels of assembly bias. In the absence of assembly bias, the inferred HODs generally describe the true underlying HODs well, validating the basic methodology. However, in all of the cases with assembly bias, the inferred HODs have systematic errors that are statistically significant. In most cases, these systematic errors cannot be identified using void statistics as auxiliary observables. We conclude that the galaxy-halo relationship inferred from galaxy clustering should be subject to a non-negligible systematic error induced by assembly bias. Our work suggests that efforts to model and/or constrain assembly bias should be high priorities as it is a threatening source of systematic error in galaxy evolution studies as well as the precision cosmology program. △ Less

Submitted 10 February, 2015; v1 submitted 7 November, 2013; originally announced November 2013.

Comments: 28 pages, 15 figures including an appendix. v2 includes minor revisions based upon referee comments. Accepted for publication in the MNRAS. Title changed from original "Phantom Menace of Galaxy Clustering"

Journal ref: MNRAS 443, 3044-3067 (2014)

The Dark Side of Galaxy Color: evidence from new SDSS measurements of galaxy clustering and lensing

Authors: Andrew P. Hearin, Douglas F. Watson, Matthew R. Becker, Reinabelle Reyes, Andreas A. Berlind, Andrew R. Zentner

Abstract: The age matching model has recently been shown to predict correctly the luminosity L and g-r color of galaxies residing within dark matter halos. The central tenet of the model is intuitive: older halos tend to host galaxies with older stellar populations. In this paper, we demonstrate that age matching also correctly predicts the g-r color trends exhibited in a wide variety of statistics of the g… ▽ More The age matching model has recently been shown to predict correctly the luminosity L and g-r color of galaxies residing within dark matter halos. The central tenet of the model is intuitive: older halos tend to host galaxies with older stellar populations. In this paper, we demonstrate that age matching also correctly predicts the g-r color trends exhibited in a wide variety of statistics of the galaxy distribution for stellar mass M* threshold samples. In particular, we present new measurements of the galaxy two-point correlation function and the galaxy-galaxy lensing signal as a function of M* and g-r color from the Sloan Digital Sky Survey, and show that age matching exhibits remarkable agreement with these and other statistics of low-redshift galaxies. In so doing, we also demonstrate good agreement between the galaxy-galaxy lensing observed by SDSS and the signal predicted by abundance matching, a new success of this model. We describe how age matching is a specific example of a larger class of Conditional Abundance Matching models (CAM), a theoretical framework we introduce here for the first time. CAM provides a general formalism to study correlations at fixed mass between any galaxy property and any halo property. The striking success of our simple implementation of CAM provides compelling evidence that this technique has the potential to describe the same set of data as alternative models, but with a dramatic reduction in the required number of parameters. CAM achieves this reduction by exploiting the capability of contemporary N-body simulations to determine dark matter halo properties other than mass alone, which distinguishes our model from conventional approaches to the galaxy-halo connection. △ Less

Submitted 11 November, 2013; v1 submitted 24 October, 2013; originally announced October 2013.

Comments: references added, minor adjustments to text and notation

Lyman-alpha Forest Constraints on Decaying Dark Matter

Authors: Mei-Yu Wang, Rupert A. C. Croft, Annika H. G. Peter, Andrew R. Zentner, Chris W. Purcell

Abstract: We present an analysis of high-resolution N-body simulations of decaying dark matter cosmologies focusing on the statistical properties of the transmitted Lyman-alpha forest flux in the high-redshift intergalactic medium. In this type of model a dark matter particle decays into a slightly less massive stable dark matter daughter particle and a comparably light particle. The small mass splitting pr… ▽ More We present an analysis of high-resolution N-body simulations of decaying dark matter cosmologies focusing on the statistical properties of the transmitted Lyman-alpha forest flux in the high-redshift intergalactic medium. In this type of model a dark matter particle decays into a slightly less massive stable dark matter daughter particle and a comparably light particle. The small mass splitting provides a non-relativistic kick velocity V_k to the daughter particle resulting in free-streaming and subsequent damping of small-scale density fluctuations. Current Lyman-alpha forest power spectrum measurements probe comoving scales up to ~ 2-3 h^-1 Mpc at redshifts z ~ 2-4, providing one of the most robust ways to probe cosmological density fluctuations on relatively small scales. The suppression of structure growth due to the free-streaming of dark matter daughter particles also has a significant impact on the neutral hydrogen cloud distribution, which traces the underlying dark matter distribution well at high redshift. We exploit Lyman-alpha forest power spectrum measurements to constrain the amount of free-streaming of dark matter in such models and thereby place limits on decaying dark matter based only on the dynamics of cosmological perturbations without any assumptions about the interactions of the decay products. We find that SDSS 1D Lyman-alpha forest power spectrum data place a lifetime-dependent upper limit V_k < 30-70 km/s for decay lifetimes < 10 Gyr. This is the most stringent model-independent bound on invisible dark matter decays with small mass splittings. For large mass splittings (large V_k), Lyman-alpha forest data restrict the dark matter lifetime to Gamma^-1 > 40 Gyr. Forthcoming BOSS data should be able to provide more stringent constraints on exotic dark matter, mainly because the larger BOSS quasar spectrum sample will significantly reduce statistical errors. △ Less

Submitted 27 September, 2013; originally announced September 2013.

Comments: 12 pages, 4 figures, submitted to PRD

Groups of two galaxies in SDSS: implications of colours on star formation quenching time-scales

Authors: Christopher Q. Trinh, Elizabeth J. Barton, James S. Bullock, Michael C. Cooper, Andrew R. Zentner, Risa H. Wechsler

Abstract: We have devised a method to select galaxies that are isolated in their dark matter halo (N=1 systems) and galaxies that reside in a group of exactly two (N=2 systems). Our N=2 systems are widely-separated (up to $\sim$\,200\,$h^{-1}$\,kpc), where close galaxy-galaxy interactions are not dominant. We apply our selection criteria to two volume-limited samples of galaxies from SDSS DR6 with… ▽ More We have devised a method to select galaxies that are isolated in their dark matter halo (N=1 systems) and galaxies that reside in a group of exactly two (N=2 systems). Our N=2 systems are widely-separated (up to $\sim$\,200\,$h^{-1}$\,kpc), where close galaxy-galaxy interactions are not dominant. We apply our selection criteria to two volume-limited samples of galaxies from SDSS DR6 with $M_{r}-5 \log_{10} h \leq$ -19 and -20 to study the effects of the environment of very sparse groups on galaxy colour. For satellite galaxies in a group of two, we find a red excess attributed to star formation quenching of 0.15\,$\pm$\,0.01 and 0.14\,$\pm$\,0.01 for the -19 and -20 samples, respectively, relative to isolated galaxies of the same stellar mass. Assuming N=1 systems are the progenitors of N=2 systems, an immediate-rapid star formation quenching scenario is inconsistent with these observations. A delayed-then-rapid star formation quenching scenario with a delay time of 3.3 and 3.7\,Gyr for the -19 and -20 samples, respectively, yields a red excess prediction in agreement with the observations. The observations also reveal that central galaxies in a group of two have a slight blue excess of 0.06\,$\pm$\,0.02 and 0.02\,$\pm$\,0.01 for the -19 and -20 samples, respectively, relative to N=1 populations of the same stellar mass. Our results demonstrate that even the environment of very sparse groups of luminous galaxies influence galaxy evolution and in-depth studies of these simple systems are an essential step towards understanding galaxy evolution in general. △ Less

Submitted 28 September, 2013; v1 submitted 24 January, 2013; originally announced January 2013.

Comments: 17 pages, 11 figures, accepted to MNRAS

Accounting for Baryons in Cosmological Constraints from Cosmic Shear

Authors: Andrew R. Zentner, Elisabetta Semboloni, Scott Dodelson, Tim Eifler, Elisabeth Krause, Andrew P. Hearin

Abstract: One of the most pernicious theoretical systematics facing upcoming gravitational lensing surveys is the uncertainty introduced by the effects of baryons on the power spectrum of the convergence field. One method that has been proposed to account for these effects is to allow several additional parameters (that characterize dark matter halos) to vary and to fit lensing data to these halo parameters… ▽ More One of the most pernicious theoretical systematics facing upcoming gravitational lensing surveys is the uncertainty introduced by the effects of baryons on the power spectrum of the convergence field. One method that has been proposed to account for these effects is to allow several additional parameters (that characterize dark matter halos) to vary and to fit lensing data to these halo parameters concurrently with the standard set of cosmological parameters. We test this method. In particular, we use this technique to model convergence power spectrum predictions from a set of cosmological simulations. We estimate biases in dark energy equation of state parameters that would be incurred if one were to fit the spectra predicted by the simulations either with no model for baryons, or with the proposed method. We show that neglecting baryonic effect leads to biases in dark energy parameters that are several times the statistical errors for a survey like the Dark Energy Survey. The proposed method to correct for baryonic effects renders the residual biases in dark energy equation of state parameters smaller than the statistical errors. These results suggest that this mitigation method may be applied to analyze convergence spectra from a survey like the Dark Energy Survey. For significantly larger surveys, such as will be carried out by the Large Synoptic Survey Telescope, the biases introduced by baryonic effects are much more significant. We show that this mitigation technique significantly reduces the biases for such larger surveys, but that a more effective mitigation strategy will need to be developed in order ensure that the residual biases in these surveys fall below the statistical errors. △ Less

Submitted 5 December, 2012; originally announced December 2012.

Comments: 16 pages, 9 figures. Submitted to Physical Review D. Comments Welcome

SHAM Beyond Clustering: New Tests of Galaxy-Halo Abundance Matching with Galaxy Groups

Authors: Andrew P. Hearin, Andrew R. Zentner, Andreas A. Berlind, Jeffrey A. Newman

Abstract: We construct mock catalogs of galaxy groups using subhalo abundance matching (SHAM) and undertake several new tests of the SHAM prescription for the galaxy-dark matter connection. All SHAM models we studied exhibit significant tension with galaxy groups observed in the Sloan Digital Sky Survey (SDSS). The SHAM prediction for the field galaxy luminosity function (LF) is systematically too dim, and… ▽ More We construct mock catalogs of galaxy groups using subhalo abundance matching (SHAM) and undertake several new tests of the SHAM prescription for the galaxy-dark matter connection. All SHAM models we studied exhibit significant tension with galaxy groups observed in the Sloan Digital Sky Survey (SDSS). The SHAM prediction for the field galaxy luminosity function (LF) is systematically too dim, and the group galaxy LF systematically too bright, regardless of the details of the SHAM prescription. SHAM models connecting r-band luminosity, Mr, to Vacc, the maximum circular velocity of a subhalo at the time of accretion onto the host, faithfully reproduce galaxy group abundance as a function of richness, g(N). However, SHAM models connecting Mr with Vpeak, the peak value of Vmax over the entire merger history of the halo, over-predict galaxy group abundance. Our results suggest that no SHAM model can simultaneously reproduce the observed g(N) and two-point projected galaxy clustering. Nevertheless, we also report a new success of SHAM: an accurate prediction for Phi(m12), the abundance of galaxy groups as a function of magnitude gap m12, defined as the difference between the r-band absolute magnitude of the two brightest group members. We show that it may be possible to use joint measurements of g(N) and Phi(m12) to tightly constrain the details of the SHAM implementation. Additionally, we show that the hypothesis that the luminosity gap is constructed via random draws from a universal LF provides a poor description of the data, contradicting recent claims in the literature. Finally, we test a common assumption of the Conditional Luminosity Function (CLF) formalism, that the satellite LF need only be conditioned by the brightness of the central galaxy. We find this assumption to be well-supported by the observed Phi(m12). △ Less

Submitted 22 July, 2013; v1 submitted 17 October, 2012; originally announced October 2012.

Comments: minor changes to draft; conclusions unchanged; version accepted by MNRAS

Report number: FERMILAB-PUB-12-531-A

Bailing Out the Milky Way: Variation in the Properties of Massive Dwarfs Among Galaxy-Sized Systems

Authors: Chris W. Purcell, Andrew R. Zentner

Abstract: Recent kinematical constraints on the internal densities of the Milky Way's dwarf satellites have revealed a discrepancy with the subhalo populations of simulated Galaxy-scale halos in the standard CDM model of hierarchical structure formation. This has been dubbed the "too big to fail" problem, with reference to the improbability of large and invisible companions existing in the Galactic environm… ▽ More Recent kinematical constraints on the internal densities of the Milky Way's dwarf satellites have revealed a discrepancy with the subhalo populations of simulated Galaxy-scale halos in the standard CDM model of hierarchical structure formation. This has been dubbed the "too big to fail" problem, with reference to the improbability of large and invisible companions existing in the Galactic environment. In this paper, we argue that both the Milky Way observations and simulated subhalos are consistent with the predictions of the standard model for structure formation. Specifically, we show that there is significant variation in the properties of subhalos among distinct host halos of fixed mass and suggest that this can reasonably account for the deficit of dense satellites in the Milky Way. We exploit well-tested analytic techniques to predict the properties in a large sample of distinct host halos with a variety of masses spanning the range expected of the Galactic halo. The analytic model produces subhalo populations consistent with both Via Lactea II and Aquarius, and our results suggest that natural variation in subhalo properties suffices to explain the discrepancy between Milky Way satellite kinematics and these numerical simulations. At least ~10% of Milky Way-sized halos host subhalo populations for which there is no "too big to fail" problem, even when the host halo mass is as large as M_host = 10^12.2 h^-1 M_sun. Follow-up studies consisting of high-resolution simulations of a large number of Milky Way-sized hosts are necessary to confirm our predictions. In the absence of such efforts, the "too big to fail" problem does not appear to be a significant challenge to the standard model of hierarchical formation. [abridged] △ Less

Submitted 6 December, 2012; v1 submitted 22 August, 2012; originally announced August 2012.

Comments: 12 pages, 3 figures; accepted by JCAP. Replaced with published version

Mind the Gap: Tightening the Mass-Richness Relation with Magnitude Gaps

Authors: Andrew P. Hearin, Andrew R. Zentner, Jeffrey A. Newman, Andreas A. Berlind

Abstract: We investigate the potential to improve optical tracers of cluster mass by exploiting measurements of the magnitude gap, m12, defined as the difference between the r-band absolute magnitude of the two brightest cluster members. We find that in a mock sample of galaxy groups and clusters constructed from the Bolshoi simulation, the scatter about the mass-richness relation decreases by 15-20% when m… ▽ More We investigate the potential to improve optical tracers of cluster mass by exploiting measurements of the magnitude gap, m12, defined as the difference between the r-band absolute magnitude of the two brightest cluster members. We find that in a mock sample of galaxy groups and clusters constructed from the Bolshoi simulation, the scatter about the mass-richness relation decreases by 15-20% when magnitude gap information is included. A similar trend is evident in a volume-limited, spectroscopic sample of galaxy groups observed in the Sloan Digital Sky Survey (SDSS). We find that SDSS groups with small magnitude gaps are richer than large-gap groups at fixed values of the one-dimensional velocity dispersion among group members sigma_v, which we use as a mass proxy. We demonstrate explicitly that m12 contains information about cluster mass that supplements the information provided by group richness and the luminosity of the brightest cluster galaxy, L_bcg. In so doing, we show that the luminosities of the members of a group with richness N are inconsistent with the distribution of luminosities that results from N random draws from the global galaxy luminosity function. As the cosmological constraining power of galaxy clusters is limited by the precision in cluster mass determination, our findings suggest a new way to improve the cosmological constraints derived from galaxy clusters. △ Less

Submitted 30 April, 2013; v1 submitted 4 July, 2012; originally announced July 2012.

Comments: references added

Dark Matter Direct Search Rates in Simulations of the Milky Way and Sagittarius Stream

Authors: Chris W. Purcell, Andrew R. Zentner, Mei-Yu Wang

Abstract: We analyze self-consistent N-body simulations of the Milky Way disk and the ongoing disruption of the Sagittarius dwarf satellite to study the effect of Sagittarius tidal debris on dark matter detection experiments. In agreement with significant previous work, we reiterate that the standard halo model is insufficient to describe the non-Maxwellian velocity distribution of the Milky Way halo in our… ▽ More We analyze self-consistent N-body simulations of the Milky Way disk and the ongoing disruption of the Sagittarius dwarf satellite to study the effect of Sagittarius tidal debris on dark matter detection experiments. In agreement with significant previous work, we reiterate that the standard halo model is insufficient to describe the non-Maxwellian velocity distribution of the Milky Way halo in our equilibrium halo-only and halo/galaxy models, and offer suggestions for correcting for this discrepancy. More importantly, we emphasize that the dark matter component of the leading tidal arm of the Sagittarius dwarf is significantly more extended than the stellar component of the arm, since the dark matter and stellar streams are not necessarily coaxial and may be offset by several kpc at the point at which they impact the Galactic disk. This suggests that the dark matter component of the Sagittarius debris is likely to have a non-negligible influence on dark matter detection experiments even when the stellar debris is centered several kpc from the solar neighborhood. Relative to models without an infalling Sagittarius dwarf, the Sagittarius dark matter debris in our models induces an energy-dependent enhancement of direct search event rates of as much as ~20 - 45%, an energy-dependent reduction in the amplitude of the annual modulation of the event rate by as much as a factor of two, a shift in the phase of the annual modulation by as much as ~20 days, and a shift in the recoil energy at which the modulation reverses phase. These influences of Sagittarius are of general interest in the interpretation of dark matter searches, but may be particularly important in the case of relatively light (m_X < 20 GeV) dark matter because the Sagittarius stream impacts the solar system at high speed compared to the primary halo dark matter. △ Less

Submitted 29 July, 2012; v1 submitted 29 March, 2012; originally announced March 2012.

Comments: 24 pages, 5 figures; to appear in JCAP; revised to reflect referee's comments

Journal ref: JCAP08(2012)027

Effects of Unstable Dark Matter on Large-Scale Structure and Constraints from Future Surveys

Authors: Mei-Yu Wang, Andrew R. Zentner

Abstract: In this paper we explore the effect of decaying dark matter (DDM) on large-scale structure and possible constraints from galaxy imaging surveys. DDM models have been studied, in part, as a way to address apparent discrepancies between the predictions of standard cold dark matter models and observations of galactic structure. Our study is aimed at developing independent constraints on these models.… ▽ More In this paper we explore the effect of decaying dark matter (DDM) on large-scale structure and possible constraints from galaxy imaging surveys. DDM models have been studied, in part, as a way to address apparent discrepancies between the predictions of standard cold dark matter models and observations of galactic structure. Our study is aimed at developing independent constraints on these models. In such models, DDM decays into a less massive, stable dark matter (SDM) particle and a significantly lighter particle. The small mass splitting between the parent DDM and the daughter SDM provides the SDM with a recoil or "kick" velocity vk, inducing a free-streaming suppression of matter fluctuations. This suppression may be probed via weak lensing power spectra measured by a number of forthcoming imaging surveys that aim primarily to constrain dark energy. Using scales on which linear perturbation theory alone is valid (multipoles < 300), surveys like Euclid or LSST can be sensitive to vk > 90 km/s for lifetimes ~ 1-5 Gyr. To estimate more aggressive constraints, we model nonlinear corrections to lensing power using a simple halo evolution model that is in good agreement with numerical simulations. In our most ambitious forecasts, using multipoles < 3000, we find that imaging surveys can be sensitive to vk ~ 10 km/s for lifetimes < 10 Gyr. Lensing will provide a particularly interesting complement to existing constraints in that they will probe the long lifetime regime far better than contemporary techniques. A caveat to these ambitious forecasts is that the evolution of perturbations on nonlinear scales will need to be well calibrated by numerical simulations before they can be realized. This work motivates the pursuit of such a numerical simulation campaign to constrain dark matter with cosmological weak lensing. △ Less

Submitted 11 January, 2012; originally announced January 2012.

Comments: 15 pages, 7 figures. Submitted to PRD

Constraining Satellite Galaxy Stellar Mass Loss and Predicting Intrahalo Light I: Framework and Results at Low Redshift

Authors: Douglas F. Watson, Andreas A. Berlind, Andrew R. Zentner

Abstract: We introduce a new technique that uses galaxy clustering to constrain how satellite galaxies lose stellar mass and contribute to the diffuse "intrahalo light" (IHL). We implement two models that relate satellite galaxy stellar mass loss to the detailed knowledge of subhalo dark matter mass loss. Model 1 assumes that the fractional stellar mass loss of a galaxy is proportional to the fractional amo… ▽ More We introduce a new technique that uses galaxy clustering to constrain how satellite galaxies lose stellar mass and contribute to the diffuse "intrahalo light" (IHL). We implement two models that relate satellite galaxy stellar mass loss to the detailed knowledge of subhalo dark matter mass loss. Model 1 assumes that the fractional stellar mass loss of a galaxy is proportional to the fractional amount of dark matter mass loss of its subhalo. Model 2 accounts for a delay in the time that stellar mass is lost since the galaxy resides deep in the potential well of the subhalo which may experience dark matter mass loss for some time before the galaxy is affected. We use these models to predict the stellar masses of a population of galaxies and use abundance matching to predict the clustering of several r-band luminosity threshold samples from the Sloan Digital Sky Survey. Abundance matching assuming no stellar mass loss (akin to abundance matching at the time of subhalo infall) over-estimates the correlation function on small scales (<~ 1 Mpc), while allowing too much stellar mass loss leads to an under-estimate. For each sample, we are thus able to constrain the amount of stellar mass loss required to match the observed clustering. We find that less luminous satellite galaxies experience more efficient stellar mass loss than luminous satellites. From these models, we can infer the amount of stellar mass that is deposited into the IHL. We find that both of our model predictions for the mean amount of IHL as a function of halo mass are consistent with current observational measurements. However, our two models predict a different amount of scatter in the IHL from halo to halo, with Model 2 being favored by observations. This demonstrates that a comparison to IHL measurements provides independent verification of our stellar mass loss models. (Abridged) △ Less

Submitted 6 June, 2013; v1 submitted 11 January, 2012; originally announced January 2012.

Comments: 16 pages, 10 figures, Accepted to ApJ

Constraining Cluster Physics with the Shape of X-ray Clusters: Comparison of Local X-ray Clusters versus LCDM Clusters

Authors: Erwin T. Lau, Daisuke Nagai, Andrey V. Kravtsov, Alexey Vikhlinin, Andrew R. Zentner

Abstract: Simulations of cluster formation have demonstrated that condensation of baryons into central galaxies during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder, even at radii as large as half of the virial radius. However, such simulations generally predict stellar fractions within cluster virial radii that are ~2 to 3 times larger than the stell… ▽ More Simulations of cluster formation have demonstrated that condensation of baryons into central galaxies during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder, even at radii as large as half of the virial radius. However, such simulations generally predict stellar fractions within cluster virial radii that are ~2 to 3 times larger than the stellar masses deduced from observations. In this work we compare ellipticity profiles of clusters simulated with and without baryonic cooling to the cluster ellipticity profiles derived from Chandra and ROSAT observations in an effort to constrain the fraction of gas that cools and condenses into the central galaxies within clusters. We find that the observed ellipticity profiles are fairly constant with radius, with an average ellipticity of 0.18 +/- 0.05. The observed ellipticity profiles are in good agreement with the predictions of non-radiative simulations. On the other hand, the ellipticity profiles of the clusters in simulations that include radiative cooling, star formation, and supernova feedback (but no AGN feedback) deviate significantly from the observed ellipticity profiles at all radii. The simulations with cooling overpredict (underpredict) ellipticity in the inner (outer) regions of galaxy clusters. By comparing the simulations with and without cooling, we show that the cooling of gas via cooling flows in the central regions of simulated clusters causes the gas distribution to be more oblate in the central regions, but makes the outer gas distribution more spherical. We find that late-time gas cooling and star formation are responsible for the significantly oblate gas distributions in cluster cores, but the gas shapes outside of cluster cores are set primarily by baryon dissipation at high redshift z > 2. △ Less

Submitted 5 September, 2012; v1 submitted 10 January, 2012; originally announced January 2012.

Comments: 10 pages, 6 figures, matching the published version in ApJ. Corrected missing reference in the arxiv version

Journal ref: ApJ 2012, Vol 755, 116

Pixel-z: Studying Substructure and Stellar Populations in Galaxies out to z~3 using Pixel Colors I. Systematics

Authors: Niraj Welikala, Andrew M. Hopkins, Brant E. Robertson, Andrew J. Connolly, Lidia Tasca, Anton M. Koekemoer, Olivier Ilbert, Sandro Bardelli, Jean-Paul Kneib, Andrew R. Zentner

Abstract: We perform a pixel-by-pixel analysis of 467 galaxies in the GOODS-VIMOS survey to study systematic effects in extracting properties of stellar populations (age, dust, metallicity and SFR) from pixel colors using the pixel-z method. The systematics studied include the effect of the input stellar population synthesis model, passband limitations and differences between individual SED fits to pixels a… ▽ More We perform a pixel-by-pixel analysis of 467 galaxies in the GOODS-VIMOS survey to study systematic effects in extracting properties of stellar populations (age, dust, metallicity and SFR) from pixel colors using the pixel-z method. The systematics studied include the effect of the input stellar population synthesis model, passband limitations and differences between individual SED fits to pixels and global SED-fitting to a galaxy's colors. We find that with optical-only colors, the systematic errors due to differences among the models are well constrained. The largest impact on the age and SFR e-folding time estimates in the pixels arises from differences between the Maraston models and the Bruzual&Charlot models, when optical colors are used. This results in systematic differences larger than the 2σ uncertainties in over 10 percent of all pixels in the galaxy sample. The effect of restricting the available passbands is more severe. In 26 percent of pixels in the full sample, passband limitations result in systematic biases in the age estimates which are larger than the 2σ uncertainties. Systematic effects from model differences are reexamined using Near-IR colors for a subsample of 46 galaxies in the GOODS-NICMOS survey. For z > 1, the observed optical/NIR colors span the rest frame UV-optical SED, and the use of different models does not significantly bias the estimates of the stellar population parameters compared to using optical-only colors. We then illustrate how pixel-z can be applied robustly to make detailed studies of substructure in high redshift galaxies such as (a) radial gradients of age, SFR, sSFR and dust and (b) the distribution of these properties within subcomponents such as spiral arms and clumps. Finally, we show preliminary results from the CANDELS survey illustrating how the new HST/WFC3 data can be exploited to probe substructure in z~1-3 galaxies. △ Less

Submitted 12 December, 2011; originally announced December 2011.

Comments: 37 pages, 21 figures, submitted to ApJ