2025-11-18T13:43:12.547393

Asking Fast Radio Bursts for More than Reionization History

Shaw, Ghara, Beniamini et al.
We propose different estimators to probe the intergalactic medium (IGM) during epoch of reionization (EoR) using the dispersion measure (${\rm DM}$) of the fast radio bursts. We consider three different reionization histories, which we can distinguish with a total of $\lesssim 1000\,{\rm DM}$ measurements during EoR if their redshifts are known. We note that the redshift derivatives of ${\rm DM}$ are also directly sensitive to the reionization history. The major point of this work is to explore the variance in the ${\rm DM}$ measurements and the information encoded in them. We find that the all-sky average $\overline{\rm DM}(z)$ gets biased from the line-of-sight (LoS) fluctuations in the ${\rm DM}$ measurements introduced by the ionization of IGM during EoR. We find that the ratio $σ_{\rm DM}/\overline{\rm DM}$ depends directly on the ionization bubble sizes as well as the reionization history. On the other hand, we also find that angular variance (coined as $\textit{structure function}$) of ${\rm DM}$ encodes the information about the duration of reionization and the typical bubble sizes as well. We establish the usefulness of variances in ${\rm DM}$ using toy models of reionization and later verify it with the realistic reionization simulations.
academic

Asking Fast Radio Bursts for More than Reionization History

Basic Information

  • Paper ID: 2409.03255
  • Title: Asking Fast Radio Bursts for More than Reionization History
  • Authors: Abinash Kumar Shaw, Raghunath Ghara, Paz Beniamini, Saleem Zaroubi, Pawan Kumar
  • Classification: astro-ph.CO (Cosmology and Extragalactic Astrophysics), astro-ph.HE (High Energy Astrophysical Phenomena)
  • Publication Date: October 17, 2025 (draft version)
  • Paper Link: https://arxiv.org/abs/2409.03255

Abstract

This study proposes utilizing dispersion measures (DM) of fast radio bursts (FRBs) to probe various estimators of the intergalactic medium (IGM) during the epoch of reionization (EoR). The research considers three distinct reionization histories and demonstrates that approximately 1000 DM measurements are sufficient to distinguish between these histories if redshifts are known. The study finds that the redshift derivative of DM is directly sensitive to reionization history. The work emphasizes exploring variance in DM measurements and the information it encodes, revealing that the all-sky averaged DM(z) exhibits bias due to line-of-sight (LoS) fluctuations introduced by IGM ionization during the EoR. The research demonstrates that the σ_DM/DM ratio directly depends on ionization bubble size and reionization history, while the angular variance of DM (structure function) encodes information about reionization duration and typical bubble size.

Research Background and Motivation

Scientific Background

The epoch of reionization (EoR) represents a critical phase in cosmic evolution, marking the transition of the universe from a cold neutral state to a hot ionized state. This process is driven by UV radiation from the first generation of celestial objects and is crucial for understanding the formation and evolution of early cosmic structures.

Limitations of Existing Methods

  1. 21cm Signal Detection: Faces enormous challenges from foreground radiation (~10⁴ times), thermal noise, and radio frequency interference
  2. Indirect Observations: Such as CMB Thomson scattering optical depth and high-redshift quasar spectra, suffer from model dependence and statistical variance issues
  3. Observational Constraints: Current experiments provide only loose constraints on the EoR timeline

Research Motivation

Fast radio bursts (FRBs), as high-energy short-duration radio pulses, have dispersion measures that directly reflect the free electron content along the propagation path, providing a novel and unique perspective for probing the IGM during the EoR. This research aims to develop new detection methods based on FRB DM measurements.

Core Contributions

  1. Proposed Multiple Novel Estimators: Including global average DM, redshift derivatives, global dispersion, and angular dispersion (structure function)
  2. Identified Important Systematic Biases: Demonstrated that LoS fluctuations cause significant deviations in average DM estimates compared to calculations using mean ionization fractions
  3. Established Variance-Physics Associations: The σ_DM/DM ratio is directly correlated with ionization bubble size and reionization history
  4. Developed Structure Function Methods: Angular variance encodes information about reionization duration and typical bubble size
  5. Provided Observational Predictions: Under known redshift conditions, approximately 1000 FRBs can distinguish between different reionization histories

Methodology Details

Task Definition

Utilizing FRB dispersion measures to probe IGM properties during the reionization epoch in the redshift range z=6-15, including reionization history, ionization bubble size distribution, and morphological characteristics.

Core Methodological Framework

1. Dispersion Measure Modeling

The dispersion measure contribution from the IGM is defined as:

DM_IGM(θ,z) = (13/16) * (3cH₀Ωb₀)/(8πGmH) * ∫[z_end to z] (1+z')Δ(θ,z')x_HII(θ,z')/√[Ωm₀(1+z')³+ΩΛ₀] dz'

where Δ is the matter overdensity and x_HII is the ionization fraction.

2. Structure Function Definition

To probe the characteristic scale of ionization bubbles, the DM structure function is defined as:

Ξ(δθ,z) = ⟨[DM_EoR(θ+δθ,z) - DM_EoR(θ,z)]²⟩_{θ,δ̂θ}

3. Key Estimators

  • Average DM: ⟨DM_EoR(θ,z)⟩_θ
  • DM Variance: σ²_DM(z) = ⟨DM_EoR(θ,z)⟩²_θ - DM_EoR²(z)
  • Redshift Derivative: dDM_EoR/dz
  • Structure Function Derivative: ∂²Ξ/∂z∂δθ

Technical Innovations

  1. Systematic Bias Identification: First quantification of bias effects from LoS fluctuations on average DM estimates
  2. Multi-Scale Information Extraction: Combining global average quantities with local fluctuation information
  3. Structure Function Method: Innovative use of second-order derivatives ∂²Ξ/∂z∂δθ to extract instantaneous and local structural information
  4. Toy Model Validation: Verification of method effectiveness through simplified binary ionization field models

Experimental Setup

Simulation Data

  1. Toy Model: Binary ionization field with spherical ionization bubbles uniformly randomly distributed
  2. Realistic Simulations: EoR light cones generated using the GRIZZLY 1D radiative transfer code
  3. Simulation Volume: 500 h⁻¹ Mpc³ with 600³ voxel resolution
  4. Redshift Range: z=6-15

Reionization Histories

Three distinct reionization scenarios are considered:

  • Faster: Rapid reionization (zmid=7.5, δz=1.0)
  • Fiducial: Baseline reionization (zmid=8.5, δz=2.0)
  • Slower: Slow reionization (zmid=9.0, δz=2.5)

Evaluation Metrics

  • Number of FRBs required for 3σ significance discrimination between different reionization histories
  • Bias and variance in DM measurements
  • Shape characteristics of structure function contours

Experimental Results

Main Findings

1. Systematic Bias Effects

Realistic simulations reveal significant differences between average DM accounting for LoS fluctuations and DM calculated using global mean ionization fractions:

  • Slow reionization: bias >2σ
  • Fast reionization: bias ≲1σ
  • Maximum bias at high redshifts (z≳12)

2. Variance-Physics Associations

The σ_DM/DM ratio exhibits clear dependence on reionization history:

  • Toy model: ratio scales as R^(1/2) with ionization bubble radius
  • Realistic simulations: ratio decreases with redshift, opposite to toy model behavior

3. Structure Function Characteristics

The ∂²Ξ/∂z∂δθ contours clearly reflect:

  • Reionization History: Peak redshift tracks zmid, contour compression along z-axis reflects reionization window
  • Ionization Bubble Size: Peak angular position corresponds to typical bubble scale

4. Observational Requirements Estimation

Number of FRBs required to distinguish different reionization histories:

Optimistic Case (perfect removal of low-redshift contributions):

  • Using DM(zmid): ≲100 FRBs
  • Using marginalized DM: 40-90 FRBs

Pessimistic Case (including complete low-redshift contributions):

  • Using DM(zmid): 80-600 FRBs
  • Using marginalized DM: 220-1000 FRBs

Ablation Studies

Ionization Bubble Size Effects

Toy model variations with different bubble radii (R=5,10,20 grid points) demonstrate:

  • Theoretical expectation σDM ∝ R^(1/2) is verified
  • Structure function peak position moves linearly with bubble size
  • Large-scale (δθ>θR) correlations decay rapidly

Reionization History Sensitivity

The slope of dDM_EoR/dz in log-log plots directly encodes reionization rate:

  • Slow history: shallowest slope
  • Fast history: steepest slope

Related Work

Traditional EoR Detection Methods

  1. 21cm Signal: Direct detection but with enormous technical challenges
  2. CMB Polarization: Provides integrated constraints but limited resolution
  3. High-Redshift Quasars: Gunn-Peterson trough and Lyα damping wing

FRB Cosmological Applications

  • Beniamini et al. (2021): First proposal of FRB feasibility for EoR detection
  • Hashimoto et al. (2021): dDM/dz sensitivity to reionization history
  • Heimersheim et al. (2022): Average DM constraints on τTh and reionization midpoint

Novelty of This Work

Compared to existing work, this research:

  1. Systematically studies variance information in DM measurements for the first time
  2. Identifies and quantifies systematic bias from LoS fluctuations
  3. Develops structure function-based methods for ionization bubble scale detection

Conclusions and Discussion

Main Conclusions

  1. Variance Contains Rich Information: DM measurement fluctuations directly encode ionization bubble size and reionization history
  2. Systematic Bias Cannot Be Ignored: LoS fluctuations cause significant bias in average DM estimates, which must be considered in data analysis
  3. Structure Function is Effective: Angular correlation analysis can probe typical ionization bubble scales
  4. Observational Feasibility: Future instruments like SKA can achieve scientific goals through ~1000 high-redshift FRBs

Limitations

  1. Redshift Determination: Method depends on precise FRB redshift measurement (10% uncertainty)
  2. Host Galaxy Contribution: Requires accurate modeling and removal of host galaxy DM contribution
  3. Low-Redshift Contamination: Imperfect removal of post-reionization IGM contributions significantly affects results
  4. FRB Origin: Based on magnetar origin assumption, but FRB physical mechanisms remain incompletely understood

Future Directions

  1. Improved Host Galaxy Modeling: Consider effects of galaxy type, inclination, and FRB location
  2. Redshift Uncertainty Handling: Develop analysis methods robust to redshift errors
  3. FRB Distribution Modeling: More precise modeling of FRB redshift distribution and spatial clustering
  4. Multi-Wavelength Joint Analysis: Combine with other EoR probes such as 21cm and CMB

In-Depth Evaluation

Strengths

  1. Strong Methodological Innovation: First systematic utilization of DM fluctuation information, opening new research directions
  2. Solid Theoretical Foundation: Estimators established from physical principles with rigorous mathematical derivations
  3. Comprehensive Verification: Complete verification chain from simple toy models to complex realistic simulations
  4. High Practical Value: Provides concrete analytical frameworks for forthcoming large FRB sample observations

Weaknesses

  1. Simplified Model Assumptions: Toy models are overly idealized; real situations are more complex
  2. Optimistic Systematic Error Treatment: Treatment of host galaxy and low-redshift contributions may be overly optimistic
  3. Limited Statistical Analysis Depth: Lacks deeper error propagation and correlation analysis
  4. Computational Complexity: Structure function calculations require large FRB samples, potentially challenging for practical applications

Impact

  1. Academic Contribution: Provides new theoretical tools and analytical methods for FRB cosmology research
  2. Practical Value: Offers important reference for scientific planning of next-generation radio instruments like SKA
  3. Interdisciplinary Significance: Bridges high-energy astrophysics and cosmology research fields

Applicable Scenarios

  1. Large-Scale FRB Surveys: High-redshift FRB observations from SKA, CHIME, and similar instruments
  2. Multi-Messenger Astrophysics: Cosmological studies combining gravitational waves, neutrinos, and other messengers
  3. Early Universe Exploration: Auxiliary means for studying first-generation galaxies and black hole formation

References

This paper cites 153 relevant references, primarily including:

  • EoR theory and observational reviews (Zaroubi 2013; Pritchard & Loeb 2012)
  • FRB physics and observations (Petroff et al. 2022; CHIME/FRB Collaboration 2021)
  • Pioneering works in FRB cosmological applications (Beniamini et al. 2021; Hashimoto et al. 2021)
  • EoR numerical simulation methods (Ghara et al. 2015, 2018)

Overall Assessment: This is a high-quality theoretical astrophysics paper that makes important contributions to FRB cosmology applications. The methodology is highly innovative, the analysis is comprehensive and in-depth, and it provides valuable theoretical guidance for future observational research. Despite some model simplifications and assumption limitations, the overall scientific value and practical significance are substantial.