Search

We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $\Delta\!\Sigma(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.150.75$ for the $\Delta\!\Sigma$ measurements, selected based on their photometric redshifts. For theoretical template, we use the "minimal bias" model for the cosmological clustering observables for the flat $\Lambda$CDM cosmological model. We compare the model predictions with the measurements in each redshift bin on large scales, $R>12$ and $8~h^{-1}\mathrm{Mpc}$ for $\Delta\!\Sigma(R)$ and $w_\mathrm{p}(R)$, respectively, where the perturbation theory-inspired model is valid. When we employ weak priors on cosmological parameters, without CMB information, we find $S_8=0.936^{+0.092}_{-0.086}$, $\sigma_8=0.85^{+0.16}_{-0.11}$, and $\Omega_\mathrm{m}=0.283^{+0.12}_{-0.035}$ for the flat $\Lambda$CDM model. Although the central value of $S_8$ appears to be larger than those inferred from other cosmological experiments, we find that the difference is consistent with expected differences due to sample variance, and our results are consistent with the other results to within the statistical uncertainties. (abriged), Comment: 24 pages, 19 figures, 4 tables, to be submitted to Phys. Rev. D

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

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

Galaxy clusters are a promising probe of late-time structure growth, but constraints on cosmology from cluster abundances are currently limited by systematics in their inferred masses. One unmitigated systematic effect in weak-lensing mass inference is ignoring the presence of baryons and treating the entire cluster as a dark matter halo. In this work we present a new flexible model for cluster densities that captures both the baryonic and dark matter profiles, a new general technique for calculating the lensing signal of an arbitrary density profile, and a methodology for stacking those lensing signal to appropriately model stacked weak-lensing measurements of galaxy cluster catalogues. We test this model on 1400 simulated clusters. Similarly to previous studies, we find that a dark matter-only model overestimates the average mass by $7.5\%$, but including our baryonic term reduces that to $0.7\%$. Additionally, to mitigate the computational complexity of our model, we construct an emulator (surrogate model) which accurately interpolates our model for parameter inference, while being much faster to use than the raw model. We also provide an open-source software framework for our model and emulator, called maszcal, which will serve as a platform for continued efforts to improve these mass-calibration techniques. In this work, we detail our model, the construction of the emulator, and the tests which we used to validate that our model does mitigate bias. Lastly, we describe tests of the emulator's accuracy, Comment: 27 pages, 16 figures. To be submitted to JCAP

We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST., Comment: comments welcome

We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS))., Comment: added references, comments welcome

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., Comment: 27 pages, 17 figures, submitted to APJS

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., Comment: Submitted to MNRAS; 21 pages, 12 figures

With current and upcoming experiments such as WFIRST, Euclid and LSST, we can observe up to billions of galaxies. While such surveys cannot obtain spectra for all observed galaxies, they produce galaxy magnitudes in color filters. This data set behaves like a high-dimensional nonlinear surface, an excellent target for machine learning. In this work, we use a lightcone of semianalytic galaxies tuned to match CANDELS observations from Lu et al. (2014) to train a set of neural networks on a set of galaxy physical properties. We add realistic photometric noise and use trained neural networks to predict stellar masses and average star formation rates on real CANDELS galaxies, comparing our predictions to SED fitting results. On semianalytic galaxies, we are nearly competitive with template-fitting methods, with biases of $0.01$ dex for stellar mass, $0.09$ dex for star formation rate, and $0.04$ dex for metallicity. For the observed CANDELS data, our results are consistent with template fits on the same data at $0.15$ dex bias in $M_{\rm star}$ and $0.61$ dex bias in star formation rate. Some of the bias is driven by SED-fitting limitations, rather than limitations on the training set, and some is intrinsic to the neural network method. Further errors are likely caused by differences in noise properties between the semianalytic catalogs and data. Our results show that galaxy physical properties can in principle be measured with neural networks at a competitive degree of accuracy and precision to template-fitting methods., Comment: 19 pages, 10 figures, 6 tables. Accepted for publication in 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<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints., Comment: 43 pages, 21 figures, accepted for publication in PASJ

In the past few years, several independent collaborations have presented cosmological constraints from tomographic cosmic shear analyses. These analyses differ in many aspects: the datasets, the shear and photometric redshift estimation algorithms, the theory model assumptions, and the inference pipelines. To assess the robustness of the existing cosmic shear results, we present in this paper a unified analysis of four of the recent cosmic shear surveys: the Deep Lens Survey (DLS), the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Science Verification data from the Dark Energy Survey (DES-SV), and the 450 deg$^{2}$ release of the Kilo-Degree Survey (KiDS-450). By using a unified pipeline, we show how the cosmological constraints are sensitive to the various details of the pipeline. We identify several analysis choices that can shift the cosmological constraints by a significant fraction of the uncertainties. For our fiducial analysis choice, considering a Gaussian covariance, conservative scale cuts, assuming no baryonic feedback contamination, identical cosmological parameter priors and intrinsic alignment treatments, we find the constraints (mean, 16% and 84% confidence intervals) on the parameter $S_{8}\equiv \sigma_{8}(\Omega_{\rm m}/0.3)^{0.5}$ to be $S_{8}=0.94_{-0.045}^{+0.046}$ (DLS), $0.66_{-0.071}^{+0.070}$ (CFHTLenS), $0.84_{-0.061}^{+0.062}$ (DES-SV) and $0.76_{-0.049}^{+0.048}$ (KiDS-450). From the goodness-of-fit and the Bayesian evidence ratio, we determine that amongst the four surveys, the two more recent surveys, DES-SV and KiDS-450, have acceptable goodness-of-fit and are consistent with each other. The combined constraints are $S_{8}=0.79^{+0.042}_{-0.041}$, which is in good agreement with the first year of DES cosmic shear results and recent CMB constraints from the Planck satellite., Comment: 22 pages, 15 figures, 7 tables; submitted to MNRAS

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

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 optimized source galaxy selection schemes for measuring cluster weak lensing (WL) mass profiles unaffected by cluster member dilution from the Subaru Hyper Suprime-Cam Strategic Survey Program (HSC-SSP). The ongoing HSC-SSP survey will uncover thousands of galaxy clusters to $z\lesssim1.5$. In deriving cluster masses via WL, a critical source of systematics is contamination and dilution of the lensing signal by cluster {members, and by foreground galaxies whose photometric redshifts are biased}. Using the first-year CAMIRA catalog of $\sim$900 clusters with richness larger than 20 found in $\sim$140 deg$^2$ of HSC-SSP data, we devise and compare several source selection methods, including selection in color-color space (CC-cut), and selection of robust photometric redshifts by applying constraints on their cumulative probability distribution function (PDF; P-cut). We examine the dependence of the contamination on the chosen limits adopted for each method. Using the proper limits, these methods give mass profiles with minimal dilution in agreement with one another. We find that not adopting either the CC-cut or P-cut methods results in an underestimation of the total cluster mass ($13\pm4\%$) and the concentration of the profile ($24\pm11\%$). The level of cluster contamination can reach as high as $\sim10\%$ at $R\approx 0.24$ Mpc/$h$ for low-z clusters without cuts, while employing either the P-cut or CC-cut results in cluster contamination consistent with zero to within the 0.5% uncertainties. Our robust methods yield a $\sim60\sigma$ detection of the stacked CAMIRA surface mass density profile, with a mean mass of $M_\mathrm{200c} = (1.67\pm0.05({\rm {stat}}))\times 10^{14}\,M_\odot/h$., Comment: 19 pages, 4 tables, 12 figures, accepted to PASJ special issue

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_1<3\times 10^{-7}$ (requirement: $4\times 10^{-7}$) at 0.5$^\circ$ scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star-galaxy shape correlations reveal additive systematics on $>1^\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

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/., Comment: 34 pages, 20 figures, 7 tables, moderate revision, accepted for publication in PASJ

Recently, several studies have discovered a strong discrepancy between the large-scale clustering biases of two subsamples of galaxy clusters at the same halo mass, split by their average projected membership distances $R_{\mathrm{mem}}$. The level of this discrepancy significantly exceeds the maximum halo assembly bias signal predicted by LCDM. In this study, we explore whether some of the clustering bias differences could be caused by biases in $R_{\mathrm{mem}}$ due to projection effects from other systems along the line-of-sight. We thoroughly investigate the halo assembly bias of the photometrically-detected redMaPPer clusters in SDSS, by defining a new variant of the average membership distance estimator $\tilde{R}_{\mathrm{mem}}$ that is more robust against projection effects in the cluster membership identification. Using the angular mark correlation functions of clusters, we show that the large-scale bias differences when splitting by $R_{\mathrm{mem}}$ can be largely attributed to such projection effects. After splitting by $\tilde{R}_{\mathrm{mem}}$, the anomalously large signal is reduced, giving a ratio of $1.02\pm0.14$ between the two clustering biases as measured from weak lensing. Using a realistic mock cluster catalog, we predict that the bias ratio between two $\tilde{R}_{\mathrm{mem}}$-split subsamples should be $<1.10$, which is at least 60% weaker than the maximum halo assembly bias signal (1.24) when split by halo concentration. Therefore, our results demonstrate that the level of halo assembly bias exhibited by redMaPPer clusters in SDSS is consistent with the LCDM prediction. With a ten-fold increase in cluster numbers, deeper ongoing surveys will enable a more robust detection of halo assembly bias. Our findings also have important implications for how projection effects and their impact on cluster cosmology can be quantified in photometric cluster catalogs., Comment: 11 pages, 8 figures, MNRAS published

We perform a measurement of the mass--richness relation of the redMaPPer galaxy cluster catalogue using weak lensing data from the Sloan Digital Sky Survey. We have carefully characterized a broad range of systematic uncertainties, including shear calibration errors, photo-$z$ biases, dilution by member galaxies, source obscuration, magnification bias, incorrect assumptions about cluster mass profiles, cluster centering, halo triaxiality, and projection effects. We also compare measurements of the lensing signal from two independently-produced shear and photometric redshift catalogues to characterize systematic errors in the lensing signal itself. Using a sample of 5,570 clusters from $0.1\le z\le 0.33$, the normalization of our power-law mass vs.\ $\lambda$ relation is $\log_{10}[M_{200m}/h^{-1}\ M_{\odot}]$ = $14.344 \pm 0.021$ (statistical) $\pm 0.023$ (systematic) at a richness $\lambda=40$, a 7 per cent calibration uncertainty, with a power-law index of $1.33^{+0.09}_{-0.10}$ ($1\sigma$). The detailed systematics characterization in this work renders it the definitive weak lensing mass calibration for SDSS redMaPPer clusters at this time., Comment: 18 pages, 4 figures. Cosmological dependence updated from previous version

The use of large, X-ray selected galaxy cluster catalogues for cosmological analyses requires a thorough understanding of the X-ray mass estimates. Weak gravitational lensing is an ideal method to shed light on such issues, due to its insensitivity to the cluster dynamical state. We perform a weak lensing calibration of 166 galaxy clusters from the REFLEX and BCS cluster catalogue and compare our results to the X-ray masses based on scaled luminosities from that catalogue. To interpret the weak lensing signal in terms of cluster masses, we compare the lensing signal to simple theoretical Navarro-Frenk-White models and to simulated cluster lensing profiles, including complications such as cluster substructure, projected large-scale structure, and Eddington bias. We find evidence of underestimation in the X-ray masses, as expected, with $\langle M_{\mathrm{X}}/M_{\mathrm{WL}}\rangle = 0.75 \pm 0.07$ stat. $\pm 0.05$ sys. for our best-fit model. The biases in cosmological parameters in a typical cluster abundance measurement that ignores this mass bias will typically exceed the statistical errors., Comment: 13 pages, 5 figures. Revised to address referee comments

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\'{e}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., Comment: 32 pages + 15 pages of technical appendices; 28 figures; published in MNRAS; all data and other information may be found at https://github.com/barnabytprowe/great3-public as the GREAT3 website and leaderboard are no longer live. v3 has updates to this comment only

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 ., Comment: 38 pages, 3 tables, 8 figures, accepted for publication in Astronomy & Computing

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

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., Comment: 30 pages, 13 figures, published in MNRAS. URLs in abstract no longer work, but all data are available for download as described in https://github.com/barnabytprowe/great3-public (v4 only has updates to comments)

We present details of the construction and characterization of the coaddition of the Sloan Digital Sky Survey Stripe 82 \ugriz\ imaging data. This survey consists of 275 deg$^2$ of repeated scanning by the SDSS camera of $2.5\arcdeg$ of $\delta$ over $-50\arcdeg \le \alpha \le 60\arcdeg$ centered on the Celestial Equator. Each piece of sky has $\sim 20$ runs contributing and thus reaches $\sim2$ magnitudes fainter than the SDSS single pass data, i.e. to $r\sim 23.5$ for galaxies. We discuss the image processing of the coaddition, the modeling of the PSF, the calibration, and the production of standard SDSS catalogs. The data have $r$-band median seeing of 1.1\arcsec, and are calibrated to $\le 1%$. Star color-color, number counts, and psf size vs modelled size plots show the modelling of the PSF is good enough for precision 5-band photometry. Structure in the psf-model vs magnitude plot show minor psf mis-modelling that leads to a region where stars are being mis-classified as galaxies, and this is verified using VVDS spectroscopy. As this is a wide area deep survey there are a variety of uses for the data, including galactic structure, photometric redshift computation, cluster finding and cross wavelength measurements, weak lensing cluster mass calibrations, and cosmic shear measurements., Comment: 18 page, 10 figures, submitted to ApJ. Small changes in text to be consistent with revised photo-z catalog of Reis et al. arXiv:1111.6620v2

Stripe 82 in the Sloan Digital Sky Survey was observed multiple times, allowing deeper images to be constructed by coadding the data. Here we analyze the ellipticities of background galaxies in this 275 square degree region, searching for evidence of distortions due to cosmic shear. The E-mode is detected in both real and Fourier space with $>5$-$\sigma$ significance on degree scales, while the B-mode is consistent with zero as expected. The amplitude of the signal constrains the combination of the matter density $\Omega_m$ and fluctuation amplitude $\sigma_8$ to be $\Omega_m^{0.7}\sigma_8 = 0.252^{+0.032}_{-0.052}$., Comment: 17 pages, 19 figures, submitted to ApJ. Analysis updated using revised photo-z catalog of Reis et al. arXiv:1111.6620v2. Changes in results are within the errors and the conclusions are unaffected

The shapes of distant galaxies are sheared by intervening galaxy clusters. We examine this effect in Stripe 82, a 275 square degree region observed multiple times in the Sloan Digital Sky Survey and coadded to achieve greater depth. We obtain a mass-richness calibration that is similar to other SDSS analyses, demonstrating that the coaddition process did not adversely affect the lensing signal. We also propose a new parameterization of the effect of tomography on the cluster lensing signal which does not require binning in redshift, and we show that using this parameterization we can detect tomography for stacked clusters at varying redshifts. Finally, due to the sensitivity of the tomographic detection to accurately marginalizing over the effect of the cluster mass, we show that tomography at low redshift (where dependence on exact cosmological models is weak) can be used to constrain mass profiles in clusters., Comment: 8 pages, 13 figures, submitted to ApJ. Analysis updated using revised photo-z catalog of Reis et al. arXiv:1111.6620v2. Changes in results are within the errors and the conclusions are unaffected

We present and describe a catalog of galaxy photometric redshifts (photo-z's) for the Sloan Digital Sky Survey (SDSS) Coadd Data. We use the Artificial Neural Network (ANN) technique to calculate photo-z's and the Nearest Neighbor Error (NNE) method to estimate photo-z errors for $\sim$ 13 million objects classified as galaxies in the coadd with $r < 24.5$. The photo-z and photo-z error estimators are trained and validated on a sample of $\sim 83,000$ galaxies that have SDSS photometry and spectroscopic redshifts measured by the SDSS Data Release 7 (DR7), the Canadian Network for Observational Cosmology Field Galaxy Survey (CNOC2), the Deep Extragalactic Evolutionary Probe Data Release 3(DEEP2 DR3), the VIsible imaging Multi-Object Spectrograph - Very Large Telescope Deep Survey (VVDS) and the WiggleZ Dark Energy Survey. For the best ANN methods we have tried, we find that 68% of the galaxies in the validation set have a photo-z error smaller than $\sigma_{68} =0.031$. After presenting our results and quality tests, we provide a short guide for users accessing the public data., Comment: 16 pages, 13 figures, submitted to ApJ. Analysis updated to remove proprietary BOSS data comprising small fraction (8%) of original spectroscopic training set and erroneously included. Changes in results are small compared to the errors and the conclusions are unaffected. arXiv admin note: substantial text overlap with arXiv:0708.0030

The Alpha Magnetic Spectrometer (AMS-02), which is scheduled to be deployed onboard the International Space Station later this year, will be capable of measuring the composition and spectra of GeV-TeV cosmic rays with unprecedented precision. In this paper, we study how the projected measurements from AMS-02 of stable secondary-to-primary and unstable ratios (such as boron-to-carbon and beryllium-10-to-beryllium-9) can constrain the models used to describe the propagation of cosmic rays throughout the Milky Way. We find that within the context of fairly simple propagation models, all of the model parameters can be determined with high precision from the projected AMS-02 data. Such measurements are less constraining in more complex scenarios, however, which allow for departures from a power-law form for the diffusion coefficient, for example, or for inhomogeneity or stochasticity in the distribution and chemical abundances of cosmic ray sources., Comment: 12 pages, 7 figures, 3 tables, matches published version

In recent years, a number of experiments have been conducted with the goal of studying cosmic rays at GeV to TeV energies. This is a particularly interesting regime from the perspective of indirect dark matter detection. To draw reliable conclusions regarding dark matter from cosmic ray measurements, however, it is important to first understand the propagation of cosmic rays through the magnetic and radiation fields of the Milky Way. In this paper, we constrain the characteristics of the cosmic ray propagation model through comparison with observational inputs, including recent data from the CREAM experiment, and use these constraints to estimate the corresponding uncertainties in the spectrum of cosmic ray electrons and positrons from dark matter particles annihilating in the halo of the Milky Way., Comment: 21 pages, 9 figures

Recent preliminary results from the PAMELA experiment indicate the presence of an excess of cosmic ray positrons above 10 GeV. In this letter, we consider possibility that this signal is the result of dark matter annihilations taking place in the halo of the Milky Way. Rather than focusing on a specific particle physics model, we take a phenomenological approach and consider a variety of masses and two-body annihilation modes, including W+W-, ZZ, b bbar, tau+ tau-, mu+ mu-, and e+e. We also consider a range of diffusion parameters consistent with current cosmic ray data. We find that a significant upturn in the positron fraction above 10 GeV is compatible with a wide range of dark matter annihilation modes, although very large annihilation cross sections and/or boost factors arising from inhomogeneities in the local dark matter distribution are required to produce the observed intensity of the signal. We comment on constraints from gamma rays, synchrotron emission, and cosmic ray antiproton measurements., Comment: 4 pages, 1 figure

Very high energy gamma-rays are expected to be absorbed by the extragalactic background light over cosmological distances via the process of electron-positron pair production. Recent observations of cosmologically distant gamma-ray emitters by ground based gamma-ray telescopes have, however, revealed a surprising degree of transparency of the universe to very high energy photons. One possible mechanism to explain this observation is the oscillation between photons and axion-like-particles (ALPs). Here we explore this possibility further, focusing on photon-ALP conversion in the magnetic fields in and around gamma-ray sources and in the magnetic field of the Milky Way, where some fraction of the ALP flux is converted back into photons. We show that this mechanism can be efficient in allowed regions of the ALP parameter space, as well as in typical configurations of the Galactic Magnetic Field. As case examples, we consider the spectrum observed from two HESS sources: 1ES1101-232 at redshift z=0.186 and H 2356-309 at z=0.165. We also discuss features of this scenario which could be used to distinguish it from standard or other exotic models., Comment: 7 pages, 4 figures. Matches published version

In 240 days of X-ray monitoring of M82, we have discovered an X-ray periodicity at $62.0 \pm 2.5$ days with a peak to peak amplitude corresponding to an isotropic luminosity of $2.4 \times 10^{40} \rm erg s^{-1}$ in M82 and an X-ray flare reaching a peak luminosity of $9.8 \times 10^{40} \rm erg s^{-1}$. The periodicity and flare likely originate from the ultraluminous X-ray source (ULX) in M82 which has been identified as a possible intermediate mass black hole. We suggest that the 62~day modulation is due to orbital motion within an X-ray binary with a Roche-lobe overflowing companion star which would imply that the average density of the companion star is near $5 \times 10^{-5} \rm g cm^{-3}$ and is therefore a giant or supergiant. Chandra observations just after the flare show an energy spectrum that is consistent with a power-law with no evidence of a thermal component or line emission. Radio observations made with the VLA during the flare allow us to rule out a blazar identification for the source and place strong constraints on relativistically beamed models of the X-ray emission. The Chandra observations reveal a second X-ray source reached a flux of $4.4 \times 10^{-12} \rm erg cm^{-2} s^{-1}$ in the 0.3-7 keV band which is dramatically higher than any flux previously seen from this source and corresponds to an isotropic luminosity of $1.1 \times 10^{40} \rm erg s^{-1}$. This source is a second ultraluminous X-ray source in M82 and may give rise to the QPOs detected from the central region of M82., Comment: 15 pages, one color figure, to appear in ApJ

The ultraluminous x-ray source (ULX) in the galaxy M82 has been identified as a possible intermediate-mass black hole. We have found that the x-ray flux from M82 is modulated with a peak-to-peak amplitude corresponding to an isotropic luminosity of 2.4x10^40 erg/s in M82 and a period of 62.0 +/- 2.5 days, which we interpret as the orbital period of the ULX binary. This orbital period implies that the mass-donor star must be a giant or supergiant. Large mass-transfer rates, sufficient to fuel the ULX, are expected for a giant-phase mass donor in an x-ray binary. The giant phase has a short lifetime, indicating that we see the ULX in M82 in a brief and unusual period of its evolution., Comment: 3 pages, appeared in Science

UNKNOWN

UNKNOWN