Search

We present a novel simulation-based cosmological analysis of galaxy-galaxy lensing and galaxy redshift-space clustering. Compared to analysis methods based on perturbation theory, our simulation-based approach allows us to probe a much wider range of scales, $0.4 \, h^{-1} \, \mathrm{Mpc}$ to $63 \, h^{-1} \, \mathrm{Mpc}$, including highly non-linear scales, and marginalises over astrophysical effects such as assembly bias. We apply this framework to data from the Baryon Oscillation Spectroscopic Survey LOWZ sample cross-correlated with state-of-the-art gravitational lensing catalogues from the Kilo Degree Survey and the Dark Energy Survey. We show that gravitational lensing and redshift-space clustering when analysed over a large range of scales place tight constraints on the growth-of-structure parameter $S_8 = \sigma_8 \sqrt{\Omega_{\rm m} / 0.3}$. Overall, we infer $S_8 = 0.792 \pm 0.022$ when analysing the combination of galaxy-galaxy lensing and projected galaxy clustering and $S_8 = 0.771 \pm 0.027$ for galaxy redshift-space clustering. These findings highlight the potential constraining power of full-scale studies over studies analysing only large scales, and also showcase the benefits of analysing multiple large-scale structure surveys jointly. Our inferred values for $S_8$ fall below the value inferred from the CMB, $S_8 = 0.834 \pm 0.016$. While this difference is not statistically significant by itself, our results mirror other findings in the literature whereby low-redshift large scale structure probes infer lower values for $S_8$ than the CMB, the so-called $S_8$-tension., Comment: 22 pages, 16 figures, submitted to MNRAS, comments welcome

The Merian survey is mapping $\sim$ 850 degrees$^2$ of the Hyper Suprime-Cam Strategic Survey Program (HSC-SSP) wide layer with two medium-band filters on the 4-meter Victor M. Blanco telescope at the Cerro Tololo Inter-American Observatory, with the goal of carrying the first high signal-to-noise (S/N) measurements of weak gravitational lensing around dwarf galaxies. This paper presents the design of the Merian filter set: N708 ($\lambda_c = 7080 \unicode{x212B}$, $\Delta\lambda = 275\unicode{x212B}$) and N540 ($\lambda_c = 5400\unicode{x212B}$, $\Delta\lambda = 210\unicode{x212B}$). The central wavelengths and filter widths of N708 and N540 were designed to detect the $\rm H\alpha$ and $\rm [OIII]$ emission lines of galaxies in the mass range $8, Comment: 18 pages, 14 figures, submitted to MNRAS

The Dark Energy Spectroscopic Instrument (DESI) was designed to conduct a survey covering 14,000 deg$^2$ over five years to constrain the cosmic expansion history through precise measurements of Baryon Acoustic Oscillations (BAO). The scientific program for DESI was evaluated during a five month Survey Validation (SV) campaign before beginning full operations. This program produced deep spectra of tens of thousands of objects from each of the stellar (MWS), bright galaxy (BGS), luminous red galaxy (LRG), emission line galaxy (ELG), and quasar target classes. These SV spectra were used to optimize redshift distributions, characterize exposure times, determine calibration procedures, and assess observational overheads for the five-year program. In this paper, we present the final target selection algorithms, redshift distributions, and projected cosmology constraints resulting from those studies. We also present a `One-Percent survey' conducted at the conclusion of Survey Validation covering 140 deg$^2$ using the final target selection algorithms with exposures of a depth typical of the main survey. The Survey Validation indicates that DESI will be able to complete the full 14,000 deg$^2$ program with spectroscopically-confirmed targets from the MWS, BGS, LRG, ELG, and quasar programs with total sample sizes of 7.2, 13.8, 7.46, 15.7, and 2.87 million, respectively. These samples will allow exploration of the Milky Way halo, clustering on all scales, and BAO measurements with a statistical precision of 0.28% over the redshift interval $z, Comment: 42 pages, 16 figures, submitted to AJ

We construct accurate emulators for the projected and redshift space galaxy correlation functions and excess surface density as measured by galaxy-galaxy lensing, based on Halo Occupation Distribution (HOD) modeling. Using the complete Mira-Titan suite of 111 $N$-body simulations, our emulators vary over eight cosmological parameters and include the effects of neutrino mass and dynamical dark energy. We demonstrate that our emulators are sufficiently accurate for the analysis of the BOSS DR12 CMASS galaxy sample over the range 0.5 < r < 50 Mpc/h. Furthermore, we show that our emulators are capable of recovering unbiased cosmological constraints from realistic mock catalogs over the same range. Our mock catalog tests show the efficacy of combining small scale galaxy-galaxy lensing with redshift space clustering and that we can constrain the growth rate and σ_8 to 7% and 4.5% respectively for a CMASS-like sample using only the measurements covered by our emulator. With the inclusion of a CMB prior on H_0, this reduces to a 2% measurement on the growth rate., 31 pages, 13 figures, revised to match published version

We propose that observations of super-massive galaxies contain cosmological constraining power similar to conventional cluster cosmology, and we provide promising indications that the associated systematic errors are comparably easier to control. We consider a fiducial spectroscopic and stellar mass complete sample of galaxies drawn from the Dark Energy Spectroscopic Survey (DESI) and forecast how constraints on Omega_m-sigma_8 from this sample will compare with those from number counts of clusters based on richness. At fixed number density, we find that massive galaxies offer similar constraints to galaxy clusters. However, a mass-complete galaxy sample from DESI has the potential to probe lower halo masses than standard optical cluster samples (which are typically limited to richness above 20 and halo mass above 10^13.5); additionally, it is straightforward to cleanly measure projected galaxy clustering for such a DESI sample, which we show can substantially improve the constraining power on Omega_m. We also compare the constraining power of stellar mass-limited samples to those from larger but mass-incomplete samples (e.g., the DESI Bright Galaxy Survey, BGS, Sample); relative to a lower number density stellar mass-limited samples, we find that a BGS-like sample improves statistical constraints by 60% for Omega_m and 40% for sigma_8, but this uses small scale information which will be harder to model for BGS. Our initial assessment of the systematics associated with supermassive galaxy cosmology yields promising results. The proposed samples have a 10% satellite fraction, but we show that cosmological constraints may be robust to the impact of satellites. These findings motivate future work to realize the potential of super-massive galaxies to probe lower halo masses than richness-based clusters and to avoid persistent systematics associated with optical cluster finding.

A. Leauthaud et al., Lensing without borders is a cross-survey collaboration created to assess the consistency of galaxy–galaxy lensing signals (ΔΣ) across different data sets and to carry out end-to-end tests of systematic errors. We perform a blind comparison of the amplitude of ΔΣ using lens samples from BOSS and six independent lensing surveys. We find good agreement between empirically estimated and reported systematic errors which agree to better than 2.3σ in four lens bins and three radial ranges. For lenses with zL > 0.43 and considering statistical errors, we detect a 3–4σ correlation between lensing amplitude and survey depth. This correlation could arise from the increasing impact at higher redshift of unrecognized galaxy blends on shear calibration and imperfections in photometric redshift calibration. At zL > 0.54, amplitudes may additionally correlate with foreground stellar density. The amplitude of these trends is within survey-defined systematic error budgets that are designed to include known shear and redshift calibration uncertainty. Using a fully empirical and conservative method, we do not find evidence for large unknown systematics. Systematic errors greater than 15 per cent (25 per cent) ruled out in three lens bins at 68 per cent (95 per cent) confidence at z < 0.54. Differences with respect to predictions based on clustering are observed to be at the 20–30 per cent level. Our results therefore suggest that lensing systematics alone are unlikely to fully explain the ‘lensing is low’ effect at z < 0.54. This analysis demonstrates the power of cross-survey comparisons and provides a promising path for identifying and reducing systematics in future lensing analyses., We acknowledge use of the lux supercomputer at UC Santa Cruz, funded by NSF MRI grant AST 1828315. This material is based on work supported by the U.D Department of Energy, Office of Science, Office of High Energy Physics under Award Number DE-SC0019301. AL acknowledges support from the David and Lucille Packard foundation, and from the Alfred.P Sloan foundation. CH acknowledges support from the European Research Council under grant number 647112, and support from the Max Planck Society and the Alexander von Humboldt Foundation in the framework of the Max Planck-Humboldt Research Award endowed by the Federal Ministry of Education and Research. KK acknowledges support from the Royal Society and Imperial College. HH is supported by a Heisenberg grant of the Deutsche Forschungsgemeinschaft (Hi 1495/5-1), as well as an ERC Consolidator Grant (no. 770935). AHW is supported by an European Research Council Consolidator Grant (no. 770935). Funding for the DES Projects has been provided by the U.S. Department of Energy; the U.S. National Science Foundation; the Ministry of Science and Education of Spain; the Science and Technology Facilities Council of the United Kingdom; the Higher Education Funding Council for England; the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign; the Kavli Institute of Cosmological Physics at the University of Chicago; the Center for Cosmology and Astro-Particle Physics at the Ohio State University; the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University; Financiadora de Estudos e Projetos; Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro; Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência; Tecnologia e Inovação; the Deutsche Forschungsgemeinschaft; and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory; the University of California at Santa Cruz; the University of Cambridge; Centro de Investigaciones Energéticas; Medioambientales y Tecnológicas-Madrid; the University of Chicago; University College London; the DES-Brazil Consortium; the University of Edinburgh; the Eidgenössische Technische Hochschule (ETH) Zürich; Fermi National Accelerator Laboratory; the University of Illinois at Urbana-Champaign; the Institut de Ciències de l’Espai (IEEC/CSIC); the Institut de Física d’Altes Energies; Lawrence Berkeley National Laboratory; the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe; the University of Michigan; NFS’s NOIRLab; the University of Nottingham; The Ohio State University; the University of Pennsylvania; the University of Portsmouth; SLAC National Accelerator Laboratory; Stanford University; the University of Sussex; Texas A&M University; and the OzDES Membership Consortium. Based in part on observations at Cerro Tololo Inter-American Observatory at NSF’s NOIRLab (NOIRLab Prop. ID 2012B-0001; PI: J. Frieman), which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The DES data management system is supported by the National Science Foundation under grant numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MICINN under grants ESP2017-89838, PGC2018-094773, PGC2018-102021, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciência e Tecnologia (INCT) do e-Universo (CNPq grant 465376/2014-2). This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics. The Hyper Suprime-Cam (HSC) collaboration includes the astronomical communities of Japan and Taiwan, and Princeton University. The HSC instrumentation and software were developed by the National Astronomical Observatory of Japan (NAOJ), the Kavli Institute for the Physics and Mathematics of the Universe (Kavli IPMU), the University of Tokyo, the High Energy Accelerator Research Organisation (KEK), the Academia Sinica Institute for Astronomy and Astrophysics in Taiwan (ASIAA), and Princeton University. Funding was contributed by the FIRST program from the Japanese Cabinet Office, the Ministry of Education, Culture, Sports, Science, and Technology (MEXT), the Japan Society for the Promotion of Science (JSPS), Japan Science and Technology Agency (JST), the Toray Science Foundation, NAOJ, Kavli IPMU, KEK, ASIAA, and Princeton University. This paper makes use of software developed for the Large Synoptic Survey Telescope. We thank the LSST Project for making their code available as free software at http://dm.lsst.org This paper is based [in part] on data collected at the Subaru Telescope and retrieved from the HSC data archive system, which is operated by Subaru Telescope and Astronomy Data Center (ADC) at NAOJ. Data analysis was in part carried out with the cooperation of Center for Computational Astrophysics (CfCA), NAOJ. The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg, and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation grant no. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation. Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme IDs 177.A-3016, 177.A-3017, 177.A-3018, and 179.A-2004, and on data products produced by the KiDS consortium. The KiDS production team acknowledges support from: Deutsche Forschungsgemeinschaft, ERC, NOVA, and NWO-M grants; Target; the University of Padova, and the University Federico II (Naples). This work was supported by the Department of Energy, Laboratory Directed Research and Development program at SLAC National Accelerator Laboratory, under contract DE-AC02-76SF00515 and as a part of the Panofsky Fellowship awarded to DG. MM is partially funded by FAPERJ, CNPq, and CONICET. BM acknowledges support from the Brazilian funding agency FAPERJ.

The evolution of the structural and kinematic properties of early-type galaxies (ETGs), their scaling relations, as well as their stellar metallicity and age contain precious information on the assembly history of these systems. We present results on the evolution of the stellar mass-velocity dispersion relation of ETGs, focusing in particular on the effects of some selection criteria used to define ETGs. We also try to shed light on the role that in-situ and ex-situ stellar populations have in massive ETGs, providing a possible explanation of the observed metallicity distributions., Comment: 5 pages, 2 figures. To appear in Proceedings of IAU Symposium 359 (T. Storchi-Bergmann, R. Overzier, W. Forman & R. Riffel, eds.)

The outer light (stellar halos) of massive galaxies has recently emerged as a possible low scatter tracer of dark matter halo mass. To test the robustness of outer light measurements across different data sets, we compare the surface brightness profiles of massive galaxies using four independent data sets: the Hyper Suprime-Cam survey (HSC), the Dark Energy Camera Legacy Survey (DECaLS), the Sloan Digital Sky Survey (SDSS), and the Dragonfly Wide Field Survey (Dragonfly). We use customized pipelines for HSC and DECaLS to achieve better sky background subtraction. For galaxies at $z 200$ kpc without significant bias. At $0.19, Comment: 19 pages, 10 figures. Submitted to MNRAS. Main results are shown in Fig 5-10

We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($\Delta\!\Sigma$) and projected galaxy clustering ($w_{\rm p}$) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range $0.15, Comment: 44 pages, 30 figures, 6 tables, submitted to Phys. Rev. D, YouTube video is available at https://www.youtube.com/watch?v=jRhQvuee97Q

Secondary halo properties beyond mass, such as the mass accretion rate (MAR), concentration, and the half mass scale, are essential in understanding the formation of large-scale structure and dark matter halos. In this paper, we study the impact of secondary halo properties on the galaxy-galaxy lensing observable, $\Delta\Sigma$. We build an emulator trained on N-body simulations to model $\Delta\Sigma$ and quantify the impact of different secondary parameters on the $\Delta\Sigma$ profile. We focus on the impact of MAR on $\Delta\Sigma$. We show that a 3$\sigma$ detection of variations in MAR at fixed halo mass could be achieved with the Hyper Suprime Cam survey in combination with a proxy for MAR with scatter $\sigma_{\Gamma_\mathrm{dyn}|\mathrm{obs}}

The canonical $\Lambda$CDM cosmological model makes precise predictions for the clustering and lensing properties of galaxies. It has been shown that the lensing amplitude of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) is lower than expected given their clustering properties. We present new measurements and modelling of galaxies in the BOSS LOWZ sample. We focus on the radial and stellar mass dependence of the lensing amplitude mis-match. We find an amplitude mis-match of around $35\%$ when assuming $\Lambda$CDM with Planck Cosmological Microwave Background (CMB) constraints. This offset is independent of halo mass and radial scale in the range $M_{\rm halo}\sim 10^{13.3} - 10^{13.9} h^{-1} M_\odot$ and $r=0.1 - 60 \, h^{-1} \mathrm{Mpc}$ ($k \approx 0.05 - 20 \, h \, {\rm Mpc}^{-1}$). The observation that the offset is both mass and scale independent places important constraints on the degree to which astrophysical processes (baryonic effects, assembly bias) can fully explain the effect. This scale independence also suggests that the "lensing is low" effect on small and large radial scales probably have the same physical origin. Resolutions based on new physics require a nearly uniform suppression, relative to $\Lambda$CDM predictions, of the amplitude of matter fluctuations on these scales. The possible causes of this are tightly constrained by measurements of the CMB and of the low-redshift expansion history., Comment: 15 pages, 7 figures, matches version accepted for publication in MNRAS

We present an anomaly detection method using Wasserstein generative adversarial networks (WGANs) on optical galaxy images from the wide-field survey conducted with the Hyper Suprime-Cam (HSC) on the Subaru Telescope in Hawai'i. The WGAN is trained on the entire sample, and learns to generate realistic HSC-like images that follow the distribution of the training data. We identify images which are less well-represented in the generator's latent space, and which the discriminator flags as less realistic; these are thus anomalous with respect to the rest of the data. We propose a new approach to characterize these anomalies based on a convolutional autoencoder (CAE) to reduce the dimensionality of the residual differences between the real and WGAN-reconstructed images. We construct a subsample of ~9,000 highly anomalous images from our nearly million object sample, and further identify interesting anomalies within these; these include galaxy mergers, tidal features, and extreme star-forming galaxies. The proposed approach could boost unsupervised discovery in the era of big data astrophysics., Comment: Accepted to the 2020 NeurIPS Machine Learning and the Physical Science Workshop

We recommend a conceptual design study for a spectroscopic facility in the southern hemisphere comprising a large diameter telescope, fiber system, and spectrographs collectively optimized for massively-multiplexed spectroscopy. As a baseline, we propose an 11.4-meter aperture, optical spectroscopic survey telescope with a five square degree field of view. Using current technologies, the facility could be equipped with 15,000 robotically-controlled fibers feeding spectrographs over 360, Comment: Submitted to the Astro2020 Decadal Survey as a facilities white paper

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2

We use the UniverseMachine to analyze the source of scatter between the central galaxy mass, the total stellar mass in the halo, and the dark matter halo mass. We also propose a new halo mass estimator, the cen+N mass: the sum of the stellar mass of the central and the N most massive satellites. We show that, when real space positions are perfectly known, the cen+N mass has scatter competitive with that of richness-based estimators. However, in redshift space, the cen+N mass suffers less from projection effects in the UniverseMachine model. The cen+N mass is therefore a viable low scatter halo mass estimator, and should be considered an important tool to constrain cosmology with upcoming spectroscopic data from DESI. We analyze the scatter in stellar mass at fixed halo mass and show that the total stellar mass in a halo is uncorrelated with secondary halo properties, but that the central stellar mass is a function of both halo mass and halo age. This is because central galaxies in older halos have had more time to grow via accretion. If the UniverseMachine model is correct, accurate galaxy-halo modeling of mass selected samples therefore needs to consider halo age in addition to mass., Comment: 13 pages, 11 figures, submitted to MNRAS

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov., Comment: 9-page APC White Paper submitted to the Astro2020 Decadal Survey. To be published in BAAS. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov

The advent of new deep+wide photometric lensing surveys will open up the possibility of direct measurements of the dark matter halos of dwarf galaxies. The HSC wide survey will be the first with the statistical capability of measuring the lensing signal with high signal-to-noise at log(M*)=8. At this same mass scale, LSST will have the most overall constraining power with a predicted signal-to-noise for the galaxy-galaxy lensing signal around dwarfs of SN=200. WFIRST and LSST will have the greatest potential to push below the log(M*) = 7 mass scale thanks to the depth of their imaging data. Studies of the dark matter halos of dwarf galaxies at z=0.1 with gravitational lensing are soon within reach. However, further work will be required to develop optimized strategies for extracting dwarfs samples from these surveys, determining redshifts, and accurately measuring lensing on small radial scales. Dwarf lensing will be a new and powerful tool to constrain the halo masses and inner density slopes of dwarf galaxies and to distinguish between baryonic feedback and modified dark matter scenarios., Comment: 6 pages, 2 figures, submitted to MNRAS

We present the clustering properties of low-$z$ $(z\leq1.4)$ galaxies selected by the Hyper Suprime-Cam Subaru Strategic Program Wide layer over $145$ deg$^{2}$. The wide-field and multi-wavelength observation yields $5,064,770$ galaxies at $0.3\leq z\leq1.4$ with photometric redshifts and physical properties. This enables the accurate measurement of angular correlation functions and subsequent halo occupation distribution (HOD) analysis allows the connection between baryonic properties and dark halo properties. The fraction of less-massive satellite galaxies at $z\lesssim1$ is found to be almost constant at $\sim20\%$, but it gradually decreases beyond $M_{\star} \sim 10^{10.4}h^{-2}M_{\odot}$. However, the abundance of satellite galaxies at $z>1$ is quite small even for less-massive galaxies due to the rarity of massive centrals at high-$z$. This decreasing trend is connected to the small satellite fraction of Lyman break galaxies at $z>3$. The stellar-to-halo mass ratios at $0.3\leq z\leq1.4$ are almost consistent with the predictions obtained using the latest empirical model; however, we identify small excesses from the theoretical model at the massive end. The pivot halo mass is found to be unchanged at $10^{11.9-12.1}h^{-1}M_{\odot}$ at $0.3\leq z\leq1.4$, and we systematically show that $10^{12}h^{-1}M_{\odot}$ is a universal pivot halo mass up to $z\sim5$ that is derived using only the clustering/HOD analyses. Nevertheless, halo masses with peaked instantaneous baryon conversion efficiencies are much smaller than the pivot halo mass regardless of a redshift, and the most efficient stellar-mass assembly is thought to be in progress in $10^{11.0-11.5}h^{-1}M_{\odot}$ dark haloes., 33 pages, 16 figures, submitted to ApJ

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i, Comment: 43 pages, 21 figures, accepted for publication in PASJ

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements., Comment: This document is an annual report. A higher version of this document can be found here (http://www.wfirst-hls-cosmology.org/products/). This document does not constitute an official WFIRST Project document

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations., Comment: 39 pages, 21 figures, 2 tables. Submitted to Publications of the Astronomical Society of Japan

Galaxy cluster centring is a key issue for precision cosmology studies using galaxy surveys. Mis-identification of central galaxies causes systematics in various studies such as cluster lensing, satellite kinematics, and galaxy clustering. The red-sequence Matched-filter Probabilistic Percolation (redMaPPer) estimates the probability that each member galaxy is central from photometric information rather than specifying one central galaxy. The redMaPPer estimates can be used for calibrating the off-centring effect, however, the centring algorithm has not previously been well-tested. We test the centring probabilities of redMaPPer cluster catalog using the projected cross correlation between redMaPPer clusters with photometric red galaxies and galaxy-galaxy lensing. We focus on the subsample of redMaPPer clusters in which the redMaPPer central galaxies (RMCGs) are not the brightest member galaxies (BMEM) and both of them have spectroscopic redshift. This subsample represents nearly 10% of the whole cluster sample. We find a clear difference in the cross-correlation measurements between RMCGs and BMEMs, and the estimated centring probability is 74$\pm$10% for RMCGs and 13$\pm$4% for BMEMs in the Gaussian offset model and 78$\pm$9% for RMCGs and 5$\pm$5% for BMEMs in the NFW offset model. These values are in agreement with the centring probability values reported by redMaPPer (75% for RMCG and 10% for BMEMs) within 1$\sigma$. Our analysis provides a strong consistency test of the redMaPPer centring probabilities. Our results suggest that redMaPPer centring probabilities are reliably estimated. We confirm that the brightest galaxy in the cluster is not always the central galaxy as has been shown in previous works., Comment: 9 pages, 4 figures, 3 tables

The Subaru Strategic Program (SSP) is an ambitious multi-band survey using the Hyper Suprime-Cam (HSC) on the Subaru telescope. The Wide layer of the SSP is both wide and deep, reaching a detection limit of i~26.0 mag. At these depths, it is challenging to achieve accurate, unbiased, and consistent photometry across all five bands. The HSC data are reduced using a pipeline that builds on the prototype pipeline for the Large Synoptic Survey Telescope. We have developed a Python-based, flexible framework to inject synthetic galaxies into real HSC images called SynPipe. Here we explain the design and implementation of SynPipe and generate a sample of synthetic galaxies to examine the photometric performance of the HSC pipeline. For stars, we achieve 1% photometric precision at i~19.0 mag and 6% precision at i~25.0 in the i-band. For synthetic galaxies with single-Sersic profiles, forced CModel photometry achieves 13% photometric precision at i~20.0 mag and 18% precision at i~25.0 in the i-band. We show that both forced PSF and CModel photometry yield unbiased color estimates that are robust to seeing conditions. We identify several caveats that apply to the version of HSC pipeline used for the first public HSC data release (DR1) that need to be taking into consideration. First, the degree to which an object is blended with other objects impacts the overall photometric performance. This is especially true for point sources. Highly blended objects tend to have larger photometric uncertainties, systematically underestimated fluxes and slightly biased colors. Second, >20% of stars at 22.5< i < 25.0 mag can be misclassified as extended objects. Third, the current CModel algorithm tends to strongly underestimate the half-light radius and ellipticity of galaxy with i>21.5 mag., Submitted to the PASJ special issue for Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP); 26 pages, 17 figures, 4 tables. The SynPipe code can be found here: https://github.com/dr-guangtou/synpipe

We report the discovery of a diffuse stellar cloud with an angular extent $\gtrsim30^{\prime\prime}$, which we term "Sumo Puff", in data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). While we do not have a redshift for this object, it is in close angular proximity to a post-merger galaxy at redshift $z=0.0431$ and is projected within a few virial radii (assuming similar redshifts) of two other ${\sim}L_\star$ galaxies, which we use to bracket a potential redshift range of $0.0055 < z < 0.0431$. The object's light distribution is flat, as characterized by a low Sersic index ($n\sim0.3$). It has a low central $g$-band surface brightness of ${\sim}26.4$ mag arcsec$^{-2}$, large effective radius of ${\sim}13^{\prime\prime}$ (${\sim}11$ kpc at $z=0.0431$ and ${\sim}1.5$ kpc at $z=0.0055$), and an elongated morphology ($b/a\sim0.4$). Its red color ($g-i\sim1$) is consistent with a passively evolving stellar population and similar to the nearby post-merger galaxy, and we may see tidal material connecting Sumo Puff with this galaxy. We offer two possible interpretations for the nature of this object: (1) it is an extreme, galaxy-size tidal feature associated with a recent merger event, or (2) it is a foreground dwarf galaxy with properties consistent with a quenched, disturbed ultra-diffuse galaxy. We present a qualitative comparison with simulations that demonstrates the feasibility of forming a structure similar to this object in a merger event. Follow-up spectroscopy and/or deeper imaging to confirm the presence of the bridge of tidal material will be necessary to reveal the true nature of this object., 10 pages, 5 figures, submitted to PASJ for the HSC-SSP special issue

We present wide-field (167 deg$^2$) weak lensing mass maps from the Hyper Supreme-Cam Subaru Strategic Program (HSC-SSP). We compare these weak lensing based dark matter maps with maps of the distribution of the stellar mass associated with luminous red galaxies. We find a strong correlation between these two maps with a correlation coefficient of $\rho=0.54\pm0.03$ (for a smoothing size of $8'$). This correlation is detected even with a smaller smoothing scale of $2'$ ($\rho=0.34\pm 0.01$). This detection is made uniquely possible because of the high source density of the HSC-SSP weak lensing survey ($\bar{n}\sim 25$ arcmin$^{-2}$). We also present a variety of tests to demonstrate that our maps are not significantly affected by systematic effects. By using the photometric redshift information associated with source galaxies, we reconstruct a three-dimensional mass map. This three-dimensional mass map is also found to correlate with the three-dimensional galaxy mass map. Cross-correlation tests presented in this paper demonstrate that the HSC-SSP weak lensing mass maps are ready for further science analyses., Comment: 14 pages, 13 figures, accepted for publication in PASJ

Massive galaxies display extended light profiles that can reach several hundreds of kilo parsecs. These stellar halos provide a fossil record of galaxy assembly histories. Using data that is both wide (~100 square degree) and deep (i>28.5 mag/arcsec^2 in i-band), we present a systematic study of the stellar halos of a sample of more than 3000 galaxies at 0.3 < z < 0.5 with $\log M_{\star}/M_{\odot} > 11.4$. Our study is based on high-quality (0.6 arcsec seeing) imaging data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP), which enables us to individually estimate surface mass density profiles to 100 kpc without stacking. As in previous work, we find that more massive galaxies exhibit more extended outer profiles. When this extended light is not properly accounted for as a result of shallow imaging or inadequate profile modeling, the derived stellar mass function can be significantly underestimated at the highest masses. Across our sample, the ellipticity of outer light profiles increases substantially as we probe larger radii. We show for the first time that these ellipticity gradients steepen dramatically as a function of galaxy mass, but we detect no mass-dependence in outer color gradients. Our results support the two-phase formation scenario for massive galaxies in which outer envelopes are built up at late times from a series of merging events. We provide surface mass surface mass density profiles in a convenient tabulated format to facilitate comparisons with predictions from numerical simulations of galaxy formation., Submitted to MNRAS; 23 pages, 8 figures, 2 appendix; Data will be made available here: http://massivegalaxies.com/ once the paper is published

We present results from a set of simulations designed to constrain the weak lensing shear calibration for the Hyper Suprime-Cam (HSC) survey. These simulations include HSC observing conditions and galaxy images from the Hubble Space Telescope (HST), with fully realistic galaxy morphologies and the impact of nearby galaxies included. We find that the inclusion of nearby galaxies in the images is critical to reproducing the observed distributions of galaxy sizes and magnitudes, due to the non-negligible fraction of unrecognized blends in ground-based data, even with the excellent typical seeing of the HSC survey (0.58" in the $i$-band). Using these simulations, we detect and remove the impact of selection biases due to the correlation of weights and the quantities used to define the sample (S/N and apparent size) with the lensing shear. We quantify and remove galaxy property-dependent multiplicative and additive shear biases that are intrinsic to our shear estimation method, including a $\sim 10$ per cent-level multiplicative bias due to the impact of nearby galaxies and unrecognized blends. Finally, we check the sensitivity of our shear calibration estimates to other cuts made on the simulated samples, and find that the changes in shear calibration are well within the requirements for HSC weak lensing analysis. Overall, the simulations suggest that the weak lensing multiplicative biases in the first-year HSC shear catalog are controlled at the 1 per cent level., Comment: 29 pages, 18 figures, 3 tables. Matches MNRAS accepted version. The Hubble Space Telescope postage stamp images used as the inputs to the simulations were publicly released at https://hsc-release.mtk.nao.ac.jp/doc/index.php/weak-lensing-simulation-catalog-pdr1/

We present initial results from the Subaru Strategic Program (SSP) with Hyper Suprime-Cam (HSC) on a comprehensive survey of emission-line galaxies at z5\times10^5$ Mpc$^3$ is essential to overcome cosmic variance. Since the current data have not reached the full depth expected for the HSC-SSP, the color cut in i-NB816 or z-NB921 induces a bias towards star-forming galaxies with large equivalent widths, primarily seen in the stellar mass functions for the H$��$ emitters at z ~ 0.25--0.40. Even so, the emission-line galaxies clearly cover a wide range of luminosity, stellar mass, and environment, thus demonstrating the usefulness of the NB data from the HSC-SSP to investigate star-forming galaxies at z, 28 pages, 22 figures, 5 tables, accepted for publication in PASJ (Subaru HSC-SSP special issue), The catalogs will be available at https://hsc-release.mtk.nao.ac.jp/

We present and characterize the catalog of galaxy shape measurements that will be used for cosmological weak lensing measurements in the Wide layer of the first year of the Hyper Suprime-Cam (HSC) survey. The catalog covers an area of 136.9 deg$^2$ split into six fields, with a mean $i$-band seeing of $0.58$ arcsec and $5\sigma$ point-source depth of $i\sim 26$. Given conservative galaxy selection criteria for first year science, the depth and excellent image quality results in unweighted and weighted source number densities of 24.6 and 21.8 arcmin$^{-2}$, respectively. We define the requirements for cosmological weak lensing science with this catalog, then focus on characterizing potential systematics in the catalog using a series of internal null tests for problems with point-spread function (PSF) modeling, shear estimation, and other aspects of the image processing. We find that the PSF models narrowly meet requirements for weak lensing science with this catalog, with fractional PSF model size residuals of approximately $0.003$ (requirement: 0.004) and the PSF model shape correlation function $\rho_11^\circ$ scales that are sufficiently large as to require mitigation in cosmic shear measurements. Finally, we discuss the dominant systematics and the planned algorithmic changes to reduce them in future data reductions., Comment: 42 figures, 4 tables, v3 matches accepted version that will be published in PASJ (minor changes from v2). For high-resolution figures and cross-references with other HSC articles that will be in the same PASJ issue, please see the published version of the article

We present an optically-selected cluster catalog from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The HSC images are sufficiently deep to detect cluster member galaxies down to $M_*\sim 10^{10.2}M_\odot$ even at $z\sim 1$, allowing a reliable cluster detection at such high redshifts. We apply the CAMIRA algorithm to the HSC Wide S16A dataset covering $\sim 232$ deg$^2$ to construct a catalog of 1921 clusters at redshift $0.115$ that roughly corresponds to $M_{\rm 200m}\gtrsim 10^{14}h^{-1}M_\odot$. We confirm good cluster photometric redshift performance, with the bias and scatter in $��z/(1+z)$ being better than 0.005 and 0.01 over most of the redshift range, respectively. We compare our cluster catalog with large X-ray cluster catalogs from XXL and XMM-LSS surveys and find good correlation between richness and X-ray properties. We also study the miscentering effect from the distribution of offsets between optical and X-ray cluster centers. We confirm the high ($>0.9$) completeness and purity for high mass clusters by analyzing mock galaxy catalogs., 18 pages, 18 figures, accepted for publication in PASJ; cluster catalogs are available at http://www.slac.stanford.edu/~oguri/cluster/

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $��_{\rm int}(g-r), 18pages, 10 figures, accepted as PASJ special issue

We present the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $��^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$��$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $��_{��z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $��=0.5$%. At $3, Accepted for publication in ApJ

For more than two decades we have known that galaxy morphological segregation is present in the Local Universe. It is important to see how this relation evolves with cosmic time. To investigate how galaxy assembly took place with cosmic time, we explore the evolution of the morphology-density relation up to redshift z~1 using about 10000 galaxies drawn from the zCOSMOS Galaxy Redshift Survey. Taking advantage of accurate HST/ACS morphologies from the COSMOS survey, of the well-characterised zCOSMOS 3D environment, and of a large sample of galaxies with spectroscopic redshift, we want to study here the evolution of the morphology-density relation up to z~1 and its dependence on galaxy luminosity and stellar mass. The multi-wavelength coverage of the field also allows a first study of the galaxy morphological segregation dependence on colour. We further attempt to disentangle between processes that occurred early in the history of the Universe or late in the life of galaxies. The zCOSMOS field benefits of high-resolution imaging in the F814W filter from the Advanced Camera for Survey (ACS). We use standard morphology classifiers, optimised for being robust against band-shifting and surface brightness dimming, and a new, objective, and automated method to convert morphological parameters into early, spiral, and irregular types. We use about 10000 galaxies down to I_AB=22.5 with a spectroscopic sampling rate of 33% to characterise the environment of galaxies up to z~1 from the 100 kpc scales of galaxy groups up to the 100 Mpc scales of the cosmic web. ABRIDGED, Comment: 23 pages, 12 figures, accepted for publication in Astronomy and Astrophysics

[ABRIGED] We present the optical identification of a sample of 695 X-ray sources detected in the first 1.3 deg^2 of the XMM-COSMOS survey, down to a 0.5-2 keV (2-10 keV) limiting flux of ~10^-15 erg cm-2 s-1 (~5x10^-15 erg cm^-2 s-1). We were able to associate a candidate optical counterpart to ~90% (626) of the X-ray sources, while for the remaining ~10% of the sources we were not able to provide a unique optical association due to the faintness of the possible optical counterparts (I_AB>25) or to the presence of multiple optical sources within the XMM-Newton error circles. We also cross-correlated the candidate optical counterparts with the Subaru multicolor and ACS catalogs and with the Magellan/IMACS, zCOSMOS and literature spectroscopic data; the spectroscopic sample comprises 248 objects (~40% of the full sample). Our analysis reveals that for ~80% of the counterparts there is a very good agreement between the spectroscopic classification, the morphological parameters as derived from ACS data, and the optical to near infrared colors. About 20% of the sources show an apparent mismatch between the morphological and spectroscopic classifications. All the ``extended'' BL AGN lie at redshift, 17 pages, ApJS COSMOS Special Issue, 2007 in press. Full resolution version is available at http://www.mpe.mpg.de/XMMCosmos/PAPERS/

Deviations from general relativity, such as could be responsible for the cosmic acceleration, would influence the growth of large scale structure and the deflection of light by that structure. We clarify the relations between several different model independent approaches to deviations from general relativity appearing in the literature, devising a translation table. We examine current constraints on such deviations, using weak gravitational lensing data of the CFHTLS and COSMOS surveys, cosmic microwave background radiation data of WMAP5, and supernova distance data of Union2. Markov Chain Monte Carlo likelihood analysis of the parameters over various redshift ranges yields consistency with general relativity at the 95% confidence level., Comment: 11 pages; 7 figures; typographical errors corrected; this is the published version

We investigate how the shape of the galaxy two-point correlation function as measured in the zCOSMOS survey depends on local environment, quantified in terms of the density contrast on scales of 5 Mpc/h. We show that the flat shape previously observed at redshifts between z=0.6 and z=1 can be explained by this volume being simply 10% over-abundant in high-density environments, with respect to a Universal density probability distribution function. When galaxies corresponding to the top 10% tail of the distribution are excluded, the measured w_p(r_p) steepens and becomes indistinguishable from LCDM predictions on all scales. This is the same effect recognised by Abbas & Sheth in the SDSS data at z~0 and explained as a natural consequence of halo-environment correlations in a hierarchical scenario. Galaxies living in high-density regions trace dark matter halos with typically higher masses, which are more correlated. If the density probability distribution function of the sample is particularly rich in high-density regions because of the variance introduced by its finite size, this produces a distorted two-point correlation function. We argue that this is the dominant effect responsible for the observed "peculiar" clustering in the COSMOS field., 6 pages, 3 figures, accepted for publication in MNRAS

We investigate the (large-scale) bar fraction in a mass-complete sample of M > 10^10.5 Msun disk galaxies at 0.2 < z < 0.6 in the COSMOS field. The fraction of barred disks strongly depends on mass, disk morphology, and specific star formation rate (SSFR). At intermediate stellar mass (10^10.5 < M < 10^11 Msun) the bar fraction in early-type disks is much higher, at all redshifts, by a factor ~2, than that in late-type disks. This trend is reversed at higher stellar mass (M > 10^11 Msun), where the fraction of bars in early-type disks becomes significantly lower, at all redshifts, than that in late-type disks. The bar fractions for galaxies with low and high SSFRs closely follow those of the morphologically-selected early-type and late-type populations, respectively. This indicates a close correspondence between morphology and SSFR in disk galaxies at these earlier epochs. Interestingly, the total bar fraction in 10^10.5 < M < 10^11 Msun disks is built up by a factor of ~2 over the redshift interval explored, while for M > 10^11 Msun disks it remains roughly constant. This indicates that, already by z ~ 0.6, spectral and morphological transformations in the most massive disk galaxies have largely converged to the familiar Hubble sequence that we observe in the local Universe, while for intermediate mass disks this convergence is ongoing until at least z ~ 0.2. Moreover, these results highlight the importance of employing mass-limited samples for quantifying the evolution of barred galaxies. Finally, the evolution of the barred galaxy populations investigated does not depend on the large-scale environmental density (at least, on the scales which can be probed with the available photometric redshifts)., 10 pages, 4 figures, updated to reflect version accepted by MNRAS

We present three bright z+ dropout candidates selected from deep Near-Infrared (NIR) imaging of the COSMOS 2 square degree field. All three objects match the 0.8-8um colors of other published z>7 candidates but are three magnitudes brighter, facilitating further study. Deep spectroscopy of two of the candidates covering 0.64-1.02um with Keck-DEIMOS and all three covering 0.94-1.10um and 1.52-1.80um with Keck-NIRSPEC detects weak spectral features tentatively identified as Ly-alpha at z=6.95 and z=7.69 in two of the objects. The third object is placed at z~1.6 based on a 24um and weak optical detection. A comparison with the spectral energy distributions of known z1um properties of all three objects can be matched to optically detected sources with photometric redshifts at z~1.8, so the non-detection in the i+ and z+ bands are the primary factors which favors a z>7 solution. If any of these objects are at z~7 the bright end of the luminosity function is significantly higher at z>7 than suggested by previous studies, but consistent within the statistical uncertainty and the dark matter halo distribution. If these objects are at low redshift, the Lyman-Break selection must be contaminated by a previously unknown population of low redshift objects with very strong breaks in their broad band spectral energy distributions and blue NIR colors. The implications of this result on luminosity function evolution at high redshift is discussed. We show that the primary limitation of z>7 galaxy searches with broad filters is the depth of the available optical data., 15 Pages, 15 figures, accepted to ApJ

To derive a new H$��$ luminosity function and to understand the clustering properties of star-forming galaxies at $z \approx 0.24$, we have made a narrow-band imaging survey for H$��$ emitting galaxies in the HST COSMOS 2 square degree field. We used the narrow-band filter NB816 ($��_c = 8150$ ��, $����= 120$ ��) and sampled H$��$ emitters with $EW_{\rm obs}(\rm H��+ [N\textsc{ii}]) > 12$ ��in a redshift range between $z=0.233$ and $z=0.251$ corresponding to a depth of 70 Mpc. We obtained 980 H$��$ emitting galaxies in a sky area of 5540 arcmin$^2$, corresponding to a survey volume of $3.1 \times 10^4 {\rm Mpc^3}$. We derive a H$��$ luminosity function with a best-fit Schechter function parameter set of $��= -1.35^{+0.11}_{-0.13}$, $\log��_* = -2.65^{+0.27}_{-0.38}$, and $\log L_* ({\rm erg s^{-1}}) = 41.94^{+0.38}_{-0.23}$. The H$��$ luminosity density is $2.7^{+0.7}_{-0.6} \times 10^{39}$ ergs s$^{-1}$ Mpc$^{-3}$. After subtracting the AGN contribution (15 %) to the H$��$ luminosity density, the star formation rate density is evaluated as $1.8^{+0.7}_{-0.4} \times 10^{-2}$ $M_{\sun}$ yr$^{-1}$ Mpc$^{-3}$. The angular two-point correlation function of H$��$ emitting galaxies of $\log L({\rm H��}) > 39.8$ is well fit by a power law form of $w(��) = 0.013^{+0.002}_{-0.001} ��^{-0.88 \pm 0.03}$, corresponding to the correlation function of $��(r) = (r/1.9{\rm Mpc})^{-1.88}$. We also find that the H$��$ emitters with higher H$��$ luminosity are more strongly clustered than those with lower luminosity., 13 pages, 9 fifgure, to apper in ApJS