BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Risa H. Wechsler

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMastro-phhep-phhep-exComputationhep-thMachine LearningArtificial Intelligenceastro-ph.HEgr-qc

Catalog footprint

What is connected

102works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

102 published item(s)

preprint2022arXiv

A Spectroscopic Road Map for Cosmic Frontier: DESI, DESI-II, Stage-5

In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Bounds on velocity-dependent dark matter-proton scattering from Milky Way satellite abundance

We use the latest measurements of the Milky Way satellite population from the Dark Energy Survey and Pan-STARRS1 to infer the most stringent astrophysical bound to date on velocity-dependent interactions between dark matter particles and protons. We model the momentum-transfer cross section as a power law of the relative particle velocity $v$ with a free normalizing amplitude, $σ_\text{MT}=σ_0 v^n$, to broadly capture the interactions arising within the non-relativistic effective theory of dark matter-proton scattering. The scattering leads to a momentum and heat transfer between the baryon and dark matter fluids in the early Universe, ultimately erasing structure on small physical scales and reducing the abundance of low-mass halos that host dwarf galaxies today. From the consistency of observations with the cold collisionless dark matter paradigm, using a new method that relies on the most robust predictions of the linear perturbation theory, we infer an upper limit on $σ_0$ of $1.4\times10^{-23}$, $2.1\times10^{-19}$, and $1.0\times10^{-12}\ \mathrm{cm}^2$, for interaction models with $n=2,4,6$, respectively, for a dark matter particle mass of $10\ \mathrm{MeV}$. These results improve observational limits on dark matter--proton scattering by orders of magnitude and thus provide an important guide for viable sub-GeV dark matter candidates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

COMAP Early Science: V. Constraints and Forecasts at $z \sim 3$

We present the current state of models for the $z\sim3$ carbon monoxide (CO) line-intensity signal targeted by the CO Mapping Array Project (COMAP) Pathfinder in the context of its early science results. Our fiducial model, relating dark matter halo properties to CO luminosities, informs parameter priors with empirical models of the galaxy-halo connection and previous CO(1-0) observations. The Pathfinder early science data spanning wavenumbers $k=0.051$-$0.62\,$Mpc$^{-1}$ represent the first direct 3D constraint on the clustering component of the CO(1-0) power spectrum. Our 95% upper limit on the redshift-space clustering amplitude $A_{\rm clust}\lesssim70\,μ$K$^2$ greatly improves on the indirect upper limit of $420\,μ$K$^2$ reported from the CO Power Spectrum Survey (COPSS) measurement at $k\sim1\,$Mpc$^{-1}$. The COMAP limit excludes a subset of models from previous literature, and constrains interpretation of the COPSS results, demonstrating the complementary nature of COMAP and interferometric CO surveys. Using line bias expectations from our priors, we also constrain the squared mean line intensity-bias product, $\langle{Tb}\rangle^2\lesssim50\,μ$K$^2$, and the cosmic molecular gas density, $ρ_\text{H2}<2.5\times10^8\,M_\odot\,$Mpc$^{-3}$ (95% upper limits). Based on early instrument performance and our current CO signal estimates, we forecast that the five-year Pathfinder campaign will detect the CO power spectrum with overall signal-to-noise of 9-17. Between then and now, we also expect to detect the CO-galaxy cross-spectrum using overlapping galaxy survey data, enabling enhanced inferences of cosmic star-formation and galaxy-evolution history.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extending the SAGA Survey (xSAGA) I: Satellite Radial Profiles as a Function of Host Galaxy Properties

We present "Extending the Satellites Around Galactic Analogs Survey" (xSAGA), a method for identifying low-$z$ galaxies on the basis of optical imaging, and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify $z < 0.03$ galaxies from more distant objects using image cutouts from the DESI Legacy Imaging Surveys. From the sample of $> 100,000$ CNN-selected low-$z$ galaxies, we identify $>20,000$ probable satellites located between 36-300 projected kpc from NASA-Sloan Atlas central galaxies in the stellar mass range $9.5 < \log(M_\star/M_\odot) < 11$. We characterize the incompleteness and contamination for CNN-selected samples, and apply corrections in order to estimate the true number of satellites as a function of projected radial distance from their hosts. Satellite richness depends strongly on host stellar mass, such that more massive host galaxies have more satellites, and on host morphology, such that elliptical hosts have more satellites than disky hosts with comparable stellar masses. We also find a strong inverse correlation between satellite richness and the magnitude gap between a host and its brightest satellite. The normalized satellite radial distribution between 36-300 kpc does not depend strongly on host stellar mass, morphology, or magnitude gap. The satellite abundances and radial distributions we measure are in reasonable agreement with predictions from hydrodynamic simulations. Our results deliver unprecedented statistical power for studying satellite galaxy populations, and highlight the promise of using machine learning for extending galaxy samples of wide-area surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Insights into the origin of halo mass profiles from machine learning

The mass distribution of dark matter haloes is the result of the hierarchical growth of initial density perturbations through mass accretion and mergers. We use an interpretable machine-learning framework to provide physical insights into the origin of the spherically-averaged mass profile of dark matter haloes. We train a gradient-boosted-trees algorithm to predict the final mass profiles of cluster-sized haloes, and measure the importance of the different inputs provided to the algorithm. We find two primary scales in the initial conditions (ICs) that impact the final mass profile: the density at approximately the scale of the haloes' Lagrangian patch $R_L$ ($R\sim 0.7\, R_L$) and that in the large-scale environment ($R\sim 1.7~R_L$). The model also identifies three primary time-scales in the halo assembly history that affect the final profile: (i) the formation time of the virialized, collapsed material inside the halo, (ii) the dynamical time, which captures the dynamically unrelaxed, infalling component of the halo over its first orbit, (iii) a third, most recent time-scale, which captures the impact on the outer profile of recent massive merger events. While the inner profile retains memory of the ICs, this information alone is insufficient to yield accurate predictions for the outer profile. As we add information about the haloes' mass accretion history, we find a significant improvement in the predicted profiles at all radii. Our machine-learning framework provides novel insights into the role of the ICs and the mass assembly history in determining the final mass profile of cluster-sized haloes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optical selection bias and projection effects in stacked galaxy cluster weak lensing

Cosmological constraints from current and upcoming galaxy cluster surveys are limited by the accuracy of cluster mass calibration. In particular, optically identified galaxy clusters are prone to selection effects that can bias the weak lensing mass calibration. We investigate the selection bias of the stacked cluster lensing signal associated with optically selected clusters, using clusters identified by the redMaPPer algorithm in the Buzzard simulations as a case study. We find that at a given cluster halo mass, the residuals of redMaPPer richness and weak lensing signal are positively correlated. As a result, for a given richness selection, the stacked lensing signal is biased high compared with what we would expect from the underlying halo mass probability distribution. The cluster lensing selection bias can thus lead to overestimated mean cluster mass and biased cosmology results. We show that the lensing selection bias exhibits a strong scale-dependence and is approximately 20 to 60 percent for $ΔΣ$ at large scales. This selection bias largely originates from spurious member galaxies within +/- 20 to 60 Mpc/h along the line of sight, highlighting the importance of quantifying projection effects associated with the broad redshift distribution of member galaxies in photometric cluster surveys. While our results qualitatively agree with those in the literature, accurate quantitative modelling of the selection bias is needed to achieve the goals of cluster lensing cosmology and will require synthetic catalogues covering a wide range of galaxy-halo connection models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Priors on red galaxy stochasticity from hybrid effective field theory

We investigate the stochastic properties of typical red galaxy samples in a controlled numerical environment. We use Halo Occupation Distribution (HOD) modelling to create mock realizations of three separate bright red galaxy samples consistent with datasets used for clustering and lensing analyses in modern galaxy surveys. Second-order Hybrid Effective Field Theory (HEFT) is used as a field-level forward model to describe the full statistical distribution of these tracer samples, and their stochastic power spectra are directly measured and compared to the Poisson shot-noise prediction. While all of the galaxy samples we consider are hosted within haloes with sub-Poisson stochasticity, we observe that the galaxy samples themselves possess stochasticities that range from sub-Poisson to super-Poisson, in agreement with predictions from the halo model. As an application of our methodology, we place priors on the expected degree of non-Poisson stochasticity in cosmological analyses using such samples. We expect these priors will be useful in reducing the complexity of the full parameter space for future analyses using second-order Lagrangian bias models. More generally, the techniques outlined here present the first application of hybrid EFT methods to characterize models of the galaxy--halo connection at the field level, revealing new connections between once-disparate modelling frameworks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021: Vera C. Rubin Observatory as a Flagship Dark Matter Experiment

Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requires astronomical observations, which still provide the only empirical evidence for dark matter to date. We have a unique opportunity right now to create a dark matter experiment with Rubin Observatory Legacy Survey of Space and Time (LSST). This experiment will be a coordinated effort to perform dark matter research, and provide a large collaborative team of scientists with the necessary organizational and funding supports. This approach leverages existing investments in Rubin. Studies of dark matter with Rubin LSST will also guide the design of, and confirm the results from, other dark matter experiments. Supporting a collaborative team to carry out a dark matter experiment with Rubin LSST is the key to achieving the dark matter science goals that have already been identified as high priority by the high-energy physics and astronomy communities.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stringent $σ_8$ constraints from small-scale galaxy clustering using a hybrid MCMC+emulator framework

We present a novel simulation-based hybrid emulator approach that maximally derives cosmological and Halo Occupation Distribution (HOD) information from non-linear galaxy clustering, with sufficient precision for DESI Year 1 (Y1) analysis. Our hybrid approach first samples the HOD space on a fixed cosmological simulation grid to constrain the high-likelihood region of cosmology+HOD parameter space, and then constructs the emulator within this constrained region. This approach significantly reduces the parameter volume emulated over, thus achieving much smaller emulator errors with fixed number of training points. We demonstrate that this combined with state-of-the-art simulations result in tight emulator errors comparable to expected DESI Y1 LRG sample variance. We leverage the new AbacusSummit simulations and apply our hybrid approach to CMASS non-linear galaxy clustering data. We infer constraints on $σ_8 = 0.762\pm0.024$ and $fσ_8 (z_{eff} = 0.52) = 0.444\pm0.016$, the tightest among contemporary galaxy clustering studies. We also demonstrate that our $fσ_8$ constraint is robust against secondary biases and other HOD model choices, a critical first step towards showcasing the robust cosmology information accessible in non-linear scales. We speculate that the additional statistical power of DESI Y1 should tighten the growth rate constraints by at least another 50-60%, significantly elucidating any potential tension with Planck. We also address the "lensing is low" tension, where we find that the combined effect of a lower $fσ_8$ and environment-based bias lowers the predicted lensing signal by 15%, accounting for approximately 50% of the discrepancy between the lensing measurement and clustering-based predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Supernova-driven turbulent metal mixing in high redshift galactic disks: metallicity fluctuations in the interstellar medium and its imprints on metal poor stars in the Milky Way

The extent to which turbulence mixes gas in the face of recurrent infusions of fresh metals by supernovae (SN) could help provide important constraints on the local star formation conditions. This includes predictions of the metallicity dispersion amongst metal poor stars, which suggests that the interstellar medium was not very well mixed at these early times. The purpose of this {\it Letter} is to help isolate, via a series of numerical experiments, some of the key processes that regulate turbulent mixing of SN elements in galactic disks. We study the gas interactions in small simulated patches of a galaxy disk with the goal of resolving the small-scale mixing effects of metals at pc scales, which enables us to measure the turbulent diffusion coefficient in various galaxy environments. By investigating the statistics of variations of $α$ elements in these simulations, we are able to derive constraints not only on the allowed range of intrinsic yield variations in SN explosions but also on the star formation history of the Milky Way. We argue that the observed dispersion of [Mg/Fe] in metal poor halo stars is compatible with the star-forming conditions expected in dwarf satellites or in an early low star-forming Milky Way progenitor. In particular, metal variations in stars that have not been phase-mixed can be used to infer the star-forming conditions of disrupted dwarf satellites.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Tidal disruption of solitons in self-interacting ultralight axion dark matter

Ultralight axions (ULAs) are promising dark matter candidates that can have a distinct impact on the formation and evolution of structure on nonlinear scales relative to the cold, collisionless dark matter (CDM) paradigm. However, most studies of structure formation in ULA models do not include the effects of self-interactions, which are expected to arise generically. Here, we study how the tidal evolution of solitons is affected by ULA self-interaction strength and sign. Specifically, using the pseudospectral solver UltraDark.jl, we simulate the tidal disruption of self-interacting solitonic cores as they orbit a $10^{11}~M_{\mathrm{\odot}}$ Navarro-Frenk-White CDM host halo potential for a range of orbital parameters, assuming a fiducial ULA particle mass of $10^{-22}\mathrm{eV}$. We find that repulsive (attractive) self-interactions significantly accelerate (decelerate) soliton tidal disruption. We also identify a degeneracy between the self-interaction strength and soliton mass that determines the efficiency of tidal disruption, such that disruption timescales are affected at the $\sim 50\%$ level for variations in the dimensionless ULA self-coupling from $λ=-10^{-92}$ to $λ=10^{-92}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

ADDGALS: Simulated Sky Catalogs for Wide Field Galaxy Surveys

We present a method for creating simulated galaxy catalogs with realistic galaxy luminosities, broad-band colors, and projected clustering over large cosmic volumes. The technique, denoted ADDGALS (Adding Density Dependent GAlaxies to Lightcone Simulations), uses an empirical approach to place galaxies within lightcone outputs of cosmological simulations. It can be applied to significantly lower-resolution simulations than those required for commonly used methods such as halo occupation distributions, subhalo abundance matching, and semi-analytic models, while still accurately reproducing projected galaxy clustering statistics down to scales of r ~ 100 kpc/h. We show that \addgals\ catalogs reproduce several statistical properties of the galaxy distribution as measured by the Sloan Digital Sky Survey (SDSS) main galaxy sample, including galaxy number densities, observed magnitude and color distributions, as well as luminosity- and color-dependent clustering. We also compare to cluster-galaxy cross correlations, where we find significant discrepancies with measurements from SDSS that are likely linked to artificial subhalo disruption in the simulations. Applications of this model to simulations of deep wide-area photometric surveys, including modeling weak-lensing statistics, photometric redshifts, and galaxy cluster finding are presented in DeRose et al (2019), and an application to a full cosmology analysis of Dark Energy Survey (DES) Year 3 like data is presented in DeRose etl al (2021). We plan to publicly release a 10,313 square degree catalog constructed using ADDGALS with magnitudes appropriate for several existing and planned surveys, including SDSS, DES, VISTA, WISE, and LSST.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}<-17$ galaxies, with more than $10^4$ at redshifts $z>7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Mock Lightcones and Theory Friendly Catalogs for the CANDELS Survey

We present mock catalogs created to support the interpretation of the CANDELS survey. We extract halos along past lightcones from the Bolshoi Planck dissipationless N-body simulations and populate these halos with galaxies using two different independently developed semi-analytic models of galaxy formation and the empirical model UniverseMachine. Our mock catalogs have geometries that encompass the footprints of observations associated with the five CANDELS fields. In order to allow field-to-field variance to be explored, we have created eight realizations of each field. In this paper, we present comparisons with observable global galaxy properties, including counts in observed frame bands, luminosity functions, color-magnitude distributions and color-color distributions. We additionally present comparisons with physical galaxy parameters derived from SED fitting for the CANDELS observations, such as stellar masses and star formation rates. We find relatively good agreement between the model predictions and CANDELS observations for luminosity and stellar mass functions. We find poorer agreement for colors and star formation rate distributions. All of the mock lightcones as well as curated "theory friendly" versions of the observational CANDELS catalogs are made available through a web-based data hub.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimization of the Observing Cadence for the Rubin Observatory Legacy Survey of Space and Time: a pioneering process of community-focused experimental design

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraints on Dark Matter Microphysics from the Milky Way Satellite Population

Alternatives to the cold, collisionless dark matter (DM) paradigm in which DM behaves as a collisional fluid generically suppress small-scale structure. Herein we use the observed population of Milky Way (MW) satellite galaxies to constrain the collisional nature of DM, focusing on DM-baryon scattering. We first derive analytic upper limits on the velocity-independent DM-baryon scattering cross section by translating the upper bound on the lowest mass of halos inferred to host satellites into a characteristic cutoff scale in the linear matter power spectrum. We then confirm and improve these results through a detailed probabilistic inference of the MW satellite population that marginalizes over relevant astrophysical uncertainties. This yields $95\%$ confidence upper limits on the DM-baryon scattering cross section of $2\times10^{-29}\ \rm{cm}^2$ ($6\times 10^{-27}\ \rm{cm}^2$) for DM particle masses $m_χ$ of~$10\ \rm{keV}$ ($10\ \rm{GeV}$); these limits scale as $m_χ^{1/4}$ for $m_χ\ll 1\ \rm{GeV}$ and $m_χ$ for~$m_χ\gg 1\ \rm{GeV}$. This analysis improves upon cosmological bounds derived from cosmic-microwave-background anisotropy measurements by multiple orders of magnitude over a wide range of DM masses, excluding regions of parameter space previously unexplored by other methods, including direct-detection experiments. Our work reveals a mapping between DM-baryon scattering and other alternative DM models, and we discuss the implications of our results for warm and fuzzy DM scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Signatures of Velocity-Dependent Dark Matter Self-Interactions in Milky Way-mass Halos

We explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, $w$, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of $w$ that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of $w$. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of $\sim 200\ \rm{km\ s}^{-1}$, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The SAGA Survey. II. Building a Statistical Sample of Satellite Systems around Milky Way-like Galaxies

We present the Stage II results from the ongoing Satellites Around Galactic Analogs (SAGA) Survey. Upon completion, the SAGA Survey will spectroscopically identify satellite galaxies brighter than $ M_{r,o} = -12.3 $ around 100 Milky Way (MW) analogs at $ z \sim 0.01 $. In Stage II, we have more than quadrupled the sample size of Stage I, delivering results from 127 satellites around 36 MW analogs with an improved target selection strategy and deep photometric imaging catalogs from the Dark Energy Survey and the Legacy Surveys. We have obtained 25,372 galaxy redshifts, peaking around $ z = 0.2 $. These data significantly increase spectroscopic coverage for very low redshift objects in $ 17 < r_o < 20.75 $ around SAGA hosts, creating a unique data set that places the Local Group in a wider context. The number of confirmed satellites per system ranges from zero to nine, and correlates with host galaxy and brightest satellite luminosities. We find that the number and the luminosities of MW satellites are consistent with being drawn from the same underlying distribution as SAGA systems. The majority of confirmed SAGA satellites are star forming, and the quenched fraction increases as satellite stellar mass and projected radius from the host galaxy decrease. Overall, the satellite quenched fraction among SAGA systems is lower than that in the Local Group. We compare the luminosity functions and radial distributions of SAGA satellites with theoretical predictions based on cold dark matter simulations and an empirical galaxy-halo connection model and find that the results are broadly in agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Constraining the scatter in the galaxy-halo connection at Milky Way masses

We develop and implement two new methods for constraining the scatter in the relationship between galaxies and dark matter halos. These new techniques are sensitive to the scatter at low halo masses, making them complementary to previous constraints that are dependent on clustering amplitudes or rich galaxy groups, both of which are only sensitive to more massive halos. In both of our methods, we use a galaxy group finder to locate central galaxies in the SDSS main galaxy sample. Our first technique uses the small-scale cross-correlation of central galaxies with all lower mass galaxies. This quantity is sensitive to the satellite fraction of low-mass galaxies, which is in turn driven by the scatter between halos and galaxies. The second technique uses the kurtosis of the distribution of line-of-sight velocities between central galaxies and neighboring galaxies. This quantity is sensitive to the distribution of halo masses that contain the central galaxies at fixed stellar mass. Theoretical models are constructed using peak halo circular velocity, $V_{\rm peak}$, as our property to connect galaxies to halos. The cross-correlation technique yields a constraint of $σ[ M_\ast|V_{\rm peak}]=0.27\pm 0.05$ dex, corresponding to a scatter in $\log M_\ast$ at fixed $M_h$ of $σ[ M_\ast|M_h]=0.38\pm 0.06$ dex at $M_h=10^{11.8}$ Msun. The kurtosis technique yields $σ[ M_\ast|V_{\rm peak}]=0.30\pm0.03$, corresponding to $σ[ M_\ast|M_h]=0.34\pm 0.04$ at $M_h=10^{12.2}$ Msun. The values of $σ[ M_\ast|M_h]$ are significantly larger than the constraints at higher masses, in agreement with the results of hydrodynamic simulations. This increase is only partly due to the scatter between $V_{\rm peak}$ and $M_h$, and it represents an increase of nearly a factor of two relative to the values inferred from clustering and group studies at high masses.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Impact of the Fiducial Cosmology Assumption on BAO Cosmological Parameter Inference

Standard analysis pipelines for measurements of Baryon Acoustic Oscillations (BAO) in galaxy surveys make use of a fiducial cosmological model to guide the data compression required to transform from observed redshifts and angles to the measured angular and radial BAO peak positions. In order to remove any dependence on the fiducial cosmology from the results, all models compared to the data should mimic the compression and its dependence on the fiducial model. In practice, approximations are made when testing models: (1) There is assumed to be no residual dependence on the fiducial cosmology after reconstruction, (2) differences in the distance--redshift relationship are assumed to match a linear scaling, and (3) differences in clustering between true and fiducial models are assumed to be removed by the free parameters used to null the non-BAO signal. We test these approximations using the current standard measurement procedure with a set of halo catalogs from the {\sc Aemulus} suite of $N$-body simulations, which span a range of $w\mathrm{CDM}$ cosmological models. We focus on reconstruction of the primordial BAO and locating the BAO. For the range of $w\mathrm{CDM}$ cosmologies covered by the {\sc Aemulus} suite, we find no evidence for systematic errors in the measured BAO shift parameters $α_{\parallel}$ and $α_{\bot}$ to $< 0.1\%$. However, the measured errors $σ_{α_{\parallel}}$ and $σ_{α_{\bot}}$ show a notable absolute increase by up to $+0.001$ and $+0.002$ respectively in the case that the fiducial cosmology does not match the truth. These effects on the inferred BAO scale will be important given the precision of measurements expected from future surveys including DESI, Euclid, and WFIRST.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Biased Tracers in Redshift Space in the EFT of Large-Scale Structure

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Connecting CO Intensity Mapping to Molecular Gas and Star Formation in the Epoch of Galaxy Assembly

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraining the Mass-Richness Relationship of redMaPPer Clusters with Angular Clustering

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detection of lensing substructure using ALMA observations of the dusty galaxy SDP.81

We study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations, we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a $M=10^{8.96\pm 0.12} M_{\odot}$ subhalo near one of the images, with a significance of $6.9σ$ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter subhalos down to $M\sim 2\times 10^7 M_{\odot}$, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted dark matter subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 dataset (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of $Λ$CDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Disentangling redshift-space distortions and nonlinear bias using the 2D power spectrum

We present the nonlinear 2D galaxy power spectrum, $P(k,μ)$, in redshift space, measured from the Dark Sky simulations, using galaxy catalogs constructed with both halo occupation distribution and subhalo abundance matching methods, chosen to represent an intermediate redshift sample of luminous red galaxies. We find that the information content in individual $μ$ (cosine of the angle to the line of sight) bins is substantially richer then multipole moments, and show that this can be used to isolate the impact of nonlinear growth and redshift space distortion (RSD) effects. Using the $μ<0.2$ simulation data, which we show is not impacted by RSD effects, we can successfully measure the nonlinear bias to an accuracy of $\sim 5$% at $k<0.6 h$Mpc$^{-1}$. This use of individual $μ$ bins to extract the nonlinear bias successfully removes a large parameter degeneracy when constraining the linear growth rate of structure. We carry out a joint parameter estimation, using the low $μ$ simulation data to constrain the nonlinear bias, and $μ\ge0.2$ to constrain the growth rate and show that $f$ can be constrained to $\sim 26\, (22)$% to a $k_{\rm max}< 0.4\, (0.6) h$Mpc$^{-1}$ from clustering alone using a simple dispersion model, for a range of galaxy models. Our analysis of individual $μ$ bins also reveals interesting physical effects which arise simply from different methods of populating halos with galaxies. We find a prominent turnaround scale, at which RSD damping effects are greater then the nonlinear growth, which differs not only for each $μ$ bin but also for each galaxy model. These features may provide unique signatures which could be used to shed light on the galaxy-dark matter connection.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy Cluster Mass Estimation from Stacked Spectroscopic Analysis

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_σ(λ)$, derived from stacked pairwise spectroscopy of clusters with richness $λ$. 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})|λ\rangle = \ln(M_{30}) + α_m \ln(λ/30)$ with $M_{30} = 1.56 \pm 0.35 \times 10^{14} M_\odot$ and $α_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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

RHAPSODY-G simulations II - Baryonic growth and metal enrichment in massive galaxy clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Concentration Dependence of the Galaxy-Halo Connection: Modeling Assembly Bias with Abundance Matching

Empirical methods for connecting galaxies to their dark matter halos have become essential for interpreting measurements of the spatial statistics of galaxies. In this work, we present a novel approach for parameterizing the degree of concentration dependence in the abundance matching method. This new parameterization provides a smooth interpolation between two commonly used matching proxies: the peak halo mass and the peak halo maximal circular velocity. This parameterization controls the amount of dependence of galaxy luminosity on halo concentration at a fixed halo mass. Effectively this interpolation scheme enables abundance matching models to have adjustable assembly bias in the resulting galaxy catalogs. With the new 400 Mpc/h DarkSky Simulation, whose larger volume provides lower sample variance, we further show that low-redshift two-point clustering and satellite fraction measurements from SDSS can already provide a joint constraint on this concentration dependence and the scatter within the abundance matching framework.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

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-$α$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The DESI Experiment Part II: Instrument Design

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= λ/Δλ$ 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Eating Habits of Milky Way Mass Halos: Destroyed Dwarf Satellites and the Metallicity Distribution of Accreted Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Faber-Jackson relation and Fundamental Plane from halo abundance matching

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Critical Assessment of Stellar Mass Measurement Methods

In this paper we perform a comprehensive study of the main sources of random and systematic errors in stellar mass measurement for galaxies using their Spectral Energy Distributions (SEDs). We use mock galaxy catalogs with simulated multi-waveband photometry (from U-band to mid-infrared) and known redshift, stellar mass, age and extinction for individual galaxies. Given different parameters affecting stellar mass measurement (photometric S/N ratios, SED fitting errors, systematic effects, the inherent degeneracies and correlated errors), we formulated different simulated galaxy catalogs to quantify these effects individually. We studied the sensitivity of stellar mass estimates to the codes/methods used, population synthesis models, star formation histories, nebular emission line contributions, photometric uncertainties, extinction and age. For each simulated galaxy, the difference between the input stellar masses and those estimated using different simulation catalogs, $Δ\log(M)$, was calculated and used to identify the most fundamental parameters affecting stellar masses. We measured different components of the error budget, with the results listed as follows: (1). no significant bias was found among different codes/methods, with all having comparable scatter; (2). A source of error is found to be due to photometric uncertainties and low resolution in age and extinction grids; (3). The median of stellar masses among different methods provides a stable measure of the mass associated with any given galaxy; (4). The deviations in stellar mass strongly correlate with those in age, with a weaker correlation with extinction; (5). the scatter in the stellar masses due to free parameters are quantified, with the sensitivity of the stellar mass to both the population synthesis codes and inclusion of nebular emission lines studied.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discovery of two gravitationally lensed quasars in the Dark Energy Survey

We present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rhapsody-G simulations: galaxy clusters as baryonic closed boxes and the covariance between hot gas and galaxies

Within a sufficiently large cosmic volume, conservation of baryons implies a simple `closed box' view in which the sum of the baryonic components must equal a constant fraction of the total enclosed mass. We present evidence from Rhapsody-G hydrodynamic simulations of massive galaxy clusters that the closed-box expectation may hold to a surprising degree within the interior, non-linear regions of haloes. At a fixed halo mass, we find a significant anti-correlation between hot gas mass fraction and galaxy mass fraction (cold gas + stars), with a rank correlation coefficient of -0.69 within $R_{500c}$. Because of this anti-correlation, the total baryon mass serves as a low-scatter proxy for total cluster mass. The fractional scatter of total baryon fraction scales approximately as $0.02 (Δ_c/100)^{0.6}$, while the scatter of either gas mass or stellar mass is larger in magnitude and declines more slowly with increasing radius. We discuss potential observational tests using cluster samples selected by optical and hot gas properties; the simulations suggest that joint selection on stellar and hot gas has potential to achieve 5% scatter in total halo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Dependence of Subhalo Abundance on Halo Concentration

Hierarchical structure formation implies that the number of subhalos within a dark matter halo depends not only on halo mass, but also on the formation history of the halo. This dependence on the formation history, which is highly correlated with halo concentration, can account for the super-Poissonian scatter in subhalo occupation at a fixed halo mass that has been previously measured in simulations. Here we propose a model to predict the subhalo abundance function for individual host halos, that incorporates both halo mass and concentration. We combine results of cosmological simulations with a new suite of zoom-in simulations of Milky Way-mass halos to calibrate our model. We show the model can successfully reproduce the mean and the scatter of subhalo occupation in these simulations. The implications of this correlation between subhalo abundance and halo concentration are further investigated. We also discuss cases in which inferences about halo properties can be affected if this correlation between subhalo abundance and halo concentration is ignored; in these cases our model would give a more accurate inference. We propose that with future deep surveys, satellite occupation in the low-mass regime can be used to verify the existence of halo assembly bias.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Relation Between SFR and Stellar Mass for Galaxies at 3.5 $\le z\le$ 6.5 in CANDELS

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $σ(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Tully-Fisher and mass-size relations from halo abundance matching

The Tully-Fisher relation (TFR) expresses the connection between rotating galaxies and the dark matter haloes they inhabit, and therefore contains a wealth of information about galaxy formation. We construct a general framework to investigate whether models based on halo abundance matching are able to reproduce the observed stellar mass TFR and mass-size relation (MSR), and use the data to constrain galaxy formation parameters. Our model tests a range of plausible scenarios, differing in the response of haloes to disc formation, the relative angular momentum of baryons and dark matter, the impact of selection effects, and the abundance matching parameters. We show that agreement with the observed TFR puts an upper limit on the scatter between galaxy and halo properties, requires weak or reversed halo contraction, and favours selection effects that preferentially eliminate fast-rotating galaxies. The MSR constrains the ratio of the disc to halo specific angular momentum to be approximately in the range 0.6-1.2. We identify and quantify two problems that models of this nature face. (1) They predict too large an intrinsic scatter for the MSR, and (2) they predict too strong an anticorrelation between the TFR and MSR residuals. We argue that resolving these problems requires introducing a correlation between stellar surface density and enclosed dark matter mass. Finally, we explore the expected difference between the TFRs of central and satellite galaxies, finding that in the favoured models this difference should be detectable in a sample of ~700 galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Connecting Direct Dark Matter Detection Experiments to Cosmologically Motivated Halo Models

Several direct detection experiments, including recently CDMS-II, have reported signals consistent with 5 to 10 GeV dark matter (DM) that appear to be in tension with null results from XENON and LUX experiments; these indicate a careful review of the theoretical basis, including the galactic DM velocity distribution function (VDF). We establish a VDF parameter space from DM-only cosmological simulations and illustrate that seemingly contradictory experimental results can be made consistent within this parameter space. Future experimental limits should be reported after they are marginalized over a range of VDF parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dark Sky Simulations: Early Data Release

The Dark Sky Simulations are an ongoing series of cosmological N-body simulations designed to provide a quantitative and accessible model of the evolution of the large-scale Universe. Such models are essential for many aspects of the study of dark matter and dark energy, since we lack a sufficiently accurate analytic model of non-linear gravitational clustering. In July 2014, we made available to the general community our early data release, consisting of over 55 Terabytes of simulation data products, including our largest simulation to date, which used $1.07 \times 10^{12}~(10240^3)$ particles in a volume $8h^{-1}\mathrm{Gpc}$ across. Our simulations were performed with 2HOT, a purely tree-based adaptive N-body method, running on 200,000 processors of the Titan supercomputer, with data analysis enabled by yt. We provide an overview of the derived halo catalogs, mass function, power spectra and light cone data. We show self-consistency in the mass function and mass power spectrum at the 1% level over a range of more than 1000 in particle mass. We also present a novel method to distribute and access very large datasets, based on an abstraction of the World Wide Web (WWW) as a file system, remote memory-mapped file access semantics, and a space-filling curve index. This method has been implemented for our data release, and provides a means to not only query stored results such as halo catalogs, but also to design and deploy new analysis techniques on large distributed datasets.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Formation of disk galaxies in preheated media: a preventative feedback model

We introduce a semi-analytic galaxy formation model implementing a self-consistent treatment for the hot halo gas configuration and the assembly of central disks. Using the model, we explore a preventative feedback model, in which the circum-halo medium is assumed to be preheated up to a certain entropy level by early starbursts or other processes, and compare it with an ejective feedback model, in which baryons are first accreted into dark matter halos and subsequently ejected out by feedback. The model demonstrates that when the medium is preheated to an entropy comparable to the halo virial entropy the baryon accretion can be largely reduced and delayed. In addition, the preheated medium can establish an extended low density gaseous halo when it accretes into the dark matter halos, and result in a specific angular momentum of the cooling gas large enough to form central disks as extended as those observed. Combined with simulated halo assembly histories, the preventative feedback model can reproduce remarkably well a number of observational scaling relations. These include the cold baryon (stellar plus cold gas) mass fraction-halo mass relations, star formation histories, disk size-stellar mass relation and its evolution, and the number density of low-mass galaxies as a function of redshift. In contrast, the conventional ejective feedback model fails to reproduce these observational trends. Using the model, we demonstrate that the properties of disk galaxies are closely tied to the thermal state of hot halo gas and even possibly the circum-halo medium, which suggests that observational data for the disk properties and circum-galactic hot/warm medium may jointly provide interesting constraints for galaxy formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

redMaPPer IV: Photometric Membership Identification of Cluster Galaxies with 1% Precision

In order to study the galaxy population of galaxy clusters with photometric data one must be able to accurately discriminate between cluster members and non-members. The redMaPPer cluster finding algorithm treats this problem probabilistically. Here, we utilize SDSS and GAMA spectroscopic membership rates to validate the redMaPPer membership probability estimates for clusters with $z\in[0.1,0.3]$. We find small - but correctable - biases, sourced by three different systematics. The first two were expected a priori, namely blue cluster galaxies and correlated structure along the line of sight. The third systematic is new: the redMaPPer template fitting exhibits a non-trivial dependence on photometric noise, which biases the original redMaPPer probabilities when utilizing noisy data. After correcting for these effects, we find exquisite agreement ($\approx 1\%$) between the photometric probability estimates and the spectroscopic membership rates, demonstrating that we can robustly recover cluster membership estimates from photometric data alone. As a byproduct of our analysis we find that on average unavoidable projection effects from correlated structure contribute $\approx 6\%$ of the richness of a redMaPPer galaxy cluster. This work also marks the second public release of the SDSS redMaPPer cluster catalog.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spectroscopic failures in photometric redshift calibration: cosmological biases and survey requirements

We use N-body-spectro-photometric simulations to investigate the impact of incompleteness and incorrect redshifts in spectroscopic surveys to photometric redshift training and calibration and the resulting effects on cosmological parameter estimation from weak lensing shear-shear correlations. The photometry of the simulations is modeled after the upcoming Dark Energy Survey and the spectroscopy is based on a low/intermediate resolution spectrograph with wavelength coverage of 5500Å < λ < 9500Å. The principal systematic errors that such a spectroscopic follow-up encounters are incompleteness (inability to obtain spectroscopic redshifts for certain galaxies) and wrong redshifts. Encouragingly, we find that a neural network-based approach can effectively describe the spectroscopic incompleteness in terms of the galaxies' colors, so that the spectroscopic selection can be applied to the photometric sample. Hence, we find that spectroscopic incompleteness yields no appreciable biases to cosmology, although the statistical constraints degrade somewhat because the photometric survey has to be culled to match the spectroscopic selection. Unfortunately, wrong redshifts have a more severe impact: the cosmological biases are intolerable if more than a percent of the spectroscopic redshifts are incorrect. Moreover, we find that incorrect redshifts can also substantially degrade the accuracy of training set based photo-z estimators. The main problem is the difficulty of obtaining redshifts, either spectroscopically or photometrically, for objects at z > 1.3. We discuss several approaches for reducing the cosmological biases, in particular finding that photo-z error estimators can reduce biases appreciably.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Comparison between Semi-Analytic Model Predictions for the CANDELS Survey

We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

CANDELS Observations of the Environmental Dependence of the Color-Mass-Morphology Relation at z = 1.6

We study the environmental dependence of color, stellar mass, and morphology by comparing galaxies in a forming cluster to those in the field at z = 1:6 with Hubble Space Telescope near-infrared imaging in the CANDELS/UDS field. We quantify the morphology of the galaxies using the effective radius, reff, and Sérsic index, n. In both the cluster and field, approximately half of the bulge-dominated galaxies (n > 2) reside on the red sequence of the color-magnitude diagram, and most disk-dominated galaxies (n < 2) have colors expected for star-forming galaxies. There is weak evidence that cluster galaxies have redder rest-frame U - B colors and higher stellar masses compared to the field. Star-forming galaxies in both the cluster and field show no significant differences in their morphologies. In contrast, there is evidence that quiescent galaxies in the cluster have larger median effective radii and smaller Sérsic indices compared to the field with a significance of 2?. These differences are most pronounced for galaxies at clustercentric distances 1 Mpc < Rproj < 1.5 Mpc, which have low Sérsic indices and possibly larger effective radii, more consistent with star-forming galaxies at this epoch and in contrast to other quiescent galaxies. We argue that star-forming galaxies are processed under the influence of the cluster environment at distances greater than the cluster-halo virial radius. Our results are consistent with models where gas accretion onto these galaxies is suppressed from processes associated with the cluster environment.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Galaxies in X-ray Groups. III. Satellite Color and Morphology Transformations

While the star formation rates and morphologies of galaxies have long been known to correlate with their local environment, the process by which these correlations are generated is not well understood. Galaxy groups are thought to play an important role in shaping the physical properties of galaxies before entering massive clusters at low redshift, and transformations of satellite galaxies likely dominate the buildup of local environmental correlations. To illuminate the physical processes that shape galaxy evolution in dense environments, we study a sample of 116 X-ray selected galaxy groups at z=0.2-1 with halo masses of 10^13-10^14 M_sun and centroids determined with weak lensing. We analyze morphologies based on HST imaging and colors determined from 31 photometric bands for a stellar mass-limited population of 923 satellite galaxies and a comparison sample of 16644 field galaxies. Controlling for variations in stellar mass across environments, we find significant trends in the colors and morphologies of satellite galaxies with group-centric distance and across cosmic time. Specifically at low stellar mass (log(M_stellar/M_sun) = 9.8-10.3), the fraction of disk-dominated star-forming galaxies declines from >50% among field galaxies to <20% among satellites near the centers of groups. This decline is accompanied by a rise in quenched galaxies with intermediate bulge+disk morphologies, and only a weak increase in red bulge-dominated systems. These results show that both color and morphology are influenced by a galaxy's location within a group halo. We suggest that strangulation and disk fading alone are insufficient to explain the observed morphological dependence on environment, and that galaxy mergers or close tidal encounters must play a role in building up the population of quenched galaxies with bulges seen in dense environments at low redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitationally Consistent Halo Catalogs and Merger Trees for Precision Cosmology

We present a new algorithm for generating merger trees and halo catalogs which explicitly ensures consistency of halo properties (mass, position, and velocity) across timesteps. Our algorithm has demonstrated the ability to improve both the completeness (through detecting and inserting otherwise missing halos) and purity (through detecting and removing spurious objects) of both merger trees and halo catalogs. In addition, our method is able to robustly measure the self-consistency of halo finders; it is the first to directly measure the uncertainties in halo positions, halo velocities, and the halo mass function for a given halo finder based on consistency between snapshots in cosmological simulations. We use this algorithm to generate merger trees for two large simulations (Bolshoi and Consuelo) and evaluate two halo finders (ROCKSTAR and BDM). We find that both the ROCKSTAR and BDM halo finders track halos extremely well; in both, the number of halos which do not have physically consistent progenitors is at the 1-2% level across all halo masses. Our code is publicly available at http://code.google.com/p/consistent-trees . Our trees and catalogs are publicly available at http://hipacc.ucsc.edu/Bolshoi/ .

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Improved Mock Galaxy Catalogs for the DEEP2 Galaxy Redshift Survey from Subhalo Abundance and Environment Matching

We develop empirical methods for modeling the galaxy population and populating cosmological N-body simulations with mock galaxies according to the observed properties of galaxies in survey data. We use these techniques to produce a new set of mock catalogs for the DEEP2 Galaxy Redshift Survey based on the output of the high-resolution Bolshoi simulation, as well as two other simulations with different cosmological parameters, all of which we release for public use. The mock-catalog creation technique uses subhalo abundance matching to assign galaxy luminosities to simulated dark-matter halos. It then adds color information to the resulting mock galaxies in a manner that depends on the local galaxy density, in order to reproduce the measured color-environment relation in the data. In the course of constructing the catalogs, we test various models for including scatter in the relation between halo mass and galaxy luminosity, within the abundance-matching framework. We find that there is no constant-scatter model that can simultaneously reproduce both the luminosity function and the autocorrelation function of DEEP2. This result has implications for galaxy-formation theory, and it restricts the range of contexts in which the mocks can be usefully applied. Nevertheless, careful comparisons show that our new mocks accurately reproduce a wide range of the other properties of the DEEP2 catalog, suggesting that they can be used to gain a detailed understanding of various selection effects in DEEP2.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mergers and Mass Accretion for Infalling Halos Both End Well Outside Cluster Virial Radii

We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reionization Histories of Milky Way Mass Halos

We investigate the connection between the epoch of reionization and the present day universe, by examining the extended mass reionization histories of dark matter halos identified at z=0. We combine an N-body dark matter simulation of a 600 Mpc volume with a three-dimensional, seminumerical reionization model. This provides reionization redshifts for each particle, which can then be connected with the properties of their halos at the present time. We find that the vast majority of present-day halos with masses larger than ~ few x 10^11 Msun reionize earlier than the rest of the universe. We also find significant halo-to-halo diversity in mass reionization histories, and find that in realistic inhomogenous models, the material within a given halo is not expected to reionize at the same time. In particular, the scatter in reionization times within individual halos is typically larger than the scatter among halos. From our fiducial reionization model, we find that the typical 68% scatter in reionization times within halos is ~ 115 Myr for 10^(12 \pm 0.25) Msun halos, decreasing slightly to ~ 95 Myr for 10^(15 \pm 0.25) Msun halos. We find a mild correlation between reionization history and environment: halos with shorter reionization histories are typically in more clustered environments, with the strongest trend on a scale of ~ 20 Mpc. Material in Milky Way mass halos with short reionization histories is preferentially reionized in relatively large HII regions, implying reionization mostly by sources external to the progenitors of the present-day halo. We investigate the impact on our results of varying the reionization model parameters, which span a range of reionization scenarios with varying timing and morphology.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rhapsody. I. Structural Properties and Formation History From a Statistical Sample of Re-simulated Cluster-size Halos

We present the first results from the Rhapsody cluster re-simulation project: a sample of 96 "zoom-in" simulations of dark matter halos of 10^14.8 +- 0.05 Msun/h, selected from a 1 (Gpc/h)^3 volume. This simulation suite is the first to resolve this many halos with ~5x10^6 particles per halo in the cluster-mass regime, allowing us to statistically characterize the distribution of and correlation between halo properties at fixed mass. We focus on the properties of the main halos and how they are affected by formation history, which we track back to z=12, over five decades in mass. We give particular attention to the impact of the formation history on the density profiles of the halos. We find that the deviations from the Navarro-Frenk-White (NFW) model and the Einasto model depend on formation time. Late-forming halos tend to have considerable deviations from both models, partly due to the presence of massive subhalos, while early-forming halos deviate less but still significantly from the NFW model and are better described by the Einasto model. We find that the halo shapes depend only moderately on formation time. Departure from spherical symmetry impacts the density profiles through the anisotropic distribution of massive subhalos. Further evidence of the impact of subhalos is provided by analyzing the phase-space structure. A detailed analysis of the properties of the subhalo population in Rhapsody is presented in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rhapsody. II. Subhalo Properties and the Impact of Tidal Stripping From a Statistical Sample of Cluster-Size Halos

We discuss the properties of subhalos in cluster-size halos, using a high-resolution statistical sample: the Rhapsody simulations introduced in Wu et al. (2012). We demonstrate that the criteria applied to select subhalos have significant impact on the inferred properties of the sample, including the scatter in the number of subhalos, the correlation between the subhalo number and formation time, and the shape of subhalos' spatial distribution and velocity structure. We find that the number of subhalos, when selected using the peak maximum circular velocity in their histories (a property expected to be closely related to the galaxy luminosity), is uncorrelated with the formation time of the main halo. This is in contrast to the previously reported correlation from studies where subhalos are selected by the current maximum circular velocity; we show that this difference is a result of the tidal stripping of the subhalos. We also find that the dominance of the main halo and the subhalo mass fraction are strongly correlated with halo concentration and formation history. These correlations are important to take into account when interpreting results from cluster samples selected with different criteria. Our sample also includes a fossil cluster, which is presented separately and placed in the context of the rest of the sample.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spatially unassociated galaxies contribute significantly to the blended submillimetre galaxy population: predictions for follow-up observations of ALMA sources

There is anecdotal evidence that spatially and physically unassociated galaxies blended into a single submillimetre (submm) source contribute to the submm galaxy (SMG) population. This work is the first to theoretically predict the number counts of such sources. We generate mock SMG catalogues using lightcones derived from the Bolshoi cosmological simulation; to assign submm flux densities to the mock galaxies, we use a fitting function previously derived from the results of dust radiative transfer performed on hydrodynamical simulations of isolated disc and merging galaxies. We then calculate submm number counts for different beam sizes and without blending. We predict that > ~50 per cent of blended SMGs have at least one spatially unassociated component with S_850 > 1 mJy. For a 15-arcsec beam, blends of >2 galaxies in which at least one component is spatially unassociated dominate the blended sources with total S_850 > ~3 mJy. The distribution of the redshift separations amongst the components is strongly bimodal. The typical redshift separation of spatially unassociated blended sources is ~1. Our predictions for the contributions of spatially unassociated components and the distribution of redshift separations are not testable with currently available data, but they will be easily tested once sufficiently accurate redshifts for the individual subcomponents (resolved by, e.g., ALMA) of a sufficient number of single-dish-detected blended SMGs are available.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Average Star Formation Histories of Galaxies in Dark Matter Halos from z=0-8

We present a robust method to constrain average galaxy star formation rates, star formation histories, and the intracluster light as a function of halo mass. Our results are consistent with observed galaxy stellar mass functions, specific star formation rates, and cosmic star formation rates from z=0 to z=8. We consider the effects of a wide range of uncertainties on our results, including those affecting stellar masses, star formation rates, and the halo mass function at the heart of our analysis. As they are relevant to our method, we also present new calibrations of the dark matter halo mass function, halo mass accretion histories, and halo-subhalo merger rates out to z=8. We also provide new compilations of cosmic and specific star formation rates; more recent measurements are now consistent with the buildup of the cosmic stellar mass density at all redshifts. Implications of our work include: halos near 10^12 Msun are the most efficient at forming stars at all redshifts, the baryon conversion efficiency of massive halos drops markedly after z ~ 2.5 (consistent with theories of cold-mode accretion), the ICL for massive galaxies is expected to be significant out to at least z ~ 1-1.5, and dwarf galaxies at low redshifts have higher stellar mass to halo mass ratios than previous expectations and form later than in most theoretical models. Finally, we provide new fitting formulae for star formation histories that are more accurate than the standard declining tau model. Our approach places a wide variety of observations relating to the star formation history of galaxies into a self-consistent framework based on the modern understanding of structure formation in LCDM. Constraints on the stellar mass-halo mass relationship and star formation rates are available for download at http://www.peterbehroozi.com/data.html .

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Connection between Galaxies and Dark Matter Structures in the Local Universe

We provide new constraints on the connection between galaxies in the local universe, identified by the Sloan Digital Sky Survey (SDSS), and dark matter halos and their constituent substructures in the $Λ$CDM model using WMAP7 cosmological parameters. Predictions for the abundance and clustering properties of dark matter halos, and the relationship between dark matter hosts and substructures, are based on a high-resolution cosmological simulation, the Bolshoi simulation. We associate galaxies with halos and subhalos using subhalo abundance matching, performing a comprehensive analysis which investigates the underlying assumptions of this technique including (a) which halo property is most closely associated with galaxy stellar masses and luminosities, (b) how much scatter is in this relationship, and (c) how much subhalos can be stripped before their galaxies are destroyed. The models are jointly constrained by new measurements of the projected two-point galaxy clustering and the observed conditional stellar mass function of galaxies in groups. The data put tight constraints on the satellite fraction of galaxies as a function of galaxy stellar mass, on the scatter between halo and galaxy properties, and on the underlying conditional stellar mass function. These data rule out several halo properties commonly used in abundance matching, largely because the satellite fractions in the models disagree with those data. We show that an abundance matching model that associates galaxies with the peak circular velocity of their halos is in good agreement with the data, when scatter of $0.20 \pm 0.03$ dex in stellar mass at a given peak velocity is included. This will yield important constraints for galaxy formation models, and also provides encouraging indications that the galaxy--halo connection can be modeled with sufficient fidelity for future precision studies of the dark Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Rockstar Phase-Space Temporal Halo Finder and the Velocity Offsets of Cluster Cores

We present a new algorithm for identifying dark matter halos, substructure, and tidal features. The approach is based on adaptive hierarchical refinement of friends-of-friends groups in six phase-space dimensions and one time dimension, which allows for robust (grid-independent, shape-independent, and noise-resilient) tracking of substructure; as such, it is named Rockstar (Robust Overdensity Calculation using K-Space Topologically Adaptive Refinement). Our method is massively parallel (up to 10^5 CPUs) and runs on the largest current simulations (>10^10 particles) with high efficiency (10 CPU hours and 60 gigabytes of memory required per billion particles analyzed). A previous paper (Knebe et al 2011) has shown Rockstar to have class-leading recovery of halo properties; we expand on these comparisons with more tests and higher-resolution simulations. We show a significant improvement in substructure recovery as compared to several other halo finders and discuss the theoretical and practical limits of simulations in this regard. Finally, we present results which demonstrate conclusively that dark matter halo cores are not at rest relative to the halo bulk or satellite average velocities and have coherent velocity offsets across a wide range of halo masses and redshifts. For massive clusters, these offsets can be up to 350 km/s at z=0 and even higher at high redshifts. Our implementation is publicly available at http://code.google.com/p/rockstar .

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unbound Particles in Dark Matter Halos

We investigate unbound dark matter particles in halos by tracing particle trajectories in a simulation run to the far future (a = 100). We find that the traditional sum of kinetic and potential energies is a very poor predictor of which dark matter particles will eventually become unbound from halos. We also study the mass fraction of unbound particles, which increases strongly towards the edges of halos, and decreases significantly at higher redshifts. We discuss implications for dark matter detection experiments, precision calibrations of the halo mass function, the use of baryon fractions to constrain dark energy, and searches for intergalactic supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Virial scaling of galaxies in clusters: bright to faint is cool to hot

By combining galaxy tracers from high-resolution N-body and hydrodynamical simulations, we present a consistent picture of the behaviour of galaxy velocities in massive clusters. In haloes above ~ 10^14 Msun, the brightest satellite galaxies are slightly cooler compared to the dark matter, while fainter satellites are hotter. Within the virial radius of a cluster, the mean velocity dispersion based on the 100 brightest galaxies is a factor of 1.065 +/- 0.005 (stat) +/- 0.027 (sys) higher than that of the dark matter (corresponding to a ~10-15 per cent bias in the dynamical mass estimate) while that based on only the five brightest galaxies is 0.868 +/- 0.039 (stat) +/- 0.035 (sys). These trends are approximately independent of redshift. The velocity structure is sensitive to the modelling of galaxies in clusters, indicative of the complex interplay of tidal stripping, dynamical friction, and merging. Velocity dispersions derived from instantaneous subhalo properties are larger than those employing either peak subhalo properties or hydrodynamical galaxy tracers. The latter two methods are consistent, implying that stacked spectroscopic analysis of cluster samples should, after correction for projection, show a trend towards slightly higher velocities when fainter galaxies are included, with an unbiased measure of dark matter velocity dispersion coming from approximately 30 galaxies per cluster. We show evidence that the velocity distribution function of bright galaxies near the cluster centre has a low-velocity tail due to strong dynamical friction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A High Throughput Workflow Environment for Cosmological Simulations

The next generation of wide-area sky surveys offer the power to place extremely precise constraints on cosmological parameters and to test the source of cosmic acceleration. These observational programs will employ multiple techniques based on a variety of statistical signatures of galaxies and large-scale structure. These techniques have sources of systematic error that need to be understood at the percent-level in order to fully leverage the power of next-generation catalogs. Simulations of large-scale structure provide the means to characterize these uncertainties. We are using XSEDE resources to produce multiple synthetic sky surveys of galaxies and large-scale structure in support of science analysis for the Dark Energy Survey. In order to scale up our production to the level of fifty 10^10-particle simulations, we are working to embed production control within the Apache Airavata workflow environment. We explain our methods and report how the workflow has reduced production time by 40% compared to manual management.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Galaxies in X-ray Groups. II. A Weak Lensing Study of Halo Centering

Locating the centers of dark matter halos is critical for understanding the mass profiles of halos as well as the formation and evolution of the massive galaxies that they host. The task is observationally challenging because we cannot observe halos directly, and tracers such as bright galaxies or X-ray emission from hot plasma are imperfect. In this paper we quantify the consequences of miscentering on the weak lensing signal from a sample of 129 X-ray selected galaxy groups in the COSMOS field with redshifts 0<z<1 and halo masses in the range 10^13 - 10^14 M_sun. By measuring the stacked lensing signal around eight different candidate centers (such as the brightest member galaxy, the mean position of all member galaxies, or the X-ray centroid), we determine which candidates best trace the center of mass in halos. In this sample of groups, we find that massive galaxies near the X-ray centroids trace the center of mass to <~75 kpc, while the X-ray position and centroids based on the mean position of member galaxies have larger offsets primarily due to the statistical uncertainties in their positions (typically ~50-150 kpc). Approximately 30% of groups in our sample have ambiguous centers with multiple bright or massive galaxies, and these groups show disturbed mass profiles that are not well fit by standard models, suggesting that they are merging systems. We find that halo mass estimates from stacked weak lensing can be biased low by 5-30% if inaccurate centers are used and the issue of miscentering is not addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Halo-to-Halo Similarity and Scatter in the Velocity Distribution of Dark Matter

We examine the Velocity Distribution Function (VDF) in dark matter halos from Milky Way to cluster mass scales. We identify an empirical model for the VDF with a wider peak and a steeper tail than a Maxwell--Boltzmann distribution, and discuss physical explanations. We quantify sources of scatter in the VDF of cosmological halos and their implication for direct detection of dark matter. Given modern simulations and observations, we find that the most significant uncertainty in the VDF of the Milky Way arises from the unknown radial position of the solar system relative to the dark matter halo scale radius.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Lack of Evolution in Galaxy Star Formation Efficiency

Using reconstructed galaxy star formation histories, we calculate the instantaneous efficiency of galaxy star formation (i.e., the star formation rate divided by the baryon accretion rate) from $z=8$ to the present day. This efficiency exhibits a clear peak near a characteristic halo mass of 10^11.7 Msun, which coincides with longstanding theoretical predictions for the mass scale relevant to virial shock heating of accreted gas. Above the characteristic halo mass, the efficiency falls off as the mass to the minus four-thirds power; below the characteristic mass, the efficiency falls off at an average scaling of mass to the two-thirds power. By comparison, the shape and normalization of the efficiency change very little since z=4. We show that a time-independent star formation efficiency simply explains the shape of the cosmic star formation rate since z=4 in terms of dark matter accretion rates. The rise in the cosmic star formation from early times until z=2 is especially sensitive to galaxy formation efficiency. The mass dependence of the efficiency strongly limits where most star formation occurs, with the result that two-thirds of all star formation has occurred inside halos within a factor of three of the characteristic mass, a range that includes the mass of the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sample variance in photometric redshift calibration: cosmological biases and survey requirements

We use N-body/photometric galaxy simulations to examine the impact of sample variance of spectroscopic redshift samples on the accuracy of photometric redshift (photo-z) determination and calibration of photo-z errors. We estimate the biases in the cosmological parameter constraints from weak lensing and derive requirements on the spectroscopic follow-up for three different photo-z algorithms chosen to broadly span the range of algorithms available. We find that sample variance is much more relevant for the photo-z error calibration than for photo-z training, implying that follow-up requirements are similar for different algorithms. We demonstrate that the spectroscopic sample can be used for training of photo-zs and error calibration without incurring additional bias in the cosmological parameters. We provide a guide for observing proposals for the spectroscopic follow-up to ensure that redshift calibration biases do not dominate the cosmological parameter error budget. For example, assuming optimistically (pessimistically) that the weak lensing shear measurements from the Dark Energy Survey could obtain 1-sigma constraints on the dark energy equation of state w of 0.035 (0.055), implies a follow-up requirement of 150 (40) patches of sky with a telescope such as Magellan, assuming a 1/8^2 deg effective field of view and 400 galaxies per patch. Assuming (optimistically) a VVDS-like spectroscopic completeness with purely random failures, this could be accomplished with about 75 (20) nights of observation. For more realistic assumptions regarding spectroscopic completeness, or in the presence of other sources of systematics not considered here, further degradations to dark energy constraints are possible. We test several approaches for reducing the requirements. Abridged

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Correlated Formation Histories of Massive Galaxies and Their Dark Matter Halos

Using observations in the COSMOS field, we report an intriguing correlation between the star formation activity of massive (~10^{11.4}\msol) central galaxies, their stellar masses, and the large-scale (~10 Mpc) environments of their group-mass (~10^{13.6}\msol) dark matter halos. Probing the redshift range z=[0.2,1.0], our measurements come from two independent sources: an X-ray detected group catalog and constraints on the stellar-to-halo mass relation derived from a combination of clustering and weak lensing statistics. At z=1, we find that the stellar mass in star-forming centrals is a factor of two less than in passive centrals at the same halo mass. This implies that the presence or lack of star formation in group-scale centrals cannot be a stochastic process. By z=0, the offset reverses, probably as a result of the different growth rates of these objects. A similar but weaker trend is observed when dividing the sample by morphology rather than star formation. Remarkably, we find that star-forming centrals at z~1 live in groups that are significantly more clustered on 10 Mpc scales than similar mass groups hosting passive centrals. We discuss this signal in the context of halo assembly and recent simulations, suggesting that star-forming centrals prefer halos with higher angular momentum and/or formation histories with more recent growth; such halos are known to evolve in denser large-scale environments. If confirmed, this would be evidence of an early established link between the assembly history of halos on large scales and the future properties of the galaxies that form inside them.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Cosmic Abundance of Classical Milky Way Satellites

We study the abundance of satellites akin to the brightest, classical dwarf spheroidals around galaxies similar in magnitude and isolation to the Milky Way and M31 in the Sloan Digital Sky Survey. From a combination of photometric and spectroscopic redshifts, we bound the mean and the intrinsic scatter in the number of satellites down to ten magnitudes fainter than the Milky Way. Restricting to magnitudes brighter than Sagittarius, we show that the Milky Way is not a significant statistical outlier in its population of classical dwarf spheroidals. At fainter magnitudes, we find an upper limit of 13 on the mean number of satellites brighter than the Fornax dwarf spheroidal. Methods to improve these limits that utilize full photometric redshift distributions hold promise, but are currently limited by incompleteness at the very lowest redshifts. Theoretical models are left to explain why the majority of dark matter subhalos that orbit Milky Way-like galaxies are inefficient at making galaxies at the luminosity scale of the brightest dwarf spheroidals, or why these subhalos predicted by Lambda-CDM do not exist.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5σ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5σ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cosmological Constraints from Galaxy Clustering and the Mass-to-Number Ratio of Galaxy Clusters

We place constraints on the average density (Omega_m) and clustering amplitude (sigma_8) of matter using a combination of two measurements from the Sloan Digital Sky Survey: the galaxy two-point correlation function, w_p, and the mass-to-galaxy-number ratio within galaxy clusters, M/N, analogous to cluster M/L ratios. Our w_p measurements are obtained from DR7 while the sample of clusters is the maxBCG sample, with cluster masses derived from weak gravitational lensing. We construct non-linear galaxy bias models using the Halo Occupation Distribution (HOD) to fit both w_p and M/N for different cosmological parameters. HOD models that match the same two-point clustering predict different numbers of galaxies in massive halos when Omega_m or sigma_8 is varied, thereby breaking the degeneracy between cosmology and bias. We demonstrate that this technique yields constraints that are consistent and competitive with current results from cluster abundance studies, even though this technique does not use abundance information. Using w_p and M/N alone, we find Omega_m^0.5*sigma_8=0.465+/-0.026, with individual constraints of Omega_m=0.29+/-0.03 and sigma_8=0.85+/-0.06. Combined with current CMB data, these constraints are Omega_m=0.290+/-0.016 and sigma_8=0.826+/-0.020. All errors are 1-sigma. The systematic uncertainties that the M/N technique are most sensitive to are the amplitude of the bias function of dark matter halos and the possibility of redshift evolution between the SDSS Main sample and the maxBCG sample. Our derived constraints are insensitive to the current level of uncertainties in the halo mass function and in the mass-richness relation of clusters and its scatter, making the M/N technique complementary to cluster abundances as a method for constraining cosmology with future galaxy surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Galaxies in X-ray Groups I: Robust Membership Assignment and the Impact of Group Environments on Quenching

Understanding the mechanisms that lead dense environments to host galaxies with redder colors, more spheroidal morphologies, and lower star formation rates than field populations remains an important problem. As most candidate processes ultimately depend on host halo mass, accurate characterizations of the local environment, ideally tied to halo mass estimates and spanning a range in halo mass and redshift are needed. In this work, we present and test a rigorous, probabalistic method for assigning galaxies to groups based on precise photometric redshifts and X-ray selected groups drawn from the COSMOS field. The groups have masses in the range 10^13 < M_200c/M_sun < 10^14 and span redshifts 0<z<1. We characterize our selection algorithm via tests on spectroscopic subsamples, including new data obtained at the VLT, and by applying our method to detailed mock catalogs. We find that our group member galaxy sample has a purity of 84% and completeness of 92% within 0.5 R200c. We measure the impact of uncertainties in redshifts and group centering on the quality of the member selection with simulations based on current data as well as future imaging and spectroscopic surveys. As a first application of our new group member catalog which will be made publicly available, we show that member galaxies exhibit a higher quenched fraction compared to the field at fixed stellar mass out to z~1, indicating a significant relationship between star formation and environment at group scales. We also address the suggestion that dusty star forming galaxies in such groups may impact the high-l power spectrum of the cosmic microwave background and find that such a population cannot explain the low power seen in recent SZ measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

How Common are the Magellanic Clouds?

We introduce a probabilistic approach to the problem of counting dwarf satellites around host galaxies in databases with limited redshift information. This technique is used to investigate the occurrence of satellites with luminosities similar to the Magellanic Clouds around hosts with properties similar to the Milky Way in the object catalog of the Sloan Digital Sky Survey. Our analysis uses data from SDSS Data Release 7, selecting candidate Milky-Way-like hosts from the spectroscopic catalog and candidate analogs of the Magellanic Clouds from the photometric catalog. Our principal result is the probability for a Milky-Way-like galaxy to host N_{sat} close satellites with luminosities similar to the Magellanic Clouds. We find that 81 percent of galaxies like the Milky Way are have no such satellites within a radius of 150 kpc, 11 percent have one, and only 3.5 percent of hosts have two. The probabilities are robust to changes in host and satellite selection criteria, background-estimation technique, and survey depth. These results demonstrate that the Milky Way has significantly more satellites than a typical galaxy of its luminosity; this fact is useful for understanding the larger cosmological context of our home galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Statistics of Satellite Galaxies Around Milky Way-Like Hosts

We calculate the probability that a Milky-Way-like halo in the standard cosmological model has the observed number of Magellanic Clouds (MCs). The statistics of the number of MCs in the LCDM model are in good agreement with observations of a large sample of SDSS galaxies. Under the sub-halo abundance matching assumption of a relationship with small scatter between galaxy r-band luminosities and halo internal velocities v_max, we make detailed comparisons to similar measurements using SDSS DR7 data by Liu et al. (2010). Models and observational data give very similar probabilities for having zero, one, and two MC-like satellites. In both cases, Milky Way-luminosity hosts have just a \sim 10% chance of hosting two satellites similar to the Magellanic Clouds. In addition, we present a prediction for the probability for a host galaxy to have Nsats satellite galaxies as a function of the magnitudes of both the host and satellite. This probability and its scaling with host properties is significantly different from that of mass-selected objects because of scatter in the mass- luminosity relation and because of variations in the star formation efficiency with halo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The integrated stellar content of dark matter halos

Measurements of the total amount of stars locked up in galaxies as a function of host halo mass contain key clues about the efficiency of processes that regulate star formation. We derive the total stellar mass fraction f_star as a function of halo mass M500c from z=0.2 to z=1 using two complementary methods. First, we derive f_star using a statistical Halo Occupation Distribution model jointly constrained by data from lensing, clustering, and the stellar mass function. This method enables us to probe f_star over a much wider halo mass range than with group or cluster catalogs. Second, we derive f_star at group scales using a COSMOS X-ray group catalog and we show that the two methods agree to within 30%. We quantify the systematic uncertainty on f_star using abundance matching methods and we show that the statistical uncertainty on f_star (~10%) is dwarfed by systematic uncertainties associated with stellar mass measurements (~45% excluding IMF uncertainties). Assuming a Chabrier IMF, we find 0.012<f_star<0.025 at M500c=10^13 Msun and 0.0057<f_star<0.015 at M500c=10^14 Msun. These values are significantly lower than previously published estimates. We investigate the cause of this difference and find that previous work has overestimated f_star due to a combination of inaccurate stellar mass estimators and/or because they have assumed that all galaxies in groups are early type galaxies with a constant mass-to-light ratio. Contrary to previous claims, our results suggest that the mean value of f_star is always significantly lower than f_gas for halos above 10^13 Msun. Combining our results with recently published gas mas fractions, we find a shortfall in f_star+f_gas at R500c compared to the cosmic mean. This shortfall varies with halo mass and becomes larger towards lower halos masses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Mass Distribution and Assembly of the Milky Way from the Properties of the Magellanic Clouds

We present a new measurement of the mass of the Milky Way (MW) based on observed properties of its largest satellite galaxies, the Magellanic Clouds (MCs), and an assumed prior of a ΛCDM universe. The large, high-resolution Bolshoi cosmological simulation of this universe provides a means to statistically sample the dynamical properties of bright satellite galaxies in a large population of dark matter halos. The observed properties of the MCs, including their circular velocity, distance from the center of the MW, and velocity within the MW halo, are used to evaluate the likelihood that a given halo would have each or all of these properties; the posterior PDF for any property of the MW system can thus be constructed. This method provides a constraint on the MW virial mass, 1.2 +0.7 -0.3(stat.) +0.3 -0.4 (sys.) x 10^12 M\odot (68% confidence), which is consistent with recent determinations that involve very different assumptions. In addition, we calculate the posterior PDF for the density profile of the MW and its satellite accretion history. Although typical satellites of 10^12 M\odot halos are accreted over a wide range of epochs over the last 10 Gyr, we find a \sim72% probability that the Magellanic Clouds were accreted within the last Gyr, and a 50% probability that they were accreted together.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Comprehensive Analysis of Uncertainties Affecting the Stellar Mass - Halo Mass Relation for 0<z<4

We conduct a comprehensive analysis of the relationship between central galaxies and their host dark matter halos, as characterized by the stellar mass-halo mass (SM-HM) relation, with rigorous consideration of uncertainties. Our analysis focuses on results from the abundance matching technique, which assumes that every dark matter halo or subhalo above a specific mass threshold hosts one galaxy. We discuss the quantitative effects of uncertainties in observed galaxy stellar mass functions (GSMFs) (including stellar mass estimates and counting uncertainties), halo mass functions (including cosmology and uncertainties from substructure), and the abundance matching technique used to link galaxies to halos (including scatter in this connection). Our analysis results in a robust estimate of the SM-HM relation and its evolution from z=0 to z=4. The shape and evolution are well constrained for z < 1. The largest uncertainties at these redshifts are due to stellar mass estimates; however, failure to account for scatter in stellar masses at fixed halo mass can lead to errors of similar magnitude in the SM-HM relation for central galaxies in massive halos. We also investigate the SM-HM relation to z=4, although the shape of the relation at higher redshifts remains fairly unconstrained when uncertainties are taken into account. These results will provide a powerful tool to inform galaxy evolution models. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A GMBCG Galaxy Cluster Catalog of 55,424 Rich Clusters from SDSS DR7

We present a large catalog of optically selected galaxy clusters from the application of a new Gaussian Mixture Brightest Cluster Galaxy (GMBCG) algorithm to SDSS Data Release 7 data. The algorithm detects clusters by identifying the red sequence plus Brightest Cluster Galaxy (BCG) feature, which is unique for galaxy clusters and does not exist among field galaxies. Red sequence clustering in color space is detected using an Error Corrected Gaussian Mixture Model. We run GMBCG on 8240 square degrees of photometric data from SDSS DR7 to assemble the largest ever optical galaxy cluster catalog, consisting of over 55,000 rich clusters across the redshift range from 0.1 < z < 0.55. We present Monte Carlo tests of completeness and purity and perform cross-matching with X-ray clusters and with the maxBCG sample at low redshift. These tests indicate high completeness and purity across the full redshift range for clusters with 15 or more members.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Annealing a Follow-up Program: Improvement of the Dark Energy Figure of Merit for Optical Galaxy Cluster Surveys

The precision of cosmological parameters derived from galaxy cluster surveys is limited by uncertainty in relating observable signals to cluster mass. We demonstrate that a small mass-calibration follow-up program can significantly reduce this uncertainty and improve parameter constraints, particularly when the follow-up targets are judiciously chosen. To this end, we apply a simulated annealing algorithm to maximize the dark energy information at fixed observational cost, and find that optimal follow-up strategies can reduce the observational cost required to achieve a specified precision by up to an order of magnitude. Considering clusters selected from optical imaging in the Dark Energy Survey, we find that approximately 200 low-redshift X-ray clusters or massive Sunyaev-Zel'dovich clusters can improve the dark energy figure of merit by 50%, provided that the follow-up mass measurements involve no systematic error. In practice, the actual improvement depends on (1) the uncertainty in the systematic error in follow-up mass measurements, which needs to be controlled at the 5% level to avoid severe degradation of the results; and (2) the scatter in the optical richness-mass distribution, which needs to be made as tight as possible to improve the efficacy of follow-up observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Counts-in-Cylinders in the Sloan Digital Sky Survey with Comparisons to N-body Simulations

Environmental statistics provide a necessary means of comparing the properties of galaxies in different environments and a vital test of models of galaxy formation within the prevailing, hierarchical cosmological model. We explore counts-in-cylinders, a common statistic defined as the number of companions of a particular galaxy found within a given projected radius and redshift interval. Galaxy distributions with the same two-point correlation functions do not necessarily have the same companion count distributions. We use this statistic to examine the environments of galaxies in the Sloan Digital Sky Survey, Data Release 4. We also make preliminary comparisons to four models for the spatial distributions of galaxies, based on N-body simulations, and data from SDSS DR4 to study the utility of the counts-in-cylinders statistic. There is a very large scatter between the number of companions a galaxy has and the mass of its parent dark matter halo and the halo occupation, limiting the utility of this statistic for certain kinds of environmental studies. We also show that prevalent, empirical models of galaxy clustering that match observed two- and three-point clustering statistics well fail to reproduce some aspects of the observed distribution of counts-in-cylinders on 1, 3 and 6-Mpc/h scales. All models that we explore underpredict the fraction of galaxies with few or no companions in 3 and 6-Mpc/h cylinders. Roughly 7% of galaxies in the real universe are significantly more isolated within a 6 Mpc/h cylinder than the galaxies in any of the models we use. Simple, phenomenological models that map galaxies to dark matter halos fail to reproduce high-order clustering statistics in low-density environments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precision Measurements of the Cluster Red Sequence using an Error Corrected Gaussian Mixture Model

The red sequence is an important feature of galaxy clusters and plays a crucial role in optical cluster detection. Measurement of the slope and scatter of the red sequence are affected both by selection of red sequence galaxies and measurement errors. In this paper, we describe a new error corrected Gaussian Mixture Model for red sequence galaxy identification. Using this technique, we can remove the effects of measurement error and extract unbiased information about the intrinsic properties of the red sequence. We use this method to select red sequence galaxies in each of the 13,823 clusters in the maxBCG catalog, and measure the red sequence ridgeline location and scatter of each. These measurements provide precise constraints on the variation of the average red galaxy populations in the observed frame with redshift. We find that the scatter of the red sequence ridgeline increases mildly with redshift, and that the slope decreases with redshift. We also observe that the slope does not strongly depend on cluster richness. Using similar methods, we show that this behavior is mirrored in a spectroscopic sample of field galaxies, further emphasizing that ridgeline properties are independent of environment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Dark Matter at the End of the Galaxy

Dark matter density profiles based upon Lambda-CDM cosmology motivate an ansatz velocity distribution function with fewer high velocity particles than the Maxwell-Boltzmann distribution or proposed variants. The high velocity tail of the distribution is determined by the outer slope of the dark matter halo, the large radius behavior of the Galactic dark matter density. N-body simulations of Galactic halos reproduce the high velocity behavior of this ansatz. Predictions for direct detection rates are dramatically affected for models where the threshold scattering velocity is within 30% of the escape velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Impact of Theoretical Uncertainties in the Halo Mass Function and Halo Bias on Precision Cosmology

We study the impact of theoretical uncertainty in the dark matter halo mass function and halo bias on dark energy constraints from imminent galaxy cluster surveys. We find that for an optical cluster survey like the Dark Energy Survey, the accuracy required on the predicted halo mass function to make it an insignificant source of error on dark energy parameters is ~ 1%. The analogous requirement on the predicted halo bias is less stringent (~ 5%), particularly if the observable-mass distribution can be well constrained by other means. These requirements depend upon survey area but are relatively insensitive to survey depth. The most stringent requirements are likely to come from a survey over a significant fraction of the sky that aims to observe clusters down to relatively low mass, Mth ~ 10^13.7 Msun/h; for such a survey, the mass function and halo bias must be predicted to accuracies of ~ 0.5% and ~ 1%, respectively. These accuracies represent a limit on the practical need to calibrate ever more accurate halo mass and bias functions. We find that improving predictions for the mass function in the low-redshift and low-mass regimes is the most effective way to improve dark energy constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Voronoi Tessellation cluster finder in 2+1 dimensions

We present a detailed description of the Voronoi Tessellation (VT) cluster finder algorithm in 2+1 dimensions, which improves on past implementations of this technique. The need for cluster finder algorithms able to produce reliable cluster catalogs up to redshift 1 or beyond and down to $10^{13.5}$ solar masses is paramount especially in light of upcoming surveys aiming at cosmological constraints from galaxy cluster number counts. We build the VT in photometric redshift shells and use the two-point correlation function of the galaxies in the field to both determine the density threshold for detection of cluster candidates and to establish their significance. This allows us to detect clusters in a self consistent way without any assumptions about their astrophysical properties. We apply the VT to mock catalogs which extend to redshift 1.4 reproducing the $Λ$CDM cosmology and the clustering properties observed in the SDSS data. An objective estimate of the cluster selection function in terms of the completeness and purity as a function of mass and redshift is as important as having a reliable cluster finder. We measure these quantities by matching the VT cluster catalog with the mock truth table. We show that the VT can produce a cluster catalog with completeness and purity $>80%$ for the redshift range up to $\sim 1$ and mass range down to $\sim 10^{13.5}$ solar masses.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

ArborZ: Photometric Redshifts Using Boosted Decision Trees

Precision photometric redshifts will be essential for extracting cosmological parameters from the next generation of wide-area imaging surveys. In this paper we introduce a photometric redshift algorithm, ArborZ, based on the machine-learning technique of Boosted Decision Trees. We study the algorithm using galaxies from the Sloan Digital Sky Survey and from mock catalogs intended to simulate both the SDSS and the upcoming Dark Energy Survey. We show that it improves upon the performance of existing algorithms. Moreover, the method naturally leads to the reconstruction of a full probability density function (PDF) for the photometric redshift of each galaxy, not merely a single "best estimate" and error, and also provides a photo-z quality figure-of-merit for each galaxy that can be used to reject outliers. We show that the stacked PDFs yield a more accurate reconstruction of the redshift distribution N(z). We discuss limitations of the current algorithm and ideas for future work.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Connecting Galaxies, Halos, and Star Formation Rates Across Cosmic Time

A simple, observationally-motivated model is presented for understanding how halo masses, galaxy stellar masses, and star formation rates are related, and how these relations evolve with time. The relation between halo mass and galaxy stellar mass is determined by matching the observed spatial abundance of galaxies to the expected spatial abundance of halos at multiple epochs -- i.e. more massive galaxies are assigned to more massive halos at each epoch. Halos at different epochs are connected by halo mass accretion histories estimated from N-body simulations. The halo--galaxy connection at fixed epochs in conjunction with the connection between halos across time provides a connection between observed galaxies across time. With approximations for the impact of merging and accretion on the growth of galaxies, one can then directly infer the star formation histories of galaxies as a function of stellar and halo mass. This model is tuned to match both the observed evolution of the stellar mass function and the normalization of the observed star formation rate -- stellar mass relation to z~1. The data demands, for example, that the star formation rate density is dominated by galaxies with Mstar~10^{10.0-10.5} Msun from 0<z<1, and that such galaxies over these epochs reside in halos with M~10^{11.5-12.5} Msun. The star formation rate -- halo mass relation is approximately Gaussian over the range 0<z<1 with a mildly evolving mean and normalization. This model is then used to shed light on a number of issues, including 1) a clarification of ``downsizing'', 2) the lack of a sharp characteristic halo mass at which star formation is truncated, and 3) the dominance of star formation over merging to the stellar build-up of galaxies with Mstar<10^11 Msun at z<1.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cosmological Constraints from the SDSS maxBCG Cluster Catalog

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog to simultaneously constrain cosmology and the richness--mass relation of the clusters. Assuming a flat \LambdaCDM cosmology, we find σ_8(Ω_m/0.25)^{0.41} = 0.832\pm 0.033 after marginalization over all systematics. In common with previous studies, our error budget is dominated by systematic uncertainties, the primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and uncertainty in the scatter of the richness--mass relation. Our constraints are fully consistent with the WMAP five-year data, and in a joint analysis we find σ_8=0.807\pm 0.020 and Ω_m=0.265\pm 0.016, an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster abundances. The remarkable consistency among these results demonstrates that cluster abundance constraints are not only tight but also robust, and highlight the power of optically-selected cluster samples to produce precision constraints on cosmological parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Interpreting the Clustering of Distant Red Galaxies

We analyze the angular clustering of z~2.3 distant red galaxies (DRGs) measured by Quadri et al 2008. We find that, with robust estimates of the measurement errors and realistic halo occupation distribution modeling, the measured clustering can be well fit within standard halo occupation models, in contrast to previous results. However, in order to fit the strong break in w(theta) at theta=10 arcsec, nearly all satellite galaxies in the DRG luminosity range are required to be DRGs. Within this luminosity-threshold sample, the fraction of galaxies that are DRGs is ~44%, implying that the formation of DRGs is more efficient for satellite galaxies than for central galaxies. Despite the evolved stellar populations contained within DRGs at z=2.3, 90% of satellite galaxies in the DRG luminosity range have been accreted within 500 Myr. Thus, satellite DRGs must have known they would become satellites well before the time of their accretion. This implies that the formation of DRGs correlates with large-scale environment at fixed halo mass, although the large-scale bias of DRGs can be well fit without such assumptions. Further data are required to resolve this issue. Using the observational estimate that ~30% of DRGs have no ongoing star formation, we infer a timescale for star formation quenching for satellite galaxies of 450 Myr, although the uncertainty on this number is large. However, unless all non-star forming satellite DRGs were quenched before accretion, the quenching timescale is significantly shorter than z~0 estimates. Down to the completeness limit of the Quadri et al sample, we find that the halo masses of central DRGs are ~50% higher than non-DRGs in the same luminosity range, but at the highest halo masses the central galaxies are DRGs only ~2/3 of the time.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Impact of Inhomogeneous Reionization on the Satellite Galaxy Population of the Milky Way

We use the publicly available subhalo catalogs from the Via Lactea simulation along with a Gpc-scale N-body simulation to understand the impact of inhomogeneous reionization on the satellite galaxy population of the Milky Way. The large-volume simulation is combined with a model for reionization that allows us to predict the distribution of reionization times for Milky Way mass halos. Motivated by this distribution, we identify candidate satellite galaxies in the simulation by requiring that any subhalo must grow above a specified mass threshold before it is reionized; after this time the photoionizing background will suppress both the formation of stars and the accretion of gas. We show that varying the reionization time over the range expected for Milky Way mass halos can change the number of satellite galaxies by roughly two orders of magnitude. This conclusion is in contradiction with a number of studies in the literature, and we conclude that this is a result of inconsistent application of the results of Gnedin (2000). We compare our satellite galaxies to observations using both abundance matching and stellar population synthesis methods to assign luminosities to our subhalos and account for observational completeness effects. Additionally, if we assume that the mass threshold is set by the virial temperature Tvir = 8e3K we find that our model accurately matches the vmax distribution, radial distribution, and luminosity function of observed Milky Way satellites for a reionization time zreion = 9.6^{1.0}_{-2.1}, assuming that the Via Lacteasubhalo distribution is representative of the Milky Way. This results in the presence of 119^{+202}_{-50} satellite galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Merger Histories of Galaxy Halos and Implications for Disk Survival

We study the merger histories of galaxy dark matter halos using a high resolution LCDM N-body simulation. Our merger trees follow ~17,000 halos with masses M_0 = (10^11--10^13) Msun at z=0 and track accretion events involving objects as small as m = 10^10 Msun. We find that mass assembly is remarkably self-similar in m/M_0, and dominated by mergers that are ~10% of the final halo mass. While very large mergers, m > 0.4 M_0, are quite rare, sizeable accretion events, m ~ 0.1 M_0, are common. Over the last 10 Gyr, an overwhelming majority (~95%) of Milky Way-sized halos with M_0 = 10^12 Msun have accreted at least one object with greater total mass than the Milky Way disk (m > 5x10^10 Msun), and approximately 70% have accreted an object with more than twice that mass (m > 10^11 Msun). Our results raise serious concerns about the survival of thin-disk dominated galaxies within the current paradigm for galaxy formation in a CDM universe. In order to achieve a ~70% disk-dominated fraction in Milky Way-sized CDM halos, mergers involving m ~ 2x10^11 Msun objects must not destroy disks. Considering that most thick disks and bulges contain old stellar populations, the situation is even more restrictive: these mergers must not heat disks or drive gas into their centers to create young bulges.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Isolating Triggered Star Formation

Galaxy pairs provide a potentially powerful means of studying triggered star formation from galaxy interactions. We use a large cosmological N-body simulation coupled with a well-tested semi-analytic substructure model to demonstrate that the majority of galaxies in close pairs reside within cluster or group-size halos and therefore represent a biased population, poorly suited for direct comparison to ``field'' galaxies. Thus, the frequent observation that some types of galaxies in pairs have redder colors than ``field'' galaxies is primarily a selection effect. We select galaxy pairs that are isolated in their dark matter halos with respect to other massive subhalos (N=2 halos) and a control sample of isolated galaxies (N=1 halos) for comparison. We then apply these selection criteria to a volume-limited subset of the 2dF Galaxy Redshift Survey with M_Bj <= -19 and obtain the first clean measure of the typical fraction of galaxies affected by triggered star formation and the average elevation in the star formation rate. We find that 24% (30.5%) of these L^\star and sub-L^{\star} galaxies in isolated 50 (30) kpc/h pairs exhibit star formation that is boosted by a factor of >~ 5 above their average past value, while only 10% of isolated galaxies in the control sample show this level of enhancement. Thus, 14% (20 %) of the galaxies in these close pairs show clear triggered star formation. The isolation criteria we develop provide a means to constrain star formation and feedback prescriptions in hydrodynamic simulations and a very general method of understanding the importance of triggered star formation in a cosmological context. (Abridged.)

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

MaxBCG: A Red Sequence Galaxy Cluster Finder

Measurements of galaxy cluster abundances, clustering properties, and mass to- light ratios in current and future surveys can provide important cosmological constraints. Digital wide-field imaging surveys, the recently-demonstrated fidelity of red-sequence cluster detection techniques, and a new generation of realistic mock galaxy surveys provide the means for construction of large, cosmologicallyinteresting cluster samples, whose selection and properties can be understood in unprecedented depth. We present the details of the "maxBCG" algorithm, a cluster-detection technique tailored to multi-band CCD-imaging data. MaxBCG primarily relies on an observational cornerstone of massive galaxy clusters: they are marked by an overdensity of bright, uniformly red galaxies. This detection scheme also exploits classical brightest cluster galaxies (BCGs), which are often found at the center of these same massive clusters. (ABRIDGED)

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Close Galaxy Counts as a Probe of Hierarchical Structure Formation

Standard LCDM predicts that the major merger rate of galaxy-size dark matter halos rises rapidly with redshift. The average number of close companions per galaxy, Nc, is often used to infer the galaxy merger rate, however, recent observational studies suggest that Nc evolves very little with redshift. Here we use a "hybrid" N- body simulation plus analytic substructure model to predict Nc directly. We identify dark matter subhalos with galaxies and show that the observed lack of close pair count evolution arises because the high merger rate per halo at early times is counteracted by a decrease in the number of halos massive enough to host a galaxy pair. We compare our results to data compiled from the DEEP2, SSRS2, and the UZC redshift surveys. Observed pair counts match our predictions if we assume a monotonic mapping between galaxy luminosity and the maximum circular velocity that each subhalo had when it was first accreted onto its host halo. This suggests that satellite galaxies are significantly more resilient to destruction than are dissipationless dark matter subhalos. We argue that while Nc does not provide a direct measure of the halo merger rate, it offers a powerful means to constrain the Halo Occupation Distribution and the spatial distribution of galaxies within halos. Interpreted in this way, close pair counts provide a useful test of galaxy formation processes on < 100 kpc scales.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Galaxy-galaxy Lensing: Dissipationless Simulations Versus the Halo Model

Galaxy-galaxy lensing is a powerful probe of the relation between galaxies and dark matter halos, but its theoretical interpretation requires a careful modeling of various contributions, such as the contribution from central and satellite galaxies. For this purpose, a phenomenological approach based on the halo model has been developed, allowing for fast exploration of the parameter space of models. In this paper, we investigate the ability of the halo model to extract information from the g-g weak lensing signal by comparing it to high-resolution dissipationless simulations that resolve subhalos. We find that the halo model reliably determines parameters such as the host halo mass of central galaxies, the fraction of galaxies that are satellites, and their radial distribution inside larger halos. If there is a significant scatter present in the central galaxy host halo mass distribution, then the mean and median mass of that distribution can differ significantly from one another, and the halo model mass determination lies between the two. This result suggests that when analyzing the data, galaxy subsamples with a narrow central galaxy halo mass distribution, such as those based on stellar mass, should be chosen for a simpler interpretation of the results.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Measurement of Galaxy Cluster Sizes, Radial Profiles, and Luminosity Functions from SDSS Photometric Data

Imaging data from the Sloan Digital Sky Survey is used to measure the empirical size-richness relation for a large sample of galaxy clusters. Using population subtraction methods, we determine the radius at which the cluster galaxy number density is 200/Omega_m times the mean galaxy density, without assuming a model for the radial distribution of galaxies in clusters. If these galaxies are unbiased on Mpc scales, this galaxy-density-based R_200 reflects the characteristic radii of clusters. We measure the scaling of this characteristic radius with richness over an order of magnitude in cluster richness, from rich clusters to poor groups. We use this information to examine the radial profiles of galaxies in clusters as a function of cluster richness, finding that the concentration of the galaxy distribution decreases with richness and is systematically lower than the concentrations measured for dark matter profiles in N-body simulations. Using these scaled radii, we investigate the behavior of the cluster luminosity function, and find that it is well matched by a Schechter function for galaxies brighter than M_r = -18 only after the central galaxy has been removed. We find that the luminosity function varies with richness and with distance from the cluster center, underscoring the importance of using an aperture that scales with cluster mass to compare physically equivalent regions of these different systems. We note that the lowest richness systems in our catalog have properties consistent with those expected of the earliest-forming halos; our cluster-finding algorithm, in addition to reliably finding clusters, may be efficient at finding fossil groups.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Modeling Galaxy-Mass Correlations in Dissipationless Simulations

We use high-resolution dissipationless simulations of the concordance flat LCDM model to make predictions for the galaxy--mass (GM) correlations and compare them to the recent SDSS weak lensing measurements.We use a simple observationally motivated scheme to assign luminosities and colors to the halos.This allows us to closely match the selection criteria used to define observational samples.The simulations reproduce the observed GM correlation function and its observed dependencies on luminosity and color.The luminosity dependence of the correlation function is primarily determined by the changing relative contribution of central and satellite galaxies at different luminosities. The color dependence of the GM correlations reflects the difference in the typical environments of blue and red galaxies. We also find agreement between the predicted and observed cross-bias, b_x=b/r,at all probed scales.The GM correlation coefficient, r, is close to unity on scales >1/h Mpc.The cross bias is thus expected to measure the actual bias of galaxy clustering on these scales.The aperture mass-to-light ratio is independent of galaxy color.The aperture mass scales approximately linearly with luminosity at L_r>10^{10}h^{-2} Lsun, while at lower luminosities the scaling is shallower: L_r^{0.5}. We show that most of the luminous galaxies (M_r<-21) are near the centers of their halos and their GM correlation function at r<100/h kpc can therefore be interpreted as the average dark matter density profile of these galaxies. We find that for galaxies in a given narrow luminosity range, there is a broad and possibly non-gaussian distribution of halo virial masses. Therefore, the average relation between mass and luminosity derived from the weak lensing analyses should be interpreted with caution.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The Ultimate Halo Mass in a LCDM Universe

In the far future of an accelerating LCDM cosmology, the cosmic web of large-scale structure consists of a set of increasingly isolated halos in dynamical equilibrium. We examine the approach of collisionless dark matter to hydrostatic equilibrium using a large N-body simulation evolved to scale factor a = 100, well beyond the vacuum--matter equality epoch, a_eq ~ 0.75, and 53/h Gyr into the future for a concordance model universe (Omega_m ~ 0.3, Omega_Lambda ~ 0.7). The radial phase-space structure of halos -- characterized at a < a_eq by a pair of zero-velocity surfaces that bracket a dynamically active accretion region -- simplifies at a > 10 a_eq when these surfaces merge to create a single zero-velocity surface, clearly defining the halo outer boundary, rhalo, and its enclosed mass, mhalo. This boundary approaches a fixed physical size encompassing a mean interior density ~ 5 times the critical density, similar to the turnaround value in a classical Einstein-deSitter model. We relate mhalo to other scales currently used to define halo mass (m200, mvir, m180b) and find that m200 is approximately half of the total asymptotic cluster mass, while m180b follows the evolution of the inner zero velocity surface for a < 2 but becomes much larger than the total bound mass for a > 3. The radial density profile of all bound halo material is well fit by a truncated Hernquist profile. An NFW profile provides a somewhat better fit interior to r200 but is much too shallow in the range r200 < r < rhalo.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Future Evolution of Structure in an Accelerating Universe

Current cosmological data indicate that our universe contains a substantial component of dark vacuum energy that is driving the cosmos to accelerate. We examine the immediate and longer term consequences of this dark energy (assumed here to have a constant density). Using analytic calculations and supporting numerical simulations, we present criteria for test bodies to remain bound to existing structures. We show that collapsed halos become spatially isolated and dynamically relax to a particular density profile with logarithmic slope steeper than -3 at radii beyond r_200. The asymptotic form of the space-time metric is then specified. We develop this scenario further by determining the effects of the accelerating expansion on the background radiation fields and individual particles. In an appendix, we generalize these results to include quintessence.

0 citations0 reviewsNo code linkedOpen access