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$.

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

The Etherington distance duality relation, which relates the luminosity distance, the angular diameter distance and the redshift of objects, depends only upon a conservation law for light that traces back directly to the Lorentzian spacetime geometry. We show that this duality relation indeed survives transition to the most general linear electrodynamics without birefringence, which rests on a spacetime geometry constituted by a Lorentzian metric and two scalar fields, a dilaton and an axion. By computing the Poynting vector and optical scalar transport in the geometrical optics limit of this framework, we derive the modification of the light flux in the presence of a dilaton field and present improved constraints on the gradient of the dilaton field from observations of the Cosmic Microwave Background spectrum. Although this flux modification may seem applicable also to fundamental modifications of the Etherington relation, we show that the distance duality relation still holds true. Thus any deviations within this classical theory would imply non-metricities, once astrophysical sources for attenuation, such as dust, are accounted for. Moreover, using the most up-to-date measurements of the luminosity distances to Supernovae of Type Ia, and the inferred angular diameter distances from the baryon acoustic feature measurements, we perform a scale-free and nearly model-independent test of the Etherington distance duality relation between redshifts of $0.38$ and $0.61$. We find consistency with the standard distance duality relation and constrain the optical depth of light between these two redshifts to $\Delta\tau=-0.006\pm0.046$.

We develop a novel method of measuring the lensing distortion profiles of clusters with stacking the scaled amplitudes of background galaxy ellipticities as a function of the scaled centric radius according to the NFW prediction of each cluster, based on the assumption that the different clusters in a sample follow the universal NFW profile. First we demonstrate the feasibility of this method using both the analytical NFW model and simulated halos in high-resolution $N$-body simulations. We then apply, as a proof of concept, this method to the Subaru weak lensing data and the XMM/Chandra X-ray observables for a sample of 50 massive clusters in the redshift range $0.15\le z\le 0.3$, where their halo masses range over an order of magnitude. To estimate the NFW parameters of each cluster, we use the halo mass proxy relation of X-ray observables, based on either the hydrostatic equilibrium or the gas mass, and then infer the halo concentration from the model $c(M)$ relation. We evaluate a performance of the NFW scaling analysis by measuring the scatters of 50 cluster lensing profiles relative to the NFW predictions over a range of radii, $0.14\le R/[h^{-1}{\rm Mpc}]\le 2.8$. We found a 4 - 6$\sigma$ level evidence of the universal NFW profile in 50 clusters, for both the X-ray halo mass proxy relations, although the gas mass appears to be a better proxy of the underlying true mass. By comparing the measurements with the simulations of cluster lensing taking into account the statistical errors of intrinsic galaxy shapes in the Subaru data, we argue that additional halo mass errors or intrinsic scatters of $\sigma_{\ln M_{500c}}\sim 0.2$ - $0.3$ could reconcile a difference between the measurements and the simulations.