The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The Lyman-$\alpha$ forest is a highly non-linear field with a lot of information available in the data beyond the power spectrum. The flux probability distribution function (PDF) has been used as a successful probe of small-scale physics. In this paper we argue that measuring coefficients of the Legendre polyonomial expansion of the PDF offers several advantages over measuring the binned values as is commonly done. In particular, $n$-th coefficient can be expressed as a linear combination of the first $n$ moments, allowing these coefficients to be measured in the presence of noise and allowing a clear route for marginalisation over mean flux. Moreover, in the presence of noise, our numerical work shows that a finite number of coefficients are well measured with a very sharp transition into noise dominance. This compresses the available information into a small number of well-measured quantities. We find that the amount of recoverable information is a very non-linear function of spectral noise that strongly favors fewer quasars measured at better signal to noise.

The squeezed-limit bispectrum, which is generated by nonlinear gravitational evolution as well as inflationary physics, measures the correlation of three wavenumbers, in the configuration where one wavenumber is much smaller than the other two. Since the squeezed-limit bispectrum encodes the impact of a large-scale fluctuation on the small-scale power spectrum, it can be understood as how the small-scale power spectrum "responds" to the large-scale fluctuation. Viewed in this way, the squeezed-limit bispectrum can be calculated using the response approach even in the cases which do not submit to perturbative treatment. To illustrate this point, we apply this approach to the cross-correlation between the large-scale quasar density field and small-scale Lyman-$\alpha$ forest flux power spectrum. In particular, using separate universe simulations which implement changes in the large-scale density, velocity gradient, and primordial power spectrum amplitude, we measure how the Lyman-$\alpha$ forest flux power spectrum responds to the local, long-wavelength quasar overdensity, and equivalently their squeezed-limit bispectrum. We perform a Fisher forecast for the ability of future experiments to constrain local non-Gaussianity using the bispectrum of quasars and the Lyman-$\alpha$ forest. Combining with quasar and Lyman-$\alpha$ forest power spectra to constrain the biases, we find that for DESI the expected $1-\sigma$ constraint is ${\rm err}[f_{\rm NL}]\sim60$. Ability for DESI to measure $f_{\rm NL}$ through this channel is limited primarily by the aliasing and instrumental noise of the Lyman-$\alpha$ forest flux power spectrum. The combination of response approach and separate universe simulations provides a novel technique to explore the constraints from the squeezed-limit bispectrum between different observables.

Using a series of carefully constructed numerical experiments based on hydrodynamic cosmological SPH simulations, we attempt to build an intuition for the relevant physics behind the large scale density ($b_\delta$) and velocity gradient ($b_\eta$) biases of the Lyman-$\alpha$ forest. Starting with the fluctuating Gunn-Peterson approximation applied to the smoothed total density field in real-space, and progressing through redshift-space with no thermal broadening, redshift-space with thermal broadening and hydrodynamicaly simulated baryon fields, we investigate how approximations found in the literature fare. We find that Seljak's 2012 analytical formulae for these bias parameters work surprisingly well in the limit of no thermal broadening and linear redshift-space distortions. We also show that his $b_\eta$ formula is exact in the limit of no thermal broadening. Since introduction of thermal broadening significantly affects its value, we speculate that a combination of large-scale measurements of $b_\eta$ and the small scale flux PDF might be a sensitive probe of the thermal state of the IGM. We find that large-scale biases derived from the smoothed total matter field are within 10-20\% to those based on hydrodynamical quantities, in line with other measurements in the literature.

We present the analysis on our new limits of the dark matter (DM) halo consisting of primordial black holes (PBHs) or massive compact halo objects (MACHOs). We present a search of the first two years of publicly available Kepler mission data for potential signatures of gravitational microlensing caused by these objects, as well as an extensive analysis of the astrophysical sources of background error. These include variable stars, flare events, and comets or asteroids which are moving through the Kepler field. We discuss the potential of detecting comets using the Kepler lightcurves, presenting measurements of two known comets and one unidentified object, most likely an asteroid or comet. After removing the background events with statistical cuts, we find no microlensing candidates. We therefore present our Monte Carlo efficiency calculation in order to constrain the PBH DM with masses in the range of 2 x 10^-9 solar masses to 10^-7 solar masses. We find that PBHs in this mass range cannot make up the entirety of the DM, thus closing a full order of magnitude in the allowed mass range for PBH DM.

Primordial Black Holes (PBHs) remain a Dark Matter (DM) candidate of the Standard Model of Particle Physics. Previously, we proposed a new method of constraining the remaining PBH DM mass range using microlensing of stars monitored by NASA's Kepler mission. We improve this analysis using a more accurate treatment of the population of the Kepler source stars, their variability and limb-darkening. We extend the theoretically detectable PBH DM mass range down to $2\times10^{-10} M_\sun$, two orders of magnitude below current limits and one third order of magnitude below our previous estimate. We address how to extract the DM properties such as mass and spatial distribution if PBH microlensing events were detected. We correct an error in a well-known finite-source limb-darkening microlensing formula and also examine the effects of varying the light curve cadence on PBH DM detectability. We also introduce an approximation for estimating the predicted rate of detection per star as a function of the star's properties, thus allowing for selection of source stars in future missions, and extend our analysis to planned surveys, such as WFIRST.

If the Dark Matter consists of primordial black holes (PBHs), we show that gravitational lensing of stars being monitored by NASA's Kepler search for extra-solar planets can cause significant numbers of detectable microlensing events. A search through the roughly 150,000 lightcurves would result in large numbers of detectable events for PBHs in the mass range $5 \ten{-10}\msun$ to $\aten{-4}\msun$. Non-detection of these events would close almost two orders of magnitude of the mass window for PBH dark matter. The microlensing rate is higher than previously noticed due to a combination of the exceptional photometric precision of the Kepler mission and the increase in cross section due to the large angular sizes of the relatively nearby Kepler field stars. We also present a new formalism for calculating optical depth and microlensing rates in the presence of large finite-source effects.