Search

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., Comment: submitted to MNRAS

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.

Within the $\Lambda$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 ($\lambda_{\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) $\lambda_{\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., Comment: Accepted to MNRAS 09/26/2023

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., Comment: 16+5 pages, 9+3 figures, 5+2 tables, Fig. 6 shows the main result. To be submitted to MNRAS, comments welcome

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.

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., Comment: 14 pages, 8 figures, submitted to Monthly Notices of the Royal Astronomical Society

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., Comment: 19 pages, 9 figures, appendices, submitted to MNRAS

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., Comment: 15 pages, 7 + 1 figures, accepted for publication in MNRAS, comments welcome

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 $\sigma_{\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., Comment: 24 pages, 25 figures, published in MNRAS

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 $\lambda_{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., Comment: Discussion on isocurvature perturbations. To appear in PRD

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.

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., Comment: 10 pages, 2 figures, submitted to MNRAS, comments welcome

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., Comment: 30 pages, 14 figures, accepted for publication in MNRAS

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 $\Delta \Sigma(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}})+\Delta \Sigma(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 $\Delta \Sigma(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., Comment: Figures 3 and 4 show the main results. Published in Monthly Notices of the Royal Astronomical Society

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., Comment: Science white paper submitted to the Astro2020 Decadal Survey

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., Comment: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

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/ ., Comment: 96 pages, 22 figures, 1 table

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., Comment: 18 pages, 10 figures, submitted to MNRAS, comments welcome

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., Comment: 24 pages, 15 figures; resubmitted to MNRAS after first referee report

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{\Gamma}_k(t)\,t}$ with $\widetilde{\Gamma}_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^{-\Gamma_0 t}\,(t/t_{nr})^{\Gamma_0\,t_{nr}/2}$. The Hubble time scale $\propto 1/H(t)$ introduces an energy uncertainty $\Delta E \sim H(t)$ which relaxes the constraints of kinematic thresholds. This opens new decay channels into heavier particles for $2\pi 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., Comment: more comments. published version

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., Comment: 21 pages, 11 figures, submitted to Monthly Notices of the Royal Astronomical Society

We compute the weak lensing convergence power spectrum, $C^{\kappa\kappa}(\theta)$, 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$., Comment: 9 pages, 3 figures

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., Comment: 11 pages, 4 figures; fixed calculation error in generating main figure

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., Comment: 25 pages, 22 figures, submitted to MNRAS, comments welcome

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., Comment: 24 pages, 15 figures. Accepted for publication in MNRAS

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, $\Delta$. 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 $\Delta$. For conventional halo definitions ($\Delta \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 $\Delta \sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $\Delta$ 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, $\Delta \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., Comment: 19 pages, 13 figures. Updated to published version. Main result summarized in Figure 10

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., Comment: 14+1 pages, 9+1 figures. Figures 2 and 8 highlight the main results. Accepted by MNRAS. Fig 9 updated; main result unchanged

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., Comment: 7 pages, 4 figures

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., Comment: 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

We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python., Comment: Revisions match version accepted for publication in AAS

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., Comment: 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

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., Comment: 14 pages, 11 figures, 3 tables. Updated to include acknowledgements

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., Comment: 14 pages, 12 figures, some typos are corrected, submitted to MNRAS

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., Comment: 22 pages. Appendix. B added. Matches the version accepted by MNRAS. Mock galaxy catalog and HOD table are available at http://www.massivegalaxies.com

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_{\chi} \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..., Comment: 13 pages, 7 figures

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., Comment: 17 pages, 14 figures, references added, labels added in figure 1 & 2, submitted to MNRAS

Systematic uncertainties that have been subdominant in past large-scale structure (LSS) surveys are likely to exceed statistical uncertainties of current and future LSS data sets, potentially limiting the extraction of cosmological information. Here we present a general framework (PCA marginalization) to consistently incorporate systematic effects into a likelihood analysis. This technique naturally accounts for degeneracies between nuisance parameters and can substantially reduce the dimension of the parameter space that needs to be sampled. As a practical application, we apply PCA marginalization to account for baryonic physics as an uncertainty in cosmic shear tomography. Specifically, we use CosmoLike to run simulated likelihood analyses on three independent sets of numerical simulations, each covering a wide range of baryonic scenarios differing in cooling, star formation, and feedback mechanisms. We simulate a Stage III (Dark Energy Survey) and Stage IV (Large Synoptic Survey Telescope/Euclid) survey and find a substantial bias in cosmological constraints if baryonic physics is not accounted for. We then show that PCA marginalization (employing at most 3 to 4 nuisance parameters) removes this bias. Our study demonstrates that it is possible to obtain robust, precise constraints on the dark energy equation of state even in the presence of large levels of systematic uncertainty in astrophysical processes. We conclude that the PCA marginalization technique is a powerful, general tool for addressing many of the challenges facing the precision cosmology program., Comment: Comments welcome, 22 pages, 16 figures

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 |\Delta 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., Comment: No major changes. References added, and text expanded. Published in Phys. Rev. D, and replaced to reflect the published version

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., Comment: 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

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., Comment: 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"

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., Comment: references added, minor adjustments to text and notation

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., Comment: 12 pages, 4 figures, submitted to PRD

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., Comment: 17 pages, 11 figures, accepted to MNRAS

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., Comment: 16 pages, 9 figures. Submitted to Physical Review D. Comments Welcome

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)., Comment: minor changes to draft; conclusions unchanged; version accepted by MNRAS

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], Comment: 12 pages, 3 figures; accepted by JCAP. Replaced with published version

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., Comment: references added