Showing 1–50 of 244 results for author: Lupton, R

The little coadd that could: Estimating shear from coadded images

Authors: Robert Armstrong, Erin Sheldon, Eric Huff, Jim Bosch, Eli Rykoff, Rachel Mandelbaum, Arun Kannawadi, Peter Melchior, Robert Lupton, Matthew R. Becker, Yusra Al-Sayyed, The LSST Dark Energy Science Collaboration

Abstract: Upcoming wide field surveys will have many overlapping epochs of the same region of sky. The conventional wisdom is that in order to reduce the errors sufficiently for systematics-limited measurements, like weak lensing, we must do simultaneous fitting of all the epochs. Using current algorithms this will require a significant amount of computing time and effort. In this paper, we revisit the pote… ▽ More Upcoming wide field surveys will have many overlapping epochs of the same region of sky. The conventional wisdom is that in order to reduce the errors sufficiently for systematics-limited measurements, like weak lensing, we must do simultaneous fitting of all the epochs. Using current algorithms this will require a significant amount of computing time and effort. In this paper, we revisit the potential of using coadds for shear measurements. We show on a set of image simulations that the multiplicative shear bias can be constrained below the 0.1% level on coadds, which is sufficient for future lensing surveys. We see no significant differences between simultaneous fitting and coadded approaches for two independent shear codes: Metacalibration and BFD. One caveat of our approach is the assumption of a principled coadd, i.e. the PSF is mathematically well-defined for all the input images. This requires us to reject CCD images that do not fully cover the coadd region. We estimate that the number of epochs that must be rejected for a survey like LSST is on the order of 20%, resulting in a small loss in depth of less than 0.1 magnitudes. We also put forward a cell-based coaddition scheme that meets the above requirements for unbiased weak lensing shear estimation in the context of LSST. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Instrument Signature Removal and Calibration Products for the Rubin Legacy Survey of Space and Time

Authors: Andrés A. Plazas Malagón, Chris Waters, Alex Broughton, Eli Rykoff, Agnès Ferté, Merlin Fisher-Levine, Robert Lupton

Abstract: The Vera C. Rubin Legacy Survey of Space and Time (LSST) will conduct an unprecedented optical survey of the southern sky, imaging the entire available sky every few nights for 10 years. To achieve its ambitious science goals of probing dark energy and dark matter, mapping the Milky Way, and exploring the transient optical sky, the systematic errors in the LSST data must be exquisitely controlled.… ▽ More The Vera C. Rubin Legacy Survey of Space and Time (LSST) will conduct an unprecedented optical survey of the southern sky, imaging the entire available sky every few nights for 10 years. To achieve its ambitious science goals of probing dark energy and dark matter, mapping the Milky Way, and exploring the transient optical sky, the systematic errors in the LSST data must be exquisitely controlled. Instrument signature removal (ISR) is a critical early step in LSST data processing to remove inherent camera effects from the raw images and produce accurate representations of the incoming light. This paper describes the current state of the ISR pipelines implemented in the LSST Science Pipelines software. The key steps in ISR are outlined, and the process of generating and verifying the necessary calibration products to carry out ISR is also discussed. Finally, an overview is given of how the Rubin data management system utilizes a data Butler and calibration collections to organize datasets and match images to appropriate calibrations during processing. Precise ISR will be essential to realize the potential of LSST to revolutionize astrophysics. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: Submitted for publication to the Journal of Astronomical Telescopes, Instruments, and Systems as part of the Special Section "Image Sensors for Precision Astronomy" (ISPA 2024 Conference at SLAC/KIPAC, March 12-14, 2024, https://www.bnl.gov/ispaworkshop/)

The structure of the stellar halo of the Andromeda galaxy explored with the NB515 for Subaru/HSC. I.: New Insights on the stellar halo up to 120 kpc

Authors: Itsuki Ogami, Mikito Tanaka, Yutaka Komiyama, Masashi Chiba, Puragra Guhathakurta, Evan N. Kirby, Rosemary F. G. Wyse, Carrie Filion, Karoline M. Gilbert, Ivanna Escala, Masao Mori, Takanobu Kirihara, Masayuki Tanaka, Miho N. Ishigaki, Kohei Hayashi, Myun Gyoon Lee, Sanjib Sharma, Jason S. Kalirai, Robert H. Lupton

Abstract: We analyse the M31 halo and its substructure within a projected radius of 120 kpc using a combination of Subaru/HSC NB515 and CFHT/MegaCam g- & i-bands. We succeed in separating M31's halo stars from foreground contamination with $\sim$ 90 \% accuracy by using the surface gravity sensitive NB515 filter. Based on the selected M31 halo stars, we discover three new substructures, which associate with… ▽ More We analyse the M31 halo and its substructure within a projected radius of 120 kpc using a combination of Subaru/HSC NB515 and CFHT/MegaCam g- & i-bands. We succeed in separating M31's halo stars from foreground contamination with $\sim$ 90 \% accuracy by using the surface gravity sensitive NB515 filter. Based on the selected M31 halo stars, we discover three new substructures, which associate with the Giant Southern Stream (GSS) based on their photometric metallicity estimates. We also produce the distance and photometric metallicity estimates for the known substructures. While these quantities for the GSS are reproduced in our study, we find that the North-Western stream shows a steeper distance gradient than found in an earlier study, suggesting that it is likely to have formed in an orbit closer to the Milky Way. For two streams in the eastern halo (Stream C and D), we identify distance gradients that had not been resolved. Finally, we investigate the global halo photometric metallicity distribution and surface brightness profile using the NB515-selected halo stars. We find that the surface brightness of the metal-poor and metal-rich halo populations, and the all population can be fitted to a power-law profile with an index of $α= -1.65 \pm 0.02$, $-2.82\pm0.01$, and $-2.44\pm0.01$, respectively. In contrast to the relative smoothness of the halo profile, its photometric metallicity distribution appears to be spatially non-uniform with nonmonotonic trends with radius, suggesting that the halo population had insufficient time to dynamically homogenize the accreted populations. △ Less

Submitted 1 January, 2024; originally announced January 2024.

Comments: 24 pages, 26 figures, 5 tables, submitted to MNRAS

Final Results of Search for New Milky Way Satellites in the Hyper Suprime-Cam Subaru Strategic Program Survey: Discovery of Two More Candidates

Authors: Daisuke Homma, Masashi Chiba, Yutaka Komiyama, Masayuki Tanaka, Sakurako Okamoto, Mikito Tanaka, Miho N. Ishigaki, Kohei Hayashi, Nobuo Arimoto, Robert H. Lupton, Michael A. Strauss, Satoshi Miyazaki, Shiang-Yu Wang, Hitoshi Murayama

Abstract: We present the final results of our search for new Milky Way (MW) satellites using the data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) survey over $\sim 1,140$ deg$^2$. In addition to three candidates that we already reported, we have identified two new MW satellite candidates in the constellation of Sextans at a heliocentric distance of $D_{\odot} \simeq 126$kpc, and Virgo at… ▽ More We present the final results of our search for new Milky Way (MW) satellites using the data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) survey over $\sim 1,140$ deg$^2$. In addition to three candidates that we already reported, we have identified two new MW satellite candidates in the constellation of Sextans at a heliocentric distance of $D_{\odot} \simeq 126$kpc, and Virgo at $D_{\odot} \simeq 151$kpc, named Sextans II and Virgo III, respectively. Their luminosities (Sext II:$M_V\simeq-3.9$mag; Vir III:$M_V\simeq-2.7$mag) and half-light radii (Sext II:$r_h\simeq154$ pc; Vir III:$r_h\simeq 44$ pc) place them in the region of size-luminosity space of ultra-faint dwarf galaxies (UFDs). Including four previously known satellites, there are a total of nine satellites in the HSC-SSP footprint. This discovery rate of UFDs is much higher than that predicted from the recent models for the expected population of MW satellites in the framework of cold dark matter models, thereby suggesting that we encounter a too many satellites problem. Possible solutions to settle this tension are also discussed. △ Less

Submitted 30 April, 2024; v1 submitted 9 November, 2023; originally announced November 2023.

Comments: 20 pages, 17 figures, 4 tables, accepted for publication in PASJ

The Merian Survey: Design, Construction, and Characterization of a Filter Set Optimized to Find Dwarf Galaxies and Measure their Dark Matter Halo Properties with Weak Lensing

Authors: Yifei Luo, Alexie Leauthaud, Jenny Greene, Song Huang, Erin Kado-Fong, Shany Danieli, Ting S. Li, Jiaxuan Li, Diana Blanco, Erik J. Wasleske, Joseph Wick, Abby Mintz, Runquan Guan, Annika H. G. Peter, Vivienne Baldassare, Alyson Brooks, Arka Banerjee, Joy Bhattacharyya, Zheng Cai, Xinjun Chen, Jim Gunn, Sean D. Johnson, Lee S. Kelvin, Mingyu Li, Xiaojing Lin , et al. (6 additional authors not shown)

Abstract: 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 desig… ▽ More 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 ($λ_c = 7080 \unicode{x212B}$, $Δλ= 275\unicode{x212B}$) and N540 ($λ_c = 5400\unicode{x212B}$, $Δλ= 210\unicode{x212B}$). The central wavelengths and filter widths of N708 and N540 were designed to detect the $\rm Hα$ and $\rm [OIII]$ emission lines of galaxies in the mass range $8<\rm \log M_*/M_\odot<9$ by comparing Merian fluxes with HSC broad-band fluxes. Our filter design takes into account the weak lensing S/N and photometric redshift performance. Our simulations predict that Merian will yield a sample of $\sim$ 85,000 star-forming dwarf galaxies with a photometric redshift accuracy of $σ_{Δz/(1+z)}\sim 0.01$ and an outlier fraction of $η=2.8\%$ over the redshift range $0.058<z<0.10$. With 60 full nights on the Blanco/Dark Energy Camera (DECam), the Merian survey is predicted to measure the average weak lensing profile around dwarf galaxies with lensing $\rm S/N \sim 32$ within $r<0.5$ Mpc and lensing $\rm S/N \sim 90$ within $r<1.0$ Mpc. This unprecedented sample of star-forming dwarf galaxies will allow for studies of the interplay between dark matter and stellar feedback and their roles in the evolution of dwarf galaxies. △ Less

Submitted 3 April, 2024; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: 18 pages, 14 figures, accepted for publication in MNRAS

Quasar Luminosity Function at z = 7

Authors: Yoshiki Matsuoka, Masafusa Onoue, Kazushi Iwasawa, Michael A. Strauss, Nobunari Kashikawa, Takuma Izumi, Tohru Nagao, Masatoshi Imanishi, Masayuki Akiyama, John D. Silverman, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Kohei Inayoshi, Rikako Ishimoto, Toshihiro Kawaguchi, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Chien-Hsiu Lee, Robert H. Lupton , et al. (19 additional authors not shown)

Abstract: We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We… ▽ More We present the quasar luminosity function (LF) at $z = 7$, measured with 35 spectroscopically confirmed quasars at $6.55 < z < 7.15$. The sample of 22 quasars from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, combined with 13 brighter quasars in the literature, covers an unprecedentedly wide range of rest-frame ultraviolet magnitudes over $-28 < M_{1450} < -23$. We found that the binned LF flattens significantly toward the faint end populated by the SHELLQs quasars. A maximum likelihood fit to a double power-law model has a break magnitude $M^*_{1450} = -25.60^{+0.40}_{-0.30}$, a characteristic density $Φ^* = 1.35^{+0.47}_{-0.30}$ Gpc$^{-3}$ mag$^{-1}$, and a bright-end slope $β= -3.34^{+0.49}_{-0.57}$, when the faint-end slope is fixed to $α= -1.2$ as observed at $z \le 6$. The overall LF shape remains remarkably similar from $z = 4$ to $7$, while the amplitude decreases substantially toward higher redshifts, with a clear indication of an accelerating decline at $z \ge 6$. The estimated ionizing photon density, $10^{48.2 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$, is less than 1 % of the critical rate to keep the intergalactic medium ionized at $z = 7$, and thus indicates that quasars are not a major contributor to cosmic reionization. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: The Astrophysical Journal Letters, in press

Hyper Suprime-Cam Year 3 Results: Cosmology from Galaxy Clustering and Weak Lensing with HSC and SDSS using the Minimal Bias Model

Authors: Sunao Sugiyama, Hironao Miyatake, Surhud More, Xiangchong Li, Masato Shirasaki, Masahiro Takada, Yosuke Kobayashi, Ryuichi Takahashi, Takahiro Nishimichi, Atsushi J. Nishizawa, Markus M. Rau, Tianqing Zhang, Roohi Dalal, Rachel Mandelbaum, Michael A. Strauss, Takashi Hamana, Masamune Oguri, Ken Osato, Arun Kannawadi, Robert Armstrong, Yutaka Komiyama, Robert H. Lupton, Nate B. Lust, Satoshi Miyazaki, Hitoshi Murayama , et al. (5 additional authors not shown)

Abstract: We present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering 416 deg$^2$, and the SDSS DR11 spectroscopic galaxies spanning the redshift range $[0.15, 0.70]$. We subdivide the SDSS galaxies into three volume-limited samples separated in redshift, each of which acts as a… ▽ More We present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering 416 deg$^2$, and the SDSS DR11 spectroscopic galaxies spanning the redshift range $[0.15, 0.70]$. We subdivide the SDSS galaxies into three volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, $w_{\rm p}(R)$. We also use the measurements of the galaxy-galaxy weak lensing signal $ΔΣ(R)$ for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, $ξ_{\pm}(\vartheta)$, measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat $Λ$CDM cosmology. We use conservative scale cuts, $R>12$ and $8~h^{-1}$Mpc, for $ΔΣ$ and $w_{\rm p}$, respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. Our baseline analysis yields $S_8=0.775^{+0.043}_{-0.038}$ (68% C.I.) for the $Λ$CDM model, after marginalizing over uncertainties in other parameters. Our value of $S_8$ is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model. △ Less

Submitted 27 December, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 29 pages, 17 figures, 4 tables, version accepted for publication in PRD together with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

Journal ref: Phys. Rev. D 108, 123521 (2023)

Hyper Suprime-Cam Year 3 Results: Cosmology from Galaxy Clustering and Weak Lensing with HSC and SDSS using the Emulator Based Halo Model

Authors: Hironao Miyatake, Sunao Sugiyama, Masahiro Takada, Takahiro Nishimichi, Xiangchong Li, Masato Shirasaki, Surhud More, Yosuke Kobayashi, Atsushi J. Nishizawa, Markus M. Rau, Tianqing Zhang, Ryuichi Takahashi, Roohi Dalal, Rachel Mandelbaum, Michael A. Strauss, Takashi Hamana, Masamune Oguri, Ken Osato, Wentao Luo, Arun Kannawadi, Bau-Ching Hsieh, Robert Armstrong, Yutaka Komiyama, Robert H. Lupton, Nate B. Lust , et al. (9 additional authors not shown)

Abstract: We present cosmology results from a blinded joint analysis of cosmic shear, $ξ_{\pm}(\vartheta)$, galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected galaxy clustering, $w_{\rm p}(R)$, measured from the Hyper Suprime-Cam three-year (HSC-Y3) shape catalog and the Sloan Digital Sky Survey (SDSS) DR11 spectroscopic galaxy catalog - a 3$\times$2pt cosmology analysis. We define luminosity-cut samples… ▽ More We present cosmology results from a blinded joint analysis of cosmic shear, $ξ_{\pm}(\vartheta)$, galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected galaxy clustering, $w_{\rm p}(R)$, measured from the Hyper Suprime-Cam three-year (HSC-Y3) shape catalog and the Sloan Digital Sky Survey (SDSS) DR11 spectroscopic galaxy catalog - a 3$\times$2pt cosmology analysis. We define luminosity-cut samples of SDSS galaxies to serve as the tracers of $w_{\rm p}$ and as the lens samples for $Δ\!Σ$ in three spectroscopic redshift bins spanning the range $0.15<z<0.7$. For the $ξ_{\pm}$ and $Δ\!Σ$ measurements, we use a single source sample over 416 deg$^2$, selected from HSC-Y3 based on having photometric redshifts (photo-$z$) greater than 0.75. For cosmological parameter inference, we use Dark Emulator combined with a halo occupation distribution prescription to model $w_{\rm p}$ and $Δ\!Σ$ down to quasi-nonlinear scales. In our baseline analysis we employ an uninformative flat prior of the residual photo-$z$ error to model a residual bias in the mean redshift of HSC source galaxies. We obtain a robust constraint on the cosmological parameters for the flat $Λ$CDM model: $S_8=σ_8(Ω_{\rm m}/0.3)^{0.5}=0.763^{+0.040}_{-0.036}$ (68% C.I.), or the best-constrained parameter given by $S'_8=σ_8(Ω_{\rm m}/0.3)^{0.22}=0.721\pm 0.028$, determined with about 4% fractional precision. Our HSC-Y3 data exhibits about 2.5$σ$ tension with the Planck inferred $S_8$ value for the $Λ$CDM model, and hints at a non-zero residual photo-$z$ bias implying that the true mean redshift of the HSC galaxies at $z\gtrsim 0.75$ is higher than that implied by the original photo-$z$ estimates. △ Less

Submitted 6 April, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 39 pages, 25 figures, 5 tables, for coordinated submission to PRD with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/. arXiv admin note: text overlap with arXiv:2111.02419

Hyper Suprime-Cam Year 3 Results: Measurements of Clustering of SDSS-BOSS Galaxies, Galaxy-Galaxy Lensing and Cosmic Shear

Authors: Surhud More, Sunao Sugiyama, Hironao Miyatake, Markus Michael Rau, Masato Shirasaki, Xiangchong Li, Atsushi J. Nishizawa, Ken Osato, Tianqing Zhang, Masahiro Takada, Takashi Hamana, Ryuichi Takahashi, Roohi Dalal, Rachel Mandelbaum, Michael A. Strauss, Yosuke Kobayashi, Takahiro Nishimichi, Masamune Oguri, Wentao Luo, Arun Kannawadi, Bau-Ching Hsieh, Robert Armstrong, James Bosch, Yutaka Komiyama, Robert H. Lupton , et al. (9 additional authors not shown)

Abstract: We use the Sloan Digital Sky Survey (SDSS) BOSS galaxies and their overlap with approximately 416 sq. degree of deep $grizy$-band imaging from the Subaru Hyper Suprime-Cam Survey (HSC). We measure three two-point correlations that form the basis of the cosmological inference presented in our companion papers, Miyatake et al. and Sugiyama et al. We use three approximately volume limited subsamples… ▽ More We use the Sloan Digital Sky Survey (SDSS) BOSS galaxies and their overlap with approximately 416 sq. degree of deep $grizy$-band imaging from the Subaru Hyper Suprime-Cam Survey (HSC). We measure three two-point correlations that form the basis of the cosmological inference presented in our companion papers, Miyatake et al. and Sugiyama et al. We use three approximately volume limited subsamples of spectroscopic galaxies by their $i$-band magnitude from the SDSS-BOSS: LOWZ (0.1<z<0.35), CMASS1 (0.43<z<0.55) and CMASS2 (0.55<z<0.7), respectively. We present high signal-to-noise ratio measurements of the projected correlation functions of these galaxies, which is expected to be proportional to the matter correlation function times the bias of galaxies on large scales. In order to break the degeneracy between the amplitude of the matter correlation and the bias of these galaxies, we use the distortions of the shapes of galaxies in HSC due to weak gravitational lensing, to measure the galaxy-galaxy lensing signal, which probes the galaxy-matter cross-correlation of the SDSS-BOSS galaxies. We also measure the cosmic shear correlation functions from HSC galaxies which is related to the projected matter correlation function. We demonstrate the robustness of our measurements with a variety of systematic tests. Our use of a single sample of HSC source galaxies is crucial to calibrate any residual systematic biases in the inferred redshifts of our galaxies. We also describe the construction of a suite of mocks: i) spectroscopic galaxy catalogs which obey the clustering and abundance of each of the three SDSS-BOSS subsamples, and ii) galaxy shape catalogs which obey the footprint of the HSC survey and have been appropriately sheared by the large-scale structure expected in a $Λ$-CDM model. We use these mock catalogs to compute the covariance of each of our observables. △ Less

Submitted 16 November, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 31 pages, 24 figures, version accepted for publication in PRD together with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Two-point Correlation Functions

Authors: Xiangchong Li, Tianqing Zhang, Sunao Sugiyama, Roohi Dalal, Ryo Terasawa, Markus M. Rau, Rachel Mandelbaum, Masahiro Takada, Surhud More, Michael A. Strauss, Hironao Miyatake, Masato Shirasaki, Takashi Hamana, Masamune Oguri, Wentao Luo, Atsushi J. Nishizawa, Ryuichi Takahashi, Andrina Nicola, Ken Osato, Arun Kannawadi, Tomomi Sunayama, Robert Armstrong, James Bosch, Yutaka Komiyama, Robert H. Lupton , et al. (10 additional authors not shown)

Abstract: We perform a blinded cosmology analysis with cosmic shear two-point correlation functions (2PCFs) measured from more than 25 million galaxies in the Hyper Suprime-Cam three-year shear catalog in four tomographic redshift bins ranging from 0.3 to 1.5. After conservative masking and galaxy selection, the survey covers 416 deg$^2$ of the northern sky with an effective galaxy number density of 15 arcm… ▽ More We perform a blinded cosmology analysis with cosmic shear two-point correlation functions (2PCFs) measured from more than 25 million galaxies in the Hyper Suprime-Cam three-year shear catalog in four tomographic redshift bins ranging from 0.3 to 1.5. After conservative masking and galaxy selection, the survey covers 416 deg$^2$ of the northern sky with an effective galaxy number density of 15 arcmin$^{-2}$ over the four redshift bins. The 2PCFs adopted for cosmology analysis are measured in the angular range: $7.1 < θ/{\rm arcmin} < 56.6$ for $ξ_+$ and $31.2 <θ/{\rm arcmin} < 248$ for $ξ_-$, with a total signal-to-noise ratio of 26.6. We apply a conservative, wide, flat prior on the photometric redshift errors on the last two tomographic bins, and the relative magnitudes of the cosmic shear amplitude across four redshift bins allow us to calibrate the photometric redshift errors. With this flat prior on redshift errors, we find $Ω_{\rm m}=0.256_{-0.044}^{+0.056}$ and $S_8\equiv σ_8 \sqrt{Ω_{\rm m}/0.3}=0.769_{-0.034}^{+0.031}$ (both 68\% CI) for a flat $Λ$ cold dark matter cosmology. We find, after unblinding, that our constraint on $S_8$ is consistent with the Fourier space cosmic shear and the 3$\times$2pt analyses on the same HSC dataset. We carefully study the potential systematics from astrophysical and systematic model uncertainties in our fiducial analysis using synthetic data, and report no biases (including projection bias in the posterior space) greater than $0.5σ$ in the estimation of $S_8$. Our analysis hints that the mean redshifts of the two highest tomographic bins are higher than initially estimated. In addition, a number of consistency tests are conducted to assess the robustness of our analysis. Comparing our result with Planck-2018 cosmic microwave background observations, we find a ~$2σ$ tension for the $Λ$CDM model. △ Less

Submitted 30 November, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 38 pages, 32 figures, 4 tables (PRD in press.)

Hyper Suprime-Cam Year 3 Results: Cosmology from Cosmic Shear Power Spectra

Authors: Roohi Dalal, Xiangchong Li, Andrina Nicola, Joe Zuntz, Michael A. Strauss, Sunao Sugiyama, Tianqing Zhang, Markus M. Rau, Rachel Mandelbaum, Masahiro Takada, Surhud More, Hironao Miyatake, Arun Kannawadi, Masato Shirasaki, Takanori Taniguchi, Ryuichi Takahashi, Ken Osato, Takashi Hamana, Masamune Oguri, Atsushi J. Nishizawa, Andrés A. Plazas Malagón, Tomomi Sunayama, David Alonso, Anže Slosar, Robert Armstrong , et al. (13 additional authors not shown)

Abstract: We measure weak lensing cosmic shear power spectra from the three-year galaxy shear catalog of the Hyper Suprime-Cam (HSC) Subaru Strategic Program imaging survey. The shear catalog covers $416 \ \mathrm{deg}^2$ of the northern sky, with a mean $i$-band seeing of 0.59 arcsec and an effective galaxy number density of 15 $\mathrm{arcmin}^{-2}$ within our adopted redshift range. With an $i$-band magn… ▽ More We measure weak lensing cosmic shear power spectra from the three-year galaxy shear catalog of the Hyper Suprime-Cam (HSC) Subaru Strategic Program imaging survey. The shear catalog covers $416 \ \mathrm{deg}^2$ of the northern sky, with a mean $i$-band seeing of 0.59 arcsec and an effective galaxy number density of 15 $\mathrm{arcmin}^{-2}$ within our adopted redshift range. With an $i$-band magnitude limit of 24.5 mag, and four tomographic redshift bins spanning $0.3 \leq z_{\mathrm{ph}} \leq 1.5$ based on photometric redshifts, we obtain a high-significance measurement of the cosmic shear power spectra, with a signal-to-noise ratio of approximately 26.4 in the multipole range $300<\ell<1800$. The accuracy of our power spectrum measurement is tested against realistic mock shear catalogs, and we use these catalogs to get a reliable measurement of the covariance of the power spectrum measurements. We use a robust blinding procedure to avoid confirmation bias, and model various uncertainties and sources of bias in our analysis, including point spread function systematics, redshift distribution uncertainties, the intrinsic alignment of galaxies and the modeling of the matter power spectrum. For a flat $Λ$CDM model, we find $S_8 \equiv σ_8 (Ω_m/0.3)^{0.5} =0.776^{+0.032}_{-0.033}$, which is in excellent agreement with the constraints from the other HSC Year 3 cosmology analyses, as well as those from a number of other cosmic shear experiments. This result implies a $\sim$$2σ$-level tension with the Planck 2018 cosmology. We study the effect that various systematic errors and modeling choices could have on this value, and find that they can shift the best-fit value of $S_8$ by no more than $\sim$$0.5σ$, indicating that our result is robust to such systematics. △ Less

Submitted 4 April, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: 35 pages, 18 figures, 6 tables, for coordinated submission to PRD with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

Journal ref: Physical Review D, Volume 108, Issue 12, December 2023, article id.123519

Fringing Analysis and Simulation for the Vera C. Rubin Observatory's Legacy Survey of Space and Time

Authors: Zhiyuan Guo, Chris W. Walter, Craig Lage, Robert H. Lupton, The LSST Dark Energy Science Collaboration

Abstract: The presence of fringing in astronomical CCD images will impact photometric quality and measurements. Yet its impact on the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) has not been fully studied. We present a detailed study on fringing for Charge-Coupled Devices (CCDs) already implemented on the Rubin Observatory LSST Camera's focal plane. After making physical measurements… ▽ More The presence of fringing in astronomical CCD images will impact photometric quality and measurements. Yet its impact on the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) has not been fully studied. We present a detailed study on fringing for Charge-Coupled Devices (CCDs) already implemented on the Rubin Observatory LSST Camera's focal plane. After making physical measurements and knowing the compositions, we have developed a model for the e2v CCDs. We present a method to fit for the internal height variation of the epoxy layer within the sensors based on fringing measurements in a laboratory setting. This method is generic enough that it can be easily modified to work for other CCDs. Using the derived fringing model, we successfully reproduce comparable fringing amplitudes that match the observed levels in images taken by existing telescopes with different optical designs. This model is then used to forecast the expected level of fringing in a single LSST y-band sky background exposure with Rubin telescope optics in the presence of a realistic time varying sky spectrum. The predicted fringing amplitude in LSST images ranges from $0.04\%$ to $0.2\%$ depending on the location of a CCD on the focal plane. We find that the predicted variation in surface brightness caused by fringing in LSST y-band skybackground images is about $0.6\ μ\rm{Jy}\ \rm{arcsec}^{-2}$, which is 40 times larger than the current measurement error. We conclude that it is necessary to include fringing correction in the Rubin's LSST image processing pipeline. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 17 pages, 16 figures, submitted to PASP

PSFs of coadded images

Authors: Rachel Mandelbaum, Mike Jarvis, Robert H. Lupton, James Bosch, Arun Kannawadi, Michael D. Murphy, Tianqing Zhang, the LSST Dark Energy Science Collaboration

Abstract: We provide a detailed exploration of the connection between choice of coaddition schemes and the point-spread function (PSF) of the resulting coadded images. In particular, we investigate what properties of the coaddition algorithm lead to the final coadded image having a well-defined PSF. The key elements of this discussion are as follows: 1. We provide an illustration of how linear coaddition… ▽ More We provide a detailed exploration of the connection between choice of coaddition schemes and the point-spread function (PSF) of the resulting coadded images. In particular, we investigate what properties of the coaddition algorithm lead to the final coadded image having a well-defined PSF. The key elements of this discussion are as follows: 1. We provide an illustration of how linear coaddition schemes can produce a coadd that lacks a well-defined PSF even for relatively simple scenarios and choices of weight functions. 2. We provide a more formal demonstration of the fact that a linear coadd only has a well-defined PSF in the case that either (a) each input image has the same PSF or (b) the coadd is produced with weights that are independent of the signal. 3. We discuss some reasons that two plausible nonlinear coaddition algorithms (median and clipped-mean) fail to produce a consistent PSF profile for stars. 4. We demonstrate that all nonlinear coaddition procedures fail to produce a well-defined PSF for extended objects. In the end, we conclude that, for any purpose where a well-defined PSF is desired, one should use a linear coaddition scheme with weights that do not correlate with the signal and are approximately uniform across typical objects of interest. △ Less

Submitted 2 February, 2023; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: 13 pages, 4 figures; pedagogical article; v2 accepted for publication in the Open Journal of Astrophysics

Journal ref: The Open Journal of Astrophysics, Vol. 6, 2023

Faro: A framework for measuring the scientific performance of petascale Rubin Observatory data products

Authors: Leanne P. Guy, Keith Bechtol, Jeffrey L. Carlin, Erik Dennihy, Peter S. Ferguson, K. Simon Krughoff, Robert H. Lupton, Colin T. Slater, Krzysztof Findeisen, Arun Kannawadi, Lee S. Kelvin, Nate B. Lust, Lauren A. MacArthur, Michael N. Martinez, Sophie L. Reed, Dan S. Taranu, W. Michael Wood-Vasey

Abstract: The Vera C. Rubin Observatory will advance many areas of astronomy over the next decade with its unique wide-fast-deep multi-color imaging survey, the Legacy Survey of Space and Time (LSST). The LSST will produce approximately 20TB of raw data per night, which will be automatically processed by the LSST Science Pipelines to generate science-ready data products -- processed images, catalogs and ale… ▽ More The Vera C. Rubin Observatory will advance many areas of astronomy over the next decade with its unique wide-fast-deep multi-color imaging survey, the Legacy Survey of Space and Time (LSST). The LSST will produce approximately 20TB of raw data per night, which will be automatically processed by the LSST Science Pipelines to generate science-ready data products -- processed images, catalogs and alerts. To ensure that these data products enable transformative science with LSST, stringent requirements have been placed on their quality and scientific fidelity, for example on image quality and depth, astrometric and photometric performance, and object recovery completeness. In this paper we introduce faro, a framework for automatically and efficiently computing scientific performance metrics on the LSST data products for units of data of varying granularity, ranging from single-detector to full-survey summary statistics. By measuring and monitoring metrics, we are able to evaluate trends in algorithmic performance and conduct regression testing during development, compare the performance of one algorithm against another, and verify that the LSST data products will meet performance requirements by comparing to specifications. We present initial results using faro to characterize the performance of the data products produced on simulated and precursor data sets, and discuss plans to use faro to verify the performance of the LSST commissioning data products. △ Less

Submitted 30 June, 2022; originally announced June 2022.

Comments: 12 pages, 9 figues, Proc. SPIE 12189, Software and Cyberinfrastructure for Astronomy VII, Montreal, CA, July 2022

Report number: SPIE Paper Number: 12189-24

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). XVI. 69 New Quasars at 5.8 < z < 7.0

Authors: Yoshiki Matsuoka, Kazushi Iwasawa, Masafusa Onoue, Takuma Izumi, Nobunari Kashikawa, Michael A. Strauss, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, John D. Silverman, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Rikako Ishimoto, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Chien-Hsiu Lee, Robert H. Lupton , et al. (19 additional authors not shown)

Abstract: We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z… ▽ More We present the spectroscopic discovery of 69 quasars at 5.8 < z < 7.0, drawn from the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) imaging survey data. This is the 16th publication from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, and completes identification of all but the faintest candidates (i.e., i-band dropouts with zAB < 24 and y-band detections, and z-band dropouts with yAB < 24) with Bayesian quasar probability Pq > 0.1 in the HSC-SSP third public data release (PDR3). The sample reported here also includes three quasars with Pq < 0.1 at z ~ 6.6, which we selected in an effort to completely cover the reddest point sources with simple color cuts. The number of high-z quasars discovered in SHELLQs has now grown to 162, including 23 type-II quasar candidates. This paper also presents identification of seven galaxies at 5.6 < z < 6.7, an [O III] emitter at z = 0.954, and 31 Galactic cool stars and brown dwarfs. High-z quasars and galaxies comprise 75 % and 16 % respectively of all the spectroscopic SHELLQs objects that pass our latest selection algorithm with the PDR3 photometry. That is, a total of 91 % of the objects lie at z > 5.6. This demonstrates that the algorithm has very high efficiency, even though we are probing an unprecedentedly low-luminosity population down to M1450 ~ -21 mag. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series

HSC Year 1 cosmology results with the minimal bias method: HSC$\times$BOSS galaxy-galaxy weak lensing and BOSS galaxy clustering

Authors: Sunao Sugiyama, Masahiro Takada, Hironao Miyatake, Takahiro Nishimichi, Masato Shirasaki, Yosuke Kobayashi, Surhud More, Ryuichi Takahashi, Ken Osato, Masamune Oguri, Jean Coupon, Chiaki Hikage, Bau-Ching Hsieh, Yotaka Komiyama, Alexie Leauthaud, Xiangchong Li, Wentao Luo, Robert H. Lupton, Hitoshi Murayama, Atsushi J. Nishizawa, Youngsoo Park, Paul A. Price, Melanie Simet, Joshua S. Speagle, Michael A. Strauss , et al. (1 additional authors not shown)

Abstract: We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.15<z<0.7$. We use luminosity-limited samples as lens samples for $Δ\!Σ$ and as large-scale structure tracers for $w_\mathrm{p}$ in three red… ▽ More We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $Δ\!Σ(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.15<z<0.7$. We use luminosity-limited samples as lens samples for $Δ\!Σ$ and as large-scale structure tracers for $w_\mathrm{p}$ in three redshift bins, and use the HSC-Y1 galaxy catalog to define a secure sample of source galaxies at $z_\mathrm{ph}>0.75$ for the $Δ\!Σ$ measurements, selected based on their photometric redshifts. For theoretical template, we use the "minimal bias" model for the cosmological clustering observables for the flat $Λ$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 $Δ\!Σ(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}$, $σ_8=0.85^{+0.16}_{-0.11}$, and $Ω_\mathrm{m}=0.283^{+0.12}_{-0.035}$ for the flat $Λ$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) △ Less

Submitted 21 November, 2021; originally announced November 2021.

Comments: 24 pages, 19 figures, 4 tables, to be submitted to Phys. Rev. D

Report number: IPMU21-0078, YITP-21-126

Reaching for the Edge I: Probing the Outskirts of Massive Galaxies with HSC, DECaLS, SDSS, and Dragonfly

Authors: Jiaxuan Li, Song Huang, Alexie Leauthaud, John Moustakas, Shany Danieli, Jenny E. Greene, Roberto Abraham, Felipe Ardila, Erin Kado-Fong, Deborah Lokhorst, Robert Lupton, Paul Price

Abstract: 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… ▽ More 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<0.05$, Dragonfly has the best control of systematics, reaching surface brightness levels of $μ_r \sim 30$ mag/arcsec$^{2}$. At $0.19<z<0.50$, HSC can reliably recover surface brightness profiles to $μ_{r} \sim 28.5$ mag/arcsec$^{2}$ reaching $R=100 - 150$ kpc. DECaLS surface brightness profiles show good agreement with HSC but are noisier at large radii. The median profiles of galaxy ensembles in both HSC and DECaLS reach $R > 200$ kpc without significant bias. At $0.19<z<0.50$, DECaLS and HSC measurements of the stellar mass contained within 100 kpc agree within 0.05 dex. Finally, we use weak gravitational lensing to show that measurements of outer light with DECaLS at $0.19<z<0.50$ show a similar promise as HSC as a low scatter proxy of halo mass. The tests and results from this paper represent an important step forward for accurate measurements of the outer light of massive galaxies and demonstrate that outer light measurements from DECam imaging will be a promising method for finding galaxy clusters for DES and DESI. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 19 pages, 10 figures. Submitted to MNRAS. Main results are shown in Fig 5-10

Cosmological inference from the emulator based halo model II: Joint analysis of galaxy-galaxy weak lensing and galaxy clustering from HSC-Y1 and SDSS

Authors: Hironao Miyatake, Sunao Sugiyama, Masahiro Takada, Takahiro Nishimichi, Masato Shirasaki, Yosuke Kobayashi, Rachel Mandelbaum, Surhud More, Masamune Oguri, Ken Osato, Youngsoo Park, Ryuichi Takahashi, Jean Coupon, Chiaki Hikage, Bau-Ching Hsieh, Alexie Leauthaud, Xiangchong Li, Wentao Luo, Robert H. Lupton, Satoshi Miyazaki, Hitoshi Murayama, Atsushi J. Nishizawa, Paul A. Price, Melanie Simet, Joshua S. Speagle , et al. (3 additional authors not shown)

Abstract: We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($Δ\!Σ$) 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<z<0.7$. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of… ▽ More We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($Δ\!Σ$) 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<z<0.7$. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of $w_{\rm p}$ in three spectroscopic redshift bins, and as the lens samples for $Δ\!Σ$. For the $Δ\!Σ$ measurements, we select a single sample of 4 million source galaxies over 140 deg$^2$ from HSC-Y1 with photometric redshifts (photo-$z$) greater than 0.75, enabling a better handle of photo-$z$ errors by comparing the $Δ\!Σ$ amplitudes for the three lens redshift bins. For cosmological parameter inference, we use an input galaxy-halo connection model built on the {\tt Dark Emulator} package with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the $Δ\!Σ$ and $w_{\rm p}$ measurements on scales from $R\simeq 3$ and $2\,h^{-1}{\rm Mpc}$, respectively, up to $30\,h^{-1}{\rm Mpc}$ assuming a flat $Λ$CDM cosmology. With various tests using mock catalogs described in Miyatake et al. (2021), we show that any bias in the clustering amplitude $S_8\equiv σ_8(Ω_{\rm m}/0.3)^{0.5}$ due to uncertainties in the galaxy-halo connection is less than $\sim50$\% of the statistical uncertainty of $S_8$, {\it unless} the assembly bias effect is unexpectedly large. Our best-fit models have $S_8=0.795^{+0.049}_{-0.042}$ (mode and 68\% credible interval) for the flat $Λ$CDM model; we find tighter constraints on the quantity $S_8(α=0.17)\equivσ_8(Ω_{\rm m}/0.3)^{0.17} =0.745^{+0.039}_{-0.031}$. (abriged) △ Less

Submitted 29 November, 2021; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: 44 pages, 30 figures, 6 tables, submitted to Phys. Rev. D, YouTube video is available at https://www.youtube.com/watch?v=jRhQvuee97Q

Report number: IPMU21-0069, YITP-21-125

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Yusra AlSayyad, Makoto Ando, Robert Armstrong, James Bosch, Eiichi Egami, Hisanori Furusawa, Junko Furusawa, Sumiko Harasawa, Yuichi Harikane, Bau-Ching Hsieh, Hiroyuki Ikeda, Kei Ito, Ikuru Iwata, Tadayuki Kodama, Michitaro Koike, Mitsuru Kokubo, Yutaka Komiyama, Xiangchong Li, Yongming Liang, Yen-Ting Lin, Robert H. Lupton, Nate B Lust, Lauren A. MacArthur, Ken Mawatari , et al. (42 additional authors not shown)

Abstract: The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of ob… ▽ More The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 25 pages, 19 figures, submitted to PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

DESC DC2 Data Release Note

Authors: LSST Dark Energy Science Collaboration, Bela Abolfathi, Robert Armstrong, Humna Awan, Yadu N. Babuji, Franz Erik Bauer, George Beckett, Rahul Biswas, Joanne R. Bogart, Dominique Boutigny, Kyle Chard, James Chiang, Johann Cohen-Tanugi, Andrew J. Connolly, Scott F. Daniel, Seth W. Digel, Alex Drlica-Wagner, Richard Dubois, Eric Gawiser, Thomas Glanzman, Salman Habib, Andrew P. Hearin, Katrin Heitmann, Fabio Hernandez, Renée Hložek , et al. (32 additional authors not shown)

Abstract: In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with th… ▽ More In preparation for cosmological analyses of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), the LSST Dark Energy Science Collaboration (LSST DESC) has created a 300 deg$^2$ simulated survey as part of an effort called Data Challenge 2 (DC2). The DC2 simulated sky survey, in six optical bands with observations following a reference LSST observing cadence, was processed with the LSST Science Pipelines (19.0.0). In this Note, we describe the public data release of the resulting object catalogs for the coadded images of five years of simulated observations along with associated truth catalogs. We include a brief description of the major features of the available data sets. To enable convenient access to the data products, we have developed a web portal connected to Globus data services. We describe how to access the data and provide example Jupyter Notebooks in Python to aid first interactions with the data. We welcome feedback and questions about the data release via a GitHub repository. △ Less

Submitted 13 June, 2022; v1 submitted 12 January, 2021; originally announced January 2021.

Comments: 25 pages, 3 figures; 9 tables. A detailed changelog can be found in Appendix A. To obtain data, visit the DESC Data Portal at https://data.lsstdesc.org/

The LSST DESC DC2 Simulated Sky Survey

Authors: LSST Dark Energy Science Collaboration, Bela Abolfathi, David Alonso, Robert Armstrong, Éric Aubourg, Humna Awan, Yadu N. Babuji, Franz Erik Bauer, Rachel Bean, George Beckett, Rahul Biswas, Joanne R. Bogart, Dominique Boutigny, Kyle Chard, James Chiang, Chuck F. Claver, Johann Cohen-Tanugi, Céline Combet, Andrew J. Connolly, Scott F. Daniel, Seth W. Digel, Alex Drlica-Wagner, Richard Dubois, Emmanuel Gangler, Eric Gawiser , et al. (55 additional authors not shown)

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

Submitted 26 January, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 39 pages, 19 figures, version accepted for publication in ApJS

A comparative study of satellite galaxies in Milky Way-like galaxies from HSC, DECaLS and SDSS

Authors: Wenting Wang, Masahiro Takada, Xiangchong Li, Scott G. Carlsten, Ting-Wen Lan, Jingjing Shi, Hironao Miyatake, Surhud More, Rachael L. Beaton, Robert Lupton, Yen-Ting Lin, Tian Qiu, Wentao Luo

Abstract: We conduct a comprehensive and statistical study of the luminosity functions (LFs) for satellite galaxies, by counting photometric galaxies from HSC, DECaLS and SDSS around isolated central galaxies (ICGs) and paired galaxies from the SDSS/DR7 spectroscopic sample. Results of different surveys show very good agreement. The satellite LFs can be measured down to $M_V\sim-10$, and for central primary… ▽ More We conduct a comprehensive and statistical study of the luminosity functions (LFs) for satellite galaxies, by counting photometric galaxies from HSC, DECaLS and SDSS around isolated central galaxies (ICGs) and paired galaxies from the SDSS/DR7 spectroscopic sample. Results of different surveys show very good agreement. The satellite LFs can be measured down to $M_V\sim-10$, and for central primary galaxies as small as $8.5<\log_{10}M_\ast/M_\odot<9.2$ and $9.2<\log_{10}M_\ast/M_\odot<9.9$, implying there are on average 3--8 satellites with $M_V<-10$ around LMC-mass ICGs. The bright end cutoff of satellite LFs and the satellite abundance are both sensitive to the magnitude gap between the primary and its companions, indicating galaxy systems with larger magnitude gaps are on average hosted by less massive dark matter haloes. By selecting primaries with stellar mass similar to our MW, we discovered that i) the averaged satellite LFs of ICGs with different magnitude gaps to their companions and of galaxy pairs with different colour or colour combinations all show steeper slopes than the MW satellite LF; ii) there are on average more satellites with $-15<M_V<-10$ than those in our MW; iii) there are on average 1.5 to 2.5 satellites with $M_V<-16$ around ICGs, consistent with our MW; iv) even after accounting for the large scatter predicted by numerical simulations, the MW satellite LF is uncommon at $M_V>-12$. Hence the MW and its satellite system are statistically atypical of our sample of MW-mass systems. In consequence, our MW is not a good representative of other MW-mass galaxies. Strong cosmological implications based on only MW satellites await additional discoveries of fainter satellites in extra-galactic systems. Interestingly, the MW satellite LF is typical among other MW-mass systems within 40~Mpc in the local Universe, perhaps implying the Local Volume is an under-dense region. △ Less

Submitted 6 November, 2020; v1 submitted 15 September, 2020; originally announced September 2020.

Comments: accepted by MNRAS, with more detailed discussions about the scatter

Proper motion measurements for stars up to $100$ kpc with Subaru HSC and SDSS Stripe 82

Authors: Tian Qiu, Wenting Wang, Masahiro Takada, Naoki Yasuda, Željko Ivezić, Robert H. Lupton, Masashi Chiba, Miho Ishigaki, Yutaka Komiyama

Abstract: We present proper motion measurements for more than $0.55$ million main-sequence stars, by comparing astrometric positions of matched stars between the multi-band imaging datasets from the Hyper Suprime-Cam (HSC) Survey and the SDSS Stripe 82. In doing this we use $3$ million galaxies to recalibrate the astrometry and set up a common reference frame between the two catalogues. The exquisite depth… ▽ More We present proper motion measurements for more than $0.55$ million main-sequence stars, by comparing astrometric positions of matched stars between the multi-band imaging datasets from the Hyper Suprime-Cam (HSC) Survey and the SDSS Stripe 82. In doing this we use $3$ million galaxies to recalibrate the astrometry and set up a common reference frame between the two catalogues. The exquisite depth and the nearly $12$ years of time baseline between HSC and SDSS enable high-precision measurements of statistical proper motions for stars down to $i\simeq 24$. A validation of our method is demonstrated by the agreement with the $Gaia$ proper motions, to the precision better than $0.1$ mas yr$^{-1}$. To retain the precision, we make a correction of the subtle effects due to the differential chromatic refraction in the SDSS images based on the comparison with the $Gaia$ proper motions against colour of stars, which is validated using the SDSS spectroscopic quasars. Combining with the photometric distance estimates for individual stars based on the precise HSC photometry, we show a significant detection of the net proper motions for stars in each bin of distance out to $100$ kpc. The two-component tangential velocities after subtracting the apparent motions due to our own motion display rich phase-space structures including a clear signature of the Sagittarius stream in the halo region of distance range $[10,\ 35]$ kpc. We also measure the tangential velocity dispersion in the distance range $5-20$ kpc and find that the data are consistent with a constant isotropic dispersion of $80\pm 10\ {\rm km/s}$. More distant stars appear to have random motions with respect to the Galactic centre on average. △ Less

Submitted 24 December, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: accepted by MNRAS, add more details

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

Authors: F. Javier Sánchez, Chris W. Walter, Humna Awan, James Chiang, Scott F. Daniel, Eric Gawiser, Tom Glanzman, David P. Kirkby, Rachel Mandelbaum, Anže Slosar, W. Michael Wood-Vasey, Yusra AlSayyad, Colin J. Burke, Seth W. Digel, Mike Jarvis, Tony Johnson, Heather Kelly, Simon Krughoff, Robert H. Lupton, Phil J. Marshall, John R. Peterson, Paul A. Price, Glenn Sembroski, Brian Van Klaveren, Matthew P. Wiesner , et al. (1 additional authors not shown)

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

Submitted 5 July, 2020; v1 submitted 3 January, 2020; originally announced January 2020.

Comments: 21 pages, 21 figures. Accepted at MNRAS

Report number: FERMILAB-PUB-19-618-AE

Morphological Star-Galaxy Separation

Authors: Colin T. Slater, Željko Ivezić, Robert H. Lupton

Abstract: We discuss the statistical foundations of morphological star-galaxy separation. We show that many of the star-galaxy separation metrics in common use today (e.g. by SDSS or SExtractor) are closely related both to each other, and to the model odds ratio derived in a Bayesian framework by Sebok (1979). While the scaling of these algorithms with the noise properties of the sources varies, these diffe… ▽ More We discuss the statistical foundations of morphological star-galaxy separation. We show that many of the star-galaxy separation metrics in common use today (e.g. by SDSS or SExtractor) are closely related both to each other, and to the model odds ratio derived in a Bayesian framework by Sebok (1979). While the scaling of these algorithms with the noise properties of the sources varies, these differences do not strongly differentiate their performance. We construct a model of the performance of a star-galaxy separator in a realistic survey to understand the impact of observational signal-to-noise ratio (or equivalently, 5-sigma limiting depth) and seeing on classification performance. The model quantitatively demonstrates that, assuming realistic densities and angular sizes of stars and galaxies, 10% worse seeing can be compensated for by approximately 0.4 magnitudes deeper data to achieve the same star-galaxy classification performance. We discuss how to probabilistically combine multiple measurements, either of the same type (e.g., subsequent exposures), or differing types (e.g., multiple bandpasses), or differing methodologies (e.g., morphological and color-based classification). These methods are increasingly important for observations at faint magnitudes, where the rapidly rising number density of small galaxies makes star-galaxy classification a challenging problem. However, because of the significant role that the signal-to-noise ratio plays in resolving small galaxies, surveys with large-aperture telescopes, such as LSST, will continue to see improving star-galaxy separation as they push to these fainter magnitudes. △ Less

Submitted 17 December, 2019; originally announced December 2019.

Comments: Accepted to AJ

Joint Survey Processing of LSST, Euclid and WFIRST: Enabling a broad array of astrophysics and cosmology through pixel level combinations of datasets

Authors: R. Chary, G. Brammer, P. Capak, W. Dawson, A. Faisst, S. Fajardo-Acosta, H. C. Ferguson, C. J. Grillmair, S. Hemmati, A. Koekemoer, B. Lee, R. Lupton, S. Malhotra, P. Melchior, I. Momcheva, J. Newman, J. Masiero, R. Paladini, A. Prakash, J. Rhodes, B. Rusholme, M. Schneider, N. Stickley, A. Smith, W. M. Wood-Vasey , et al. (1 additional authors not shown)

Abstract: Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper… ▽ More Joint survey processing (JSP) is the pixel level combination of LSST, Euclid, and WFIRST datasets. By combining the high spatial resolution of the space-based datasets with deep, seeing-limited, ground-based images in the optical bands, systematics like source confusion and astrometric mismatch can be addressed to derive the highest precision optical/infrared photometric catalogs. This white paper highlights the scientific motivation, computational and algorithmic needs to build joint pixel level processing capabilities, which the individual projects by themselves will not be able to support. Through this white paper, we request that the Astro2020 decadal committee recognize the JSP effort as a multi-agency project with the natural outcome being a collaborative effort among groups which are normally supported by a single agency. JSP will allow the U.S. (and international) astronomical community to manipulate the flagship data sets and undertake innovative science investigations ranging from solar system object characterization, exoplanet detections, nearby galaxy rotation rates and dark matter properties, to epoch of reionization studies. It will also result in the ultimate constraints on cosmological parameters and the nature of dark energy, with far smaller uncertainties and a better handle on systematics than by any one survey alone. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Comments: Astro2020 APC White Paper, 11 pages

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). X. Discovery of 35 Quasars and Luminous Galaxies at 5.7 $\le$ z $\le$ 7.0

Authors: Yoshiki Matsuoka, Kazushi Iwasawa, Masafusa Onoue, Nobunari Kashikawa, Michael A. Strauss, Chien-Hsiu Lee, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Takuma Izumi, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Shuhei Koyama, Robert H. Lupton, Takeo Minezaki , et al. (22 additional authors not shown)

Abstract: We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has… ▽ More We report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: ApJ, in press. arXiv admin note: text overlap with arXiv:1803.01861

Journal ref: ApJ 883:183, 2019

Future Uses of the LSST Facility: Input from the LSST Project Science Team

Authors: S. Kahn, Robert Blum, Chuck Claver, Andrew Connolly, Leanne Guy, Željko Ivezić, Robert H. Lupton, William O'Mullane, Steven Ritz, Michael Strauss, Christopher Stubbs, Sandrine J. Thomas, J. Anthony Tyson, Bo Xin

Abstract: In this white paper, we discuss future uses of the LSST facility after the planned 10-year survey is complete. We expect the LSST survey to profoundly affect the scientific landscape over the next ten years, and it is likely that unexpected discoveries may drive its future scientific program. We discuss various operations and instrument options that could be considered for an extended LSST mission… ▽ More In this white paper, we discuss future uses of the LSST facility after the planned 10-year survey is complete. We expect the LSST survey to profoundly affect the scientific landscape over the next ten years, and it is likely that unexpected discoveries may drive its future scientific program. We discuss various operations and instrument options that could be considered for an extended LSST mission beyond ten years. △ Less

Submitted 21 July, 2019; originally announced July 2019.

Comments: 10 pages, no figures. arXiv admin note: text overlap with arXiv:1905.04669

Boötes IV: A New Milky Way Satellite Discovered in the Subaru Hyper Suprime-Cam Survey and Implications for the Missing Satellite Problem

Authors: Daisuke Homma, Masashi Chiba, Yutaka Komiyama, Masayuki Tanaka, Sakurako Okamoto, Mikito Tanaka, Miho N. Ishigaki, Kohei Hayashi, Nobuo Arimoto, Scott G. Carlsten, Robert H. Lupton, Michael A. Strauss, Satoshi Miyazaki, Gabriel Torrealba, Shiang-Yu Wang, Hitoshi Murayama

Abstract: We report on the discovery of a new Milky Way (MW) satellite in Boötes based on data from the on-going Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP). This satellite, named Boötes IV, is the third ultra-faint dwarf that we have discovered in the HSC-SSP. We have identified a statistically significant (32.3$σ$) overdensity of stars having characteristics of a metal-poor, old stellar populat… ▽ More We report on the discovery of a new Milky Way (MW) satellite in Boötes based on data from the on-going Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP). This satellite, named Boötes IV, is the third ultra-faint dwarf that we have discovered in the HSC-SSP. We have identified a statistically significant (32.3$σ$) overdensity of stars having characteristics of a metal-poor, old stellar population. The distance to this stellar system is $D_{\odot}=209^{+20}_{-18}$ kpc with a $V$-band absolute magnitude of $M_V=-4.53^{+0.23}_{-0.21}$ mag. Boötes IV has a half-light radius of $r_h=462^{+98}_{-84}$ pc and an ellipticity of $0.64^{+0.05}_{-0.05}$, which clearly suggests that this is a dwarf satellite galaxy. We also found another overdensity that appears to be a faint globular cluster with $M_V=-0.20^{+0.59}_{-0.83}$ mag and $r_h=5.9^{+1.5}_{-1.3}$ pc located at $D_{\odot}=46^{+4}_{-4}$ kpc. Adopting the recent prediction for the total population of satellites in a MW-sized halo by Newton et al. (2018), which combined the characteristics of the observed satellites by SDSS and DES with the subhalos obtained in $Λ$CDM models, we estimate that there should be about two MW satellites at $M_V\le0$ in the $\sim676$ deg$^2$ covered by HSC-SSP, whereas that area includes six satellites. Thus, the observed number of satellites is larger than the theoretical prediction. On the face of it, we have a problem of too many satellites, instead of the well-known missing satellites problem whereby the $Λ$CDM theory overpredicts the number of satellites in a MW-sized halo. This may imply that the models need more refinements for the assignment of subhalos to satellites such as considering those found by the current deeper survey. [abridged] △ Less

Submitted 17 June, 2019; originally announced June 2019.

Comments: 13 pages, 11 figures, 2 tables, accepted for publication in PASJ

Cosmological constraints from cosmic shear two-point correlation functions with HSC survey first-year data

Authors: T. Hamana, M. Shirasaki, S. Miyazaki, C. Hikage, M. Oguri, S. More, R. Armstrong, A. Leauthaud, R. Mandelbaum, H. Miyatake, A. J. Nishizawa, M. Simet, M. Takada, H. Aihara, J. Bosch, Y. Komiyama, R. Lupton, H. Murayama, M. A. Strauss, M. Tanaka

Abstract: We present measurements of cosmic shear two-point correlation functions (TPCFs) from Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) first-year data, and derived cosmological constraints based on a blind analysis. The HSC first-year shape catalog is divided into four tomographic redshift bins ranging from $z=0.3$ to 1.5 with equal widths of $Δz =0.3$. The unweighted galaxy number densities in… ▽ More We present measurements of cosmic shear two-point correlation functions (TPCFs) from Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) first-year data, and derived cosmological constraints based on a blind analysis. The HSC first-year shape catalog is divided into four tomographic redshift bins ranging from $z=0.3$ to 1.5 with equal widths of $Δz =0.3$. The unweighted galaxy number densities in each tomographic bin are 5.9, 5.9, 4.3, and 2.4 arcmin$^{-2}$ from lower to higher redshifts, respectively. We adopt the standard TPCF estimators, $ξ_\pm$, for our cosmological analysis, given that we find no evidence of the significant B-mode shear. The TPCFs are detected at high significance for all ten combinations of auto- and cross-tomographic bins over a wide angular range, yielding a total signal-to-noise ratio of 19 in the angular ranges adopted in the cosmological analysis, $7'<θ<56'$ for $ξ_+$ and $28'<θ<178'$ for $ξ_-$. We perform the standard Bayesian likelihood analysis for cosmological inference from the measured cosmic shear TPCFs, including contributions from intrinsic alignment of galaxies as well as systematic effects from PSF model errors, shear calibration uncertainty, and source redshift distribution errors. We adopt a covariance matrix derived from realistic mock catalogs constructed from full-sky gravitational lensing simulations that fully account for survey geometry and measurement noise. For a flat $Λ$ cold dark matter model, we find $S_8 \equiv σ_8\sqrt{Ω_m/0.3}=0.823_{-0.028}^{+0.032}$, and $Ω_m=0.332_{-0.096}^{+0.050}$. We carefully check the robustness of the cosmological results against astrophysical modeling uncertainties and systematic uncertainties in measurements, and find that none of them has a significant impact on the cosmological constraints. △ Less

Submitted 10 January, 2022; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: Replaced with the bug-fixed version. The results presented in the original publication were found to be affected by bugs in the software used for the numerical computation. The erratum will be published in PASJ. In this version, revised results obtained from corrected computations are presented. Revised chains are available from http://th.nao.ac.jp/MEMBER/hamanatk/HSC16aCSTPCFbugfix/index.html

Second Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Yusra AlSayyad, Makoto Ando, Robert Armstrong, James Bosch, Eiichi Egami, Hisanori Furusawa, Junko Furusawa, Andy Goulding, Yuichi Harikane, Chiaki Hikage, Paul T. P. Ho, Bau-Ching Hsieh, Song Huang, Hiroyuki Ikeda, Masatoshi Imanishi, Kei Ito, Ikuru Iwata, Anton T. Jaelani, Ryota Kakuma, Kojiro Kawana, Satoshi Kikuta, Umi Kobayashi, Michitaro Koike, Yutaka Komiyama , et al. (40 additional authors not shown)

Abstract: This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially ob… ▽ More This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 22 August, 2019; v1 submitted 29 May, 2019; originally announced May 2019.

Comments: 26 pages, 23 figures, accepted for publication in PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

Discovery of the First Low-Luminosity Quasar at z > 7

Authors: Yoshiki Matsuoka, Masafusa Onoue, Nobunari Kashikawa, Michael A. Strauss, Kazushi Iwasawa, Chien-Hsiu Lee, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Takuma Izumi, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Shuhei Koyama, Robert H. Lupton, Takeo Minezaki , et al. (23 additional authors not shown)

Abstract: We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric l… ▽ More We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in ApJ Letters

Journal ref: ApJL 872:L2, 2019

An Overview of the LSST Image Processing Pipelines

Authors: James Bosch, Yusra AlSayyad, Robert Armstrong, Eric Bellm, Hsin-Fang Chiang, Siegfried Eggl, Krzysztof Findeisen, Merlin Fisher-Levine, Leanne P. Guy, Augustin Guyonnet, Željko Ivezić, Tim Jenness, Gábor Kovács, K. Simon Krughoff, Robert H. Lupton, Nate B. Lust, Lauren A. MacArthur, Joshua Meyers, Fred Moolekamp, Christopher B. Morrison, Timothy D. Morton, William O'Mullane, John K. Parejko, Andrés A. Plazas, Paul A. Price , et al. (9 additional authors not shown)

Abstract: The Large Synoptic Survey Telescope (LSST) is an ambitious astronomical survey with a similarly ambitious Data Management component. Data Management for LSST includes processing on both nightly and yearly cadences to generate transient alerts, deep catalogs of the static sky, and forced photometry light-curves for billions of objects at hundreds of epochs, spanning at least a decade. The algorithm… ▽ More The Large Synoptic Survey Telescope (LSST) is an ambitious astronomical survey with a similarly ambitious Data Management component. Data Management for LSST includes processing on both nightly and yearly cadences to generate transient alerts, deep catalogs of the static sky, and forced photometry light-curves for billions of objects at hundreds of epochs, spanning at least a decade. The algorithms running in these pipelines are individually sophisticated and interact in subtle ways. This paper provides an overview of those pipelines, focusing more on those interactions than the details of any individual algorithm. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: Submitted to proceedings for ADASS XXVIII

The stellar halo of isolated central galaxies in the Hyper Suprime-Cam imaging survey

Authors: Wenting Wang, Jiaxin Han, Alessandro Sonnenfeld, Naoki Yasuda, Xiangchong Li, Yipeng Jing, Surhud More, Paul A. Price, Robert Lupton, Eli Rykoff, David V. Stark, Ting-Wen Lan, Masahiro Takada, Song Huang, Wentao Luo, Neta A. Bahcall, Yutaka Komiyama

Abstract: We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise lev… ▽ More We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$σ$ significance. PSF also slightly flattens the measured colour profiles. △ Less

Submitted 4 May, 2019; v1 submitted 12 November, 2018; originally announced November 2018.

Comments: accepted by MNRAS - Significant changes have been made compared with the first version, including discussions on the extended PSF wings, robustness of our results to source detection and masking thresholds and more detailed investigations on the indications of positive colour gradients

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). V. Quasar Luminosity Function and Contribution to Cosmic Reionization at z = 6

Authors: Yoshiki Matsuoka, Michael A. Strauss, Nobunari Kashikawa, Masafusa Onoue, Kazushi Iwasawa, Ji-Jia Tang, Chien-Hsiu Lee, Masatoshi Imanishi, Tohru Nagao, Masayuki Akiyama, Naoko Asami, James Bosch, Hisanori Furusawa, Tomotsugu Goto, James E. Gunn, Yuichi Harikane, Hiroyuki Ikeda, Takuma Izumi, Toshihiro Kawaguchi, Nanako Kato, Satoshi Kikuta, Kotaro Kohno, Yutaka Komiyama, Robert H. Lupton, Takeo Minezaki , et al. (22 additional authors not shown)

Abstract: We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-… ▽ More We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $α= -1.23^{+0.44}_{-0.34}$, a bright-end slope $β= -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $Φ^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization. △ Less

Submitted 5 November, 2018; originally announced November 2018.

Comments: Accepted for publication in ApJ

Journal ref: ApJ 869:150, 2018

Cosmology from cosmic shear power spectra with Subaru Hyper Suprime-Cam first-year data

Authors: Chiaki Hikage, Masamune Oguri, Takashi Hamana, Surhud More, Rachel Mandelbaum, Masahiro Takada, Fabian Köhlinger, Hironao Miyatake, Atsushi J. Nishizawa, Hiroaki Aihara, Robert Armstrong, James Bosch, Jean Coupon, Anne Ducout, Paul Ho, Bau-Ching Hsieh, Yutaka Komiyama, François Lanusse, Alexie Leauthaud, Robert H. Lupton, Elinor Medezinski, Sogo Mineo, Shoken Miyama, Satoshi Miyazaki, Ryoma Murata , et al. (12 additional authors not shown)

Abstract: 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 co… ▽ More 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 $Λ$ cold dark matter ($Λ$CDM) model, we find $S_8\equiv σ_8(Ω_{\rm m}/0.3)^α=0.800^{+0.029}_{-0.028}$ for $α=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $α=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σ$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints. △ Less

Submitted 26 January, 2019; v1 submitted 24 September, 2018; originally announced September 2018.

Comments: 43 pages, 21 figures, accepted for publication in PASJ

A collimated beam projector for precise telescope calibration

Authors: Michael Coughlin, T. M. C. Abbott, Kairn Brannon, Chuck Claver, Peter Doherty, Merlin Fisher-Levine, Patrick Ingraham, Robert Lupton, Nicholas Mondrik, Christopher Stubbs

Abstract: The precise determination of the instrumental response function versus wavelength is a central ingredient in contemporary photometric calibration strategies. This typically entails propagating narrowband illumination through the system pupil, and comparing the detected photon rate across the focal plane to the amount of incident light as measured by a calibrated photodiode. However, stray light ef… ▽ More The precise determination of the instrumental response function versus wavelength is a central ingredient in contemporary photometric calibration strategies. This typically entails propagating narrowband illumination through the system pupil, and comparing the detected photon rate across the focal plane to the amount of incident light as measured by a calibrated photodiode. However, stray light effects and reflections/ghosting (especially on the edges of filter passbands) in the optical train constitute a major source of systematic uncertainty when using a flat-field screen as the illumination source. A collimated beam projector that projects a mask onto the focal plane of the instrument can distinguish focusing light paths from stray and scattered light, allowing for a precise determination of instrumental throughput. This paper describes the conceptual design of such a system, outlines its merits, and presents results from a prototype system used with the Dark Energy Camera wide field imager on the 4-meter Blanco telescope. A calibration scheme that blends results from flat-field images with collimated beam projector data to obtain the equivalent of an illumination correction at high spectral and angular resolution is also presented. In addition to providing a precise system throughput calibration, by monitoring the evolution of the intensity and behaviour of the ghosts in the optical system, the collimated beam projector can be used to track the evolution of the filter transmission properties and various anti-reflective coatings in the optical system. △ Less

Submitted 15 May, 2018; originally announced May 2018.

A Study of the Point Spread Function in SDSS Images

Authors: Bo Xin, Željko Ivezić, Robert H. Lupton, John R. Peterson, Peter Yoachim, R. Lynne Jones, Charles F. Claver, George Angeli

Abstract: We use SDSS imaging data in $ugriz$ passbands to study the shape of the point spread function (PSF) profile and the variation of its width with wavelength and time. We find that the PSF profile is well described by theoretical predictions based on von Kármán's turbulence theory. The observed PSF radial profile can be parametrized by only two parameters, the profile's full width at half maximum (FW… ▽ More We use SDSS imaging data in $ugriz$ passbands to study the shape of the point spread function (PSF) profile and the variation of its width with wavelength and time. We find that the PSF profile is well described by theoretical predictions based on von Kármán's turbulence theory. The observed PSF radial profile can be parametrized by only two parameters, the profile's full width at half maximum (FWHM) and a normalization of the contribution of an empirically determined "instrumental" PSF. The profile shape is very similar to the "double gaussian plus power-law wing" decomposition used by SDSS image processing pipeline, but here it is successfully modeled with two free model parameters, rather than six as in SDSS pipeline. The FWHM variation with wavelength follows the $λ^α$ power law, where $α\approx-0.3$ and is correlated with the FWHM itself. The observed behavior is much better described by von Kármán's turbulence theory with the outer scale parameter in the range 5$-$100 m, than by the Kolmogorov's turbulence theory. We also measure the temporal and angular structure functions for FWHM and compare them to simulations and results from literature. The angular structure function saturates at scales beyond 0.5$-$1.0 degree. The power spectrum of the temporal behavior is found to be broadly consistent with a damped random walk model with characteristic timescale in the range $\sim5-30$ minutes, though data show a shallower high-frequency behavior. The latter is well fit by a single power law with index in the range $-1.5$ to $-1.0$. A hybrid model is likely needed to fully capture both the low-frequency and high-frequency behavior of the temporal variations of atmospheric seeing. △ Less

Submitted 8 May, 2018; originally announced May 2018.

Comments: 16 pages, 9 figures, submitted to the Astronomical Journal

Wavelength Dependent PSFs and their impact on Weak Lensing Measurements

Authors: S. G. Carlsten, Michael A. Strauss, Robert H. Lupton, Joshua E. Meyers, Satoshi Miyazaki

Abstract: We measure and model the wavelength dependence of the PSF in the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) survey. We find that PSF chromaticity is present in that redder stars appear smaller than bluer stars in the $g, r,$ and $i$-bands at the 1-2 per cent level and in the $z$ and $y$-bands at the 0.1-0.2 per cent level. From the color dependence of the PSF, we fit a model between th… ▽ More We measure and model the wavelength dependence of the PSF in the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) survey. We find that PSF chromaticity is present in that redder stars appear smaller than bluer stars in the $g, r,$ and $i$-bands at the 1-2 per cent level and in the $z$ and $y$-bands at the 0.1-0.2 per cent level. From the color dependence of the PSF, we fit a model between the monochromatic PSF trace radius, $R$, and wavelength of the form $R(λ)\propto λ^{b}$. We find values of $b$ between -0.2 and -0.5, depending on the epoch and filter. This is consistent with the expectations of a turbulent atmosphere with an outer scale length of $\sim 10-100$ m, indicating that the atmosphere is dominating the chromaticity. We find evidence in the best seeing data that the optical system and detector also contribute some wavelength dependence. Meyers and Burchat (2015) showed that $b$ must be measured to an accuracy of $\sim 0.02$ not to dominate the systematic error budget of the Large Synoptic Survey Telescope (LSST) weak lensing (WL) survey. Using simple image simulations, we find that $b$ can be inferred with this accuracy in the $r$ and $i$-bands for all positions in the LSST field of view, assuming a stellar density of 1 star arcmin$^{-2}$ and that the optical PSF can be accurately modeled. Therefore, it is possible to correct for most, if not all, of the bias that the wavelength-dependent PSF will introduce into an LSST-like WL survey. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Comments: 14 pages, 10 figures. Submitted to MNRAS. Comments welcome

WFIRST Science Investigation Team "Cosmology with the High Latitude Survey" Annual Report 2017

Authors: Olivier Doré, Christopher Hirata, Yun Wang, David Weinberg, Ivano Baronchelli, Andrew Benson, Peter Capak, Ami Choi, Tim Eifler, Shoubaneh Hemmati, Shirley Ho, Albert Izard, Bhuvnesh Jain, Mike Jarvis, Alina Kiessling, Elisabeth Krause, Elena Massara, Dan Masters, Alex Merson, Hironao Miyatake, Andres Plazas Malagon, Rachel Mandelbaum, Lado Samushia, Chaz Shapiro, Melanie Simet , et al. (24 additional authors not shown)

Abstract: 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 obje… ▽ More 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. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Comments: 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

Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs). IV. Discovery of 41 Quasars and Luminous Galaxies at 5.7 < z < 6.9

Authors: Y. Matsuoka, K. Iwasawa, M. Onoue, N. Kashikawa, M. A. Strauss, C. -H. Lee, M. Imanishi, T. Nagao, M. Akiyama, N. Asami, J. Bosch, H. Furusawa, T. Goto, J. E. Gunn, Y. Harikane, H. Ikeda, T. Izumi, T. Kawaguchi, N. Kato, S. Kikuta, K. Kohno, Y. Komiyama, R. H. Lupton, T. Minezaki, S. Miyazaki , et al. (22 additional authors not shown)

Abstract: We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidate… ▽ More We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes. △ Less

Submitted 20 May, 2018; v1 submitted 5 March, 2018; originally announced March 2018.

Comments: ApJS in press. arXiv admin note: substantial text overlap with arXiv:1704.05854

Journal ref: ApJS 237 (2018) 5

SCARLET: Source separation in multi-band images by Constrained Matrix Factorization

Authors: Peter Melchior, Fred Moolekamp, Maximilian Jerdee, Robert Armstrong, Ai-Lei Sun, James Bosch, Robert Lupton

Abstract: We present the source separation framework SCARLET for multi-band images, which is based on a generalization of the Non-negative Matrix Factorization to alternative and several simultaneous constraints. Our approach describes the observed scene as a mixture of components with compact spatial support and uniform spectra over their support. We present the algorithm to perform the matrix factorizatio… ▽ More We present the source separation framework SCARLET for multi-band images, which is based on a generalization of the Non-negative Matrix Factorization to alternative and several simultaneous constraints. Our approach describes the observed scene as a mixture of components with compact spatial support and uniform spectra over their support. We present the algorithm to perform the matrix factorization and introduce constraints that are useful for optical images of stars and distinct stellar populations in galaxies, in particular symmetry and monotonicity with respect to the source peak position. We also derive the treatment of correlated noise and convolutions with band-dependent point spread functions, rendering our approach applicable to coadded images observed under variable seeing conditions. SCARLET thus yields a PSF-matched photometry measurement with an optimally chosen weight function given by the mean morphology in all available bands. We demonstrate the performance of SCARLET for deblending crowded extragalactic scenes and on an AGN jet -- host galaxy separation problem in deep 5-band imaging from the Hyper Suprime-Cam Stategic Survey Program. Using simulations with prominent crowding we show that SCARLET yields superior results to the HSC-SDSS deblender for the recovery of total fluxes, colors, and morphologies. Due to its non-parametric nature, a conceptual limitation of SCARLET is its sensitivity to undetected sources or multiple stellar population within detected sources, but an iterative strategy that adds components at the location of significant residuals appears promising. The code is implemented in Python with C++ extensions and is available at https://github.com/fred3m/scarlet △ Less

Submitted 22 July, 2018; v1 submitted 27 February, 2018; originally announced February 2018.

Comments: accepted by Astronomy & Computing

First release of high-redshift superluminous supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA). I. Photometric properties

Authors: Takashi J. Moriya, Masaomi Tanaka, Naoki Yasuda, Ji-an Jiang, Chien-Hsiu Lee, Keiichi Maeda, Tomoki Morokuma, Ken'ichi Nomoto, Robert M. Quimby, Nao Suzuki, Ichiro Takahashi, Masayuki Tanaka, Nozomu Tominaga, Masaki Yamaguchi, Stephanie R. Bernard, Jeff Cooke, Chris Curtin, Lluis Galbany, Santiago Gonzalez-Gaitan, Giuliano Pignata, Tyler Pritchard, Yutaka Komiyama, Robert H. Lupton

Abstract: We report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), the transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at the spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with the host-galaxy photometric redshifts of 3.2 (H… ▽ More We report our first discoveries of high-redshift supernovae from the Subaru HIgh-Z sUpernova CAmpaign (SHIZUCA), the transient survey using Subaru/Hyper Suprime-Cam. We report the discovery of three supernovae at the spectroscopically-confirmed redshifts of 2.399 (HSC16adga), 1.965 (HSC17auzg), and 1.851 (HSC17dbpf), and two supernova candidates with the host-galaxy photometric redshifts of 3.2 (HSC16apuo) and 4.2 (HSC17dsid), respectively. In this paper, we present their photometric properties and the spectroscopic properties of the confirmed high-redshift supernovae are presented in the accompanying paper (Curtin et al. 2019). The supernovae with the confirmed redshifts of z ~ 2 have the rest ultraviolet peak magnitudes close to -21 mag and they are likely superluminous supernovae. The discovery of three supernovae at z ~ 2 roughly corresponds to the approximate event rate of ~ 900 +/- 520 Gpc-3 yr-1 with Poisson error, which is consistent with the total superluminous supernova rate estimated by extrapolating the local rate based on the cosmic star-formation history. Adding unconfirmed superluminous supernova candidates would increase the event rate. Our superluminous supernova candidates at the redshifts of around 3 and 4 indicate the approximate superluminous supernova rates of ~ 400 +/- 400 Gpc-3 yr-1 (z ~ 3) and ~ 500 +/- 500 Gpc-3 yr-1 (z ~ 4) with Poisson errors. Our initial results demonstrate the outstanding capability of Hyper Suprime-Cam to discover high-redshift supernovae. △ Less

Submitted 19 February, 2019; v1 submitted 24 January, 2018; originally announced January 2018.

Comments: 27 pages, 15 figures, 6 tables, accepted by The Astrophysical Journal Supplements

Journal ref: The Astrophysical Journal Supplement Series, Volume 241, Issue 2, article id. 16, 19 pp. (2019)

Stellar Stream and Halo Structure in the Andromeda Galaxy From a Subaru/Hyper Suprime-Cam Survey

Authors: Yutaka Komiyama, Masashi Chiba, Mikito Tanaka, Masayuki Tanaka, Takanobu Kirihara, Yohei Miki, Masao Mori, Robert H. Lupton, Puragra Guhathakurta, Jason S. Kalirai, Karoline Gilbert, Evan Kirby, Myun Gyoon Lee, In Sung Jang, Sanjib Sharma, Kohei Hayashi

Abstract: We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the N… ▽ More We present wide and deep photometry of the northwest part of the halo of the Andromeda galaxy (M31) using Hyper Suprime-Cam on the Subaru Telescope. The survey covers 9.2 deg$^{2}$ field in the $g$, $i$, and $NB515$ bands and shows a clear red giant branch (RGB) of M31's halo stars and a pronounced red clump (RC) feature. The spatial distribution of RC stars shows a prominent stream feature, the North Western (NW) Stream, and a diffuse substructure in the south part of our survey field. We estimate the distances based on the RC method and obtain $(m-M)$ = 24.63$\pm 0.191$(random)$\pm0.057$(systematic) and 24.29$\pm 0.211$(random)$\pm0.057$(systematic) mag for the NW stream and diffuse substructure, respectively, implying that the NW Stream is located behind M31, whereas the diffuse substructure is located in front. We also estimate line-of-sight distances along the NW Stream and find that the south part of the stream is $\sim$20 kpc closer to us relative to the north part. The distance to the NW Stream inferred from the isochrone fitting to the color-magnitude diagram favors the RC-based distance, but the TRGB-based distance estimated for $NB515$-selected RGB stars does not agree with it. The surface number density distribution of RC stars across the NW Stream is found to be approximately Gaussian with a FWHM of $\sim$25 arcmin (5.7 kpc), with a slight skew to the south-west side. That along the NW Stream shows a complicated structure including variations in number density and a significant gap in the stream. △ Less

Submitted 11 December, 2017; originally announced December 2017.

Comments: 24 pages, 18 figures, accepted for publication in ApJ

Exploring the brighter fatter effect with the Hyper Suprime-Cam

Authors: William R. Coulton, Robert Armstrong, Kendrick M. Smith, Robert H. Lupton, David N. Spergel

Abstract: The brighter fatter effect has been postulated to arise due to the build up of a transverse electric field, produced as photo-charges accumulate in the pixels' potential wells. We investigate the brighter fatter effect in Hyper Suprime-Cam by examining flat fields and moments of stars. We observe deviations from the expected linear relation in the photon transfer curve, luminosity dependent correl… ▽ More The brighter fatter effect has been postulated to arise due to the build up of a transverse electric field, produced as photo-charges accumulate in the pixels' potential wells. We investigate the brighter fatter effect in Hyper Suprime-Cam by examining flat fields and moments of stars. We observe deviations from the expected linear relation in the photon transfer curve, luminosity dependent correlations between pixels in flat field images and a luminosity dependent point spread function (PSF) in stellar observations. Under the key assumptions of translation invariance and Maxwell's equations in the quasi-static limit, we give a first-principles proof that the effect can be parametrized by a translationally invariant scalar kernel. We describe how this kernel can be estimated from flat fields and discuss how this kernel has been used to remove the brighter fatter distortions in Hyper Suprime-Cam images. We find that our correction restores the expected linear relation in the photon transfer curves and significantly reduces, but does not completely remove, the luminosity dependence of the PSF over a wide range of magnitudes. △ Less

Submitted 16 November, 2017; originally announced November 2017.

Comments: 15 pages, 20 figures, submitted to the Astronomical Journal

Scientific Synergy Between LSST and Euclid

Authors: Jason Rhodes, Robert C. Nichol, Éric Aubourg, Rachel Bean, Dominique Boutigny, Malcolm N. Bremer, Peter Capak, Vincenzo Cardone, Benoît Carry, Christopher J. Conselice, Andrew J. Connolly, Jean-Charles Cuillandre, N. A. Hatch, George Helou, Shoubaneh Hemmati, Hendrik Hildebrandt, Renée Hložek, Lynne Jones, Steven Kahn, Alina Kiessling, Thomas Kitching, Robert Lupton, Rachel Mandelbaum, Katarina Markovic, Phil Marshall , et al. (12 additional authors not shown)

Abstract: Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions.… ▽ More Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating. △ Less

Submitted 29 November, 2017; v1 submitted 23 October, 2017; originally announced October 2017.

Comments: 33 pages, 7 figures, to appear in ApJS, revised with additional references and minor corrections

Illuminating Low-Surface-Brightness Galaxies with the Hyper Suprime-Cam Survey

Authors: Johnny P. Greco, Jenny E. Greene, Michael A. Strauss, Lauren A. MacArthur, Xzavier Flowers, Andy D. Goulding, Song Huang, Ji Hoon Kim, Yutaka Komiyama, Alexie Leauthaud, Lukas Leisman, Robert H. Lupton, Cristóbal Sifón, Shiang-Yu Wang

Abstract: We present a catalog of extended low-surface-brightness galaxies (LSBGs) identified in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first ${\sim}$200 deg$^2$ of the survey, we have uncovered 781 LSBGs, spanning red ($g-i\geq0.64$) and blue ($g-i<0.64$) colors and a wide range of morphologies. Since we focus on extended galaxies ($r_\mathrm{eff}=2.5$-… ▽ More We present a catalog of extended low-surface-brightness galaxies (LSBGs) identified in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first ${\sim}$200 deg$^2$ of the survey, we have uncovered 781 LSBGs, spanning red ($g-i\geq0.64$) and blue ($g-i<0.64$) colors and a wide range of morphologies. Since we focus on extended galaxies ($r_\mathrm{eff}=2.5$-$14^{\prime\prime}$), our sample is likely dominated by low-redshift objects. We define LSBGs to have mean surface brightnesses $\barμ_\mathrm{eff}(g)>24.3$ mag arcsec$^{-2}$, which allows nucleated galaxies into our sample. As a result, the central surface brightness distribution spans a wide range of $μ_0(g)=18$-$27.4$ mag arcsec$^{-2}$, with 50% and 95% of galaxies fainter than 24.3 and 22 mag arcsec$^{-2}$, respectively. Furthermore, the surface brightness distribution is a strong function of color, with the red distribution being much broader and generally fainter than that of the blue LSBGs, and this trend shows a clear correlation with galaxy morphology. Red LSBGs typically have smooth light profiles that are well-characterized by single-component Sérsic functions. In contrast, blue LSBGs tend to have irregular morphologies and show evidence for ongoing star formation. We crossmatch our sample with existing optical, HI, and ultraviolet catalogs to gain insight into the physical nature of the LSBGs. We find that our sample is diverse, ranging from dwarf spheroidals and ultra-diffuse galaxies in nearby groups to gas-rich irregulars to giant LSB spirals, demonstrating the potential of the HSC-SSP to provide a truly unprecedented view of the LSBG population. △ Less

Submitted 9 May, 2018; v1 submitted 13 September, 2017; originally announced September 2017.

Comments: 21 pages, 14 figures, 1 table; accepted for publication in ApJ

Journal ref: ApJ, 857, 104 (2018)

The Hyper Suprime-Cam Software Pipeline

Authors: James Bosch, Robert Armstrong, Steven Bickerton, Hisanori Furusawa, Hiroyuki Ikeda, Michitaro Koike, Robert Lupton, Sogo Mineo, Paul Price, Tadafumi Takata, Masayuki Tanaka, Naoki Yasuda, Yusra AlSayyad, Andrew C. Becker, William Coulton, Jean Coupon, Jose Garmilla, Song Huang, K. Simon Krughoff, Dustin Lang, Alexie Leauthaud, Kian-Tat Lim, Nate B. Lust, Lauren A. MacArthur, Rachel Mandelbaum , et al. (10 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 18 May, 2017; originally announced May 2017.

Comments: 39 pages, 21 figures, 2 tables. Submitted to Publications of the Astronomical Society of Japan

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

Authors: Rachel Mandelbaum, Hironao Miyatake, Takashi Hamana, Masamune Oguri, Melanie Simet, Robert Armstrong, James Bosch, Ryoma Murata, François Lanusse, Alexie Leauthaud, Jean Coupon, Surhud More, Masahiro Takada, Satoshi Miyazaki, Joshua S. Speagle, Masato Shirasaki, Cristóbal Sifón, Song Huang, Atsushi J. Nishizawa, Elinor Medezinski, Yuki Okura, Nobuhiro Okabe, Nicole Czakon, Ryuichi Takahashi, Will Coulton , et al. (6 additional authors not shown)

Abstract: 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σ$ point-source depth of $i\sim 26$. Given conservative galaxy selection crit… ▽ More 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σ$ 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 $ρ_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. △ Less

Submitted 26 December, 2017; v1 submitted 18 May, 2017; originally announced May 2017.

Comments: 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

Characterization and Photometric Performance of the Hyper Suprime-Cam Software Pipeline

Authors: Song Huang, Alexie Leauthaud, Ryoma Murata, James Bosch, Paul Price, Robert Lupton, Rachel Mandelbaum, Claire Lackner, Steve Bickerton, Satoshi Miyazaki, Jean Coupon, Masayuki Tanaka

Abstract: 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 proto… ▽ More 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. △ Less

Submitted 3 May, 2017; originally announced May 2017.

Comments: 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