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

Rachel Mandelbaum

Rachel Mandelbaum 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.IMastro-ph.GAastro-phhep-exApplicationsgr-qcMachine LearningArtificial Intelligencehep-ph

Catalog footprint

What is connected

99works
10topics
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

99 published item(s)

preprint2026arXiv

Diagnosing the Effects of Spectroscopic Training Set Imperfection on Photometric Redshift Performance

Most LSST extragalactic science will rely on photometric redshifts (photo-$z$) to extract distance information for the galaxies. However, an incomplete or non-representative training set can introduce bias into photo-$z$ estimation. It is necessary to understand how various forms of training set imperfection, such as incompleteness and non-trivial spectroscopic target selection, affect photo-$z$ estimation algorithms, and to identify metrics best-suited to quantify the impact. This work aims to systematically study metrics for diagnosing how various photo-$z$ methods react to certain types of training set incompleteness and non-representativeness. We use methods available through the open-source Python library Redshift Assessment Infrastructure Layers (RAIL) to systematically test the algorithms CMNN, GPz, FlexZBoost, and PZFlow on mock training data degraded in accordance with several existing spectroscopic sky surveys, as well as under conditions of inverse redshift incompleteness, which approximately mimics observed patterns of incompleteness at high redshift. We employ the algorithm TrainZ as a control. Finally, we quantify photo-$z$ algorithm performance using a variety of statistical metrics implemented externally to RAIL. We determine that the Kullback-Liebler Divergence, Wasserstein Distance, and Probability Integral Transform are particularly informative metrics with which to assess the impact of training set imperfection on algorithmic performance. We also find that inverse redshift incompleteness effects alone lack the complexity to realistically represent anticipated training data.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Hyrax: An Extensible Framework for Rapid ML Experimentation and Unsupervised Discovery in the Era of Rubin, Roman, and Euclid

The NSF-DOE Vera C. Rubin Observatory, Roman Space Telescope, Euclid, and other next-generation surveys will deliver imaging, spectroscopic, and time-domain data at scales that increasingly shift the bottleneck in astronomical machine learning (ML) projects from model design to infrastructure. We present Hyrax, an open-source, modular, GPU-enabled Python framework that supports the full ML lifecycle in astronomy: from data acquisition and training to inference and experiment comparison, with capabilities including multimodal dataset support, integrated vector databases for similarity search, and interactive two- and three-dimensional latent-space exploration for unsupervised discovery. We demonstrate Hyrax's versatility through five representative applications on real survey data: (i) unsupervised representation learning on $\sim 4\times10^5$ Rubin Legacy Survey of Space and Time (LSST) Data Preview 1 (DP1) galaxies, surfacing new merger and low-surface-brightness candidates missing from reference Euclid and Dark Energy Survey catalogs, while also isolating imaging artifacts -- all without labeled training data; (ii) hybrid density-based clustering for identifying cluster-scale gravitational lens candidates in DP1 data; (iii) multimodal early-time transient classification in the Zwicky Transient Facility leveraging light curves, spectra, images, and metadata; (iv) supervised false-positive filtering in shift-and-stack searches for distant solar system objects in the Dark Energy Camera Ecliptic Exploration Project survey; and (v) supervised detection of semi-resolved dwarf galaxies in Hyper Suprime-Cam and LSST-like imaging using synthetic source injection. Together, these results demonstrate that Hyrax provides astronomy-specific ML infrastructure that enables systematic discovery and rapid methodological iteration across next-generation astronomical surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

From Data to Software to Science with the Rubin Observatory LSST

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Galaxies and Halos on Graph Neural Networks: Deep Generative Modeling Scalar and Vector Quantities for Intrinsic Alignment

In order to prepare for the upcoming wide-field cosmological surveys, large simulations of the Universe with realistic galaxy populations are required. In particular, the tendency of galaxies to naturally align towards overdensities, an effect called intrinsic alignments (IA), can be a major source of systematics in the weak lensing analysis. As the details of galaxy formation and evolution relevant to IA cannot be simulated in practice on such volumes, we propose as an alternative a Deep Generative Model. This model is trained on the IllustrisTNG-100 simulation and is capable of sampling the orientations of a population of galaxies so as to recover the correct alignments. In our approach, we model the cosmic web as a set of graphs, where the graphs are constructed for each halo, and galaxy orientations as a signal on those graphs. The generative model is implemented on a Generative Adversarial Network architecture and uses specifically designed Graph-Convolutional Networks sensitive to the relative 3D positions of the vertices. Given (sub)halo masses and tidal fields, the model is able to learn and predict scalar features such as galaxy and dark matter subhalo shapes; and more importantly, vector features such as the 3D orientation of the major axis of the ellipsoid and the complex 2D ellipticities. For correlations of 3D orientations the model is in good quantitative agreement with the measured values from the simulation, except for at very small and transition scales. For correlations of 2D ellipticities, the model is in good quantitative agreement with the measured values from the simulation on all scales. Additionally, the model is able to capture the dependence of IA on mass, morphological type and central/satellite type.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Galaxy Distribution Incompleteness Testing Using Self-Organizing Maps

The calibration of redshift distributions for photometric samples using spectroscopic surveys is plagued by the difficulty in modelling the selection functions of spectroscopic surveys. In this work, we analyse how these selection functions impact redshift inference and quantify the induced biases using local calibration tests in photometry space. The study is carried out using simulations that mimic the radial selection function of a spectroscopic survey and an accompanying mock catalog of a photometric galaxy survey catalog. We use a self-organizing map to partition the photometry space and perform a local $χ^2$ test to study the probability calibration of redshift inferences that use the spectroscopic data for calibration. The goal of this work is to investigate the effect of uncorrected selection functions in the calibration data on redshift prediction accuracy and critically discuss mitigation methods. In particular we test culling-based bias correction techniques, that aim to remove redshift calibration biases by identifying regions in photometry with few spectroscopic calibration data, and propose avenues for future research. We found that removing regions in color-magnitude space that are underpopulated with spectroscopic calibration data does not remove all biases in redshift inference induced by the selection function.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

GREAT3 results I: systematic errors in shear estimation and the impact of real galaxy morphology

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the Sérsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Intrinsic alignments of bulges and discs

Galaxies exhibit coherent alignments with local structure in the Universe. This effect, called Intrinsic Alignments (IA), is an important contributor to the systematic uncertainties for wide-field weak lensing surveys. On cosmological distance scales, intrinsic shape alignments have been observed in red galaxies, which are usually bulge-dominated; while blue galaxies, which are mostly disc-dominated, exhibit shape alignments consistent with a null detection. However, disc-dominated galaxies typically consist of two prominent structures: disc and bulge. Since the bulge component has similar properties as elliptical galaxies and is thought to have formed in a similar fashion, naturally one could ask whether the bulge components exhibit similar alignments as ellipticals? In this paper, we investigate how different components of galaxies exhibit IA in the TNG100-1 cosmological hydrodynamical simulation, as well as the dependence of IA on the fraction of stars in rotation-dominated structures at $z=0$. The measurements were controlled for mass differences between the samples. We find that the bulges exhibit significantly higher IA signals, with a nonlinear alignment model amplitude of $A_I = 2.98^{+0.36}_{-0.37}$ compared to the amplitude for the galaxies as a whole (both components), $A_I = 1.13^{+0.37}_{-0.35}$. The results for bulges are statistically consistent with those for elliptical galaxies, which have $A_I = 3.47^{+0.57}_{-0.57}$. These results highlight the importance of studying galaxy dynamics in order to understand galaxy alignments and their cosmological implications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Enabling Flagship Dark Energy Experiments to Reach their Full Potential

A new generation of powerful dark energy experiments will open new vistas for cosmology in the next decade. However, these projects cannot reach their utmost potential without data from other telescopes. This white paper focuses in particular on the compelling benefits of ground-based spectroscopic and photometric observations to complement the Vera C. Rubin Observatory, as well as smaller programs in aid of a DESI-2 experiment and CMB-S4. These additional data sets will both improve dark energy constraints from these flagship projects beyond what would possible on their own and open completely new windows into fundamental physics. For example, additional photometry and single-object spectroscopy will provide necessary follow-up information for supernova and strong lensing cosmology, while highly-multiplexed spectroscopy both from smaller facilities over wide fields and from larger facilities over narrower regions of sky will yield more accurate photometric redshift estimates for weak lensing and galaxy clustering measurements from the Rubin Observatory, provide critical spectroscopic host galaxy redshifts for supernova Hubble diagrams, provide improved understanding of limiting astrophysical systematic effects, and enable new measurements that probe the nature of gravity. A common thread is that access to complementary data from a range of telescopes/instruments would have a substantial impact on the rate of advance of dark energy science in the coming years.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Rubin Observatory after LSST

The Vera C. Rubin Observatory will begin the Legacy Survey of Space and Time (LSST) in 2024, spanning an area of 18,000 square degrees in six bands, with more than 800 observations of each field over ten years. The unprecedented data set will enable great advances in the study of the formation and evolution of structure and exploration of physics of the dark universe. The observations will hold clues about the cause for the accelerated expansion of the universe and possibly the nature of dark matter. During the next decade, LSST will be able to confirm or dispute if tensions seen today in cosmological data are due to new physics. New and unexpected phenomena could confirm or disrupt our current understanding of the universe. Findings from LSST will guide the path forward post-LSST. The Rubin Observatory will still be a uniquely powerful facility even then, capable of revealing further insights into the physics of the dark universe. These could be obtained via innovative observing strategies, e.g., targeting new probes at shorter timescales than with LSST, or via modest instrumental changes, e.g., new filters, or through an entirely new instrument for the focal plane. This White Paper highlights some of the opportunities in each scenario from Rubin observations after LSST.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Snowmass2021: Opportunities from Cross-survey Analyses of Static Probes

Cosmological data in the next decade will be characterized by high-precision, multi-wavelength measurements of thousands of square degrees of the same patches of sky. By performing multi-survey analyses that harness the correlated nature of these datasets, we will gain access to new science, and increase the precision and robustness of science being pursued by each individual survey. However, effective application of such analyses requires a qualitatively new level of investment in cross-survey infrastructure, including simulations, associated modeling, coordination of data sharing, and survey strategy. The scientific gains from this new level of investment are multiplicative, as the benefits can be reaped by even present-day instruments, and can be applied to new instruments as they come online.

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

The Third Gravitational Lensing Accuracy Testing (GREAT3) Challenge Handbook

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third in a series of image analysis challenges, with a goal of testing and facilitating the development of methods for analyzing astronomical images that will be used to measure weak gravitational lensing. This measurement requires extremely precise estimation of very small galaxy shape distortions, in the presence of far larger intrinsic galaxy shapes and distortions due to the blurring kernel caused by the atmosphere, telescope optics, and instrumental effects. The GREAT3 challenge is posed to the astronomy, machine learning, and statistics communities, and includes tests of three specific effects that are of immediate relevance to upcoming weak lensing surveys, two of which have never been tested in a community challenge before. These effects include realistically complex galaxy models based on high-resolution imaging from space; spatially varying, physically-motivated blurring kernel; and combination of multiple different exposures. To facilitate entry by people new to the field, and for use as a diagnostic tool, the simulation software for the challenge is publicly available, though the exact parameters used for the challenge are blinded. Sample scripts to analyze the challenge data using existing methods will also be provided. See http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The three-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey

We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|δm|<9\times 10^{-3}$. The bias residual $δm$ shows no dependence on redshift in the range $0<z\leq 3$. We define the requirements for cosmological weak lensing science for this shear catalog, and quantify potential systematics in the catalog using a series of internal null tests for systematics related to point-spread function modelling and shear estimation. A variety of the null tests are statistically consistent with zero or within requirements, but (i) there is evidence for PSF model shape residual correlations; and (ii) star-galaxy shape correlations reveal additive systematics. Both effects become significant on $>1$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Impact of Image Persistence in the Roman Space Telescope High-Latitude Survey

The High Latitude Survey of the Nancy Grace Roman Space Telescope is expected to measure the positions and shapes of hundreds of millions of galaxies in an area of 2220 deg$^2$. This survey will provide high-quality weak lensing data with unprecedented systematics control. The Roman Space Telescope will survey the sky in near infrared (NIR) bands using Teledyne H4RG HgCdTe photodiode arrays. These NIR arrays exhibit an effect called persistence: charges that are trapped in the photodiodes during earlier exposures are gradually released into later exposures, leading to contamination of the images and potentially to errors in measured galaxy properties such as fluxes and shapes. In this work, we use image simulations that incorporate the persistence effect to study its impact on galaxy shape measurements and weak lensing signals. No significant spatial correlations are found between the galaxy shape changes induced by persistence. On the scales of interest for weak lensing cosmology, the effect of persistence on the weak lensing correlation function is about two orders of magnitude lower than the Roman Space Telescope additive shear error budget, indicating that the persistence effect is expected to be a subdominant contributor to the systematic error budget for weak lensing with the Roman Space Telescope given its current design.

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
preprint2021arXiv

The Impact of Light Polarization Effects on Weak Lensing Systematics

A fraction of the light observed from edge-on disk galaxies is polarized due to two physical effects: selective extinction by dust grains aligned with the magnetic field, and scattering of the anisotropic starlight field. Since the reflection and transmission coefficients of the reflecting and refracting surfaces in an optical system depend on the polarization of incoming rays, this optical polarization produces both (a) a selection bias in favor of galaxies with specific orientations and (b) a polarization-dependent PSF. In this work we build toy models to obtain for the first time an estimate for the impact of polarization on PSF shapes and the impact of the selection bias due to the polarization effect on the measurement of the ellipticity used in shear measurements. In particular, we are interested in determining if this effect will be significant for WFIRST. We show that the systematic uncertainties in the ellipticity components are $8\times 10^{-5}$ and $1.1 \times 10^{-4}$ due to the selection bias and PSF errors respectively. Compared to the overall requirements on knowledge of the WFIRST PSF ellipticity ($4.7\times 10^{-4}$ per component), both of these systematic uncertainties are sufficiently close to the WFIRST tolerance level that more detailed studies of the polarization effects or more stringent requirements on polarization-sensitive instrumentation for WFIRST are required.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Impact of Observing Strategy on Cosmological Constraints with LSST

The generation-defining Vera C. Rubin Observatory will make state-of-the-art measurements of both the static and transient universe through its Legacy Survey for Space and Time (LSST). With such capabilities, it is immensely challenging to optimize the LSST observing strategy across the survey's wide range of science drivers. Many aspects of the LSST observing strategy relevant to the LSST Dark Energy Science Collaboration, such as survey footprint definition, single visit exposure time and the cadence of repeat visits in different filters, are yet to be finalized. Here, we present metrics used to assess the impact of observing strategy on the cosmological probes considered most sensitive to survey design; these are large-scale structure, weak lensing, type Ia supernovae, kilonovae and strong lens systems (as well as photometric redshifts, which enable many of these probes). We evaluate these metrics for over 100 different simulated potential survey designs. Our results show that multiple observing strategy decisions can profoundly impact cosmological constraints with LSST; these include adjusting the survey footprint, ensuring repeat nightly visits are taken in different filters and enforcing regular cadence. We provide public code for our metrics, which makes them readily available for evaluating further modifications to the survey design. We conclude with a set of recommendations and highlight observing strategy factors that require further research.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The LSST DESC DC2 Simulated Sky Survey

We describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep (WFD) area of approximately 300 deg^2 as well as a deep drilling field (DDF) of approximately 1 deg^2. We simulate 5 years of the planned 10-year survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the dataset to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic testbed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time-domain cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Role of Machine Learning in the Next Decade of Cosmology

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cosmological constraints from galaxy-lensing cross correlations using BOSS galaxies with SDSS and CMB lensing

We present cosmological parameter constraints based on a joint modeling of galaxy-lensing cross correlations and galaxy clustering measurements in the SDSS, marginalizing over small-scale modeling uncertainties using mock galaxy catalogs, without explicit modeling of galaxy bias. We show that our modeling method is robust to the impact of different choices for how galaxies occupy dark matter halos and to the impact of baryonic physics (at the $\sim2\%$ level in cosmological parameters) and test for the impact of covariance on the likelihood analysis and of the survey window function on the theory computations. Applying our results to the measurements using galaxy samples from BOSS and lensing measurements using shear from SDSS galaxies and CMB lensing from Planck, with conservative scale cuts, we obtain $S_8\equiv\left(\frac{σ_8}{0.8228}\right)^{0.8}\left(\frac{Ω_m}{0.307}\right)^{0.6}=0.85\pm0.05$ (stat.) using LOWZ $\times$ SDSS galaxy lensing, and $S_8=0.91\pm0.1$ (stat.) using combination of LOWZ and CMASS $\times$ Planck CMB lensing. We estimate the systematic uncertainty in the galaxy-galaxy lensing measurements to be $\sim6\%$ (dominated by photometric redshift uncertainties) and in the galaxy-CMB lensing measurements to be $\sim3\%$, from small scale modeling uncertainties including baryonic physics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Estimating redshift distributions using Hierarchical Logistic Gaussian processes

This work uses hierarchical logistic Gaussian processes to infer true redshift distributions of samples of galaxies, through their cross-correlations with spatially overlapping spectroscopic samples. We demonstrate that this method can accurately estimate these redshift distributions in a fully Bayesian manner jointly with galaxy-dark matter bias models. We forecast how systematic biases in the redshift-dependent galaxy-dark matter bias model affect redshift inference. Using published galaxy-dark matter bias measurements from the Illustris simulation, we compare these systematic biases with the statistical error budget from a forecasted weak gravitational lensing measurement. If the redshift-dependent galaxy-dark matter bias model is mis-specified, redshift inference can be biased. This can propagate into relative biases in the weak lensing convergence power spectrum on the 10% - 30% level. We, therefore, showcase a methodology to detect these sources of error using Bayesian model selection techniques. Furthermore, we discuss the improvements that can be gained from incorporating prior information from Bayesian template fitting into the model, both in redshift prediction accuracy and in the detection of systematic modeling biases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimising LSST Observing Strategy for Weak Lensing Systematics

The LSST survey will provide unprecedented statistical power for measurements of dark energy. Consequently, controlling systematic uncertainties is becoming more important than ever. The LSST observing strategy will affect the statistical uncertainty and systematics control for many science cases; here, we focus on weak lensing systematics. The fact that the LSST observing strategy involves hundreds of visits to the same sky area provides new opportunities for systematics mitigation. We explore these opportunities by testing how different dithering strategies (pointing offsets and rotational angle of the camera in different exposures) affect additive weak lensing shear systematics on a baseline operational simulation, using the $ρ-$statistics formalism. Some dithering strategies improve systematics control at the end of the survey by a factor of up to $\sim 3-4$ better than others. We find that a random translational dithering strategy, applied with random rotational dithering at every filter change, is the most effective of those strategies tested in this work at averaging down systematics. Adopting this dithering algorithm, we explore the effect of varying the area of the survey footprint, exposure time, number of exposures in a visit, and exposure to the Galactic plane. We find that any change that increases the average number of exposures (in filters relevant to weak lensing) reduces the additive shear systematics. Some ways to achieve this increase may not be favorable for the weak lensing statistical constraining power or for other probes, and we explore the relative trade-offs between these options given constraints on the overall survey parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The LSST DESC Data Challenge 1: Generation and Analysis of Synthetic Images for Next Generation Surveys

Data Challenge 1 (DC1) is the first synthetic dataset produced by the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC). DC1 is designed to develop and validate data reduction and analysis and to study the impact of systematic effects that will affect the LSST dataset. DC1 is comprised of $r$-band observations of 40 deg$^{2}$ to 10-year LSST depth. We present each stage of the simulation and analysis process: a) generation, by synthesizing sources from cosmological N-body simulations in individual sensor-visit images with different observing conditions; b) reduction using a development version of the LSST Science Pipelines; and c) matching to the input cosmological catalog for validation and testing. We verify that testable LSST requirements pass within the fidelity of DC1. We establish a selection procedure that produces a sufficiently clean extragalactic sample for clustering analyses and we discuss residual sample contamination, including contributions from inefficiency in star-galaxy separation and imperfect deblending. We compute the galaxy power spectrum on the simulated field and conclude that: i) survey properties have an impact of 50\% of the statistical uncertainty for the scales and models used in DC1 ii) a selection to eliminate artifacts in the catalogs is necessary to avoid biases in the measured clustering; iii) the presence of bright objects has a significant impact (2- to 6-$σ$) in the estimated power spectra at small scales ($\ell > 1200$), highlighting the impact of blending in studies at small angular scales in LSST;

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Validation of PSF Models for HST and Other Space-Based Observations

Forthcoming space-based observations will require high-quality point-spread function (PSF) models for weak gravitational lensing measurements. One approach to generating these models is using a wavefront model based on the known telescope optics. We present an empirical framework for validating such models to confirm that they match the actual PSF to within requirements by comparing the models to the observed light distributions of isolated stars. We apply this framework to Tiny Tim, the standard tool for generating model PSFs for the Hubble Space Telescope (HST), testing its models against images taken by HST's Advanced Camera for Surveys in the Wide Field Channel. We show that Tiny Tim's models, in the default configuration, differ significantly from the observed PSFs, most notably in their sizes. We find that the quality of Tiny Tim PSFs can be improved through fitting the full set of Zernike polynomial coefficients which characterise the optics, to the point where the practical significance of the difference between model and observed PSFs is negligible for most use cases, resulting in additive and multiplicative biases both of order approximately 4e-4. We also show that most of this improvement can be retained through using an updated set of Zernike coefficients, which we provide.

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

Galaxy-Galaxy Lensing in HSC: Validation Tests and the Impact of Heterogeneous Spectroscopic Training Sets

Although photometric redshifts (photo-z's) are crucial ingredients for current and upcoming large-scale surveys, the high-quality spectroscopic redshifts currently available to train, validate, and test them are substantially non-representative in both magnitude and color. We investigate the nature and structure of this bias by tracking how objects from a heterogeneous training sample contribute to photo-z predictions as a function of magnitude and color, and illustrate that the underlying redshift distribution at fixed color can evolve strongly as a function of magnitude. We then test the robustness of the galaxy-galaxy lensing signal in 120 deg$^2$ of HSC-SSP DR1 data to spectroscopic completeness and photo-z biases, and find that their impacts are sub-dominant to current statistical uncertainties. Our methodology provides a framework to investigate how spectroscopic incompleteness can impact photo-z-based weak lensing predictions in future surveys such as LSST and WFIRST.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The evolution of galaxy intrinsic alignments in the MassiveBlack II universe

We investigate the redshift evolution of the intrinsic alignments (IA) of galaxies in the \texttt{MassiveBlackII} (MBII) simulation. We select galaxy samples above fixed subhalo mass cuts ($M_h>10^{11,12,13}~M_{\odot}/h$) at $z=0.6$ and trace their progenitors to $z=3$ along their merger trees. Dark matter components of $z=0.6$ galaxies are more spherical than their progenitors while stellar matter components tend to be less spherical than their progenitors. The distribution of the galaxy-subhalo misalignment angle peaks at $\sim10~\mathrm{deg}$ with a mild increase with time. The evolution of the ellipticity-direction~(ED) correlation amplitude $ω(r)$ of galaxies (which quantifies the tendency of galaxies to preferentially point towards surrounding matter overdensities) is governed by the evolution in the alignment of underlying dark matter~(DM) subhaloes to the matter density of field, as well as the alignment between galaxies and their DM subhaloes. At scales $\sim1~\mathrm{cMpc}/h$, the alignment between DM subhaloes and matter overdensity gets suppressed with time, whereas the alignment between galaxies and DM subhaloes is enhanced. These competing tendencies lead to a complex redshift evolution of $ω(r)$ for galaxies at $\sim1~\mathrm{cMpc}/h$. At scales $>1~\mathrm{cMpc}/h$, alignment between DM subhaloes and matter overdensity does not evolve significantly; the evolution of the galaxy-subhalo misalignment therefore leads to an increase in $ω(r)$ for galaxies by a factor of $\sim4$ from $z=3$ to $0.6$ at scales $>1~\mathrm{cMpc}/h$. The balance between competing physical effects is scale dependant, leading to different conclusions at much smaller scales($\sim0.1~\mathrm{Mpc}/h$).

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Tomographic galaxy clustering with the Subaru Hyper Suprime-Cam first year public data release

We analyze the clustering of galaxies in the first public data release of the HSC Subaru Strategic Program. Despite the relatively small footprints of the observed fields, the data are an excellent proxy for the deep photometric datasets that will be acquired by LSST, and are therefore an ideal test bed for the analysis methods being implemented by the LSST DESC. We select a magnitude limited sample with $i<24.5$ and analyze it in four redshift bins covering $0.15\lesssim z \lesssim1.5$. We carry out a Fourier-space analysis of the two-point clustering of this sample, including all auto- and cross-correlations. We demonstrate the use of map-level deprojection methods to account for fluctuations in the galaxy number density caused by observational systematics. Through an HOD analysis, we place constraints on the characteristic halo masses of this sample, finding a good fit up to scales $k_{\rm max}=1\,{\rm Mpc}^{-1}$, including both auto- and cross-correlations. Our results show monotonically decreasing average halo masses, which can be interpreted in terms of the drop-out of red galaxies at high redshifts for a flux-limited sample. In terms of photometric redshift systematics, we show that additional care is needed in order to marginalize over uncertainties in the redshift distribution in galaxy clustering, and that these uncertainties can be constrained by including cross-correlations. We are able to make a $\sim3σ$ detection of lensing magnification in the HSC data. Our results are stable to variations in $σ_8$ and $Ω_c$ and we find constraints that agree well with measurements from Planck and low-redshift probes. Finally, we use our pipeline to study the clustering of galaxies as a function of limiting flux, and provide a simple fitting function for the linear galaxy bias for magnitude limited samples as a function of limiting magnitude and redshift. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cross-correlating Planck CMB lensing with SDSS: Lensing-lensing and galaxy-lensing cross-correlations

We present results from cross-correlating Planck CMB lensing maps with the Sloan Digital Sky Survey (SDSS) galaxy lensing shape catalog and BOSS galaxy catalogs. For galaxy position vs. CMB lensing cross-correlations, we measure the convergence signal around the galaxies in configuration space, using the BOSS LOWZ ($z\sim0.30$) and CMASS ($z\sim0.57$) samples. With fixed Planck 2015 cosmology, doing a joint fit with the galaxy clustering measurement, for the LOWZ (CMASS) sample we find a galaxy bias $b_g=1.75\pm0.04$ ($1.95\pm 0.02$) and galaxy-matter cross-correlation coefficient $r_{cc}=1.0\pm0.2$ ($0.8\pm 0.1$) using $20<r_p<70h^{-1}$Mpc, consistent with results from galaxy-galaxy lensing. Using the same scales and including the galaxy-galaxy lensing measurements, we constrain $Ω_m=0.284\pm0.024$ and relative calibration bias between the CMB lensing and galaxy lensing to be $b_γ=0.82^{+0.15}_{-0.14}$. The combination of galaxy lensing and CMB lensing also allows us to measure the cosmological distance ratios (with $z_l\sim0.3$, $z_s\sim0.5$) $\mathcal R=\frac{D_s D_{l,*}}{D_* D_{l,s}}=2.68\pm0.29$, consistent with predictions from the Planck 2015 cosmology ($\mathcal R=2.35$). We detect the galaxy position-CMB convergence cross-correlation at small scales, $r_p<1h^{-1}$Mpc, and find consistency with lensing by NFW halos of mass $M_h\sim10^{13}h^{-1}M_\odot$. Finally, we measure the CMB lensing-galaxy shear cross-correlation, finding an amplitude of $A=0.76\pm0.23$ ($z_\text{eff}=0.35$, $θ<2^\circ$) with respect to Planck 2015 $Λ$CDM predictions ($1σ$-level consistency). We do not find evidence for relative systematics between the CMB and SDSS galaxy lensing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detection of the Splashback Radius and Halo Assembly bias of Massive Galaxy Clusters

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$σ$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $Λ$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Enabling Dark Energy Science with Deep Generative Models of Galaxy Images

Understanding the nature of dark energy, the mysterious force driving the accelerated expansion of the Universe, is a major challenge of modern cosmology. The next generation of cosmological surveys, specifically designed to address this issue, rely on accurate measurements of the apparent shapes of distant galaxies. However, shape measurement methods suffer from various unavoidable biases and therefore will rely on a precise calibration to meet the accuracy requirements of the science analysis. This calibration process remains an open challenge as it requires large sets of high quality galaxy images. To this end, we study the application of deep conditional generative models in generating realistic galaxy images. In particular we consider variations on conditional variational autoencoder and introduce a new adversarial objective for training of conditional generative networks. Our results suggest a reliable alternative to the acquisition of expensive high quality observations for generating the calibration data needed by the next generation of cosmological surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Evidence of Halo Assembly Bias in Massive Clusters

We present significant evidence of halo assembly bias for SDSS redMaPPer galaxy clusters in the redshift range $[0.1, 0.33]$. By dividing the 8,648 clusters into two subsamples based on the average member galaxy separation from the cluster center, we first show that the two subsamples have very similar halo mass of $M_{\rm 200m}\simeq 1.9\times 10^{14}~h^{-1}M_\odot$ based on the weak lensing signals at small radii $R<\sim 10~h^{-1}{\rm Mpc}$. However, their halo bias inferred from both the large-scale weak lensing and the projected auto-correlation functions differs by a factor of $\sim$1.5, which is a signature of assembly bias. The same bias hypothesis for the two subsamples is excluded at 2.5$σ$ in the weak lensing and 4.4$σ$ in the auto-correlation data, respectively. This result could bring a significant impact on both galaxy evolution and precision cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy alignments: An overview

The alignments between galaxies, their underlying matter structures, and the cosmic web constitute vital ingredients for a comprehensive understanding of gravity, the nature of matter, and structure formation in the Universe. We provide an overview on the state of the art in the study of these alignment processes and their observational signatures, aimed at a non-specialist audience. The development of the field over the past one hundred years is briefly reviewed. We also discuss the impact of galaxy alignments on measurements of weak gravitational lensing, and discuss avenues for making theoretical and observational progress over the coming decade.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy alignments: Observations and impact on cosmology

Galaxy shapes are not randomly oriented, rather they are statistically aligned in a way that can depend on formation environment, history and galaxy type. Studying the alignment of galaxies can therefore deliver important information about the physics of galaxy formation and evolution as well as the growth of structure in the Universe. In this review paper we summarise key measurements of galaxy alignments, divided by galaxy type, scale and environment. We also cover the statistics and formalism necessary to understand the observations in the literature. With the emergence of weak gravitational lensing as a precision probe of cosmology, galaxy alignments have taken on an added importance because they can mimic cosmic shear, the effect of gravitational lensing by large-scale structure on observed galaxy shapes. This makes galaxy alignments, commonly referred to as intrinsic alignments, an important systematic effect in weak lensing studies. We quantify the impact of intrinsic alignments on cosmic shear surveys and finish by reviewing practical mitigation techniques which attempt to remove contamination by intrinsic alignments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galaxy alignments: Theory, modelling and simulations

The shapes of galaxies are not randomly oriented on the sky. During the galaxy formation and evolution process, environment has a strong influence, as tidal gravitational fields in the large-scale structure tend to align nearby galaxies. Additionally, events such as galaxy mergers affect the relative alignments of both the shapes and angular momenta of galaxies throughout their history. These "intrinsic galaxy alignments" are known to exist, but are still poorly understood. This review will offer a pedagogical introduction to the current theories that describe intrinsic galaxy alignments, including the apparent difference in intrinsic alignment between early- and late-type galaxies and the latest efforts to model them analytically. It will then describe the ongoing efforts to simulate intrinsic alignments using both N-body and hydrodynamic simulations. Due to the relative youth of this field, there is still much to be done to understand intrinsic galaxy alignments and this review summarises the current state of the field, providing a solid basis for future work.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Intrinsic alignments in redMaPPer clusters -- I. Central galaxy alignments and angular segregation of satellites

The shapes of cluster central galaxies are not randomly oriented, but rather exhibit coherent alignments with the shapes of their parent clusters as well as with the surrounding large-scale structures. In this work, we aim to identify the galaxy and cluster quantities that most strongly predict the central galaxy alignment phenomenon among a large parameter space with a sample of 8237 clusters and 94817 members within 0.1<z<0.35, based on the redMaPPer cluster catalog constructed from the Sloan Digital Sky Survey. We first quantify the alignment between the projected central galaxy shapes and the distribution of member satellites, to understand what central galaxy and cluster properties most strongly correlate with these alignments. Next, we investigate the angular segregation of satellites with respect to their central galaxy major axis directions, to identify the satellite properties that most strongly predict their angular segregation. We find that central galaxies are more aligned with their member galaxy distributions in clusters that are more elongated and have higher richness, and for central galaxies with larger physical size, higher luminosity and centering probability, and redder color. Satellites with redder color, higher luminosity, located closer to the central galaxy, and with smaller ellipticity show a stronger angular segregation toward their central galaxy major axes. Finally, we provide physical explanations for some of the identified correlations, and discuss the connection to theories of central galaxy alignments, the impact of primordial alignments with tidal fields, and the importance of anisotropic accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Intrinsic alignments of BOSS LOWZ galaxies II: Impact of shape measurement methods

Measurements of intrinsic alignments of galaxy shapes with the large-scale density field, and the inferred intrinsic alignments model parameters, are sensitive to the shape measurement methods used. In this paper we measure the intrinsic alignments of the Sloan Digital Sky Survey-III (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) LOWZ galaxies using three different shape measurement methods (re-Gaussianization, isophotal, and de Vaucouleurs), identifying a variation in the inferred intrinsic alignments amplitude at the 40% level between these methods, independent of the galaxy luminosity or other properties. We also carry out a suite of systematics tests on the shapes and their two-point correlation functions, identifying a pronounced contribution from additive PSF systematics in the de Vaucouleurs shapes. Since different methods measure galaxy shapes at different effective radii, the trends we identify in the intrinsic alignments amplitude are consistent with the interpretation that the outer regions of galaxy shapes are more responsive to tidal fields, resulting in isophote twisting and stronger alignments for isophotal shapes. We observe environment dependence of ellipticity, with brightest galaxies in groups being rounder on average compared to satellite and field galaxies. We also study the anisotropy in intrinsic alignments measurements introduced by projected shapes, finding effects consistent with predictions of the nonlinear alignment model and hydrodynamic simulations. The large variations seen using the different shape measurement methods have important implications for intrinsic alignments forecasting and mitigation with future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Intrinsic alignments of disk and elliptical galaxies in the MassiveBlack-II and Illustris simulations

We study the shapes and intrinsic alignments of disks and elliptical galaxies in the MassiveBlack-II (MBII) and Illustris cosmological hydrodynamic simulations, with volumes of $(100h^{-1}Mpc)^{3}$ and $(75h^{-1}Mpc)^{3}$ respectively. We find that simulated disk galaxies are more oblate in shape and more misaligned with the shape of their host dark matter subhalo when compared with ellipticals. The disk major axis is found to be oriented towards the location of nearby elliptical galaxies. We also find that the disks are thinner in MBII and misalignments with dark matter halo orientations are smaller in both disks and ellipticals when compared with Illustris. As a result, the intrinsic alignment correlation functions at fixed mass have a higher amplitude in MBII than in Illustris. Finally, at scales above $\sim 0.1$\hmpc, the intrinsic alignment two-point correlation functions for disk galaxies in both simulations are consistent with a null detection, unlike those for ellipticals. Despite significant differences in the treatments of hydrodynamics and baryonic physics in the simulations, we find that the $w_{δ+}$ correlation function scales similarly with transverse separation. However, the less massive galaxies show different scale dependence in the ED correlation. This result indicates that, while hydrodynamic simulations are a promising tool to study intrinsic alignments, further study is needed to understand the impact of differences in the implementations of hydrodynamics and baryonic feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Maximizing Science in the Era of LSST: A Community-Based Study of Needed US Capabilities

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On Detecting Halo Assembly Bias with Galaxy Populations

The fact that the clustering of dark matter halos depends not only on their mass, but also the formation epoch, is a prominent, albeit subtle, feature of the cold dark matter structure formation theory, and is known as assembly bias. At low mass scales ($\sim 10^{12}\,h^{-1}M_\odot$), early-forming halos are predicted to be more strongly clustered than the late-forming ones. In this study we aim to robustly detect the signature of assembly bias observationally, making use of formation time indicators of central galaxies in low mass halos as a proxy for the halo formation history. Weak gravitational lensing is employed to ensure our early- and late-forming halo samples have similar masses, and are free of contamination of satellites from more massive halos. For the two formation time indicators used (resolved star formation history and current specific star formation rate), we do not find convincing evidence of assembly bias. For a pair of early- and late-forming galaxy samples with mean mass $M_{200c} \approx 9\times10^{11}\,h^{-1}M_\odot$, the relative bias is $1.00\pm 0.12$. We attribute the lack of detection to the possibilities that either the current measurements of these indicators are too noisy, or they do not correlate well with the halo formation history. Alternative proxies for the halo formation history that should perform better are suggested for future studies.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Strong bimodality in the host halo mass of central galaxies from galaxy-galaxy lensing

We use galaxy-galaxy lensing to study the dark matter halos surrounding a sample of Locally Brightest Galaxies (LBGs) selected from the Sloan Digital Sky Survey. We measure mean halo mass as a function of the stellar mass and colour of the central galaxy. Mock catalogues constructed from semi-analytic galaxy formation simulations demonstrate that most LBGs are the central objects of their halos, greatly reducing interpretation uncertainties due to satellite contributions to the lensing signal. Over the full stellar mass range, $10.3 < \log [M_*/M_\odot] < 11.6$, we find that passive central galaxies have halos that are at least twice as massive as those of star-forming objects of the same stellar mass. The significance of this effect exceeds $3σ$ for $\log [M_*/M_\odot] > 10.7$. Tests using the mock catalogues and on the data themselves clarify the effects of LBG selection and show that it cannot artificially induce a systematic dependence of halo mass on LBG colour. The bimodality in halo mass at fixed stellar mass is reproduced by the astrophysical model underlying our mock catalogue, but the sign of the effect is inconsistent with recent, nearly parameter-free age-matching models. The sign and magnitude of the effect can, however, be reproduced by halo occupation distribution models with a simple (few-parameter) prescription for type-dependence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing gravity on large scales by combining weak lensing with galaxy clustering using CFHTLenS and BOSS CMASS

We measure a combination of gravitational lensing, galaxy clustering, and redshift-space distortions called $E_G$. The quantity $E_G$ probes both parts of metric potential and is insensitive to galaxy bias and $σ_8$. These properties make it an attractive statistic to test $Λ$CDM, General Relativity and its alternate theories. We have combined CMASS DR11 with CFHTLenS and recent measurements of $β$ from RSD analysis, and find $E_G(z = 0.57) = 0.42 \pm 0.056$, an 13\% measurement in agreement with the prediction of general relativity $E_G(z = 0.57) = 0.396 \pm 0.011$ using the Planck 2015 cosmological parameters. We have corrected our measurement for various observational and theoretical systematics. Our measurement is consistent with the first measurement of $E_G$ using CMB lensing in place of galaxy lensing (Pullen et. al. 2015a) at small scales, but shows 2.8$σ$ tension when compared with their final results including large scales. This analysis with future surveys will provide improved statistical error and better control over systematics to test General Relativity and its alternate theories.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effect of detector nonlinearity on WFIRST PSF profiles for weak gravitational lensing measurements

Weak gravitational lensing (WL) is one of the most powerful techniques to learn about the dark sector of the universe. To extract the WL signal from astronomical observations, galaxy shapes must be measured and corrected for the point spread function (PSF) of the imaging system with extreme accuracy. Future WL missions (such as the Wide-Field Infrared Survey Telescope, WFIRST) will use a family of hybrid nearinfrared CMOS detectors (HAWAII-4RG) that are untested for accurate WL measurements. Like all image sensors, these devices are subject to conversion gain nonlinearities (voltage response to collected photo-charge) that bias the shape and size of bright objects such as reference stars that are used in PSF determination. We study this type of detector nonlinearity (NL) and show how to derive requirements on it from WFIRST PSF size and ellipticity requirements. We simulate the PSF optical profiles expected for WFIRST and measure the fractional error in the PSF size and the absolute error in the PSF ellipticity as a function of star magnitude and the NL model. For our nominal NL model (a quadratic correction), we find that, uncalibrated, NL can induce an error of 0.01 (fractional size) and 0.00175 (absolute ellipticity error) in the H158 bandpass for the brightest unsaturated stars in WFIRST. In addition, our simulations show that to limit the bias of the size and ellipticity errors in the H158 band to approximately 10% of the estimated WFIRST error budget, the parameter of our quadratic NL model must be calibrated to about 1% and 2.4%, respectively. We present a fitting formula that can be used to estimate WFIRST detector NL requirements once a true PSF error budget is established.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The impact of correlated noise on galaxy shape estimation for weak lensing

The robust estimation of the tiny distortions (shears) of galaxy shapes caused by weak gravitational lensing in the presence of much larger shape distortions due to the point-spread function (PSF) has been widely investigated. One major problem is that most galaxy shape measurement methods are subject to bias due to pixel noise in the images ("noise bias"). Noise bias is usually characterized using uncorrelated noise fields; however, real images typically have low-level noise correlations due to galaxies below the detection threshold, and some types of image processing can induce further noise correlations. We investigate the effective detection significance and its impact on noise bias in the presence of correlated noise for one method of galaxy shape estimation. For a fixed noise variance, the biases in galaxy shape estimates can differ substantially for uncorrelated versus correlated noise. However, use of an estimate of detection significance that accounts for the noise correlations can almost entirely remove these differences, leading to consistent values of noise bias as a function of detection significance for correlated and uncorrelated noise. We confirm the robustness of this finding to properties of the galaxy, the PSF, and the noise field, and quantify the impact of anisotropy in the noise correlations. Our results highlight the importance of understanding the pixel noise model and its impact on detection significances when correcting for noise bias on weak lensing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Impact of Interpixel Capacitance in CMOS Detectors on PSF shapes and Implications for WFIRST

Unlike optical CCDs, near-infrared detectors, which are based on CMOS hybrid readout technology, typically suffer from electrical crosstalk between the pixels. The interpixel capacitance (IPC) responsible for the crosstalk affects the point-spread function (PSF) of the telescope, increasing the size and modifying the shape of all objects in the images while correlating the Poisson noise. Upcoming weak lensing surveys that use these detectors, such as WFIRST, place stringent requirements on the PSF size and shape (and the level at which these are known), which in turn must be translated into requirements on IPC. To facilitate this process, we present a first study of the effect of IPC on WFIRST PSF sizes and shapes. Realistic PSFs are forward-simulated from physical principles for each WFIRST bandpass. We explore how the PSF size and shape depends on the range of IPC coupling with pixels that are connected along an edge or corner; for the expected level of IPC in WFIRST, IPC increases the PSF sizes by $\sim$5\%. We present a linear fitting formula that describes the uncertainty in the PSF size or shape due to uncertainty in the IPC, which could arise for example due to unknown time evolution of IPC as the detectors age or due to spatial variation of IPC across the detector. We also study of the effect of a small anisotropy in the IPC, which further modifies the PSF shapes. Our results are a first, critical step in determining the hardware and characterization requirements for the detectors used in the WFIRST survey.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A weak gravitational lensing recalibration of the scaling relations linking the gas properties of dark halos to their mass

We use weak gravitational lensing to measure mean mass profiles around Locally Brightest Galaxies (LBGs). These are selected from the SDSS/DR7 spectroscopic and photometric catalogues to be brighter than any neighbour projected within 1.0 Mpc and differing in redshift by $<1000$ km/s. Most ($> 83\%$) are expected to be the central galaxies of their dark matter halos. Previous stacking analyses have used this LBG sample to measure mean Sunyaev-Zeldovich flux and mean X-ray luminosity as a function of LBG stellar mass. In both cases, a simulation of the formation of the galaxy population was used to estimate effective halo mass for LBGs of given stellar mass, allowing the derivation of scaling relations between the gas properties of halos and their mass. By comparing results from a variety of simulations to our lensing data, we show that this procedure has significant model dependence reflecting: (i) the failure of any given simulation to reproduce observed galaxy abundances exactly; (ii) a dependence on the cosmology underlying the simulation; and (iii) a dependence on the details of how galaxies populate halos. We use our lensing results to recalibrate the scaling relations, eliminating most of this model dependence and explicitly accounting both for residual modelling uncertainties and for observational uncertainties in the lensing results. The resulting scaling relations link the mean gas properties of dark halos to their mass over an unprecedentedly wide range, $10^{12.5}<M_{500}/ \mathrm{M_\odot}<10^{14.5}$, and should fairly and robustly represent the full halo population.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detection of spatial correlations of fundamental plane residuals, and cosmological implications

The fundamental plane (FP) is a widely used tool to investigate the properties of early-type galaxies, and the tight relation between its parameters has spawned several cosmological applications, including its use as a distance indicator for peculiar velocity surveys and as a means to suppress intrinsic noise in cosmic size magnification measurements. Systematic trends with the large-scale structure across the FP could cause serious biases for these cosmological probes, but may also yield new insights into the early-type population. Here we report the first detection of spatial correlations among offsets in galaxy size from an FP that explicitly accounts for redshift trends, using a sample of about $95,000$ elliptical galaxies from the Sloan Digital Sky Survey. We show that these offsets correlate with the density field out to at least $10h^{-1}$Mpc at $4σ$ significance in a way that cannot be explained by systematic errors in galaxy size estimates. We propose a physical explanation for the correlations by dividing the sample into central, satellite, and field galaxies, identifying trends for each galaxy type separately. Central (satellite) galaxies lie on average above (below) the FP, which we argue could be due to a higher (lower) than average mass-to-light ratio. We fit a simple model to the correlations of FP residuals and use it to predict the impact on peculiar velocity power spectra, finding a contamination larger than $10\,\%$ for $k>0.04\,h/$Mpc. Moreover, cosmic magnification measurements based on an FP could be severely contaminated over a wide range of scales by the intrinsic FP correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Galaxy shapes and alignments in the MassiveBlack-II hydrodynamic and dark matter-only simulations

We compare the shapes and intrinsic alignments of galaxies in the MassiveBlack-II cosmological hydrodynamic simulation (MBII) to those in a dark matter-only (DMO) simulation performed with the same volume (100$h^{-1}$Mpc)$^{3}$, cosmological parameters, and initial conditions. Understanding the impact of baryonic physics on galaxy shapes and alignments and their relation to the dark matter distribution should prove useful to map the intrinsic alignments of galaxies from hydrodynamic to dark matter-only simulations. We find that dark matter subhalos are typically rounder in MBII, and the shapes of stellar matter in low mass galaxies are more misaligned with the shapes of the dark matter of the corresponding subhalos in the DMO simulation. At $z=0.06$, the fractional difference in the mean misalignment angle between MBII and DMO simulations varies from $\sim 28 \% - 12 \%$ in the mass range $10^{10.8} - 6.0 \times 10^{14} h^{-1}M_{\odot}$. We study the dark matter halo shapes and alignments as a function of radius, and find that while galaxies in MBII are more aligned with the inner parts of their dark matter subhalos, there is no radial trend in their alignments with the corresponding subhalo in the DMO simulation. This result highlights the importance of baryonic physics in determining the alignment of the galaxy with respect to the inner parts of the halo. Finally, we compare the ellipticity-direction (ED) correlation for galaxies to that for dark matter halos, finding that it is suppressed on all scales by stellar-dark matter misalignment. In the projected shape-density correlation ($w_{δ+}$), which includes ellipticity weighting, this effect is partially canceled by the higher mean ellipticities of the stellar component, but differences of order $30-40\%$ remain on scales $> 1$ Mpc over a range of subhalo masses, with scale-dependent effects below $1$ Mpc.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

GalSim: The modular galaxy image simulation toolkit

GALSIM is a collaborative, open-source project aimed at providing an image simulation tool of enduring benefit to the astronomical community. It provides a software library for generating images of astronomical objects such as stars and galaxies in a variety of ways, efficiently handling image transformations and operations such as convolution and rendering at high precision. We describe the GALSIM software and its capabilities, including necessary theoretical background. We demonstrate that the performance of GALSIM meets the stringent requirements of high precision image analysis applications such as weak gravitational lensing, for current datasets and for the Stage IV dark energy surveys of the Large Synoptic Survey Telescope, ESA's Euclid mission, and NASA's WFIRST-AFTA mission. The GALSIM project repository is public and includes the full code history, all open and closed issues, installation instructions, documentation, and wiki pages (including a Frequently Asked Questions section). The GALSIM repository can be found at https://github.com/GalSim-developers/GalSim .

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Instrumental systematics and weak gravitational lensing

We present a pedagogical review of the weak gravitational lensing measurement process and its connection to major scientific questions such as dark matter and dark energy. Then we describe common ways of parametrizing systematic errors and understanding how they affect weak lensing measurements. Finally, we discuss several instrumental systematics and how they fit into this context, and conclude with some future perspective on how progress can be made in understanding the impact of instrumental systematics on weak lensing measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intrinsic alignments of galaxies in the MassiveBlack-II simulation: analysis of two-point statistics

The intrinsic alignment of galaxies with the large-scale density field is an important astrophysical contaminant in upcoming weak lensing surveys. We present detailed measurements of the galaxy intrinsic alignments and associated ellipticity-direction (ED) and projected shape ($w_{g+}$) correlation functions for galaxies in the cosmological hydrodynamic MassiveBlack-II (MB-II) simulation. We carefully assess the effects on galaxy shapes, misalignment of the stellar component with the dark matter shape and two-point statistics of iterative weighted (by mass and luminosity) definitions of the (reduced and unreduced) inertia tensor. We find that iterative procedures must be adopted for a reliable measurement of the reduced tensor but that luminosity versus mass weighting has only negligible effects. Both ED and $w_{g+}$ correlations increase in amplitude with subhalo mass (in the range of $10^{10} - 6.0\times 10^{14}h^{-1}M_{\odot}$), with a weak redshift dependence (from $z=1$ to $z=0.06$) at fixed mass. At $z \sim 0.3$, we predict a $w_{g+}$ that is in reasonable agreement with SDSS LRG measurements and that decreases in amplitude by a factor of $\sim 5$--18 for galaxies in the LSST survey. We also compared the intrinsic alignments of centrals and satellites, with clear detection of satellite radial alignments within their host halos. Finally, we show that $w_{g+}$ (using subhalos as tracers of density) and $w_{δ+}$ (using dark matter density) predictions from the simulations agree with that of non-linear alignment models (NLA) at scales where the 2-halo term dominates in the correlations (and tabulate associated NLA fitting parameters). The 1-halo term induces a scale dependent bias at small scales which is not modeled in the NLA model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intrinsic alignments of SDSS-III BOSS LOWZ sample galaxies

Intrinsic alignments (IA) of galaxies, i.e. correlations of galaxy shapes with each other or with the density field, are a major astrophysical source of contamination for weak lensing surveys. We present the results of IA measurements of galaxies on 0.1- 200 Mpc/h scales using the SDSS-III BOSS LOWZ sample, in the redshift range 0.16<z<0.36. We extend the existing IA measurements for spectroscopic LRGs to lower luminosities, and show that the luminosity dependence of large-scale IA can be well-described by a power law. Within the limited redshift and color range of our sample, we observe no significant redshift or color dependence of IA. We measure the halo mass of LOWZ galaxies using galaxy-galaxy lensing, and show that the mass dependence of large-scale IA is also well described by a power law. We detect variation in the scale dependence of IA with mass and luminosity, which underscores the need to use flexible templates in order to remove the IA signal. We also study the environment dependence of IA by splitting the sample into field and group galaxies, which are further split into satellite and central galaxies. We show that group central galaxies are aligned with their halos at small scales and also are aligned with the tidal fields out to large scales. We also detect the radial alignments of satellite galaxies within groups, which results in a null detection of large-scale intrinsic alignments for satellites. These results can be used to construct better intrinsic alignment models for removal of this contaminant to the weak lensing signal.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Investigating Galaxy-Filament Alignments in Hydrodynamic Simulations using Density Ridges

In this paper, we study the filamentary structures and the galaxy alignment along filaments at redshift $z=0.06$ in the MassiveBlack-II simulation, a state-of-the-art, high-resolution hydrodynamical cosmological simulation which includes stellar and AGN feedback in a volume of (100 Mpc$/h$)$^3$. The filaments are constructed using the subspace constrained mean shift (SCMS; Ozertem & Erdogmus (2011) and Chen et al. (2015a)). First, we show that reconstructed filaments using galaxies and reconstructed filaments using dark matter particles are similar to each other; over $50\%$ of the points on the galaxy filaments have a corresponding point on the dark matter filaments within distance $0.13$ Mpc$/h$ (and vice versa) and this distance is even smaller at high-density regions. Second, we observe the alignment of the major principal axis of a galaxy with respect to the orientation of its nearest filament and detect a $2.5$ Mpc$/h$ critical radius for filament's influence on the alignment when the subhalo mass of this galaxy is between $10^9M_\odot/h$ and $10^{12}M_\odot/h$. Moreover, we find the alignment signal to increase significantly with the subhalo mass. Third, when a galaxy is close to filaments (less than $0.25$ Mpc$/h$), the galaxy alignment toward the nearest galaxy group depends on the galaxy subhalo mass. Finally, we find that galaxies close to filaments or groups tend to be rounder than those away from filaments or groups.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large-scale structure and the intrinsic alignment of galaxies

Coherent alignments of galaxy shapes, often called "intrinsic alignments" (IA), are the most significant source of astrophysical uncertainty in weak lensing measurements. We develop the tidal alignment model of IA and demonstrate its success in describing observational data. We also describe a technique to separate IA from galaxy-galaxy lensing measurements. Applying this technique to luminous red galaxy lenses in the Sloan Digital Sky Survey, we constrain potential IA contamination from associated sources to be below a few percent.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Luminous Red Galaxies in Clusters: Central Occupation, Spatial Distributions, and Mis-centering

Luminous Red Galaxies (LRG) from the Sloan Digital Sky Survey are considered among the best understood samples of galaxies, and they are employed in a broad range of cosmological studies. Because they form a relatively homogeneous population, with high stellar masses and red colors, they are expected to occupy halos in a relatively simple way. In this paper, we study how LRGs occupy massive halos via direct counts in clusters and we reveal several unexpected trends suggesting that the connection between LRGs and dark matter halos may not be straightforward. Using the redMaPPer cluster catalog, we derive the central occupation of LRGs as a function richness, Ncen(λ). Assuming no correlation between cluster mass and central galaxy luminosity at fixed richness, we show that clusters contain a significantly lower fraction of central LRGs than predicted from the two-point correlation function. At halo masses of 10^14.5 Msun, we find Ncen=0.73, compared to Ncen of 0.89 from correlation studies. Our central occupation function for LRGs converges to 0.95 at large halo masses. A strong anti-correlation between central luminosity and cluster mass at fixed richness is required to reconcile our results with those based on clustering studies. We also derive P_BNC, the probability that the brightest cluster member is not the central galaxy. We find P_BNC ~ 20-30% which is a factor of ~2 lower than the value found by Skibba et al. 2011. Finally, we study the radial offsets of bright non-central LRGs from cluster centers and show that bright non-central LRGs follow a different radial distribution compared to red cluster members, which follow a Navarro-Frank-White profile. This work demonstrates that even the most massive clusters do not always have an LRG at the center, and that the brightest galaxy in a cluster is not always the central galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mapping stellar content to dark matter halos using galaxy clustering and galaxy-galaxy lensing in the SDSS DR7

The mapping between the distributions of the observed galaxy stellar mass and the underlying dark matter halos provides the crucial link from theories of large-scale structure formation to interpreting the complex phenomena of galaxy formation and evolution. We develop a novel statistical method, based on the Halo Occupation Distribution model (HOD), to solve for this mapping by jointly fitting the galaxy clustering and the galaxy-galaxy lensing measured from the Sloan Digital Sky Survey (SDSS). The method, called the iHOD model, extracts maximum information from the survey by including ~80% more galaxies than the traditional HOD methods, and takes into account the incompleteness of the stellar mass samples in a statistically consistent manner. The derived stellar-to-halo mass relation not only explains the clustering and lensing of SDSS galaxies over almost four decades in stellar mass, but also successfully predicts the stellar mass functions observed in SDSS. Due to its capability of modelling significantly more galaxies, the iHOD is able to break the degeneracy between the logarithmic scatter in the stellar mass at fixed halo mass and the slope of the stellar-to-halo mass relation at high mass end, without the need to assume a strong prior on the scatter and/or use the stellar mass function as an input. We detect a decline of the scatter with halo mass, from 0.22 dex at below 10^{12} Msun to 0.18 dex at 10^{14} Msun. The model also enables stringent constraints on the satellite stellar mass functions at fixed halo mass, predicting a departure from the Schechter functional form in high mass halos. The iHOD model can be easily applied to existing and future spectroscopic datasets, greatly improving the statistical constraint on the stellar-to-halo mass relation compared to the traditional HOD methods within the same survey.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mapping stellar content to dark matter halos. II. Halo mass is the main driver of galaxy quenching

We develop a simple yet comprehensive method to distinguish the underlying drivers of galaxy quenching, using the clustering and galaxy-galaxy lensing of red and blue galaxies in SDSS. Building on the iHOD framework developed by Zu & Mandelbaum (2015a), we consider two quenching scenarios: 1) a "halo" quenching model in which halo mass is the sole driver for turning off star formation in both centrals and satellites; and 2) a "hybrid" quenching model in which the quenched fraction of galaxies depends on their stellar mass while the satellite quenching has an extra dependence on halo mass. The two best-fit models describe the red galaxy clustering and lensing equally well, but halo quenching provides significantly better fits to the blue galaxies above $10^{11} M_\odot/h^2$. The halo quenching model also correctly predicts the average halo mass of the red and blue centrals, showing excellent agreement with the direct weak lensing measurements of locally brightest galaxies. Models in which quenching is not tied to halo mass, including an age-matching model in which galaxy colour depends on halo age at fixed $M_*$, fail to reproduce the observed halo mass for massive blue centrals. We find similar critical halo masses responsible for the quenching of centrals and satellites (~$1.5\times10^{12} M\odot/h^2$), hinting at a uniform quenching mechanism for both, e.g., the virial shock-heating of infalling gas. The success of the iHOD halo quenching model provides strong evidence that the physical mechanism that quenches star formation in galaxies is tied principally to the masses of their dark matter halos rather than the properties of their stellar components.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The impact of cosmic variance on simulating weak lensing surveys

Upcoming weak lensing surveys will survey large cosmological volumes to measure the growth of cosmological structure with time and thereby constrain dark energy. One major systematic uncertainty in this process is the calibration of the weak lensing shape distortions, or shears. Most upcoming surveys plan to test several aspects of their shear estimation algorithms using sophisticated image simulations that include realistic galaxy populations based on high-resolution data from the Hubble Space Telescope (HST). However, existing datasets from the (HST) cover very small cosmological volumes, so cosmic variance could cause the galaxy populations in them to be atypical. A narrow redshift slice from such surveys could be dominated by a single large overdensity or underdensity. In that case, the morphology-density relation could alter the local galaxy populations and yield an incorrect calibration of shear estimates as a function of redshift. We directly test this scenario using the COSMOS survey, the largest-area (HST) survey to date, and show how the statistical distributions of galaxy shapes and morphological parameters (e.g., Sérsic $n$) are influenced by redshift-dependent cosmic variance. The typical variation in RMS ellipticity due to environmental effects is 5 per cent (absolute, not relative) for redshift bins of width $Δz=0.05$, which could result in uncertain shear calibration at the 1 per cent level. We conclude that the cosmic variance effects are large enough to exceed the systematic error budget of future surveys, but can be mitigated with careful choice of training dataset and sufficiently large redshift binning.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Weak Lensing Signal and the Clustering of BOSS Galaxies I: Measurements

A joint analysis of the clustering of galaxies and their weak gravitational lensing signal is well-suited to simultaneously constrain the galaxy-halo connection as well as the cosmological parameters by breaking the degeneracy between galaxy bias and the amplitude of clustering signal. In a series of two papers, we perform such an analysis at the highest redshift ($z\sim0.53$) in the literature using CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Eleventh Data Release (SDSS-III/BOSS DR11) catalog spanning 8300~deg$^2$. In this paper, we present details of the clustering and weak lensing measurements of these galaxies. We define a subsample of 400,916 CMASS galaxies based on their redshifts and stellar mass estimates so that the galaxies constitute an approximately volume-limited and similar population over the redshift range $0.47\le z\le 0.59$. We obtain a signal-to-noise ratio $S/N\simeq 56$ for the galaxy clustering measurement. We also explore the redshift and stellar mass dependence of the clustering signal. For the weak lensing measurement, we use existing deeper imaging data from the CFHTLS with publicly available shape and photometric redshift catalogs from CFHTLenS, but only in a 105~deg$^2$ area which overlaps with BOSS. This restricts the lensing measurement to only 5,084 CMASS galaxies. After careful systematic tests, we find a highly significant detection of the CMASS weak lensing signal, with total $S/N\simeq 26$. These measurements form the basis of the halo occupation distribution and cosmology analysis presented in More et al. (Paper II).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Weak Lensing Signal and the Clustering of BOSS Galaxies II: Astrophysical and Cosmological Constraints

We perform a joint analysis of the abundance, the clustering and the galaxy-galaxy lensing signal of galaxies measured from Data Release 11 of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS III-BOSS) in our companion paper, Miyatake et al. (2014). The lensing signal was obtained by using the shape catalog of background galaxies from the Canada France Hawaii Telescope Legacy Survey, which was made publicly available by the CFHTLenS collaboration, with an area overlap of about 105 deg$^2$. We analyse the data in the framework of the halo model in order to fit halo occupation parameters and cosmological parameters ($Ω_{\rm m}$ and $σ_8$) to these observables simultaneously, and thus break the degeneracy between galaxy bias and cosmology. Adopting a flat $Λ$CDM cosmology with priors on $Ω_b h^2$, $n_{\rm s}$ and $h$ from the analysis of WMAP 9-year data, we obtain constraints on the stellar mass-halo mass relation of galaxies in our sample. Marginalizing over the halo occupation distribution parameters and a number of other nuisance parameters in our model, we obtain $Ω_{\rm m}=0.310^{+0.019}_{-0.020}$ and $σ_8=0.785^{+0.044}_{-0.044}$ (68% confidence). We demonstrate the robustness of our results with respect to sample selection and a variety of systematics such as the halo off-centering effect and possible incompleteness in our sample. Our constraints are consistent, complementary and competitive with those obtained using other independent probes of these cosmological parameters. The cosmological analysis is the first of its kind to be performed at a redshift as high as $0.53$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Background sky obscuration by cluster galaxies as a source of systematic error for weak lensing

Lensing magnification and stacked shear measurements of galaxy clusters rely on measuring the density of background galaxies behind the clusters. The most common ways of measuring this quantity ignore the fact that some fraction of the sky is obscured by the cluster galaxies themselves, reducing the area in which background galaxies can be observed. We discuss the size of this effect in the Sloan Digital Sky Survey (SDSS) and the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), finding a minimum 1 per cent effect at $0.1h^{-1}$Mpc from the centers of clusters in SDSS; the effect is an order of magnitude higher in CFHTLenS. The resulting biases on cluster mass and concentration measurements are of the same order as the size of the obscuration effect, which is below the statistical errors for cluster lensing in SDSS but likely exceeds them for CFHTLenS. We also forecast the impact of this systematic error on cluster mass and magnification measurements in several upcoming surveys, and find that it typically exceeds the statistical errors. We conclude that future surveys must account for this effect in stacked lensing and magnification measurements in order to avoid being dominated by systematic error.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Galaxy and Mass Assembly (GAMA): The halo mass of galaxy groups from maximum-likelihood weak lensing

We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, $M_\star$, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium $Λ$CDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median $M_h|M_\star$ measured from weak lensing depends more sensitively on the lognormal dispersion in $M_\star$ at fixed $M_h$ than it does on the median $M_\star|M_h$. Our measurements suggest an intrinsic dispersion of $σ_{\log(M_\star)}\sim 0.15$; 4) Comparing our mass estimates with those in the catalogue, we find that the G$^3$Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Galaxy Halo Masses from Weak Gravitational Lensing

In this review, I discuss the use of galaxy-galaxy weak lensing measurements to study the masses of dark matter halos in which galaxies reside. After summarizing how weak gravitational lensing measurements can be interpreted in terms of halo mass, I review measurements that were used to derive the relationship between optical galaxy mass tracers, such as stellar mass or luminosity, and dark matter halo mass. Measurements of galaxy-galaxy lensing from the past decade have led to increasingly tight constraints on the connection between dark matter halo mass and optical mass tracers, including both the mean relationships between these quantities and the intrinsic scatter between them. I also review some of the factors that can complicate analysis, such as the choice of modeling procedure, and choices made when dividing up samples of lens galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Galaxy Shapes and Intrinsic Alignments in The MassiveBlack-II Simulation

The intrinsic alignment of galaxy shapes with the large-scale density field is a contaminant to weak lensing measurements, as well as being an interesting signature of galaxy formation and evolution (albeit one that is difficult to predict theoretically). Here we investigate the shapes and relative orientations of the stars and dark matter of halos and subhalos (central and satellite) extracted from the MassiveBlack-II simulation, a state-of-the-art high resolution hydrodynamical cosmological simulation which includes stellar and AGN feedback in a volume of $(100{h^{-1}\mathrm{Mpc}})^3$. We consider redshift evolution from $z=1$ to $0.06$ and mass evolution within the range of subhalo masses, $10^{10} -6.0 \times 10^{14.0}{h^{-1}M_{\odot}}$. The shapes of the dark matter distributions are generally more round than the shapes defined by stellar matter. The projected root-mean-square (RMS) ellipticity per component for stellar matter is measured to be $e_\text{rms} = 0.28$ at $z=0.3$ for $M_{subhalo}> 10^{12.0}{h^{-1}M_{\odot}}$, which compares favourably with observational measurements. We find that the shapes of stellar and dark matter are more round for less massive subhalos and at lower redshifts. By directly measuring the relative orientation of the stellar matter and dark matter of subgroups, we find that, on average, the misalignment between the two components is larger for less massive subhalos. The mean misalignment angle varies from $\sim 30^{\circ}-10^{\circ}$ for $M \sim 10^{10} - 10^{14} {h^{-1}M_{\odot}}$ and shows a weak dependence on redshift. We also compare the misalignment angles in central and satellite subhalos at fixed subhalo mass, and find that centrals are more misaligned than satellites. We present fitting formulae for the shapes of dark and stellar matter in subhalos and also the probability distributions of misalignment angles.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Intrinsic alignments of group and cluster galaxies in photometric surveys

Intrinsic alignments of galaxies have been shown to contaminate weak gravitational lensing observables on linear scales, $r>$ 10 $h^{-1}$Mpc, but studies of alignments in the non-linear regime have thus far been inconclusive. We present an estimator for extracting the intrinsic alignment signal of galaxies around stacked clusters of galaxies from multiband imaging data. Our estimator removes the contamination caused by galaxies that are gravitationally lensed by the clusters and scattered in redshift space due to photometric redshift uncertainties. It uses posterior probability distributions for the redshifts of the galaxies in the sample and it is easily extended to obtain the weak gravitational lensing signal while removing the intrinsic alignment contamination. We apply this algorithm to groups and clusters of galaxies identified in the Sloan Digital Sky Survey `Stripe 82' coadded imaging data over $\sim 150$ deg$^2$. We find that the intrinsic alignment signal around stacked clusters in the redshift range $0.1<z<0.4$ is consistent with zero. In terms of the tidal alignment model of Catelan et al. (2001), we set joint constraints on the strength of the alignment and the bias of the lensing groups and clusters on scales between 0.1 and $10\,h^{-1}$ Mpc, $b_LC_1ρ_{\rm crit} = -2_{-14}^{+14} \times 10^{-4}$. This constrains the contamination fraction of alignment to lensing signal to the range between $[-18,23]$ per cent below scales of 1 $h^{-1}$ Mpc at 95 per cent confidence level, and this result depends on our photometric redshift quality and selection criteria used to identify background galaxies. Our results are robust to the choice of photometric band in which the shapes are measured ($i$ and $r$) and to centring on the Brightest Cluster Galaxy or on the geometrical centre of the clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Weighing "El Gordo" with a Precision Scale: Hubble Space Telescope Weak-lensing Analysis of the Merging Galaxy Cluster ACT-CL J0102-4915 at z=0.87

(Abridged) We present a HST weak-lensing study of the merging galaxy cluster "El Gordo" (ACT-CL J0102-4915) at z=0.87 discovered by the Atacama Cosmology Telescope collaboration as the strongest SZ decrement in its ~1000 sq. deg survey. Our weak-lensing analysis confirms that ACT-CL J0102-4915 is indeed an extreme system consisting of two massive (~10^15 Msun each) subclusters with a projected separation of ~0.7 Mpc. This binary mass structure revealed by our lensing study is consistent with the cluster galaxy distribution and the dynamical study carried out with 89 spectroscopic members. We estimate the mass of ACT-CL J0102-4915 by simultaneously fitting two axisymmetric NFW profiles allowing their centers to vary. Our MCMC analysis shows that the masses of the northwestern (NW) and the southeastern (SE) components are M200c=(1.38+-0.22) x 10^15 Msun and (0.78+-0.20) x 10^15 Msun, respectively. The lensing-based velocity dispersions are consistent with their spectroscopic measurements. The centroids of both components are tightly constrained (~4") and close to the optical luminosity centers. The X-ray and mass peaks are spatially offset by ~8" (~62 kpc), which is significant at the ~2 sigma confidence level and confirms that the baryonic and dark matter in this cluster are disassociated. The dark matter peak, however, does not lead the gas peak in the direction expected if we are viewing the cluster soon after first core passage during a high speed merger. Under the assumption that the merger is happening in the plane of the sky, extrapolation of the two NFW halos to a radius r200a=2.4 Mpc yields a combined mass of M200a=(3.13+-0.56) x 10^15 Msun. This extrapolated total mass is consistent with our two-component-based dynamical analysis and previous X-ray measurements, projecting ACT-CL J0102-4915 to be the most massive cluster at z>0.6 known to date.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A highly elongated prominent lens at z=0.87: first strong lensing analysis of El Gordo

We present the first strong-lensing (SL) analysis of the galaxy cluster ACT-CL J0102-4915 (\emph{El Gordo}), in recent \emph{HST}/ACS images, revealing a prominent strong lens at a redshift of $z=0.87$. This finding adds to the already-established unique properties of \emph{El Gordo}: it is the most massive, hot, X-ray luminous, and bright Sunyaev-Zeldovich effect cluster at $z\gtrsim0.6$, and the only `bullet'-like merging cluster known at these redshifts. The lens consists of two merging massive clumps, where for a source redshift of $z_{s}\sim2$ each clump exhibits only a small, separate critical area, with a total area of $0.69\pm0.11\sq\arcmin$ over the two clumps. For a higher source redshift, $z_{s}\sim4$, the critical curves of the two clumps merge together into one bigger and very elongated lens (axis ratio $\simeq5.5$), enclosing an effective area of $1.44\pm0.22\sq\arcmin$. The critical curves continue expanding with increasing redshift so that for high-redshift sources ($z_{s}\gtrsim9$) they enclose an area of $\sim1.91\pm0.30\sq\arcmin$ (effective $θ_{e}\simeq46.8\pm3.7\arcsec$) and a mass of $6.09\pm1.04\times10^{14}M_{\odot}$. According to our model, the area of high magnification ($μ>10$) for such high redshift sources is $\simeq1.2\sq\arcmin$, and the area with $μ>5$ is $\simeq2.3\sq\arcmin$, making \emph{El Gordo} a compelling target for studying the high-redshift Universe. We obtain a strong lower limit on the total mass of \emph{El Gordo}, $\gtrsim1.7\times10^{15}M_{\odot}$ from the SL regime alone, suggesting a total mass of, roughly, $M_{200}\sim2.3\times10^{15}M_{\odot}$. Our results should be revisited when additional spectroscopic and \emph{HST} imaging data are available.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cosmological parameter constraints from galaxy-galaxy lensing and galaxy clustering with the SDSS DR7

Recent studies have shown that the cross-correlation coefficient between galaxies and dark matter is very close to unity on scales outside a few virial radii of galaxy halos, independent of the details of how galaxies populate dark matter halos. This finding makes it possible to determine the dark matter clustering from measurements of galaxy-galaxy weak lensing and galaxy clustering. We present new cosmological parameter constraints based on large-scale measurements of spectroscopic galaxy samples from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We generalise the approach of Baldauf et al. (2010) to remove small scale information (below 2 and 4 Mpc/h for lensing and clustering measurements, respectively), where the cross-correlation coefficient differs from unity. We derive constraints for three galaxy samples covering 7131 sq. deg., containing 69150, 62150, and 35088 galaxies with mean redshifts of 0.11, 0.28, and 0.40. We clearly detect scale-dependent galaxy bias for the more luminous galaxy samples, at a level consistent with theoretical expectations. When we vary both σ_8 and Ω_m (and marginalise over non-linear galaxy bias) in a flat LCDM model, the best-constrained quantity is σ_8 (Ω_m/0.25)^{0.57}=0.80 +/- 0.05 (1-sigma, stat. + sys.), where statistical and systematic errors have comparable contributions, and we fixed n_s=0.96 and h=0.7. These strong constraints on the matter clustering suggest that this method is competitive with cosmic shear in current data, while having very complementary and in some ways less serious systematics. We therefore expect that this method will play a prominent role in future weak lensing surveys. When we combine these data with WMAP7 CMB data, constraints on σ_8, Ω_m, H_0, w_{de} and \sum m_ν become 30--80 per cent tighter than with CMB data alone, since our data break several parameter degeneracies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Galaxy and Mass Assembly (GAMA): galaxy radial alignments in GAMA groups

We constrain the distributions of projected radial alignment angles of satellite galaxy shapes within the Galaxy And Mass Assembly survey group catalogue. We identify the galaxy groups using spectroscopic redshifts and measure galaxy projected ellipticities from Sloan Digital Sky Survey imaging. With a sample of 3,850 groups with 13,655 satellite galaxies with high quality shape measurements, we find a less than 2-sigma signal of radial alignments in the mean projected ellipticity components and the projected position angle when using galaxy shape estimates optimized for weak lensing measurements. Our radial alignment measurement increases to greater than 3-sigma significance relative to the expectation for no alignments if we use 2-D Sérsic model fits to define galaxy orientations. Our weak measurement of radial alignments is in conflict with predictions from dark matter N-body simulations, which we interpret as evidence for large mis-alignments of baryons and dark matter in group and cluster satellites. Within our uncertainties, that are dominated by our small sample size, we find only weak and marginally significant trends of the radial alignment angle distributions on projected distance from the group centre, host halo mass, and redshift that could be consistent with a tidal torquing mechanism for radial alignments. Using our lensing optimized shape estimators, we estimate that intrinsic alignments of galaxy group members may contribute a systematic error to the mean differential projected surface mass density of groups inferred from weak lensing observations by -1 +/- 20% at scales around 300 kpc/h from the group centre assuming a photometric redshift r.m.s. error of 10%, and given our group sample with median redshift of 0.17 and median virial masses ~10^{13} h^{-1}M_{sun}.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Growth of Cosmic Structure: Probing Dark Energy Beyond Expansion

The quantity and quality of cosmic structure observations have greatly accelerated in recent years. Further leaps forward will be facilitated by imminent projects, which will enable us to map the evolution of dark and baryonic matter density fluctuations over cosmic history. The way that these fluctuations vary over space and time is sensitive to the nature of dark matter and dark energy. Dark energy and gravity both affect how rapidly structure grows; the greater the acceleration, the more suppressed the growth of structure, while the greater the gravity, the more enhanced the growth. While distance measurements also constrain dark energy, the comparison of growth and distance data tests whether General Relativity describes the laws of physics accurately on large scales. Modified gravity models are able to reproduce the distance measurements but at the cost of altering the growth of structure (these signatures are described in more detail in the accompanying paper on Novel Probes of Gravity and Dark Energy). Upcoming surveys will exploit these differences to determine whether the acceleration of the Universe is due to dark energy or to modified gravity. To realize this potential, both wide field imaging and spectroscopic redshift surveys play crucial roles. Projects including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map more than a 1000 cubic-billion-light-year volume of the Universe. These will map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over the last 3/4 of the age of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Where are the Luminous Red Galaxies (LRGs)? Using correlation measurements and lensing to relate LRGs to dark matter halos

Nonlinear redshift-space distortions, the Finger-of-God (FoG) effect, can complicate the interpretation of the galaxy power spectrum. Here, we demonstrate the method proposed by Hikage et al. (2012) to use complimentary observations to directly constrain this effect on the data. We use catalogs of Luminous Red Galaxies (LRGs) and photometric galaxies from the SDSS DR7 to measure the redshift-space power spectrum of LRGs, the cross-correlation of LRGs with the shapes of background photometric galaxies (galaxy-galaxy weak lensing), and the projected cross-correlation of LRGs with photometric galaxies having similar photometric redshifts to the LRG spectroscopic redshift. All of these measurements use a reconstructed halo field. While we use the position of each LRG for single LRG systems, we compare the measurements using different halo-center proxies for multiple-LRG systems (4.5 per cent of all the halos): the brightest LRG position (BLRG), the faintest LRG position (FLRG) and their arithmetical mean position (Mean), respectively, in each system. We find significant differences in the measured correlations of different centers, showing consistent off-centering effects in the three observables. By comparing the measurements with a halo model that treats the satellite photometric galaxies as being distributed according to a generalized NFW profile, we find that about 40 (70) per cent of BLRGs (FLRGs) are off-centered satellite galaxies in the multiple-LRG systems. The satellite LRGs have typical off-centering radius of about 400 kpc/h, and velocity dispersion of about 500 km/s in host halos with a mean mass of 1.6x10^14 Ms/h. We show that, if LRGs in the single LRG systems have similar offsets, the residual FoG contamination in the LRG power spectrum can be significant at k>0.1 h/Mpc, which may cause a bias in cosmological parameters such as the neutrino mass.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cluster Density Profiles as a Test of Modified Gravity

We present a new test of gravitational interactions at the r\sim(0.2-20)Mpc scale, around the virial radius of dark matter halos measured through cluster-galaxy lensing of maxBCG clusters from the Sloan Digital Sky Survey (SDSS). We employ predictions from self-consistent simulations of f(R) gravity to find an upper bound on the background field amplitude of f_R0<3.5x10^-3 at the 1D-marginalized 95% confidence level. As a model-independent assessment of the constraining power of cluster profiles measured through weak gravitational lensing, we also constrain the amplitude F_0 of a phenomenological modification based on the profile enhancement induced by f(R) gravity when not including effects from the increased cluster abundance in f(R). In both scenarios, dark-matter-only simulations of the concordance model corresponding to f_R0=0 and F_0=0 are consistent with the lensing measurements, i.e., at the 68% and 95% confidence level, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Separating intrinsic alignment and galaxy-galaxy lensing

The coherent physical alignment of galaxies is an important systematic for gravitational lensing studies as well as a probe of the physical mechanisms involved in galaxy formation and evolution. We develop a formalism for treating this intrinsic alignment (IA) in the context of galaxy-galaxy lensing and present an improved method for measuring IA contamination, which can arise when sources physically associated with the lens are placed behind the lens due to photometric redshift scatter. We apply the technique to recent Sloan Digital Sky Survey (SDSS) measurements of Luminous Red Galaxy lenses and typical (L*) source galaxies with photometric redshifts selected from the SDSS imaging data. Compared to previous measurements, this method has the advantage of being fully self-consistent in its treatment of the IA and lensing signals, solving for the two simultaneously. We find an IA signal consistent with zero, placing tight constraints on both the magnitude of the IA effect and its potential contamination to the lensing signal. While these constraints depend on source selection and redshift quality, the method can be applied to any measurement that uses photometric redshifts. We obtain a model-independent upper-limit of roughly 10% IA contamination for projected separations of approximately 0.1-100 Mpc/h. With more stringent photo-z cuts and reasonable assumptions about the physics of intrinsic alignments, this upper limit is reduced to 1-2%. These limits are well below the statistical error of the current lensing measurements. Our results suggest that IA will not present intractable challenges to the next generation of galaxy-galaxy lensing experiments, and the methods presented here should continue to aid in our understanding of alignment processes and in the removal of IA from the lensing signal.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Subaru weak-lensing measurement of a z = 0.81 cluster discovered by the Atacama Cosmology Telescope Survey

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\barρ_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the Λ-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Baryon Oscillation Spectroscopic Survey of SDSS-III

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 9 Spectroscopic Galaxy Sample

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5σin both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7σ. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Ninth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-III Baryon Oscillation Spectroscopic Survey

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The SDSS-III Baryon Oscillation Spectroscopic Survey: The Quasar Luminosity Function from Data Release Nine

We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational Lensing Accuracy Testing 2010 (GREAT10) Challenge Handbook

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image analysis challenge aimed at the development of algorithms to analyze astronomical images. Specifically, the challenge is to measure varying image distortions in the presence of a variable convolution kernel, pixelization and noise. This is the second in a series of challenges set to the astronomy, computer science and statistics communities, providing a structured environment in which methods can be improved and tested in preparation for planned astronomical surveys. GREAT10 extends upon previous work by introducing variable fields into the challenge. The "Galaxy Challenge" involves the precise measurement of galaxy shape distortions, quantified locally by two parameters called shear, in the presence of a known convolution kernel. Crucially, the convolution kernel and the simulated gravitational lensing shape distortion both now vary as a function of position within the images, as is the case for real data. In addition, we introduce the "Star Challenge" that concerns the reconstruction of a variable convolution kernel, similar to that in a typical astronomical observation. This document details the GREAT10 Challenge for potential participants. Continually updated information is also available from http://www.greatchallenges.info.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Photometric Redshift Probability Distributions for Galaxies in the SDSS DR8

We present redshift probability distributions for galaxies in the SDSS DR8 imaging data. We used the nearest-neighbor weighting algorithm presented in Lima et al. 2008 and Cunha et al. 2009 to derive the ensemble redshift distribution N(z), and individual redshift probability distributions P(z) for galaxies with r < 21.8. As part of this technique, we calculated weights for a set of training galaxies with known redshifts such that their density distribution in five dimensional color-magnitude space was proportional to that of the photometry-only sample, producing a nearly fair sample in that space. We then estimated the ensemble N(z) of the photometric sample by constructing a weighted histogram of the training set redshifts. We derived P(z) s for individual objects using the same technique, but limiting to training set objects from the local color-magnitude space around each photometric object. Using the P(z) for each galaxy, rather than an ensemble N(z), can reduce the statistical error in measurements that depend on the redshifts of individual galaxies. The spectroscopic training sample is substantially larger than that used for the DR7 release, and the newly added PRIMUS catalog is now the most important training set used in this analysis by a wide margin. We expect the primary source of error in the N(z) reconstruction is sample variance: the training sets are drawn from relatively small volumes of space. Using simulations we estimated the uncertainty in N(z) at a given redshift is 10-15%. The uncertainty on calculations incorporating N(z) or P(z) depends on how they are used; we discuss the case of weak lensing measurements. The P(z) catalog is publicly available from the SDSS website.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Precision simulation of ground-based lensing data using observations from space

Current and upcoming wide-field, ground-based, broad-band imaging surveys promise to address a wide range of outstanding problems in galaxy formation and cosmology. Several such uses of ground-based data, especially weak gravitational lensing, require highly precise measurements of galaxy image statistics with careful correction for the effects of the point-spread function (PSF). In this paper, we introduce the SHERA (SHEar Reconvolution Analysis) software to simulate ground-based imaging data with realistic galaxy morphologies and observing conditions, starting from space-based data (from COSMOS, the Cosmological Evolution Survey) and accounting for the effects of the space-based PSF. This code simulates ground-based data, optionally with a weak lensing shear applied, in a model-independent way using a general Fourier space formalism. The utility of this pipeline is that it allows for a precise, realistic assessment of systematic errors due to the method of data processing, for example in extracting weak lensing galaxy shape measurements or galaxy radial profiles, given user-supplied observational conditions and real galaxy morphologies. Moreover, the simulations allow for the empirical test of error estimates and determination of parameter degeneracies, via generation of many noise maps. The public release of this software, along with a large sample of cleaned COSMOS galaxy images (corrected for charge transfer inefficiency), should enable upcoming ground-based imaging surveys to achieve their potential in the areas of precision weak lensing analysis, galaxy profile measurement, and other applications involving detailed image analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

SDSS-III: Massive Spectroscopic Surveys of the Distant Universe, the Milky Way Galaxy, and Extra-Solar Planetary Systems

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Seeing in the dark -- I. Multi-epoch alchemy

Weak lensing by large-scale structure is an invaluable cosmological tool given that most of the energy density of the concordance cosmology is invisible. Several large ground-based imaging surveys will attempt to measure this effect over the coming decade, but reliable control of the spurious lensing signal introduced by atmospheric turbulence and telescope optics remains a challenging problem. We address this challenge with a demonstration that point-spread function (PSF) effects on measured galaxy shapes in current ground-based surveys can be corrected with existing analysis techniques. In this work, we co-add existing Sloan Digital Sky Survey imaging on the equatorial stripe in order to build a data set with the statistical power to measure cosmic shear, while using a rounding kernel method to null out the effects of the anisotropic PSF. We build a galaxy catalogue from the combined imaging, characterise its photometric properties, and show that the spurious shear remaining in this catalogue after the PSF correction is negligible compared to the expected cosmic shear signal. We identify a new source of systematic error in the shear-shear auto-correlations arising from selection biases related to masking. Finally, we discuss the circumstances in which this method is expected to be useful for upcoming ground-based surveys that have lensing as one of the science goals, and identify the systematic errors that can reduce its efficacy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Seeing in the dark -- II. Cosmic shear in the Sloan Digital Sky Survey

Statistical weak lensing by large-scale structure -- cosmic shear -- is a promising cosmological tool, which has motivated the design of several large upcoming surveys. Here, we present a measurement of cosmic shear using coadded Sloan Digital Sky Survey (SDSS) imaging in 168 square degrees of the equatorial region, with r<23.5 and i<22.5, a source number density of 2.2 galaxies per square arcminute and median redshift of 0.52. These coadds were generated using a new method described in the companion Paper I that was intended to minimise systematic errors in the lensing measurement due to coherent PSF anisotropies that are otherwise prevalent in the SDSS imaging data. We present measurements of cosmic shear out to angular separations of 2 degrees, along with systematics tests that (combined with those from Paper I on the catalogue generation) demonstrate that our results are dominated by statistical rather than systematic errors. Assuming a cosmological model corresponding to WMAP7 and allowing only the amplitude of matter fluctuations to vary, we find a best-fit value of sigma_8=0.636 +0.109 -0.154 (1-sigma); without systematic errors this would be sigma_8=0.636 +0.099 -0.137 (1-sigma). Assuming a flat LCDM model, the combined constraints with WMAP7 are sigma_8=0.784 +0.028 -0.026 (1-sigma), +0.055 -0.054 (2-sigma) and Omega_m h^2=0.1303 +0.0047 -0.0048 (1-sigma)+0.009 -0.009 (2-sigma); the 2-sigma error ranges are respectively 14 and 17 per cent smaller than WMAP7 alone. Aside from the intrinsic value of such cosmological constraints from the growth of structure, we identify some important lessons for upcoming surveys that may face similar issues when combining multi-epoch data to measure cosmic shear.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Eighth Data Release of the Sloan Digital Sky Survey: First Data from SDSS-III

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Morphology of Galaxies in the Baryon Oscillation Spectroscopic Survey

We study the morphology of luminous and massive galaxies at 0.3<z<0.7 targeted in the Baryon Oscillation Spectroscopic Survey (BOSS) using publicly available Hubble Space Telescope imaging from COSMOS. Our sample (240 objects) provides a unique opportunity to check the visual morphology of these galaxies which were targeted based solely on stellar population modelling. We find that the majority (74+/-6%) possess an early-type morphology (elliptical or S0), while the remainder have a late-type morphology. This is as expected from the goals of the BOSS target selection which aimed to predominantly select slowly evolving galaxies, for use as cosmological probes, while still obtaining a fair fraction of actively star forming galaxies for galaxy evolution studies. We show that a colour cut of (g-i)>2.35 selects a sub-sample of BOSS galaxies with 90% early-type morphology - more comparable to the earlier Luminous Red Galaxy (LRG) samples of SDSS-I/II. The remaining 10% of galaxies above this cut have a late-type morphology and may be analogous to the "passive spirals" found at lower redshift. We find that 23+/-4% of the early-type galaxies are unresolved multiple systems in the SDSS imaging. We estimate that at least 50% of these are real associations (not projection effects) and may represent a significant "dry merger" fraction. We study the SDSS pipeline sizes of BOSS galaxies which we find to be systematically larger (by 40%) than those measured from HST images, and provide a statistical correction for the difference. These details of the BOSS galaxies will help users of the data fine-tune their selection criteria, dependent on their science applications. For example, the main goal of BOSS is to measure the cosmic distance scale and expansion rate of the Universe to percent-level precision - a point where systematic effects due to the details of target selection may become important.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An algorithm for the direct reconstruction of the dark matter correlation function from weak lensing and galaxy clustering

The clustering of matter on cosmological scales is an essential probe for studying the physical origin and composition of our Universe. To date, most of the direct studies have focused on shear-shear weak lensing correlations, but it is also possible to extract the dark matter clustering by combining galaxy-clustering and galaxy-galaxy-lensing measurements. In this study we develop a method that can constrain the dark matter correlation function from galaxy clustering and galaxy-galaxy-lensing measurements, by focusing on the correlation coefficient between the galaxy and matter overdensity fields. To generate a mock galaxy catalogue for testing purposes, we use the Halo Occupation Distribution approach applied to a large ensemble of N-body simulations to model pre-existing SDSS Luminous Red Galaxy sample observations. Using this mock catalogue, we show that a direct comparison between the excess surface mass density measured by lensing and its corresponding galaxy clustering quantity is not optimal. We develop a new statistic that suppresses the small-scale contributions to these observations and show that this new statistic leads to a cross-correlation coefficient that is within a few percent of unity down to 5 Mpc/h. Furthermore, the residual incoherence between the galaxy and matter fields can be explained using a theoretical model for scale-dependent bias, giving us a final estimator that is unbiased to within 1%. We also perform a comprehensive study of other physical effects that can affect the analysis, such as redshift space distortions and differences in radial windows between galaxy clustering and weak lensing observations. We apply the method to a range of cosmological models and show the viability of our new statistic to distinguish between cosmological models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Confirmation of general relativity on large scales from weak lensing and galaxy velocities

Although general relativity underlies modern cosmology, its applicability on cosmological length scales has yet to be stringently tested. Such a test has recently been proposed, using a quantity, EG, that combines measures of large-scale gravitational lensing, galaxy clustering and structure growth rate. The combination is insensitive to 'galaxy bias' (the difference between the clustering of visible galaxies and invisible dark matter) and is thus robust to the uncertainty in this parameter. Modified theories of gravity generally predict values of EG different from the general relativistic prediction because, in these theories, the 'gravitational slip' (the difference between the two potentials that describe perturbations in the gravitational metric) is non-zero, which leads to changes in the growth of structure and the strength of the gravitational lensing effect3. Here we report that EG = 0.39 +/- 0.06 on length scales of tens of megaparsecs, in agreement with the general relativistic prediction of EG $\approx$ 0.4. The measured value excludes a model within the tensor-vector-scalar gravity theory, which modifies both Newtonian and Einstein gravity. However, the relatively large uncertainty still permits models within f(R) theory, which is an extension of general relativity. A fivefold decrease in uncertainty is needed to rule out these models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precision cluster mass determination from weak lensing

Weak gravitational lensing has been used extensively in the past decade to constrain the masses of galaxy clusters, and is the most promising observational technique for providing the mass calibration necessary for precision cosmology with clusters. There are several challenges in estimating cluster masses, particularly (a) the sensitivity to astrophysical effects and observational systematics that modify the signal relative to the theoretical expectations, and (b) biases that can arise due to assumptions in the mass estimation method, such as the assumed radial profile of the cluster. All of these challenges are more problematic in the inner regions of the cluster, suggesting that their influence would ideally be suppressed for the purpose of mass estimation. However, at any given radius the differential surface density measured by lensing is sensitive to all mass within that radius, and the corrupted signal from the inner parts is spread out to all scales. We develop a new statistic that is ideal for estimation of cluster masses because it completely eliminates mass contributions below a chosen scale (which we suggest should be about 20 per cent of the virial radius), and thus reduces sensitivity to systematic and astrophysical effects. We use simulated and analytical profiles to quantify systematic biases on the estimated masses for several standard methods of mass estimation, finding that these can lead to significant mass biases that range from ten to over fifty per cent. The mass uncertainties when using our new statistic are reduced by up to a factor of ten relative to the standard methods, while only moderately increasing the statistical errors. This new method of mass estimation will enable a higher level of precision in future science work with weak lensing mass estimates for galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Testing adiabatic contraction with SDSS elliptical galaxies

We study the profiles of 75 086 elliptical galaxies from the Sloan Digital Sky Survey (SDSS) at both large (50-500 kpc/h) and small (~3 kpc/h) scales. Weak lensing observations in the outskirts of the halo are combined with measurements of the stellar velocity dispersion in the interior regions of the galaxy for stacked galaxy samples. The weak lensing measurements are well characterized by a Navarro, Frenk and White (NFW) profile. The dynamical mass measurements exceed the extrapolated NFW profile even after the estimated stellar masses are subtracted, providing evidence for the modification of the dark matter profile by the baryons. This excess mass is quantitatively consistent with the predictions of the adiabatic contraction (AC) hypothesis. Our finding suggests that the effects of AC during galaxy formation are stable to subsequent bombardment from major and minor mergers. We explore several theoretical and observational systematics and conclude that they cannot account for the inferred mass excess. The most significant source of systematic error is in the IMF, which would have to increase the stellar mass estimates by a factor of two relative to masses from the Kroupa IMF to fully explain the mass excess without AC. Such an increase would create tension with results from SAURON (Cappellari et al. 2006). We demonstrate a connection between the level of contraction of the dark matter halo profile and scatter in the size-luminosity relation, which is a projection of the fundamental plane. Whether or not AC is the mechanism supplying the excess mass, models of galaxy formation and evolution must reconcile the observed halo masses from weak lensing with the comparatively large dynamical masses at the half light radii of the galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The WiggleZ Dark Energy Survey: Direct constraints on blue galaxy intrinsic alignments at intermediate redshifts

Correlations between the intrinsic shapes of galaxy pairs, and between the intrinsic shapes of galaxies and the large-scale density field, may be induced by tidal fields. These correlations, which have been detected at low redshifts (z<0.35) for bright red galaxies in the Sloan Digital Sky Survey (SDSS), and for which upper limits exist for blue galaxies at z~0.1, provide a window into galaxy formation and evolution, and are also an important contaminant for current and future weak lensing surveys. Measurements of these alignments at intermediate redshifts (z~0.6) that are more relevant for cosmic shear observations are very important for understanding the origin and redshift evolution of these alignments, and for minimising their impact on weak lensing measurements. We present the first such intermediate-redshift measurement for blue galaxies, using galaxy shape measurements from SDSS and spectroscopic redshifts from the WiggleZ Dark Energy Survey. Our null detection allows us to place upper limits on the contamination of weak lensing measurements by blue galaxy intrinsic alignments that, for the first time, do not require significant model-dependent extrapolation from the z~0.1 SDSS observations. Also, combining the SDSS and WiggleZ constraints gives us a long redshift baseline with which to constrain intrinsic alignment models and contamination of the cosmic shear power spectrum. Assuming that the alignments can be explained by linear alignment with the smoothed local density field, we find that a measurement of σ_8 in a blue-galaxy dominated, CFHTLS-like survey would be contaminated by at most +/-0.02 (95% confidence level, SDSS and WiggleZ) or +/-0.03 (WiggleZ alone) due to intrinsic alignments. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Galaxy density profiles and shapes -- I. simulation pipeline for lensing by realistic galaxy models

Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, whch we aim to investigate in a consistent way. In this first paper in a series we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Galaxy density profiles and shapes -- II. selection biases in strong lensing surveys

[Abridged] Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio, and image separation distribution. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sersic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.

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

Results of the GREAT08 Challenge: An image analysis competition for cosmological lensing

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Halo mass - concentration relation from weak lensing

We perform a statistical weak lensing analysis of dark matter profiles around tracers of halo mass from galactic- to cluster-size halos. In this analysis we use 170,640 isolated ~L* galaxies split into ellipticals and spirals, 38,236 groups traced by isolated spectroscopic Luminous Red Galaxies (LRGs) and 13,823 MaxBCG clusters from the Sloan Digital Sky Survey (SDSS) covering a wide range of richness. Together these three samples allow a determination of the density profiles of dark matter halos over three orders of magnitude in mass, from 10^{12} M_{sun} to 10^{15} M_{sun}. The resulting lensing signal is consistent with an NFW or Einasto profile on scales outside the central region. We find that the NFW concentration parameter c_{200b} decreases with halo mass, from around 10 for galactic halos to 4 for cluster halos. Assuming its dependence on halo mass in the form of c_{200b} = c_0 [M/(10^{14}M_{sun}/h)]^β, we find c_0=4.6 +/- 0.7 (at z=0.22) and β=0.13 +/- 0.07, with very similar results for the Einasto profile. The slope (β) is in agreement with theoretical predictions, while the amplitude is about two standard deviations below the predictions for this mass and redshift, but we note that the published values in the literature differ at a level of 10-20% and that for a proper comparison our analysis should be repeated in simulations. We discuss the implications of our results for the baryonic effects on the shear power spectrum: since these are expected to increase the halo concentration, the fact that we see no evidence of high concentrations on scales above 20% of the virial radius suggests that baryonic effects are limited to small scales, and are not a significant source of uncertainty for the current weak lensing measurements of the dark matter power spectrum. [ABRIDGED]

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A robust lower limit on the amplitude of matter fluctuations in the universe from cluster abundance and weak lensing

Cluster abundance measurements are among the most sensitive probes of the amplitude of matter fluctuations in the universe, which in turn can help constrain other cosmological parameters, like the dark energy equation of state or neutrino mass. However, difficulties in calibrating the relation between the cluster observable and halo mass, and the lack of completeness information, make this technique particularly susceptible to systematic errors. Here we argue that a cluster abundance analysis using statistical weak lensing on the stacked clusters leads to a robust lower limit on the amplitude of fluctuations. The method compares the average weak lensing signal measured around the whole cluster sample to a theoretical prediction, assuming that the clusters occupy the centers of all of the most massive halos above some minimum mass threshold. If the amplitude of fluctuations is below a certain limiting value, there are too few massive clusters in this model and the theoretical prediction falls below the observations. Since any effects that modify the model assumptions can only decrease the theoretical prediction, the limiting amplitude becomes a robust lower limit. Here, we apply it to a volume limited sample of 16,000 group/cluster candidates identified from isolated luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS). We find $σ_8 (Ω_m/0.25)^{0.5}>0.62$ at the 95% c.l. after taking into account observational errors in the lensing analysis. While this is a relatively weak constraint, both the scatter in the LRG luminosity-halo mass relation and the lensing errors are large; the constraints could improve considerably in the future with more sophisticated cluster identification algorithms and smaller errors in the lensing analysis. [Abridged]

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