Geometric Search for Hawking Radiation from Nearby Primordial Black Holes

1. Problem Statement

High-energy transient timing analyses, particularly for Gamma-Ray Bursts (GRBs), universally rely on the plane-wave approximation. This assumption holds because the vast majority of GRBs are cosmological (billions of light-years away), making the curvature of their wavefronts negligible over interplanetary baselines.

However, this approximation creates a blind spot for nearby, short-duration gamma-ray flashes, specifically the final evaporation bursts of Primordial Black Holes (PBHs).

• The Challenge: A nearby burst (within the Solar System or nearby stellar neighborhood) produces a curved wavefront.
• The Degeneracy: Timing-only reconstructions suffer from a fundamental degeneracy where changes in the source's direction can mimic changes in distance. Without an independent measurement of the sky direction, one cannot distinguish between a distant source with a specific arrival time and a nearby source with a different arrival time.
• The Goal: To develop a method that breaks this degeneracy to directly measure the distance of a transient, thereby identifying potential PBH evaporation events which would otherwise be indistinguishable from cosmological GRBs.

2. Methodology

The authors propose a purely geometric method that combines imaging localizations (sky direction) with multi-spacecraft timing (arrival time differences).

The Geometric Model

• Plane-Wave Limit: For a source at infinite distance ($D \to \infty$) from direction $\hat{n}$, the time delay between two detectors $i$ and $j$ is linear:
  $$\Delta t^\infty_{ij} = \frac{(\mathbf{r}_i - \mathbf{r}_j) \cdot \hat{n}}{c}$$
• Finite-Distance Model: For a source at a finite distance $D$ along $\hat{n}$, the delay becomes non-linear due to wavefront curvature:
  $$\Delta t_{ij}(D) = \frac{|\mathbf{r}_i - D\hat{n}| - |\mathbf{r}_j - D\hat{n}|}{c}$$
• Solution: Once the imaging instrument fixes the direction $\hat{n}$, the problem reduces to solving for a single parameter, $D$. The intersection of the "delay hyperboloid" (defined by timing) and the "localization ray" (defined by imaging) yields the distance.

Analytical Framework

The authors derived an analytical inversion of the delay equation to solve for $D$ directly. They defined a residual $\delta t_{ij}$ (observed delay minus plane-wave prediction) and quantified consistency with the plane-wave hypothesis using a statistical metric $x_{ij}$.

• Curvature Scaling: The deviation from the plane-wave limit scales as $|\delta t| \propto B_\perp^2 / D$, where $B_\perp$ is the baseline component perpendicular to the source direction.
• Uncertainty Propagation: The total uncertainty combines timing errors (via Modified Cross-Correlation Function, MCCF) and localization errors ($\sigma_\theta$). The method determines a lower limit $D_{min}$ if no curvature is detected, or a specific distance $D$ if curvature is resolved.

Data Sample

• Selection: 12 short GRBs ($T_{90} < 3$ s) localized by imaging instruments (Swift/BAT) and jointly observed by Konus-Wind.
• Criteria: Events without measured redshift or secure host associations (to avoid pre-assumed cosmological distances).
• Benchmarks: Cross-checked with near-Earth networks (Fermi/GBM, Insight-HXMT, GECAM).

3. Key Contributions

1. First Quantitative Framework: The paper presents the first systematic framework for event-by-event direct distance inference for gamma-ray transients by combining imaging and timing, specifically targeting the PBH search regime.
2. Breaking the Degeneracy: It demonstrates how fixing the sky direction via imaging resolves the distance-direction degeneracy inherent in pure timing triangulation.
3. Feasibility Analysis: It establishes that current technology is already capable of probing distances up to $10^3$ AU, and future missions could extend this to $10^5$ AU.
4. Physical Inference: It links a measured distance $D$ directly to the physical properties of a PBH (Mass and Lifetime) under the Hawking radiation interpretation, moving beyond rate-based statistical limits.

4. Results

• Current Sample Analysis: Applied to the 12 Swift/Konus-Wind short GRBs, no significant deviation from the plane-wave expectation was found. All events were consistent with cosmological distances within $3\sigma$.
• Constraints: The most constraining event yielded a lower limit of $D_{min} \approx 1.2$ AU.
  • Significance: While no PBH was found, this proves the method has entered the Solar-System-scale regime. Current non-detections reflect the limitations of the specific sample and detector configurations (timing/localization precision), not the method itself.
• Distance Reach:
  • Current Capabilities: Can measure distances from AU scales up to $\sim 10^3$ AU.
  • Future Roadmap:
    • Wind-Earth Baseline ($\sim 10^6$ km) + Arcmin Localization: Reach $10 - 10^2$ AU.
    • Sun-Earth Baseline + Arcmin Localization: Reach $\sim 10^3$ AU.
    • Mars-Earth Baseline + 0.1 arcmin Localization: Reach $\sim 10^4$ AU.
    • Sun-Earth Baseline + Arcsec Localization: Reach $\sim 10^5$ AU.

5. Significance and Implications

• PBH Detection: A confirmed finite distance for a short GRB would immediately rule out a cosmological origin and strongly favor an exotic nearby source, such as a PBH evaporation.
• Direct Physical Parameters: Unlike statistical rate limits, a geometric distance measurement allows for the direct calculation of the PBH's mass ($M$) and lifetime ($\tau$) for a specific event:
  $$M \propto D \cdot F^{-1/2}$$
  $$\tau \propto D^3 \cdot F^{-3/2}$$
  (Where $F$ is the observed flux).
• Technological Drivers: The paper identifies that the primary bottlenecks are no longer clock accuracy (which is sufficient) but rather:
  1. Wide-field localization accuracy (moving from arcminutes to arcseconds).
  2. High-time-resolution data products (sub-millisecond binning).
  3. Long baselines (e.g., Mars-Earth or Sun-Earth configurations).
• Strategic Roadmap: The authors provide a clear path for future missions (e.g., piggybacking gamma-ray detectors on Mars missions or utilizing the Einstein Probe) to systematically search for PBHs in the "nearby" universe ($10^2 - 10^5$ AU), a regime previously inaccessible to direct geometric measurement.

In summary, this work transforms the search for nearby Primordial Black Holes from a statistical exercise into a direct geometric measurement problem, demonstrating that current and near-future space assets are on the verge of being able to detect and characterize these exotic objects if they exist within our solar neighborhood.