Search

Strongly lensed type Ia supernovae (SNe Ia) provide a unique cosmological probe to address the Hubble tension problem in cosmology. In addition to the sensitivity of the time delays to the value of the Hubble constant, the transient and standard candle nature of SNe Ia also enable valuable joint constraints on the model of the lens and the cosmological parameters. The upcoming Legacy Survey of Space and Time (LSST) with the Vera C. Rubin Observatory is expected to increase the number of observed SNe Ia by an order of magnitude in ten years of its lifetime. However, finding such systems in the LSST data is a challenge. In this work, we revisit the color-magnitude (CM) diagram used previously as a means to identify lensed SNe Ia and extend the work further as follows. We simulate LSST-like photometric data ($rizy$~bands) of lensed SNe Ia and analyze it in the CM parameter space. We find that a subset of lensed SNe Ia are redder compared to unlensed SNe Ia at a given magnitude, both in the rising and falling phases of their light curves and for SNe up to $z=3$. We propose a modified selection criterion based on these new results. We show that the contamination coming from the unlensed core-collapse (CC) SNe is negligible, whereas a small fraction of lensed CC SNe types Ib and Ic may get selected by this criterion as potential lensed SNe. Finally, we demonstrate that our criterion works well on a wide sample of observed unlensed SNe Ia, a handful of known multiply-imaged lensed SNe systems, and a representative sample of observed super-luminous supernovae.

We present the luminosity-halo mass relations of satellite (sLHMRs) galaxies in the SDSS redMaPPer cluster catalogue and the effects of the dense cluster environment on subhalo mass evolution. We use data from the Subaru Hyper Suprime-Cam survey Year-3 catalogue of galaxy shapes to measure the weak lensing signal around these satellites. This signal serves as a probe of the matter distribution around the satellites, thereby providing the masses of their associated subhalos. We bin our satellites based on physical observable quantities such as their luminosity or the host cluster's richness, combined with their cluster-centric radial separations. Our results indicate that although more luminous satellites tend to reside in more massive halos, the sLHMRs depend on the distance of the satellite from the cluster centre. Subhalos near the cluster centre (within $<0.3 h^{-1}Mpc$) are stripped of mass. Consequently, the ratio of subhalo mass to luminosity decreases near the cluster centre. For low luminosity galaxies ($L < 10^{10} h^{-2}L_{\odot}$), the lack of evidence of increasing subhalo masses with luminosity shows the impact of tidal stripping. We also present stellar-to-subhalo mass relations (sSHMRs) for our satellite sample evolving at different cluster-centric separations. Inferred sSHMRs in the outer radial bin appear to match that observed for the field galaxies. We show that the sSHMRs from the mock-redMaPPer run on galaxy catalogues generated by the empirical UniverseMachine galaxy formation model are in good agreement with our observational results. Satellites, when binned based on the host cluster's richness, show very little dependence of the subhalo mass on the richness.

Next-generation ground-based gravitational-wave (GW) detectors are expected to detect millions of binary black hole mergers during their operation period. A small fraction ($\sim 0.1 - 1\%$) of them will be strongly lensed by intervening galaxies and clusters, producing multiple copies of the GW signals. The expected number of lensed events and the distribution of the time delay between lensed images will depend on the mass distribution of the lenses at different redshifts. Warm dark matter or fuzzy dark matter models predict lower abundances of small mass dark matter halos as compared to the standard cold dark matter. This will result in a reduction in the number of strongly lensed GW events, especially at small time delays. Using the number of lensed events and the lensing time delay distribution, we can put a lower bound on the mass of the warm/fuzzy dark matter particle from a catalog of lensed GW events. The expected bounds from GW strong lensing from next-generation detectors are significantly better than the current constraints., Comment: 8 pages, 7 figures

The uncertainties in photometric redshifts and stellar masses from imaging surveys affect galaxy sample selection, their abundance measurements, as well as the measured weak lensing signals. We develop a framework to assess the systematic effects arising from the use of redshifts and stellar masses derived from photometric data, and explore their impact on the inferred galaxy-dark matter connection. We use galaxy catalogues from the UniverseMachine (UM) galaxy formation model to create Pz-mock galaxy samples that approximately follow the redshift errors in the Subaru HSC survey. We focus on galaxy stellar-mass thresholds ranging from $\log\left[M_*/(h^{-2}M_\odot)\right]$ from $8.6$ to $11.2$ in steps of 0.2 dex within two redshift bins $0.30-0.55$ and $0.55-0.80$. A comparison of the Pz-mock samples to true galaxy samples in UM shows a relatively mild sample contamination for thresholds with $\log\left[M_{*,\rm limit}/(h^{-2}M_\odot)\right]<10.6$, while an increasing contamination towards the more massive end. We show how such contamination affects the measured abundance and the lensing signal. A joint HOD modelling of the observables from the Pz-mock compared to the truth in the UM informs the systematic biases on the average halo masses of central galaxies in the HSC survey. Even with a reasonably conservative choice of photo-$z$ errors in Pz-mock, we show that the inferred halo masses deduced from the HSC galaxies for low-mass thresholds will have a systematic bias smaller than 0.05 dex. Beyond $\log\left[M_{*,\rm limit}/(h^{-2}M_\odot)\right]=10.6$, the inferred halo masses show an increasing systematic bias with stellar mass, reaching values of order $0.2$ dex, larger than the statistical error., Comment: 20 pages, 20 figures, Submitted to MNRAS Journal

We present the results of a nation-wide baseline survey, conducted by us, for the status of Astronomy education among secondary school students in India. The survey was administered in 10 different languages to over 2000 students from diverse backgrounds, and it explored multiple facets of their perspectives on astronomy. The topics included students' views on the incorporation of astronomy in curricula, their grasp of fundamental astronomical concepts, access to educational resources, cultural connections to astronomy, and their levels of interest and aspirations in the subject. We find notable deficiencies in students' knowledge of basic astronomical principles, with only a minority demonstrating proficiency in key areas such as celestial sizes, distances, and lunar phases. Furthermore, access to resources such as telescopes and planetariums remain limited across the country. Despite these challenges, a significant majority of students expressed a keen interest in astronomy. We further analyze the data along socioeconomic and gender lines. Particularly striking were the socioeconomic disparities, with students from resource-poor backgrounds often having lower levels of access and proficiency. Some differences were observed between genders, although not very pronounced. The insights gleaned from this study hold valuable implications for the development of a more robust astronomy curriculum and the design of effective teacher training programs in the future., Comment: 15 pages, 19 figures

We present cosmological constraints using the abundance of weak-lensing shear-selected galaxy clusters in the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The clusters are selected on the mass maps constructed using the three-year (Y3) weak-lensing data with an area of $\approx500~$deg$^2$, resulting in a sample size of $129$ clusters with high signal-to-noise ratios $\nu$ of $\nu\geq4.7$. Owing to the deep, wide-field, and uniform imaging of the HSC survey, this is by far the largest sample of shear-selected clusters, in which the selection solely depends on gravity and is free from any assumptions about the dynamical state. Informed by the optical counterparts, the shear-selected clusters span a redshift range of $z\lesssim0.7$ with a median of $z\approx0.3$. The lensing sources are securely selected at $z\gtrsim0.7$ with a median of $z\approx1.3$, leading to nearly zero cluster member contamination. We carefully account for (1) the bias in the photometric redshift of sources, (2) the bias and scatter in the weak-lensing mass using a simulation-based calibration, and (3) the measurement uncertainty that is directly estimated on the mass maps using an injection-based method developed in a companion paper (Chen et al. submitted). In a blind analysis, the fully marginalized posteriors of the cosmological parameters are obtained as $\Omega_{\mathrm{m}} = 0.50^{+0.28}_{-0.24}$, $\sigma_8 = 0.685^{+0.161}_{-0.088}$, $\hat{S}_{8}\equiv\sigma_8\left(\Omega_{\mathrm{m}}/0.3\right)^{0.25} = 0.835^{+0.041}_{-0.044}$, and $\sigma_8\left(\Omega_{\mathrm{m}}/0.3\right)^{0.5} = 0.993^{+0.084}_{-0.126}$ in a flat $\Lambda$CDM model. We compare our cosmological constraints with other studies, including those based on cluster abundances, galaxy-galaxy lensing and clustering, and Cosmic Microwave Background observed by $Planck$, and find good agreement at levels of $\lesssim2\sigma$. [abridged], Comment: Accepted for publication in the Open Journal of Astrophysics. Main results are Figures 11 and 12. The cluster catalog and the cosmological chain are accessible at https://github.com/inonchiu/hsc_shear_selected_clusters

We present the first step towards deriving cosmological constraints through the abundances of galaxy clusters selected in a $510\,\mathrm{deg}^2$ weak-lensing aperture mass map, constructed with the Year-Three shear catalog from the Hyper Suprime-Cam Subaru Strategic Program. We adopt a conservative source galaxy selection to construct a sample of $129$ weak-lensing peaks with a signal-to-noise ratio above $4.7$. We use semi-analytical injection simulations to derive the selection function and the mass--observable relation of our sample. These results take into account complicated uncertainties associated with weak-lensing measurements, such as the non-uniform survey depth and the complex survey geometry, projection effects from uncorrelated large-scale structures, and the intrinsic alignment of source galaxies. We also propose a novel modeling framework to make parts of the mass--observable relation insensitive to assumed cosmological parameters. Such a framework not only offers a great computational advantage to cosmological studies, but can also benefit future astrophysical studies using shear-selected clusters. Our results are an important step towards utilizing these cluster samples that are constructed nearly independent of any baryonic assumptions in upcoming deep-and-wide lensing surveys from the Vera Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope., Comment: 17 pages, 9 figures, submitted to OJA. The cluster catalog of the cosmological sample (signal-to-noise ratio above $4.7$) is accessible at https://github.com/inonchiu/hsc_shear_selected_clusters

We present a detailed exposition of a statistical method for estimating cosmological parameters from the observation of a large number of strongly lensed binary-black-hole (BBH) mergers observable by next (third) generation (XG) gravitational-wave (GW) detectors. This method, first presented in Jana (2023 Phys. Rev. Lett. 130 261401), compares the observed number of strongly lensed GW events and their time delay distribution (between lensed images) with observed events to infer cosmological parameters. We show that the precision of the estimation of the cosmological parameters does not have a strong dependance on the assumed BBH redshift distribution model. Using the large number of unlensed mergers, XG detectors are expected to measure the BBH redshift distribution with sufficient precision for the cosmological inference. However, a biased inference of the BBH redshift distribution will bias the estimation of cosmological parameters. An incorrect model for the distribution of lens properties can also lead to a biased cosmological inference. However, Bayesian model selection can assist in selecting the right model from a set of available parametric models for the lens distribution. We also present a way to incorporate the effect of contamination in the data due to the limited efficiency of lensing identification methods, so that it will not bias the cosmological inference., Comment: 17 pages, 5 figures

Galaxies that cause the strong gravitational lensing of background galaxies provide us crucial information about the distribution of matter around them. Traditional modelling methods that analyse such strong lenses are both time and resource consuming, require sophisticated lensing codes and modelling expertise. To study the large lens population expected from imaging surveys such as LSST, we need fast and automated analysis methods. In this work, we build and train a simple convolutional neural network with an aim to rapidly predict model parameters of gravitational lenses. We focus on the most important lens mass model parameters, namely, the Einstein radius, the axis ratio and the position angle of the major axis of the mass distribution. The network is trained on a variety of simulated data with an increasing degree of realism and shows satisfactory performance on simulated test data. The trained network is then applied to the real sample of galaxy-scale candidate lenses from the Subaru HSC, a precursor survey to LSST. Unlike the simulated lenses, we do not have the ground truth for the real lenses. Therefore, we have compared our predictions with those from YattaLens, a lens modelling pipeline. Additionally, we also compare the parameter predictions for 10 HSC lenses that were also studied by other conventional modelling methods. These comparisons show a fair quantitative agreement on the Einstein radius, although the axis ratio and the position angle from the network as well as the individual modelling methods, seem to have systematic uncertainties beyond the quoted errors., Comment: 10 pages, 9 figures

The baryonic feedback effect is considered as a possible solution to the so-called $S_8$ tension indicated in cosmic shear cosmology. The baryonic effect is more significant on smaller scales, and affects the cosmic shear two-point correlation functions (2PCFs) with different scale- and redshift-dependencies from those of the cosmological parameters. In this paper, we use the Hyper Suprime-Cam Year 3 (HSC-Y3) data to measure the cosmic shear 2PCFs ($\xi_{\pm}$) down to 0.28 arcminutes, taking full advantage of the high number density of source galaxies in the deep HSC data, to explore a possible signature of the baryonic effect. While the published HSC analysis used the cosmic shear 2PCFs on angular scales, which are sensitive to the matter power spectrum at $k\lesssim 1~h{\rm Mpc}^{-1}$, the smaller scale HSC cosmic shear signal allows us to probe the signature of matter power spectrum up to $k\simeq 20~h{\rm Mpc}^{-1}$. Using the accurate emulator of the nonlinear matter power spectrum, DarkEmulator2, we show that the dark matter-only model can provide an acceptable fit to the HSC-Y3 2PCFs down to the smallest scales. In other words, we do not find any clear signature of the baryonic effects or do not find a systematic shift in the $S_8$ value with the inclusion of the smaller-scale information as would be expected if the baryonic effect is significant. Alternatively, we use a flexible 6-parameter model of the baryonic effects, which can lead to both enhancement and suppression in the matter power spectrum compared to the dark matter-only model, to perform the parameter inference of the HSC-Y3 2PCFs. We find that the small-scale HSC data allow only a fractional suppression of up to 5 percent in the matter power spectrum at $k\sim 1~h{\rm Mpc}^{-1}$, which is not sufficient to reconcile the $S_8$ tension., Comment: 30 pages, 16 figures

Gravitational waves (GW) from the inspiral of binary compact objects offers a one-step measurement of the luminosity distance to the event, which is essential for the measurement of the Hubble constant, $H_0$, that characterizes the expansion rate of the Universe. However, unlike binary neutron stars, the inspiral of binary black holes is not expected to be accompanied by electromagnetic radiation and a subsequent determination of its redshift. Consequently, independent redshift measurements of such GW events are necessary to measure $H_0$. In this study, we present a novel Bayesian approach to infer $H_0$ from the cross-correlation between galaxies with known redshifts and individual binary black hole merger events. We demonstrate the efficacy of our method with $250$ simulated GW events distributed within $1$ Gpc in colored Gaussian noise of Advanced LIGO and Advanced Virgo detectors operating at O4 sensitivity. We show that such measurements can constrain the Hubble constant with a precision of $\lesssim 15 \%$ ($90\%$ highest density interval). We highlight the potential improvements that need to be accounted for in further studies before the method can be applied to real data., Comment: 12 pages, 5 figures, 1 table

Orphan galaxies that have lost a large fraction of the dark matter subhaloes have often been invoked in semi-analytical as well as empirical models of galaxy formation. We run a mock cluster finder that mimics the optical cluster finding technique of the redMaPPer algorithm on a catalogue of galaxies with quenched star formation from one such empirical model, the UniverseMachine, and obtain the prevalence of orphan galaxies in these clusters as a function of their cluster-centric distance. We compare the fraction of orphan galaxies with the upper limits derived based on our prior observations of the weak lensing signals around satellite galaxies from SDSS redMaPPer clusters. Although the orphan fraction from the UniverseMachine is marginally consistent with the upper limits in the innermost regions of galaxy clusters spanning [0.1, 0.3] $h^{-1}$ Mpc, we observe that the orphan fractions substantially violate the upper limits in the outer regions of galaxy clusters beyond 0.3 $h^{-1}$ Mpc. We discuss the reasons, plausible improvements to the model and how observations can be used to constrain such models further., Comment: 6 pages, 5 figures

We present cosmology results obtained from a blind joint analysis of the abundance, projected clustering, and weak lensing of galaxy clusters measured from the Sloan Digital Sky Survey (SDSS) redMaPPer cluster catalog and the Hyper-Suprime Cam (HSC) Year3 shape catalog. We present a full-forward model for the cluster observables, which includes empirical modeling for the anisotropic boosts on the lensing and clustering signals of optical clusters. We validate our analysis via mock cluster catalogs which include observational systematics, such as the projection effect and the effect of baryonic feedback, and find that our analysis can robustly constrain cosmological parameters in an unbiased manner without any informative priors on our model parameters. The joint analysis of our observables in the context of the flat $\Lambda$CDM model results in cosmological constraints for $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm m} / 0.3}=0.816^{+0.041}_{-0.039}$. Our result is consistent with the $S_8$ inference from other cosmic microwave background- and large scale structure-based cosmology analyses, including the result from the \emph{Planck} 2018 primary CMB analysis., Comment: v1: 22 pages, 15 figures, Comments welcome

Superclusters are the largest massive structures in the cosmic web on tens to hundreds of megaparsecs (Mpc) scales. They are the largest assembly of galaxy clusters in the Universe. Apart from a few detailed studies of such structures, their evolutionary mechanism is still an open question. In order to address and answer the relevant questions, a statistically significant, large catalog of superclusters covering a wide range of redshifts and sky areas is essential. Here, we present a large catalog of 662 superclusters identified using a modified $\textit{ Friends of Friends}$ algorithm applied on the WHL (Wen-Han-Liu) cluster catalog within a redshift range of $0.05 \le z \le 0.42$. We name the most massive supercluster at $z \sim 0.25$ as $\textit{Einasto Supercluster}$. We find that the median mass of superclusters is $\sim 5.8 \times 10^{15}$ M$_{\odot}$ and median size $\sim 65$ Mpc. We find that the supercluster environment slightly affects the growth of clusters. We compare the properties of the observed superclusters with the mock superclusters extracted from the Horizon Run 4 cosmological simulation. The properties of superclusters in mocks and observations are in broad agreement. We find that the density contrast of a supercluster is correlated with its maximum extent with a power law index, $\alpha \sim -2$. The phase-space distribution of mock superclusters shows that, on average, $\sim 90\%$ part of a supercluster has a gravitational influence on its constituents. We also show mock halos' average number density and peculiar velocity profiles in and around the superclusters., Comment: 23 pages, 16 figures, and 2 tables. Accepted for publication in ApJ

In the absence of numerous gravitational-wave detections with confirmed electromagnetic counterparts, the "dark siren" method has emerged as a leading technique of gravitational-wave cosmology. The method allows redshift information of such events to be inferred statistically from a catalogue of potential host galaxies. Due to selection effects, dark siren analyses necessarily depend on the mass distribution of compact objects and the evolution of their merger rate with redshift. Informative priors on these quantities will impact the inferred posterior constraints on the Hubble constant ($H_0$). It is thus crucial to vary these unknown distributions during an $H_0$ inference. This was not possible in earlier analyses due to the high computational cost, restricting them to either excluding galaxy catalogue information, or fixing the gravitational-wave population mass distribution and risking introducing bias to the $H_0$ measurement. This paper introduces a significantly enhanced version of the Python package GWCOSMO, which allows joint estimation of cosmological and compact binary population parameters. This thereby ensures the analysis is now robust to a major source of potential bias. The gravitational-wave events from the Third Gravitational-Wave Transient Catalogue are reanalysed with the GLADE+ galaxy catalogue, and an updated, more reliable measurement of $H_0=69^{+12}_{-7}$ km s$^{-1}$ Mpc$^{-1}$ is found (maximum a posteriori probability and 68% highest density interval). This improved method will enable cosmological analyses with future gravitational-wave detections to make full use of the information available (both from galaxy catalogues and the compact binary population itself), leading to promising new independent bounds on the Hubble constant., Comment: 31 pages, 11 figures, accepted for publication in JCAP

Primordial black holes, which could have formed during the early Universe through overdensities in primordial density fluctuations during inflation, are potential candidates for dark matter. We explore the use of lensing parallax of Gamma ray bursts (GRBs), which results in different fluxes being observed from two different vantage points, in order to probe the abundance of primordial black holes in the unexplored window within the mass range $[10^{-15}-10^{-11}]M_\odot$. We derive the optical depth for the lensing of GRBs with a distribution of source properties and realistic detector sensitivities. We comment on the ability of the proposed Indian twin satellite mission Daksha in its low earth orbit to conduct this experiment. If the two Daksha satellites observe 10000 GRBs simultaneously and the entirety of dark matter is made up of $[10^{-15}-10^{-12}]M_\odot$ black holes, Daksha will detect non-zero lensing events with a probability ranging from 80 to 50 per cent at the bin edges, respectively. Non-detections will not conclusively rule out primordial black holes as dark matter in this mass range. However, we show that meaningful constraints can be obtained in such a case if the two satellites are separated by at least the Earth-Moon distance., Comment: 10 pages, 7 figures, updated according to referee's comments, accepted for publication in MNRAS

We infer the connection between the stellar mass of galaxies from the Subaru Hyper Suprime-Cam (HSC) survey, and their dark matter halo masses and its evolution in two bins of redshifts between $[0.3, 0.8]$. We use the measurements of the weak lensing signal of galaxies using background sources from the Year 1 shape catalog from the HSC survey. We bin galaxies in stellar mass with varying thresholds ranging from $8.6 \leq \log [ M_*/(h^{-2} {M_\odot})] \leq 11.2$ and use stringent cuts in the selection of source galaxies to measure the weak lensing signal. We model these measurements of the weak lensing signal together with the abundance of galaxies in the halo occupation distribution framework. We obtain constraints on the halo occupation parameters of central galaxies $M_{\rm min}$ and $\sigma_{\log M}$, which correspond to the halo mass at which central galaxies for each threshold sample reach half occupancy, and its scatter, respectively, along with parameters that describe the occupation of the satellite galaxies. The measurements of abundance and weak lensing individually constrain different degeneracy directions in the $M_{\rm min}$ and $\sigma_{\log M}$ plane, thus breaking the degeneracy in these parameters. We demonstrate that the weak lensing measurements are best able to constrain the average central halo masses, $\langle M_{\rm cen} \rangle$. We compare our measurements to those obtained using the abundance and clustering of these galaxies as well as the subhalo abundance matching measurements and demonstrate qualitative agreement. We find that the galaxy-dark matter connection does not vary significantly between redshift bins we explore in this study. Uncertainties in the photometric redshift of the lens galaxies imply that more efforts are required to understand the true underlying stellar mass-halo mass relation of galaxies and its evolution over cosmic epoch.

Galaxy clusters, being the most massive objects in the Universe, exhibit the strongest alignment with the large-scale structure. However, mis-identification of members due to projection effects from the large scale structure can occur. We studied the impact of projection effects on the measurement of the intrinsic alignment of galaxy clusters, using galaxy cluster mock catalogs. Our findings showed that projection effects result in a decrease of the large scale intrinsic alignment signal of the cluster and produce a bump at $r_p\sim 1h^{-1}/Mpc$, most likely due to interlopers and missed member galaxies. This decrease in signal explains the observed similar alignment strength between bright central galaxies and clusters in the SDSS redMaPPer cluster catalog. The projection effect and cluster intrinsic alignment signal are coupled, with clusters having lower fractions of missing members or having higher fraction of interlopers exhibiting higher alignment signals in their projected shapes. We aim to use these findings to determine the impact of projection effects on galaxy cluster cosmology in future studies., Comment: 15 pages, 14 figures, Accepted by MNRAS

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 $\Delta \Sigma(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, $\xi_{\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 $\Lambda$CDM cosmology. We use conservative scale cuts, $R>12$ and $8~h^{-1}$Mpc, for $\Delta\Sigma$ 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 $\Lambda$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., Comment: 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/

We present cosmology results from a blinded joint analysis of cosmic shear, $\xi_{\pm}(\vartheta)$, galaxy-galaxy weak lensing, $\Delta\!\Sigma(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 $\Delta\!\Sigma$ in three spectroscopic redshift bins spanning the range $0.15

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

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 $\Lambda$CDM model, we find $S_8 \equiv \sigma_8 (\Omega_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\sigma$-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\sigma$, indicating that our result is robust to such systematics., Comment: 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/

We investigate the stellar mass-dependence of the galaxy size-dark matter halo radius relation for low redshift galaxies using weak gravitational lensing measurements. Our sample consists of $\sim$38,000 galaxies more massive than $10^{8}{\rm M}_{\odot}h^{-2}$ and within $z<0.3$ drawn from the overlap of GAMA survey DR4 and HSC-SSP PDR2. We divide our sample into a number of stellar mass bins and measure stacked weak lensing signals. We model the signals using a conditional stellar mass function to infer the stellar mass-halo mass relation. We fit a single S\'ersic model to HSC $i$-band images of our galaxies and obtain their three-dimensional half-light radii. We use these measurements to construct a median galaxy size-mass relation. We then combine these relations to infer the galaxy size-halo radius relation. We confirm that this relation appears linear given the statistical errors, i.e. the ratio of galaxy size to halo radius remains constant over two orders of magnitudes in stellar mass above $\sim 10^{9} {\rm M}_{\odot}h^{-2}$. Extrapolating the stellar mass-halo mass relation below this limit, we see an indication of a decreasing galaxy size-halo radius ratio with the decline in stellar mass. At stellar mass $\sim 10^{8} {\rm M}_{\odot}h^{-2}$ the ratio becomes 30% smaller than its value in linear regime. The possible existence of a such trend in dwarf galaxy sectors calls for either modification in models employing a constant fraction of halo angular momentum transferred to explain sizes of dwarfs or else points towards our lack of knowledge about dark matter haloes of low-mass galaxies., Comment: 20 pages, 12 figures, accepted in MNRAS

We present the splashback radius measurements around the SRG/eROSITA eFEDS X-ray selected galaxy clusters by cross-correlating them with HSC S19A photometric galaxies. The X-ray selection is expected to be less affected by systematics related to projection that affects optical cluster finder algorithms. We use a nearly volume-limited sample of 109 galaxy clusters selected in 0.5-2.0 keV band having luminosity $L_X > 10^{43.5}\,{\rm erg s^{-1} h^{-2}}$ within the redshift $z<0.75$ and obtain measurements of the projected cross-correlation with a signal-to-noise of $17.43$. We model our measurements to infer a three-dimensional profile and find that the steepest slope is sharper than $-3$ and associate the location with the splashback radius. We infer the value of the 3D splashback radius $r_{\rm sp} = 1.45^{+0.30}_{-0.26}\,{\rm h^{-1} Mpc}$. We also measure the weak lensing signal of the galaxy clusters and obtain halo mass $\log[M_{\rm 200m}/{\rm h^{-1}M_\odot}] = 14.52 \pm 0.06$ using the HSC-S16A shape catalogue data at the median redshift $z=0.46$ of our cluster sample. We compare our $r_{\rm sp}$ values with the spherical overdensity boundary $r_{\rm 200m} = 1.75 \pm 0.08\,{\rm h^{-1} Mpc}$ based on the halo mass which is consistent within $1.2\sigma$ with the $\Lambda$CDM predictions. Our constraints on the splashback radius, although broad, are the best measurements thus far obtained for an X-ray selected galaxy cluster sample., Comment: 15 pages, 10 figures

Which galaxies in the general population turn into active galactic nuclei (AGNs) is a keystone of galaxy formation and evolution. Thanks to SRG/eROSITA's contiguous 140 square degree pilot survey field, we constructed a large, complete, and unbiased soft X-ray flux-limited ($F_X>6.5\times 10^{-15}$ erg s$^{-1}$ cm$^{-2}$) AGN sample at low redshift, $0.05

Cosmological weak lensing measurements rely on a precise measurement of the shear two-point correlation function (2PCF) along with a deep understanding of systematics that affect it. In this work, we demonstrate a general framework for detecting and modeling the impact of PSF systematics on the cosmic shear 2PCF, and mitigating its impact on cosmological analysis. Our framework can describe leakage and modeling error from all spin-2 quantities contributed by the PSF second and higher moments, rather than just the second moments, using the cross-correlations between galaxy shapes and PSF moments. We interpret null tests using the HSC Year 3 (Y3) catalogs with this formalism, and find that leakage from the spin-2 combination of PSF fourth moments is the leading contributor to additive shear systematics, with total contamination that is an order of magnitude higher than that contributed by PSF second moments alone. We conducted a mock cosmic shear analysis for HSC Y3, and find that, if uncorrected, PSF systematics can bias the cosmological parameters $\Omega_m$ and $S_8$ by $\sim$0.3$\sigma$. The traditional second moment-based model can only correct for a 0.1$\sigma$ bias, leaving the contamination largely uncorrected. We conclude it is necessary to model both PSF second and fourth moment contamination for HSC Y3 cosmic shear analysis. We also reanalyze the HSC Y1 cosmic shear analysis with our updated systematics model, and identify a 0.07$\sigma$ bias on $\Omega_m$ when using the more restricted second moment model from the original analysis. We demonstrate how to self-consistently use the method in both real space and Fourier space, assess shear systematics in tomographic bins, and test for PSF model overfitting., Comment: 33 pages, 28 figures, match the refereed version in MNRAS

We present posterior sample redshift distributions for the Hyper Suprime-Cam Subaru Strategic Program Weak Lensing three-year (HSC Y3) analysis. Using the galaxies' photometry and spatial cross-correlations, we conduct a combined Bayesian Hierarchical Inference of the sample redshift distributions. The spatial cross-correlations are derived using a subsample of Luminous Red Galaxies (LRGs) with accurate redshift information available up to a photometric redshift of $z < 1.2$. We derive the photometry-based constraints using a combination of two empirical techniques calibrated on spectroscopic- and multiband photometric data that covers a spatial subset of the shear catalog. The limited spatial coverage induces a cosmic variance error budget that we include in the inference. Our cross-correlation analysis models the photometric redshift error of the LRGs to correct for systematic biases and statistical uncertainties. We demonstrate consistency between the sample redshift distributions derived using the spatial cross-correlations, the photometry, and the posterior of the combined analysis. Based on this assessment, we recommend conservative priors for sample redshift distributions of tomographic bins used in the three-year cosmological Weak Lensing analyses., Comment: 23 pages, 11 figures, 1 table, updated to match the accepted version in the MNRAS

We present the science case for the proposed Daksha high energy transients mission. Daksha will comprise of two satellites covering the entire sky from 1~keV to $>1$~MeV. The primary objectives of the mission are to discover and characterize electromagnetic counterparts to gravitational wave source; and to study Gamma Ray Bursts (GRBs). Daksha is a versatile all-sky monitor that can address a wide variety of science cases. With its broadband spectral response, high sensitivity, and continuous all-sky coverage, it will discover fainter and rarer sources than any other existing or proposed mission. Daksha can make key strides in GRB research with polarization studies, prompt soft spectroscopy, and fine time-resolved spectral studies. Daksha will provide continuous monitoring of X-ray pulsars. It will detect magnetar outbursts and high energy counterparts to Fast Radio Bursts. Using Earth occultation to measure source fluxes, the two satellites together will obtain daily flux measurements of bright hard X-ray sources including active galactic nuclei, X-ray binaries, and slow transients like Novae. Correlation studies between the two satellites can be used to probe primordial black holes through lensing. Daksha will have a set of detectors continuously pointing towards the Sun, providing excellent hard X-ray monitoring data. Closer to home, the high sensitivity and time resolution of Daksha can be leveraged for the characterization of Terrestrial Gamma-ray Flashes., Comment: 19 pages, 7 figures. Accepted in Experimental Astronomy. More details about the mission at https://www.dakshasat.in/

We perform a rigorous cosmology analysis on simulated type Ia supernovae (SN~Ia) and evaluate the improvement from including photometric host-galaxy redshifts compared to using only the "zspec" subset with spectroscopic redshifts from the host or SN. We use the Deep Drilling Fields (~50 deg^2) from the Photometric LSST Astronomical Time-Series Classification Challenge (PLaSTiCC), in combination with a low-z sample based on Data Challenge2 (DC2). The analysis includes light curve fitting to standardize the SN brightness, a high-statistics simulation to obtain a bias-corrected Hubble diagram, a statistical+systematics covariance matrix including calibration and photo-z uncertainties, and cosmology fitting with a prior from the cosmic microwave background. Compared to using the zspec subset, including events with SN+host photo-z results in i) more precise distances for z>0.5, ii) a Hubble diagram that extends 0.3 further in redshift, and iii) a 50 % increase in the Dark Energy Task Force figure of merit (FoM) based on the w0-wa CDM model. Analyzing 25 simulated data samples, the average bias on w0 and wa is consistent with zero. The host photo-z systematic of 0.01 reduces FoM by only 2 % because i) most z<0.5 events are in the zspec subset, ii) the combined SN+host photo-z has X 2 smaller bias, and iii) the anti-correlation between fitted redshift and color self corrects distance errors. To prepare for analysing real data, the next SNIa-cosmology analysis with photo-z's should include non SN-Ia contamination and host galaxy mis-associations., Comment: 13 pages, 9 figures, 7 tables. Submitted to ApJ

The clustering of dark matter halos depends on the assembly history of halos at fixed halo mass; a phenomenon referred to as \textit{halo assembly bias}. Halo assembly bias is readily observed in cosmological simulations of dark matter. However, it is difficult to detect it in observations. The identification of galaxy or cluster properties correlated with the formation time of the halo at fixed halo mass and the ability to select galaxy clusters free from projection effects are the two most significant hurdles in the observational detection of halo assembly bias. The latter, in particular, can cause a misleading detection of halo assembly bias by boosting the amplitude of lensing and clustering on large scales. This study uses twelve different properties of central galaxies of SDSS redMaPPer clusters derived from spectroscopy to divide the clusters into sub-samples. We test the dependence of the clustering amplitude on these properties at fixed richness. We first infer halo mass and bias using weak lensing signals around the clusters using shapes of galaxies from the SDSS survey. We validate the bias difference between the two subsamples using cluster-galaxy cross-correlations. This methodology allows us to decouple the contamination due to the projection effects from the halo assembly bias signals. We do not find any significant evidence of a difference in the clustering amplitudes correlated with any of our explored properties. Our results indicate that central galaxy properties may not correlate significantly with the halo assembly histories at fixed richness., Comment: 21 pages, 7 figures

Weak gravitational lensing directly probes the matter distribution surrounding satellite galaxies in galaxy clusters. We measure the weak lensing signal induced on the shapes of background galaxies around SDSS redMaPPer cluster satellite galaxies, which have their central galaxies assigned with a probability $P_{\rm cen}>0.95$ in the redshift range, $0.1\leq z\leq 0.33$. We use the galaxy shapes from the Subaru Hyper Suprime-Cam (HSC) survey for this purpose. We bin satellite galaxies by their distance from the cluster centre and compare it to the signal around a control sample of galaxies which do not reside in clusters but have similar colours and magnitudes. We explore the effect of environmental processes on the dark matter mass around satellites. We see hints of a difference in the mass of the subhalo of the satellite compared to the halo masses of galaxies in our control sample, especially in the innermost cluster-centric radial bin ($0.1

We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that the lensing mass profile is consistent with that of massive galaxies ($M_{*}\sim 10.8~M_\odot$). Surprisingly, the lensing signal remains unchanged when the AGN sample is split into low and high stellar mass hosts. Specifically, we find that the excess surface density (ESD) of AGNs, residing in galaxies with high stellar masses, significantly differs from that of the control sample. We further fit a halo occupation distribution model to the data to infer the posterior distribution of parameters including the average halo mass. We find that the characteristic halo mass of the full AGN population lies near the knee ($\rm log(M_h/h^{-1}M_{\odot})=12.0$) of the stellar-to-halo mass relation (SHMR). Illustrative of the results given above, the halo masses of AGNs residing in host galaxies with high stellar masses (i.e., above the knee of the SHMR) falls below the calibrated SHMR while the halo mass of the low stellar mass sample is more consistent with the established SHMR. These results indicate that massive halos with higher clustering bias tend to suppress AGN activity, probably due to the lack of available gas., Comment: 17 pages,10 figures, accepted for publication in the Astrophysical Journal

We perform a cosmic shear analysis of HSC survey first-year data (HSC-Y1) using Complete Orthogonal Sets of E/B-Integrals (COSEBIs) to derive cosmological constraints. We compute E/B-mode COSEBIs from cosmic shear two-point correlation functions measured on an angular range of $4\arcmin<\theta<180\arcmin$. We perform the standard Bayesian likelihood analysis for cosmological inference from the measured E-mode COSEBIs, including contributions from intrinsic alignments of galaxies as well as systematic effects from point spread function model errors, shear calibration uncertainties, and source redshift distribution errors. We adopt a covariance matrix derived from realistic mock catalogs constructed from full-sky gravitational lensing simulations that fully take account of the survey geometry and measurement noise. For a flat $\Lambda$ cold dark matter model, we find $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}=0.809_{-0.026}^{+0.036}$. 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. We also find that the measured B-mode COSEBIs are consistent with zero. We examine, using mock HSC-Y1 data, the consistency of our $S_8$ constraints with those derived from the other cosmic shear two-point statistics, the power spectrum analysis by Hikage et al (2019) and the two-point correlation function analysis by Hamana et al (2020), which adopt the same HSC-Y1 shape catalog, and find that all the $S_8$ constraints are consistent with each other, although expected correlations between derived $S_8$ constraints are weak., Comment: 23 pages, 19 figures, accepted for PASJ. arXiv admin note: substantial text overlap with arXiv:1906.06041

We utilize the Hyper Suprime-Cam (HSC) Wide Survey to explore the properties of galaxies located in the voids identified from the Baryon Oscillation Spectroscopic Survey (BOSS) up to z~0.7. The HSC reaches i~25, allowing us to characterize the void galaxies down to 10$^{9.2}$ solar mass. We find that the revised void galaxy densities, when including faint galaxies in voids defined by bright galaxies, are still underdense compared to the mean density from the entire field. In addition, we classify galaxies into star-forming, quiescent, and green valley populations, and find that void galaxies tend to have slightly higher fractions of star-forming galaxies under the mass and redshift control, although the significance of this result is only moderate (2$\sigma$). However, when we focus on the star-forming population, the distribution of the specific star formation rate (sSFR) of void galaxies shows little difference from that of the control galaxies. Similarly, the median sSFR of star-forming void galaxies is also in good agreement with that of the star-forming control galaxies. Moreover, the effective green valley fraction of void galaxies, defined as the number of green valley galaxies over the number of nonquiescent galaxies, is comparable to that of the control ones, supporting the suggestion that void and control galaxies evolve under similar physical processes and quenching frequencies. Our results thus favor a scenario of the galaxy assembly bias., Comment: 17 pages, 10 figures, 3 tables, ApJ accepted

Although it is generally accepted that massive galaxies form in a two-phased fashion, beginning with a rapid mass buildup through intense starburst activities, followed by primarily dry mergers that mainly deposit stellar mass at outskirts, the late time stellar mass growth of brightest cluster galaxies (BCGs), the most massive galaxies in the universe, is still not well understood. Several independent measurements have indicated a slower mass growth rate than predictions from theoretical models. We attempt to resolve the discrepancy by measuring the frequency of BCGs with multiple-cores, which serve as a proxy of the merger rates in the central region and facilitate a more direct comparison with theoretical predictions. Using 79 BCGs at $z=0.06-0.15$ with integral field spectroscopic (IFS) data from the Mapping Nearby Galaxies at APO (MaNGA) project, we obtain a multiple-core fraction of $0.11 \pm 0.04$ at $z\approx 0.1$ within a 18 kpc radius from the center, which is comparable to the value of $0.08 \pm 0.04$ derived from mock observations of 218 simulated BCGs from the cosmological hydrodynamical simulation IllustrisTNG. We find that most of cores that appear close to the BCGs from imaging data turn out to be physically associated systems. Anchoring on the similarity in the multiple-core frequency between the MaNGA and IllustrisTNG, we discuss the mass growth rate of BCGs over the past 4.5 Gyr., Comment: 38 pages, 31 figures, 5 tables; accepted by ApJ

We present a novel analysis for cluster cosmology that fully forward models the abundances, weak lensing, and the clustering of galaxy clusters. Our analysis notably includes an empirical model for the anisotropic boosts impacting the lensing and clustering signals of optical clusters. These boosts arise from a preferential selection of clusters surrounded by anisotropic large scale structure, a consequence of the limited discrimination between line-of-sight interlopers and true cluster members offered by photometric surveys. We validate our analysis via a blind cosmology challenge on mocks, and find that we can obtain tight and unbiased cosmological constraints without informative priors or external calibrations on any of our model parameters. We then apply our analysis on the SDSS redMaPPer clusters, and find results favoring low $\Omega_\mathrm{m}$ and high $\sigma_8$, combining to yield the lensing strength constraint $S_8 = 0.718_{-0.021}^{+0.024}$. We investigate potential drivers behind these results through a series of post-unblinding tests, noting that our results are consistent with existing cluster cosmology constraints but clearly inconsistent with other CMB/LSS based cosmology results. From these tests, we find hints that a suppression in the cluster lensing signal may be driving our results., Comment: v1: 19 pages, 14 figures. Comments welcome. v2: Added clarifications in text and reflected a minor bug fix for Section 3; version submitted to MNRAS. v3: version accepted and published at MNRAS

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

As the number of detected gravitational wave sources increase with increasing sensitivity of the gravitational wave observatories, observing strongly lensed pair of events will become a real possibility. Lensed GW events will have very accurately measured time delays, as well as magnification ratios. If the corresponding lens system can be identified electromagnetically and the redshifts of the lens and the host galaxy can be measured, such events can be used to constrain important astrophysical parameters of the lens system. As most lensing events will have image separations that will be significantly smaller than the GW event localization uncertainties, it is important that we develop diagnostics that will aid in the robust identification of such lensed events. We define a new statistic based on the joint probability of lensing observables that can be used to discriminate lensed pair of events from the unlensed ones. To this end, we carry out simulations of lensed GW events to infer the distribution of the relative time delays, relative magnifications sub-divided by the type of lensed images. We compare this distribution to a similar one obtained for random unlensed event pairs. Our improved statistic can be used to improve the existing ranking approach adopted by the search pipelines in order to down-select event pairs for joint parameter estimates. The distributions we obtain can further be used to define more informative priors in joint parameter estimation analyses for candidate lensed events., Comment: 8 pages, 6 figures

We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($\Delta\!\Sigma$) and projected galaxy clustering ($w_{\rm p}$) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range $0.15

Low mass star formation inside massive clusters is crucial to understand the effect of cluster environment on processes like circumstellar disk evolution, planet and brown dwarf formation. The young massive association of Cygnus OB2, with a strong feedback from massive stars, is an ideal target to study the effect of extreme environmental conditions on its extensive low-mass population. We aim to perform deep multi-wavelength studies to understand the role of stellar feedback on the IMF, brown dwarf fraction and circumstellar disk properties in the region. We introduce here, the deepest and widest optical photometry of 1.5$^\circ$ diameter region centred at Cygnus OB2 in r$_{2}$, i$_{2}$, z and Y-filters using Subaru Hyper Suprime-Cam (HSC). This work presents the data reduction, source catalog generation, data quality checks and preliminary results about the pre-main sequence sources. We obtain 713,529 sources in total, with detection down to $\sim$ 28 mag, 27 mag, 25.5 mag and 24.5 mag in r$_{2}$, i$_{2}$, z and Y-band respectively, which is $\sim$ 3 - 5 mag deeper than the existing Pan-STARRS and GTC/OSIRIS photometry. We confirm the presence of a distinct pre-main sequence branch by statistical field subtraction of the central 18$^\prime$ region. We find the median age of the region as $\sim$ 5 $\pm$ 2 Myrs with an average disk fraction of $\sim$ 9$\%$. At this age, combined with A$_V$ $\sim$ 6 - 8 mag, we detect sources down to a mass range $\sim$ 0.01 - 0.17 M$_\odot$. The deep HSC catalog will serve as the groundwork for further studies on this prominent active young cluster., Comment: Accepted for publication in Monthly Notices of Royal Astronomical Society (MNRAS), 20 pages, 21 figures, 4 Tables, typos corrected

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\sigma$) 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/., Comment: 25 pages, 19 figures, submitted to PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

Rotation periods of 53 small (diameters $2 < D < 40$ km) Jupiter Trojans (JTs) were derived using the high-cadence light curves obtained by the FOSSIL phase I survey, a Subaru/Hyper Suprime-Cam intensive program. These are the first reported periods measured for JTs with $D < 10$ km. We found a lower limit of the rotation period near 4 hr, instead of the previously published result of 5 hr (Ryan et al. 2017; Szabo et al. 2017, 2020) found for larger JTs. Assuming a rubble-pile structure for JTs, a bulk density of 0.9 gcm$^{-3}$ is required to withstand this spin rate limit, consistent with the value $0.8-1.0$ gcm$^{-3}$ (Marchis et al. 2006; Mueller et al. 2010; Buie et al. 2015; Berthier et al. 2020) derived from the binary JT system, (617) Patroclus-Menoetius system., Comment: Accepted by PSJ on July 9th, 2021. 15 pages, 7 figures, and 3 tables

We utilize the galaxy shape catalogue from the first-year data release of the Subaru Hyper Suprime-cam Survey (HSC) to study the dark matter content of galaxy groups in the Universe using weak lensing. We use galaxy groups from the Galaxy Mass and Assembly galaxy survey in approximately $100$ sq. degrees of the sky that overlap with the HSC survey as lenses. We restrict our analysis to the $1587$ groups with at least five members. We divide these groups into six bins each of group luminosity and group member velocity dispersion and measure the lensing signal with a signal-to-noise ratio of $55$ and $51$ for these two different selections, respectively. We use a Bayesian halo model framework to infer the halo mass distribution of our groups binned in the two different observable properties and constrain the power-law scaling relation, and the scatter between mean halo masses and the two group observable properties. We obtain a 5 percent constraint on the amplitude of the scaling relation between halo mass and group luminosity with $\avg{M} = (0.81\pm 0.04)\times10^{14}\hinvMsun$ for $L_{\rm grp}=10^{11.5}\hinvsqLsun$, and a power-law index of $\alpha=1.01\pm 0.07$. We constrain the amplitude of the scaling relation between halo mass and velocity dispersion to be $\avg{M}=(0.93\pm 0.05)\times10^{14}\hinvMsun$ for $\sigma=500 \kms$ and a power-law index to be $\alpha=1.52\pm0.10$. However, these scaling relations are sensitive to the exact cuts applied to the number of group members. Comparisons with similar scaling relations from the literature show that our results are consistent and have significantly reduced errors., Comment: 17 pages, 11 figures, Accepted for publication in MNRAS

We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|\delta m|<9\times 10^{-3}$. The bias residual $\delta m$ shows no dependence on redshift in the range $01$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements., Comment: 39 pages, 31 figures, 5 tables; PASJ (in press)

The abundance, distribution and inner structure of satellites of galaxy clusters can be sensitive probes of the properties of dark matter. We run 30 cosmological zoom-in simulations with self-interacting dark matter (SIDM), with a velocity-dependent cross-section, to study the properties of subhalos within cluster-mass hosts. We find that the abundance of subhalos that survive in the SIDM simulations are suppressed relative to their cold dark matter (CDM) counterparts. Once the population of disrupted subhalos -- which may host orphan galaxies -- are taken into account, satellite galaxy populations in CDM and SIDM models can be reconciled. However, even in this case, the inner structure of subhalos are significantly different in the two dark matter models. We study the feasibility of using the weak lensing signal from the subhalo density profiles to distinguish between the cold and self-interacting dark matter while accounting for the potential contribution of orphan galaxies. We find that the effects of self-interactions on the density profile of subhalos can appear degenerate with subhalo disruption in CDM, when orphans are accounted for. With current error bars from the Subaru Hyper Suprime-Cam Strategic Program, we find that subhalos in the outskirts of clusters (where disruption is less prevalent) can be used to constrain dark matter physics. In the future, the Vera C. Rubin Observatory Legacy Survey of Space and Time will give precise measurements of the weak lensing profile and can be used to constrain $\sigma_T/m$ at the $\sim 1$ cm$^2$ g$^{-1}$ level at $v\sim 2000$ km s$^{-1}$., Comment: 19 pages, 12 figures, to be submitted to ApJ