Using $N$-body simulations for the first time, we show that the anisotropic primordial non-Gaussianity (PNG) causes a scale-dependent modification, given by $1/k^2$ at small $k$ limit, in the three-dimensional power spectra of halo shapes (intrinsic alignments), whilst the conventional power spectrum of halo number density field remains unaffected. We discuss that wide-area imaging and spectroscopic surveys observing the same region of the sky allow us to constrain the quadrupole PNG coefficient $f_{\rm NL}^{s=2}$ at a precision comparable with that of the cosmic microwave background.

An optical cluster finder inevitably suffers from projection effects, where it misidentifies a superposition of galaxies in multiple halos along the line-of-sight as a single cluster. Using mock cluster catalogs built from cosmological N-body simulations, we quantify the impact of these projection effects with a particular focus on the observables of interest for cluster cosmology, namely the cluster lensing and the cluster clustering signals. We find that "observed" clusters, i.e. clusters identified by our cluster finder algorithm, exhibit lensing and clustering signals that deviate from expectations based on a statistically isotropic halo model -- while both signals agree with halo model expectations on small scales, they show unexpected boosts on large scales, by up to a factor of 1.2 or 1.4 respectively. We identify the origin of these boosts as the inherent selection bias of optical cluster finders for clusters embedded within filaments aligned with the line-of-sight, and show that a minority ($\sim 30\%$) of such clusters within the entire sample is responsible for this observed boost. We discuss the implications of our results on previous studies of optical cluster, as well as prospects for identifying and mitigating projection effects in future cluster cosmology analyses.

We use a suite of $N$-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the $\Lambda$CDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the $E$-mode field of halo shapes with the matter/halo distribution, $P_{\delta E}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those of the auto-power spectrum of the $E$ mode, $P^{(\ell)}_{EE}(k)$, based on the $E$/$B$-mode decomposition. The IA power spectra have non-vanishing amplitudes over the linear to nonlinear scales, and the large-scale amplitudes at $k\lesssim 0.1~h~{\rm Mpc}^{-1}$ are related to the matter power spectrum via a constant coefficient ($A_{\rm IA}$), similar to the linear bias parameter. We find that the cross- and auto-power spectra $P_{\delta E}$ and $P_{EE}$ at nonlinear scales, $k\gtrsim 0.1~h~{\rm Mpc}^{-1}$, show different $k$-dependences relative to the matter power spectrum, suggesting a violation of the nonlinear alignment model commonly used to model contaminations of cosmic shear signals. The IA power spectra exhibit baryon acoustic oscillations, and vary with halo samples of different masses, redshifts and cosmological parameters ($\Omega_{\rm m}, S_8$). The cumulative signal-to-noise ratio for the IA power spectra is about 60\% of that for the halo density power spectrum, where the super-sample covariance is found to give a significant contribution to the total covariance. Our results demonstrate that the IA power spectra of galaxy shapes, measured from imaging and spectroscopic surveys for an overlapping area of the sky, can be powerful tools to probe the underlying matter power spectrum, the primordial curvature perturbations, and cosmological parameters, in addition to the standard galaxy density power spectrum.

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.

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