• ### Cosmic Shear as a Probe of Galaxy Formation Physics(1605.09056)

Sept. 4, 2018 astro-ph.CO, astro-ph.GA
We evaluate the potential for current and future cosmic shear measurements from large galaxy surveys to constrain the impact of baryonic physics on the matter power spectrum. We do so using a model-independent parameterization that describes deviations of the matter power spectrum from the dark-matter-only case as a set of principal components that are localized in wavenumber and redshift. We perform forecasts for a variety of current and future datasets, and find that at least ~90% of the constraining power of these datasets is contained in no more than nine principal components. The constraining power of different surveys can be quantified using a figure of merit defined relative to currently available surveys. With this metric, we find that the final Dark Energy Survey dataset (DES Y5) and the Hyper Suprime Cam Survey will be roughly an order of magnitude more powerful than existing data in constraining baryonic effects. Upcoming Stage IV surveys (LSST, Euclid, and WFIRST) will improve upon this by a further factor of a few. We show that this conclusion is robust to marginalization over several key systematics. The ultimate power of cosmic shear to constrain galaxy formation is dependent on understanding systematics in the shear measurements at small (sub-arcminute) scales. If these systematics can be sufficiently controlled, cosmic shear measurements from DES Y5 and other future surveys have the potential to provide a very clean probe of galaxy formation and to strongly constrain a wide range of predictions from modern hydrodynamical simulations.
• ### Mapping Extragalactic Dark Matter Annihilation with Galaxy Surveys: A Systematic Study of Stacked Group Searches(1709.00416)

Dark matter in the halos surrounding galaxy groups and clusters can annihilate to high-energy photons. Recent advancements in the construction of galaxy group catalogs provide many thousands of potential extragalactic targets for dark matter. In this paper, we outline a procedure to infer the dark matter signal associated with a given galaxy group. Applying this procedure to a catalog of sources, one can create a full-sky map of the brightest extragalactic dark matter targets in the nearby Universe ($z\lesssim 0.03$), supplementing sources of dark matter annihilation from within the Local Group. As with searches for dark matter in dwarf galaxies, these extragalactic targets can be stacked together to enhance the signals associated with dark matter. We validate this procedure on mock $\textit{Fermi}$ gamma-ray data sets using a galaxy catalog constructed from the $\texttt{DarkSky}$ $N$-body cosmological simulation and demonstrate that the limits are robust, at $\mathcal{O}(1)$ levels, to systematic uncertainties on halo mass and concentration. We also quantify other sources of systematic uncertainty arising from the analysis and modeling assumptions. Our results suggest that a stacking analysis using galaxy group catalogs provides a powerful opportunity to discover extragalactic dark matter and complements existing studies of Milky Way dwarf galaxies.
• ### The Aemulus Project III: Emulation of the Galaxy Correlation Function(1804.05867)

Using the N-body simulations of the AEMULUS Project, we construct an emulator for the non-linear clustering of galaxies in real and redshift space. We construct our model of galaxy bias using the halo occupation framework, accounting for possible velocity bias. The model includes 15 parameters, including both cosmological and galaxy bias parameters. We demonstrate that our emulator achieves ~ 1% precision at the scales of interest, 0.1<r<10 h^{-1} Mpc, and recovers the true cosmology when tested against independent simulations. Our primary parameters of interest are related to the growth rate of structure, f, and its degenerate combination fsigma_8. Using this emulator, we show that the constraining power on these parameters monotonically increases as smaller scales are included in the analysis, all the way down to 0.1 h^{-1} Mpc. For a BOSS-like survey, the constraints on fsigma_8 from r<30 h^{-1} Mpc scales alone are more than a factor of two tighter than those from the fiducial BOSS analysis of redshift-space clustering using perturbation theory at larger scales. The combination of real- and redshift-space clustering allows us to break the degeneracy between f and sigma_8, yielding a 9% constraint on f alone for a BOSS-like analysis. The current AEMULUS simulations limit this model to surveys of massive galaxies. Future simulations will allow this framework to be extended to all galaxy target types, including emission-line galaxies.
• The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ~14,000 square degrees of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamic observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering an inference-based catalog which includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4um, 4.6um, 12um and 22um) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.
• ### Effects of NII and H$\alpha$ Line Blending on the WFIRST Galaxy Redshift Survey(1804.08061)

April 22, 2018 astro-ph.CO
The Wide Field Infrared Survey Telescope (WFIRST) will conduct a range of observational programs, including a galaxy redshift survey using the redshifted H$\alpha$ emission line as the primary feature for spectroscopic redshift determination. Due to the finite spectroscopic resolution of the WFIRST grism, the H$\alpha$ line and the neighboring NII lines are not separately resolved, leading to a potential bias in the redshift that depends on the NII/H$\alpha$ ratio. We investigate how this bias propagates into the redshift-space correlation functions and measurements of cosmic expansion and growth of structure obtained from the WFIRST redshift survey. Using the Buzzard simulation of mock galaxies, we explore the effect of line blending and mass (hence environment)-dependent NII/H$\alpha$ ratio on redshift-space distortion (RSD) and baryon acoustic oscillation (BAO) measurements. We measured the BAO shift parameters $\alpha_{\parallel}$ and $\alpha_{\perp}$, as well as the growth of structure parameter $f_{v}$, and calculated their errors based on the correlations between the NII/H$\alpha$ ratio and large-scale structure. We found that $\Delta\alpha_{\parallel} = 0.31 \pm 0.23 \%$ ($0.26\pm0.17\%$), $\Delta\alpha_{\perp} = -0.10\pm0.10\%$ ($-0.12 \pm 0.11 \%$), and $\Delta f_{v} = 0.17\pm0.33\%$ ($-0.20 \pm 0.30\%$) for the redshift 1.355--1.994 (0.700--1.345) bins, which would use approximately 18$\%$, 9$\%$, and 7$\%$ of the systematic error budget in a root-sum-square sense. These errors may be small enough that even without mitigation they are tolerable, and we discuss mitigations that may further reduce this error.
• ### The Aemulus Project I: Numerical Simulations for Precision Cosmology(1804.05865)

April 16, 2018 astro-ph.CO, astro-ph.GA
The rapidly growing statistical precision of galaxy surveys has lead to a need for ever-more precise predictions of the observables used to constrain cosmological and galaxy formation models. The primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper we present one such suite of simulations, forming the basis for the AEMULUS Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h^-1 Gpc)^3 simulations with mass resolution and force softening of 3.51\times 10^10 (Omega_m / 0.3) ~ h^-1 M_sun and 20 ~ h^-1 kpc respectively in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) for halos with more than 500 (200) particles respectively and scales of r>200 ~ h^-1 kpc in real space or k ~ 3 h Mpc^-1 in harmonic space for z\le 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k ~ 4 h^-1 Mpc. We provide the halo catalogs and snapshots detailed in this work to the community at https://AemulusProject.github.io.
• ### The Aemulus Project II: Emulating the Halo Mass Function(1804.05866)

Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a halo mass function emulator and demonstrate improved accuracy relative to state-of-the-art analytic models. In this work, mass is defined using an overdensity criteria of 200 relative to the mean background density. Our emulator is constructed from the AEMULUS simulations, a suite of 40 N-body simulations with snapshots from z=3 to z=0. These simulations cover the flat wCDM parameter space allowed by recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results, varying the parameters w, Omega_m, Omega_b, sigma_8, N_{eff}, n_s, and H_0. We validate our emulator using five realizations of seven different cosmologies, for a total of 35 test simulations. These test simulations were not used in constructing the emulator, and were run with fully independent initial conditions. We use our test simulations to characterize the modeling uncertainty of the emulator, and introduce a novel way of marginalizing over the associated systematic uncertainty. We confirm non-universality in our halo mass function emulator as a function of both cosmological parameters and redshift. Our emulator achieves better than 1% precision over much of the relevant parameter space, and we demonstrate that the systematic uncertainty in our emulator will remain a negligible source of error for cluster abundance studies through at least the LSST Year 1 data set.
• ### The Connection between Galaxies and their Dark Matter Halos(1804.03097)

April 9, 2018 astro-ph.CO, astro-ph.GA
In our modern understanding of galaxy formation, every galaxy forms within a dark matter halo. The formation and growth of galaxies over time is connected to the growth of the halos in which they form. The advent of large galaxy surveys as well as high-resolution cosmological simulations has provided a new window into the statistical relationship between galaxies and halos and its evolution. Here we define this galaxy-halo connection as the multi-variate distribution of galaxy and halo properties that can be derived from observations and simulations. This connection provides a key test of physical galaxy formation models; it also plays an essential role in constraints of cosmological models using galaxy surveys and in elucidating the properties of dark matter using galaxies. We review techniques for inferring the galaxy-halo connection and the insights that have arisen from these approaches. Some things we have learned are that galaxy formation efficiency is a strong function of halo mass; at its peak in halos around a pivot halo mass of 10^12 Msun, less than 20% of the available baryons have turned into stars by the present day; the intrinsic scatter in galaxy stellar mass is small, less than 0.2 dex at a given halo mass above this pivot mass; below this pivot mass galaxy stellar mass is a strong function of halo mass; the majority of stars over cosmic time were formed in a narrow region around this pivot mass. We also highlight key open questions about how galaxies and halos are connected, including understanding the correlations with secondary properties and the connection of these properties to galaxy clustering.
• ### DESCQA: An Automated Validation Framework for Synthetic Sky Catalogs(1709.09665)

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

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission, the spectroscopic follow-up to WFIRST, multiplexing its scientific return by obtaining deep 1 to 4 micron slit spectroscopy for ~90% of all galaxies imaged by the ~2200 sq deg WFIRST High Latitude Survey at z > 0.5. ATLAS spectroscopy will measure accurate and precise redshifts for ~300M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a 3D map of the Universe with ~Mpc resolution in redshift space. ATLAS will: (1) Revolutionize galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Open a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Explore the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS is a 1.5m telescope with a field of view (FoV) of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 600, a wavelength range of 1-4 microns, and a spectroscopic multiplex factor ~5,000-10,000. ATLAS is designed to fit within the NASA probe-class mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology (DMDs can reach TRL 6 within 2 years). ATLAS will lead to transformative science over the entire range of astrophysics.
• ### Spatial Clustering of Dark Matter Halos: Secondary Bias, Neighbor Bias, and the Influence of Massive Neighbors on Halo Properties(1708.08451)

Jan. 9, 2018 astro-ph.CO
We explore the phenomenon commonly known as halo assembly bias, whereby dark matter halos of the same mass are found to be more or less clustered when a second halo property is considered, for halos in the mass range $3.7 \times 10^{11} \; h^{-1} \mathrm{M_{\odot}} - 5.0 \times 10^{13} \; h^{-1} \mathrm{M_{\odot}}$. Using the Large Suite of Dark Matter Simulations (LasDamas) we consider nine commonly used halo properties and find that a clustering bias exists if halos are binned by mass or by any other halo property. This secondary bias implies that no single halo property encompasses all the spatial clustering information of the halo population. The mean values of some halo properties depend on their halo's distance to a more massive neighbor. Halo samples selected by having high values of one of these properties therefore inherit a neighbor bias such that they are much more likely to be close to a much more massive neighbor. This neighbor bias largely accounts for the secondary bias seen in halos binned by mass and split by concentration or age. However, halos binned by other mass-like properties still show a secondary bias even when the neighbor bias is removed. The secondary bias of halos selected by their spin behaves differently than that for other halo properties, suggesting that the origin of the spin bias is different than of other secondary biases.
• ### On estimation of contamination from hydrogen cyanide in carbon monoxide line intensity mapping(1706.03005)

Aug. 31, 2017 astro-ph.CO, astro-ph.GA
Line-intensity mapping surveys probe large-scale structure through spatial variations in molecular line emission from a population of unresolved cosmological sources. Future such surveys of carbon monoxide line emission, specifically the CO(1-0) line, face potential contamination from a disjoint population of sources emitting in a hydrogen cyanide emission line, HCN(1-0). This paper explores the potential range of the strength of HCN emission and its effect on the CO auto power spectrum, using simulations with an empirical model of the CO/HCN--halo connection. We find that effects on the observed CO power spectrum depend on modeling assumptions but are very small for our fiducial model based on our understanding of the galaxy--halo connection, with the bias in overall CO detection significance due to HCN expected to be less than 1%.
• ### Uncertainties in Parameters Estimated with Neural Networks: Application to Strong Gravitational Lensing(1708.08843)

Aug. 29, 2017 astro-ph.CO, astro-ph.IM
In Hezaveh et al. 2017 we showed that deep learning can be used for model parameter estimation and trained convolutional neural networks to determine the parameters of strong gravitational lensing systems. Here we demonstrate a method for obtaining the uncertainties of these parameters. We review the framework of variational inference to obtain approximate posteriors of Bayesian neural networks and apply it to a network trained to estimate the parameters of the Singular Isothermal Ellipsoid plus external shear and total flux magnification. We show that the method can capture the uncertainties due to different levels of noise in the input data, as well as training and architecture-related errors made by the network. To evaluate the accuracy of the resulting uncertainties, we calculate the coverage probabilities of marginalized distributions for each lensing parameter. By tuning a single hyperparameter, the dropout rate, we obtain coverage probabilities approximately equal to the confidence levels for which they were calculated, resulting in accurate and precise uncertainty estimates. Our results suggest that neural networks can be a fast alternative to Monte Carlo Markov Chains for parameter uncertainty estimation in many practical applications, allowing more than seven orders of magnitude improvement in speed.
• ### Weak Lensing Measurement of the Mass--Richness Relation of SDSS redMaPPer Clusters(1603.06953)

April 3, 2017 astro-ph.CO
We perform a measurement of the mass--richness relation of the redMaPPer galaxy cluster catalogue using weak lensing data from the Sloan Digital Sky Survey. We have carefully characterized a broad range of systematic uncertainties, including shear calibration errors, photo-$z$ biases, dilution by member galaxies, source obscuration, magnification bias, incorrect assumptions about cluster mass profiles, cluster centering, halo triaxiality, and projection effects. We also compare measurements of the lensing signal from two independently-produced shear and photometric redshift catalogues to characterize systematic errors in the lensing signal itself. Using a sample of 5,570 clusters from $0.1\le z\le 0.33$, the normalization of our power-law mass vs.\ $\lambda$ relation is $\log_{10}[M_{200m}/h^{-1}\ M_{\odot}]$ = $14.344 \pm 0.021$ (statistical) $\pm 0.023$ (systematic) at a richness $\lambda=40$, a 7 per cent calibration uncertainty, with a power-law index of $1.33^{+0.09}_{-0.10}$ ($1\sigma$). The detailed systematics characterization in this work renders it the definitive weak lensing mass calibration for SDSS redMaPPer clusters at this time.
• DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.
• DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.
• ### The Faber-Jackson relation and Fundamental Plane from halo abundance matching(1604.04670)

Nov. 3, 2016 astro-ph.CO, astro-ph.GA
The Fundamental Plane (FP) describes the relation between the stellar mass, size, and velocity dispersion of elliptical galaxies; the Faber-Jackson relation (FJR) is its projection onto {mass, velocity} space. In this work we redeploy and expand the framework of Desmond & Wechsler (2015) to ask whether abundance matching-based LCDM models that have shown success in matching the spatial distribution of galaxies are also capable of explaining key properties of the FJR and FP, including their scatter. Within our framework, agreement with the normalisation of the FJR requires haloes to expand in response to disc formation. We find that the tilt of the FP may be explained by a combination of the observed non-homology in galaxy structure and the variation in mass-to-light ratio produced by abundance matching with a universal initial mass function (IMF), provided that the anisotropy of stellar motions is taken into account. However, the predicted scatter around the FP is considerably increased by situating galaxies in cosmologically-motivated haloes due to variations in halo properties at fixed stellar mass, and appears to exceed that of the data. This implies that additional correlations between galaxy and halo variables may be required to fully reconcile these models with elliptical galaxy scaling relations.
• ### EmpiriciSN: Re-sampling Observed Supernova/Host Galaxy Populations using an XD Gaussian Mixture Model(1611.00363)

Nov. 1, 2016 astro-ph.IM
We describe two new open source tools written in Python for performing extreme deconvolution Gaussian mixture modeling (XDGMM) and using a conditioned model to re-sample observed supernova and host galaxy populations. XDGMM is new program for using Gaussian mixtures to do density estimation of noisy data using extreme deconvolution (XD) algorithms that has functionality not available in other XD tools. It allows the user to select between the AstroML (Vanderplas et al. 2012; Ivezic et al. 2015) and Bovy et al. (2011) fitting methods and is compatible with scikit-learn machine learning algorithms (Pedregosa et al. 2011). Most crucially, it allows the user to condition a model based on the known values of a subset of parameters. This gives the user the ability to produce a tool that can predict unknown parameters based on a model conditioned on known values of other parameters. EmpiriciSN is an example application of this functionality that can be used for fitting an XDGMM model to observed supernova/host datasets and predicting likely supernova parameters using on a model conditioned on observed host properties. It is primarily intended for simulating realistic supernovae for LSST data simulations based on empirical galaxy properties.
• ### Biased Tracers in Redshift Space in the EFT of Large-Scale Structure(1610.09321)

Oct. 28, 2016 hep-th, hep-ph, gr-qc, astro-ph.CO
The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a novel formalism that is able to accurately predict the clustering of large-scale structure (LSS) in the mildly non-linear regime. Here we provide the first computation of the power spectrum of biased tracers in redshift space at one loop order, and we make the associated code publicly available. We compare the multipoles $\ell=0,2$ of the redshift-space halo power spectrum, together with the real-space matter and halo power spectra, with data from numerical simulations at $z=0.67$. For the samples we compare to, which have a number density of $\bar n=3.8 \cdot 10^{-2}(h \ {\rm Mpc}^{-1})^3$ and $\bar n=3.9 \cdot 10^{-4}(h \ {\rm Mpc}^{-1})^3$, we find that the calculation at one-loop order matches numerical measurements to within a few percent up to $k\simeq 0.43 \ h \ {\rm Mpc}^{-1}$, a significant improvement with respect to former techniques. By performing the so-called IR-resummation, we find that the Baryon Acoustic Oscillation peak is accurately reproduced. Based on the results presented here, long-wavelength statistics that are routinely observed in LSS surveys can be finally computed in the EFTofLSS. This formalism thus is ready to start to be compared directly to observational data.
• ### The connection between the host halo and the satellite galaxies of the Milky Way(1605.02075)

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

April 14, 2016 astro-ph.CO, astro-ph.GA
We study the evolution of the stellar component and the metallicity of both the intracluster medium and of stars in massive ($M_{\rm vir}\approx 6\times 10^{14}$ M$_{\odot}/h$) simulated galaxy clusters from the Rhapsody-G suite in detail and compare them to observational results. The simulations were performed with the AMR code RAMSES and include the effect of AGN feedback at the sub-grid level. AGN feedback is required to produce realistic galaxy and cluster properties and plays a role in mixing material in the central regions and regulating star formation in the central galaxy. In both our low and high resolution runs with fiducial stellar yields, we find that stellar and ICM metallicities are a factor of two lower than in observations. We find that cool core clusters exhibit steeper metallicity gradients than non-cool core clusters, in qualitative agreement with observations. We verify that the ICM metallicities measured in the simulation can be explained by a simple "regulator" model in which the metallicity is set by a balance of stellar yield and gas accretion. It is plausible that a combination of higher resolution and higher metal yield in AMR simulation would allow the metallicity of simulated clusters to match observed values; however this hypothesis needs to be tested with future simulations. Comparison to recent literature highlights that results concerning the metallicity of clusters and cluster galaxies might depend sensitively on the scheme chosen to solve the hydrodynamics.
• ### Constraining the Mass-Richness Relationship of redMaPPer Clusters with Angular Clustering(1604.00048)

March 31, 2016 astro-ph.CO
The potential of using cluster clustering for calibrating the mass-observable relation of galaxy clusters has been recognized theoretically for over a decade. Here, we demonstrate the feasibility of this technique to achieve high precision mass calibration using redMaPPer clusters in the Sloan Digital Sky Survey North Galactic Cap. By including cross-correlations between several richness bins in our analysis we significantly improve the statistical precision of our mass constraints. The amplitude of the mass-richness relation is constrained to 7% statistical precision. However, the error budget is systematics dominated, reaching an 18% total error that is dominated by theoretical uncertainty in the bias-mass relation for dark matter halos. We perform a detailed treatment of the effects of assembly bias on our analysis, finding that the contribution of such effects to our parameter uncertainties is somewhat greater than that of measurement noise. We confirm the results from Miyatake et al. (2015) that the clustering amplitude of redMaPPer clusters depends on galaxy concentration, and provide additional evidence in support of this effect being due to some form of assembly bias. The results presented here demonstrate the power of cluster clustering for mass calibration and cosmology provided the current theoretical systematics can be ameliorated.
• ### Connecting CO Intensity Mapping to Molecular Gas and Star Formation in the Epoch of Galaxy Assembly(1503.08833)

Feb. 1, 2016 astro-ph.CO, astro-ph.GA
Intensity mapping, which images a single spectral line from unresolved galaxies across cosmological volumes, is a promising technique for probing the early universe. Here we present predictions for the intensity map and power spectrum of the CO(1-0) line from galaxies at z ~ 2.4 -2.8, based on a parameterized model for the galaxy-halo connection, and demonstrate the extent to which properties of high-redshift galaxies can be directly inferred from such observations. We find that our fiducial prediction should be detectable by a realistic experiment. Motivated by significant modeling uncertainties, we demonstrate the effect on the power spectrum of varying each parameter in our model. Using simulated observations, we infer constraints on our model parameter space with an MCMC procedure, and show corresponding constraints on the LIR-LCO relation and the CO luminosity function. These constraints would be complementary to current high-redshift galaxy observations, which can detect the brightest galaxies but not complete samples from the faint end of the luminosity function. By probing these populations in aggregate, CO intensity mapping could be a valuable tool for probing molecular gas and its relation to star formation in high-redshift galaxies.
• ### The Eating Habits of Milky Way Mass Halos: Destroyed Dwarf Satellites and the Metallicity Distribution of Accreted Stars(1601.07905)

Jan. 28, 2016 astro-ph.GA
We study the mass spectrum of destroyed dwarfs that contribute to the accreted stellar mass of Milky Way (MW) mass M_vir ~ 10^12.1 M_sun) halos using a suite of 45 zoom-in, dissipationless simulations. Empirical models are employed to relate (peak) subhalo mass to dwarf stellar mass, and we use constraints from z=0 observations and hydrodynamical simulations to estimate the metallicity distribution of the accreted stellar material. The dominant contributors to the accreted stellar mass are relatively massive dwarfs with M_star ~ 10^8-10^10 M_sun. Halos with more quiescent accretion histories tend to have lower mass progenitors (10^8-10^9 M_sun), and lower overall accreted stellar masses. Ultra-faint mass (M_star < 10^5 M_sun) dwarfs contribute a negligible amount (<< 1%) to the accreted stellar mass and, despite having low average metallicities, supply a small fraction (~2-5 %) of the very metal-poor stars with [Fe/H] < -2. Dwarfs with masses 10^5 < M_star/M_sun < 10^8 provide a substantial amount of the very metal-poor stellar material (~40-80 %), and even relatively metal-rich dwarfs with M_star > 10^8 M_sun can contribute a considerable fraction (~20-60 %) of metal-poor stars if their metallicity distributions have significant metal-poor tails. Finally, we find that the generic assumption of a quiescent assembly history for the MW halo seems to be in tension with the mass spectrum of its surviving dwarfs. We suggest that the MW could be a "transient fossil"; a quiescent halo with a recent accretion event(s) that disguises the preceding formation history of the halo.
• ### Galaxy Cluster Mass Estimation from Stacked Spectroscopic Analysis(1601.05773)

Jan. 21, 2016 astro-ph.CO, astro-ph.GA
We use simulated galaxy surveys to study: i) how galaxy membership in redMaPPer clusters maps to the underlying halo population, and ii) the accuracy of a mean dynamical cluster mass, $M_\sigma(\lambda)$, derived from stacked pairwise spectroscopy of clusters with richness $\lambda$. Using $\sim\! 130,000$ galaxy pairs patterned after the SDSS redMaPPer cluster sample study of Rozo et al. (2015 RMIV), we show that the pairwise velocity PDF of central--satellite pairs with $m_i < 19$ in the simulation matches the form seen in RMIV. Through joint membership matching, we deconstruct the main Gaussian velocity component into its halo contributions, finding that the top-ranked halo contributes $\sim 60\%$ of the stacked signal. The halo mass scale inferred by applying the virial scaling of Evrard et al. (2008) to the velocity normalization matches, to within a few percent, the log-mean halo mass derived through galaxy membership matching. We apply this approach, along with mis-centering and galaxy velocity bias corrections, to estimate the log-mean matched halo mass at $z=0.2$ of SDSS redMaPPer clusters. Employing the velocity bias constraints of Guo et al. (2015), we find $\langle \ln(M_{200c})|\lambda \rangle = \ln(M_{30}) + \alpha_m \ln(\lambda/30)$ with $M_{30} = 1.56 \pm 0.35 \times 10^{14} M_\odot$ and $\alpha_m = 1.31 \pm 0.06_{stat} \pm 0.13_{sys}$. Systematic uncertainty in the velocity bias of satellite galaxies overwhelmingly dominates the error budget.