Search

The three-point correlation function (3PCF) of a weak lensing shear field contains information that is complementary to that in the two-point correlation function (2PCF), which can help improve the cosmological parameters and calibrate astrophysical and observational systematics parameters. However, the application of the 3PCF to observed data has been limited due to the computational challenges of calculating theoretical predictions of the 3PCF from a bispectrum model. In this paper, we present a new method to compute the shear 3PCF efficiently and accurately. We employ the multipole expansion of the bispectrum to compute the shear 3PCF, and show that the method is substantially more efficient than direct numerical integration. We found that the multipole-based method can compute the shear 3PCF with 5% accuracy in 10 (40) seconds for the single (four) source redshift bin setup. The multipole-based method can be also used to compute the third-order aperture mass statistics quickly and accurately, accounting for the bin-averaging effect on the shear 3PCF. Our method provides a fast and robust tool for probing the underlying cosmological model with third-order statistics of weak lensing shear., Comment: 20 pages, 9 figures, 2 tables, code available at https://github.com/git-sunao/fastnc

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

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

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

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/

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

The extended source effect on microlensing magnification is non-negligible and must be taken into account for in an analysis of microlensing. However, the evaluation of the extended source magnification is numerically expensive because it includes the two-dimensional integral over source profile. Various studies have developed methods to reduce this integral down to the one-dimensional-integral or integral-free form, which adopt some approximations or depend on the exact form of the source profile, e.g. disk, linear/quadratic limb-darkening profile. In this paper, we develop a new method to evaluate the extended source magnification based on fast Fourier transformation (FFT), which does not adopt any approximations and is applicable to any source profiles. Our implementation of the FFT based method enables the fast evaluation of the extended source magnification as fast as $\sim1$ msec (CPU time on a laptop) and guarantees an accuracy better than 0.3%. The FFT based method can be used for the template fitting to a huge data set of light curves from the existing and upcoming surveys., Comment: 9 pages, 6 figures; Preparing for submission to ApJ; Code available at https://github.com/git-sunao/fft-extended-source

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

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

Dark matter in the form of axions is expected to form miniclusters, and their dense regions can harbor compact axion stars. Such axion stars could be discovered by microlensing events. In particular, some candidate events reported by Subaru HSC and OGLE can be explained simultaneously if the axion stars with masses of the order of the Earth mass make up about $\sim27^{+7}_{013}$ percent of dark matter. For QCD axions, this corresponds to the axion mass in the range $10^{-9}-10^{-6}$ eV, which is consistent with the experimental constraints, as well as the cosmological anthropic window of parameters., Comment: 8 pages, 2 figures; latest version for the publication in Physics Letters B

We present validation tests of emulator-based halo model method for cosmological parameter inference, assuming hypothetical measurements of the projected correlation function of galaxies, $w_{\rm p}(R)$, and the galaxy-galaxy weak lensing, $\Delta\!\Sigma(R)$, from the spectroscopic SDSS galaxies and the Hyper Suprime-Cam Year1 (HSC-Y1) galaxies. To do this, we use \textsc{Dark Emulator} developed in Nishimichi et al. based on an ensemble of $N$-body simulations, which is an emulation package enabling a fast, accurate computation of halo clustering quantities for flat-geometry $w$CDM cosmologies. Adopting the halo occupation distribution, the emulator allows us to obtain model predictions of $\Delta\!\Sigma$ and $w_{\rm p}$ for the SDSS-like galaxies at a few CPU seconds for an input set of parameters. We present performance and validation of the method by carrying out Markov Chain Monte Carlo analyses of the mock signals measured from a variety of mock catalogs that mimic the SDSS and HSC-Y1 galaxies. We show that the halo model method can recover the underlying true cosmological parameters to within the 68\% credible interval, except for the mocks including the assembly bias effect (although we consider the unrealistically-large amplitude of assembly bias effect). Even for the assembly bias mock, we demonstrate that the cosmological parameters can be recovered {\it if} the analysis is restricted to scales $R\gtrsim 10~h^{-1}{\rm Mpc}$. We also show that, by using a single population of source galaxies to infer the relative strengths of $\Delta\!\Sigma$ for multiple lens samples at different redshifts, the joint probes method allows for self-calibration of photometric redshift errors and multiplicative shear bias. Thus we conclude that the emulator-based halo model method can be safely applied to the HSC-Y1 dataset, achieving a precision of $\sigma(S_8)\simeq 0.04$., Comment: 39 pages, 22 figures and 5 tables

Primordial black holes (PBHs) can constitute the predominant fraction of dark matter (DM) if PBHs reside in the currently unconstrained "sublunar" mass range. PBHs originating from scalar perturbations generated during inflation can naturally appear with a broad spectrum in a class of models. The resulting stochastic gravitational wave (GW) background generated from such PBH production can account for the recently reported North American Nanohertz Observatory for Gravitational Waves (NANOGrav) pulsar timing array data signal, and will be testable in future GW observations by interferometer-type experiments such as Laser Interferometer Space Antenna (LISA). We show that the broad mass function of such PBH DM is already being probed by Subaru Hyper Suprime-Cam (HSC) microlensing data and is consistent with a detected candidate event. Upcoming observations of HSC will be able to provide an independent definitive test of the stochastic GW signals originating from such PBH DM production scenarios., Comment: 6 pages, 1 figure; matches published version

We assess the performance of a perturbation theory inspired method for inferring cosmological parameters from the joint measurements of galaxy-galaxy weak lensing ($\Delta\Sigma$) and the projected galaxy clustering ($w_{\rm p}$). To do this, we use a wide variety of mock galaxy catalogs constructed based on a large set of $N$-body simulations that mimic the Subaru HSC-Y1 and SDSS galaxies, and apply the method to the mock signals to address whether to recover the underlying true cosmological parameters in the mocks. We find that, as long as the appropriate scale cuts, $12$ and $8~h^{-1}{\rm Mpc}$ for $\Delta\Sigma$ and $w_{\rm p}$ respectively, are adopted, a "minimal-bias" model using the linear bias parameter $b_1$ alone and the nonlinear matter power spectrum can recover the true cosmological parameters (here focused on $\Omega_{\rm m}$ and $\sigma_8$) to within the 68% credible interval, for all the mocks we study including one in which an assembly bias effect is implemented. This is as expected if physical processes inherent in galaxy formation/evolution are confined to local, small scales below the scale cut, and thus implies that real-space observables have an advantage in filtering out the impact of small-scale nonlinear effects in parameter estimation, compared to their Fourier-space counterparts. In addition, we find that a theoretical template including the higher-order bias contributions such as nonlinear bias parameter $(b_2)$ does not improve the cosmological constraints, but rather leads to a larger parameter bias compared to the baseline $b_1$-method., Comment: Submitted to PRD, A revised version in response to the referee's comments

Primordial black holes (PBHs) are a viable candidate for dark matter if the PBH masses are in the currently unconstrained "sublunar" mass range. We revisit the possibility that PBHs were produced by nucleation of false vacuum bubbles during inflation. We show that this scenario can produce a population of PBHs that simultaneously accounts for all dark matter, explains the candidate event in Subaru Hyper Suprime-Cam (HSC) data, and contains both heavy black holes as observed by LIGO and very heavy seeds of supermassive black holes. We demonstrate with numerical studies that future observations of HSC, as well as other optical surveys, such as LSST, will be able to provide a definitive test for this generic PBH formation mechanism if it is the dominant source of dark matter., Comment: 7 pages, 3 figures; v2: published version

Microlensing of stars, e.g. in the Galactic bulge and Andromeda galaxy (M31), is among the most robust, powerful method to constrain primordial black holes (PBHs) that are a viable candidate of dark matter. If PBHs are in the mass range $M_{\rm PBH} \lower.5ex\hbox{$\; \buildrel < \over \sim \;$} 10^{-10}M_\odot$, its Schwarzschild radius ($r_{\rm Sch}$) becomes comparable with or shorter than optical wavelength ($\lambda)$ used in a microlensing search, and in this regime the wave optics effect on microlensing needs to be taken into account. For a lensing PBH with mass satisfying $r_{\rm Sch}\sim \lambda$, it causes a characteristic oscillatory feature in the microlensing light curve, and it will give a smoking gun evidence of PBH if detected, because any astrophysical object cannot have such a tiny Schwarzschild radius. Even in a statistical study, e.g. constraining the abundance of PBHs from a systematic search of microlensing events for a sample of many source stars, the wave effect needs to be taken into account. We examine the impact of wave effect on the PBH constraints obtained from the $r$-band (6210\AA) monitoring observation of M31 stars in Niikura et al. (2019), and find that a finite source size effect is dominant over the wave effect for PBHs in the mass range $M_{\rm PBH}\simeq[10^{-11},10^{-10}]M_\odot$. We also discuss that, if a denser-cadence (10~sec), $g$-band monitoring observation for a sample of white dwarfs over a year timescale is available, it would allow one to explore the wave optics effect on microlensing light curve, if it occurs, or improve the PBH constraints in $M_{\rm PBH}\lower.5ex\hbox{$\; \buildrel < \over \sim \;$} 10^{-11}M_\odot$ even from a null detection., Comment: Accepted for publication in MNRAS; 10 pages, 8 figures. Typos are corrected. Fig.4 is newly added

Primordial black holes (PBHs) have long been suggested as a viable candidate for the elusive dark matter (DM). The abundance of such PBHs has been constrained using a number of astrophysical observations, except for a hitherto unexplored mass window of $M_{\rm PBH}=[10^{-14},10^{-9}]M_\odot$. Here we carry out a dense-cadence (2~min sampling rate), 7 hour-long observation of the Andromeda galaxy (M31) with the Subaru Hyper Suprime-Cam to search for microlensing of stars in M31 by PBHs lying in the halo regions of the Milky Way (MW) and M31. Given our simultaneous monitoring of tens of millions of stars in M31, if such light PBHs make up a significant fraction of DM, we expect to find many microlensing events for the PBH DM scenario. However, we identify only a single candidate event, which translates into the most stringent upper bounds on the abundance of PBHs in the mass range $M_{\rm PBH}\simeq [10^{-11}, 10^{-6}]M_\odot$., Comment: 43 pages, 25 figures, 2 tables. Made significant revision of the microlensing event rate calculation taking into account both effects of finite source size and wave optics