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 detection of a gamma-ray burst (GRB) in the solar neighborhood would have very important implications for GRB phenomenology. The leading theories for cosmological GRBs would not be able to explain such events. The final bursts of evaporating Primordial Black Holes (PBHs), however, would be a natural explanation for local GRBs. We present a novel technique that can constrain the distance to gamma-ray bursts using detections from widely separated, non-imaging spacecraft. This method can determine the actual distance to the burst if it is local. We applied this method to constrain distances to a sample of 36 short duration GRBs detected by the Interplanetary Network (IPN) that show observational properties that are expected from PBH evaporations. These bursts have minimum possible distances in the 10^13-10^18 cm (7-10^5 AU) range, consistent with the expected PBH energetics and with a possible origin in the solar neighborhood, although none of the bursts can be unambiguously demonstrated to be local. Assuming these bursts are real PBH events, we estimate lower limits on the PBH burst evaporation rate in the solar neighborhood.

Many early universe theories predict the creation of Primordial Black Holes (PBHs). PBHs could have masses ranging from the Planck mass to 10^5 solar masses or higher depending on the size of the universe at formation. A Black Hole (BH) has a Hawking temperature which is inversely proportional to its mass. Hence a sufficiently small BH will quasi-thermally radiate particles at an ever-increasing rate as emission lowers its mass and raises its temperature. The final moments of this evaporation phase should be explosive and its description is dependent on the particle physics model. In this work we investigate the final few seconds of BH evaporation, using the Standard Model and incorporating the most recent Large Hadron Collider (LHC) results, and provide a new parameterization for the instantaneous emission spectrum. We calculate for the first time energy-dependent PBH burst light curves in the GeV/TeV energy range. Moreover, we explore PBH burst search methods and potential observational PBH burst signatures. We have found a unique signature in the PBH burst light curves that may be detectable by GeV/TeV gamma-ray observatories such as the High Altitude Water Cerenkov (HAWC) observatory. The implications of beyond the Standard Model theories on the PBH burst observational characteristics are also discussed, including potential sensitivity of the instantaneous photon detection rate to a squark threshold in the 5 -10 TeV range.

Many early universe theories predict the creation of Primordial Black Holes (PBHs). The PBHs could have masses ranging from the Planck mass to 10^5 solar masses or higher depending on the formation scenario. Hawking showed that any Black Hole (BH) has a temperature which is inversely proportional to its mass. Hence a sufficiently small BH will thermodynamically radiate particles at an ever-increasing rate, continually decreasing its mass and raising its temperature. The final moments of this evaporation phase should be explosive. In this work, we investigate the final few seconds of the BH burst using the Standard Model of particle physics and calculate the energy dependent burst time profiles in the GeV/TeV range. We use the HAWC (High Altitude Water Cherenkov) observatory as a case study and calculate PBH burst light curves which would be observed by HAWC.

Primordial Black Holes (PBHs) are of interest in many cosmological contexts. PBHs lighter than about 1012 kg are predicted to be directly detectable by their Hawking radiation. This radiation should produce both a diffuse extragalactic gamma-ray background from the cosmologically-averaged distribution of PBHs and gamma-ray burst signals from individual light black holes. The Fermi, Milagro, Veritas, HESS and HAWC observatories, in combination with new burst recognition methodologies, offer the greatest sensitivity for the detection of such black holes or placing limits on their existence.

Primordial Black Holes (PBHs) are gravitationally collapsed objects that may have been created by density fluctuations in the early universe and could have arbitrarily small masses down to the Planck scale. Hawking showed that due to quantum effects, a black hole has a temperature inversely proportional to its mass and will emit all species of fundamental particles thermally. PBHs with initial masses of ~5.0 x 10^14 g should be expiring in the present epoch with bursts of high-energy particles, including gamma radiation in the GeV - TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma-ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is well suited to perform a search for PBH bursts. Based on a search on the Milagro data, we report new PBH burst rate density upper limits over a range of PBH observation times. In addition, we report the sensitivity of the Milagro successor, the High Altitude Water Cherenkov (HAWC) observatory, to PBH evaporation events.

Report of the CF6 Working Group at Snowmass 2013. Topics addressed include ultra-high energy cosmic rays, neutrinos, gamma rays, baryogenesis, and experiments probing the fundamental nature of spacetime.

Primordial Black Holes (PBHs) created early in the universe are dark matter candidates. One method of detecting these PBHs is through their Hawking radiation. PBHs created with an initial mass of 5.0 x 10^14 g should be evaporating today with bursts of high-energy particles, including gamma radiation in the GeV - TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is ideally suited for the direct search of PBH bursts. Based on a search in Milagro data, we report PBH upper limits according to the standard model.

Several authors have claimed that the observable Hawking emission from a microscopic black hole is significantly modified by the formation of a photosphere or chromosphere around the black hole due to QED or QCD interactions between the emitted particles. Analyzing these models we identify a number of physical and geometrical effects which invalidate them. In all cases, we find that the observational signatures of a cosmic or Galactic background of black holes or an individual black hole remain essentially those of the standard Hawking model, with little change to the detection probability.

The theoretical maximum time variation in the electronic charge permitted by the Generalized Second Law of Thermodynamics applied to black holes radiating and accreting in the cosmic microwave background matches the measured cosmological variation in the fine structure constant claimed by Webb et al.. Such black holes cannot respond adiabatically to a varying fine structure constant.

Primordial Black Holes (PBHs), which may have been created in the early Universe, are predicted to be detectable by their Hawking radiation. The Fermi Gamma-ray Space Telescope observatory offers increased sensitivity to the gamma-ray bursts produced by PBHs with an initial mass of $\sim 5\times 10^{14}$ g expiring today. PBHs are candidate progenitors of unidentified Gamma-Ray Bursts (GRBs) that lack X-ray afterglow. We propose spectral lag, which is the temporal delay between the high and low energy pulses, as an efficient method to identify PBH evaporation events with the Fermi Large Area Telescope (LAT).

Primordial Black Holes (PBHs), which may have been created in the early Universe, are predicted to be detectable by their Hawking radiation. PBHs with an initial mass of 5.0 * 10^14 g should be expiring today with a burst of high energy particles. Evaporating PBHs in the solar neighborhood are candidate Gamma-Ray Bursts (GRBs) progenitors. We propose spectral lag, which is the temporal delay between the high energy photon pulse and the low energy photon pulse, as a possible method to detect PBH evaporation events with the Fermi Gamma-ray Space Telescope Observatory.

We discuss the Flambaum commentary [arXiv:0808.2518] on the application of the Generalized Second Law of Thermodynamics to derive limits on the variation of fundamental constants, for example the fine structure constant.

We discuss a variety of bremsstrahlung processes associated with charged particles emitted by evaporating black holes. We show that such particles produce a negligible number of bremsstrahlung photons from their scattering off each other, though at low frequencies inner bremsstrahlung photons dominate over the direct Hawking emission of photons. This analysis and the further analysis of the accompanying paper invalidate Heckler's claim that sufficiently hot evaporating black holes form QED photospheres.

Several authors, most notably Heckler, have claimed that the observable Hawking emission from a microscopic black hole is significantly modified by the formation of a photosphere around the black hole due to QED or QCD interactions between the emitted particles. In this paper we analyze these claims and identify a number of physical and geometrical effects which invalidate these scenarios. We point out two key problems. First, the interacting particles must be causally connected to interact, and this condition is satisfied by only a small fraction of the emitted particles close to the black hole. Second, a scattered particle requires a distance ~ E/m_e^2 for completing each bremsstrahlung interaction, with the consequence that it is improbable for there to be more than one complete bremsstrahlung interaction per particle near the black hole. These two effects have not been included in previous analyses. We conclude that the emitted particles do not interact sufficiently to form a QED photosphere. Similar arguments apply in the QCD case and prevent a QCD photosphere (chromosphere) from developing when the black hole temperature is much greater than Lambda_QCD, the threshold for QCD particle emission. Additional QCD phenomenological arguments rule out the development of a chromosphere around black hole temperatures of order Lambda_QCD. In all cases, the observational signatures of a cosmic or Galactic halo background of primordial black holes or an individual black hole remain essentially those of the standard Hawking model, with little change to the detection probability. We also consider the possibility, as proposed by Belyanin et al. and D. Cline et al., that plasma interactions between the emitted particles form a photosphere, and we conclude that this scenario too is not supported.

Here we apply the Generalized Second Law of Thermodynamics and derive upper limits on the variation in the fundamental constants. The maximum variation in the electronic charge permitted for black holes accreting and emitting in the present cosmic microwave background corresponds to a variation in the fine structure constant of ((Delta alpha) / alpha) ~ 2 x 10^-23 per second. This value matches the variation measured by Webb et al. using absorption lines in the spectra of distant quasars and suggests the variation mechanism may be a coupling between the electron and the cosmic photon background.

Here we apply the Generalized Second Law of Thermodynamics (GSL) to black holes accreting and emitting in the present Universe and derive upper limits on the variation in the gravitational constant G. The limits depend on how the gravitational mass M varies with G. Parameterizing M goes as G^n, if n > -1/2 (including n = 0), the GSL applied to the full range of black holes theoretically allowed in the present Universe does not constrain an increase in G but any decrease must be less than about |(1/G) dG/dt| = 10^-52 per second. If n < -1/2, the GSL does not constrain a decrease in G but any increase must be less than about |(1/G) dG/dt| = 10^-52 per second. At earlier redshifts, these constraints weaken as z^3. If n = -1/2, the GSL does not constrain a decrease but any increase must be less than about |(1/G) dG/dt| = (1/t). If the mass range is restricted to those black holes which have been astronomically observed, the present constraints on n > -1/2 and n < -1/2 are only weakened by a factor of about 10^8 with the tightest constraints coming from stellar mass black holes and the n = -1/2 bound does not change. The stellar mass black hole limits should constrain the variation of G in Standard Model physics and all extension models which approximate classical physics on astronomical scales.

Here it is shown that the variation in the fine structure constant measured by Webb et al. matches the theoretically derived value for the maximum variation in the electronic charge permitted by the Generalized Second Law of Thermodynamics for black holes accreting and emitting in the present cosmic microwave background. It is postulated that the constants of nature, independently or dependently, vary at the maximal rate allowed by the Generalized Second Law of Thermodynamics.

Decaying topological defects, in particular cosmic strings, can produce a significant flux of high energy neutrinos, photons and cosmic rays. According to the prevailing understanding of cosmic string dynamics in an expanding Universe, the network of long strings loses its energy first into string loops, which in turn give off most of their energy as gravitational radiation. However, it has been suggested by Vincent et al. (VHS) that particle emission may be the dominant energy loss channel for the long string network. In this case, the predicted flux of high energy particles would be much larger. Here we calculate the predicted flux of high energy gamma rays, neutrinos and cosmic ray antiprotons and protons as a function of the scale of symmetry breaking at which the strings are produced and as a function of the initial energy m_J of the particle jets which result from the string decay. Assuming the validity of the VHS scenario, we find that due to the interactions with the cosmic radiation backgrounds all fluxes but the neutrino flux are suppressed at the highest energies. This indicates that the observed events above the GZK cutoff can only be accounted for in this scenario if the primary particle is a neutrino and the scale of the symmetry breaking is somewhat less than the GUT scale. The domain of parameter space corresponding to GUT-scale symmetry breaking is excluded also by the current observations below the GZK cutoff. A new aspect of this work is the calculation of the spectrum of the tau neutrinos directly produced in the decay of the X-particles. This significantly increases the tau neutrino signal at high energies in all cosmic string scenarios.

Previous work on the neutrino spectra from high energy sources has not included the tau neutrinos directly produced by the decays in the source. Here we consider the tau neutrino component and discuss how its inclusion modifies the expected neutrino spectra. We discuss implications for interpreting any observed tau neutrino component in TeV - UHE events as evidence of nu_mu --> nu_tau oscillations.

Topological defects, particularly cosmic strings, can provide a mechanism to produce particles with energies of the order 10^21 eV and higher. Here, we report on order of magnitude calculations of fluxes from a cosmic string network which evolves according to a new scenario according to which the main channel for energy loss is the particle production rather than gravitational radiation. We compare the predicted fluxes for protons (anti-protons) and neutrinos (anti-neutrinos) with observations of extremely high energy cosmic rays.

We calculate the flux of $\gamma$-rays emitted from primordial black holes (PBHs) which are formed by a ``blue'' power-law spectrum of density fluctuations in the early universe. Gamma-ray emission from such PBHs may contribute significantly to the observed extragalactic diffuse $\gamma$-ray background (DGB). Using the observed DGB flux from the imaging Compton Telescope (COMPTEL) and the Energetic Gamma Ray Experiment Telescope (EGRET) as the upper limit of $\gamma$-ray flux from PBHs, we derive the upper limit on the spectral index $n$ of the density fluctuations. The range of initial PBH masses which can contribute to the DGB is $2\times 10^{13}\g - 5\times 10^{14}\g$, corresponding to a cosmic reheating temperature of $7\times 10^{7}\gev - 4\times 10^{8} \gev$. In this range, we find the upper limit to be $n \lsim 1.23-1.25$. This limit is stronger than those derived from the energy density in PBHs or PBH relics and matches the value of $n$ required to explain the cosmic microwave background anisotropy.

Primordial Black Holes (PBH) may have formed from the collapse of cosmic string loops. The spectral shape of the PBH mass spectrum can be determined by the scaling argument for string networks. Limits on the spectral amplitude derived from extragalactic gamma-ray and galactic gamma-ray and cosmic ray flux observations as well as constraints from the possible formation of stable black holes remnants are reanalyzed. The new constraints are remarkably close to those derived from the normalization of the cosmic string model to the cosmic microwave background anisotropies.

In theories with cosmic strings, a small fraction of string loops may collapse to form black holes. In this Letter, various constraints on such models involving black holes are considered. Hawking radiation from black holes, gamma and cosmic ray flux limits and constraints from the possible formation of stable black hole remnants are reanalyzed. The constraints which emerge from these considerations are remarkably close to those derived from the normalization of the cosmic string model to the cosmic microwave background anisotropies.

We calculate the flux of ultra high energy photons from individual ordinary (i.e. non-superconducting) cosmic strings and compare the results with the sensitivity of current and proposed TeV and EeV telescopes. Our calculations give only upper limits for the gamma ray flux, since the source of the photons, jets from particle production at cusps, may be weakened by back reaction effects. For the usual cosmic distribution of strings, the predicted bursts from strings with the value of mass per unit length associated with galaxy formation or light strings may just be detectable. A diffuse gamma ray background from light strings may also be seen by the Fly's Eye detector at above $7 \times 10^{10}$ GeV.