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.

We investigate primordial tensor non-Gaussianity in single field inflation, during a phase of non-attractor evolution when the spectrum of primordial tensor modes can be enhanced to a level detectable at interferometer scales. Making use of a tensor duality we introduced in arXiv:1808.10475, we analytically compute the full bispectrum of primordial tensor fluctuations during the non-attractor era. During this epoch the shape of the tensor bispectrum is enhanced in the squeezed limit, its amplitude can be amplified with respect to slow-roll models, and tensor non-Gaussianity can exhibit a scale dependence distinctive of our set-up. We prove that our results do not depend on the frame used for the calculations. Squeezed tensor non-Gaussianity induces a characteristic quadrupolar anisotropy on the power spectrum of the stochastic background of primordial tensor perturbations. As a step to make contact with gravitational wave experiments, we discuss the response function of a ground based Michelson interferometer to a gravitational wave background with such a feature.

We investigate whether Effective Field Theory (EFT) approaches, which have been useful in examining inflation and dark energy, can also be used to establish a systematic approach to inflationary reheating. We consider two methods. First, we extend Weinberg's background EFT to the end of inflation and reheating. We establish when parametric resonance and decay of the inflaton occurs, but also find intrinsic theoretical limitations, which make it difficult to capture some reheating models. This motivates us to next consider Cheung, et. al.'s EFT approach, which instead focuses on perturbations and the symmetry breaking induced by the cosmological background. Adapting the latter approach to reheating implies some new and important differences compared to the EFT of Inflation. In particular, there are new hierarchical scales, and we must account for inflaton oscillations during reheating, which lead to discrete symmetry breaking. Guided by the fundamental symmetries, we construct the EFT of reheating, and as an example of its usefulness we establish a new class of reheating models and the corresponding predictions for gravity wave observations. In this paper we primarily focus on the first stages of preheating. We conclude by discussing challenges for the approach and future directions. This paper builds on ideas first proposed in the note arXiv:1507.06651.

High-scale supersymmetry (SUSY) with a split spectrum has become increasingly interesting given the current experimental results. A SUSY scale above the weak scale could be naturally associated with a heavy unstable gravitino, whose decays populate the dark matter (DM) particles. In the mini-split scenario with gravitino at about the PeV scale and the lightest TeV scale neutralino being (a component of) DM, the requirement that the DM relic abundance resulting from gravitino decays does not overclose the Universe and satisfies the indirect detection constraints demand the reheating temperature to be below 10^9 - 10^{10} GeV. On the other hand, the BICEP2 result prefers a heavy inflaton with mass at around 10^{13} GeV and a reheating temperature at or above 10^9 GeV with some general assumptions. The mild tension could be alleviated if SUSY scale is even higher with the gravitino mass above the PeV scale. Intriguingly, in no-scale supergravity, gravitinos could be very heavy at about 10^{13} GeV, the inflaton mass scale, while gauginos could still be light at the TeV scale.

We examine the evolution of cosmological perturbations in a non-thermal post-inflationary history with a late-time matter domination period prior to BBN. Such a cosmology could arise naturally in the well-motivated moduli scenario in the context of supersymmetry (SUSY) -- in particular in models of Split-SUSY. Sub-horizon dark matter perturbations grow linearly during the matter dominated phase before reheating and can lead to an enhancement in the growth of substructure on small scales, even in the presence of dark matter annihilations. This suggests that a new scale (the horizon size at reheating) could be important for determining the primordial matter power spectrum. However, we find that in many non-thermal models free-streaming effects or kinetic decoupling after reheating can completely erase the enhancement leading to small-scale structures. In particular, in the moduli scenario with wino or higgsino dark matter we find that the dark matter particles produced from moduli decays would thermalize with radiation and kinetically decouple below the reheating temperature. Thus, the growth of dark matter perturbations is not sustained, and the predictions for the matter power spectrum are similar to a standard thermal history. We comment on possible exceptions, but these appear difficult to realize within standard moduli scenarios. We conclude that although enhanced structure does not provide a new probe for investigating the cosmic dark ages within these models, it does suggest that non-thermal histories offer a robust alternative to a strictly thermal post-inflationary history.

Within the Minimal Supersymmetric Standard Model (MSSM), LHC bounds suggest that scalar superpartner masses are far above the electroweak scale. Given a high superpartner mass, nonthermal dark matter is a viable alternative to WIMP dark matter generated via freezeout. In the presence of moduli fields nonthermal dark matter production is associated with a long matter dominated phase, modifying the spectral index and primordial tensor amplitude relative to those in a thermalized primordial universe. Nonthermal dark matter can have a higher self-interaction cross-section than its thermal counterpart, enhancing astrophysical bounds on its annihilation signals. We constrain the contributions to the neutralino mass from the bino, wino and higgsino using existing astrophysical bounds and direct detection experiments for models with nonthermal neutralino dark matter. Using these constraints we quantify the expected change to inflationary observables resulting from the nonthermal phase.