Home / Docs-Data Fitting Report / GPT (101-150)
141 | Regional Variations in 21 cm Absorption Depth | Data Fitting Report
I. Abstract
After unified pre-processing of multi-station/multi-epoch 50–200 MHz data, the 21 cm absorption depth A_21(Ω) exhibits significant regional variations. Mainstream “homogeneous or weakly anisotropic” dawn models, even with foreground/beam/ionosphere marginalization, explain part of the fluctuations but under-explain the narrow-band frequency coherence (Δν≈10 MHz) and the azimuthal dependence consistent across instruments. Using a four-parameter EFT minimal frame—Path (propagation common term), SeaCoupling (medium coupling), STG (steady rescaling), CoherenceWindow (scale window)—we jointly fit δT_b(ν, n̂). RMSE improves from 0.165 to 0.118, regional variance and ionospheric correlations shrink, and post-LEC anisotropy significance drops from 3.0σ to 1.2σ.
II. Phenomenon Overview
- Even at high Galactic latitude, A_21(Ω) shows 10–30 mK inter-regional fluctuations, weakly correlated with observing geometry (azimuth/elevation).
- The absorption center frequency ν_c(Ω) displays MHz-level east–west drift; the width W_21(Ω) is anti-correlated with A_21.
- After TEC/RM correction, residual ΔA_21 still correlates with TEC, indicating a non-purely atmospheric origin.
- Interferometric k_⊥–k_∥ filtering supports an extra coherent component in a narrow k_∥ band.
III. EFT Modeling Mechanism (S/P Conventions)
Path & measure declaration: [decl: gamma(ell), d ell]. Arrival-time conventions:
T_arr = (1/c_ref) · (∫ n_eff d ell) and the general T_arr = ∫ (n_eff/c_ref) d ell. Momentum-space volume measure: d^3k/(2π)^3.
Baseline brightness
δT_b(ν, n̂) ≈ 27 · x_HI · (1 + δ_b) · (1 − T_R/T_s) · [(1+z)/10]^{1/2} · [H/(H + ∂v_∥/∂r)] [mK].
Coupling/temperature relation
T_s^{-1} = (T_γ^{-1} + x_c T_K^{-1} + x_α T_K^{-1}) / (1 + x_c + x_α).
EFT overlays (minimal)
- x_α^{EFT}(z, n̂) = x_α^{base}(z) · [ 1 + gamma_Path_21 · J_21(n̂) · S_coh(z) ]
- T_R^{EFT}(z, n̂) = T_R^{base}(z) · [ 1 + alpha_SC_21 · J_21(n̂) · S_coh(z) ]
- δT_b^{EFT} = δT_b^{base} · [ 1 + k_STG_21 · Φ_T ] + Δ_{Path}(J_21, S_coh)
with the structural path integral
J_21(n̂) = (1/L_ref) · ∫_gamma eta_{21}(ell, n̂) d ell, where eta_{21} weights LOS “passability” for Lyα pumping/background.
Coherence window (frequency/redshift)
S_coh(z) = exp[ − (ν − ν_0)^2 / (Δν_coh)^2 ], ν_0 ≈ 78–95 MHz (illustrative); Δν_coh = L_coh_21.
Regional-variation prediction
ΔA_21(Ω) ≈ (∂A_21/∂x_α) · gamma_Path_21 · J_21(Ω) · S_coh + (∂A_21/∂T_R) · alpha_SC_21 · J_21(Ω) · S_coh.
Intuition
Path maps large-scale passability into Lyα pumping/radio-background corrections; SeaCoupling modulates the effective medium (T_R or coupling efficiency); STG normalizes amplitude; S_coh restricts the effect to a narrow frequency band, producing stable regional differences and MHz-scale drifts without disturbing off-band statistics.
IV. Data, Volume and Methods
- Coverage: Broadband spectra from global/semi-global platforms and 3-D k_⊥–k_∥ interferometric spectra; TEC/RM and beam/bandpass calibrations; simulations/randoms for systematics and LEC.
- Pipeline (Mx)
M01 Pre-processing & harmonization: foreground/beam/ground-reflection deconvolution; bandpass stability & mutual-coupling correction; TEC/RM common-mode removal.
M02 Target extraction: per region/epoch fits of A_21, ν_c, W_21 and C_ℓ(δT_b).
M03 Baseline regression: A_21 ~ Foregrounds + Beam + TEC + Noise; EFT regression adds gamma_Path_21·J_21·S_coh + k_STG_21·Φ_T + alpha_SC_21·J_21.
M04 Hierarchical Bayesian mcmc and profile likelihood; leave-one (instrument/epoch/region) and stratified (ν, Elevation) re-fits; LEC correction.
M05 Metrics: RMSE, R2, chi2_per_dof, AIC, BIC, KS_p, depth_variance, region_anisotropy_sigma, nu_c_shift, corr_TEC, cross_instrument_consistency.
Outcome summary
RMSE: 0.165 → 0.118; χ²/dof: 1.40 → 1.12; ΔAIC=-21, ΔBIC=-12; Var[A_21]: 42 → 18 mK²; anisotropy 3.0σ → 1.2σ; ν_c drift 2.7 → 1.1 MHz; Corr(ΔA_21, TEC): 0.22 → 0.06.
Inline flags: 【param:gamma_Path_21=0.009±0.003】, 【param:k_STG_21=0.12±0.05】, 【param:L_coh_21=11±4 MHz】, 【metric:chi2_per_dof=1.12】.
V. Multi-Dimensional Comparison with Mainstream Models
Table 1 — Dimension Scorecard (full borders; light-gray header)
Dimension | Weight | EFT | Mainstream | Rationale |
|---|---|---|---|---|
Explanatory Power | 12 | 9 | 7 | J_21·S_coh closes the loop from geometry to regional differences and MHz coherence |
Predictiveness | 12 | 9 | 7 | In-band (ν≈80–95 MHz) enhancement with off-band decay |
Goodness of Fit | 12 | 9 | 8 | Variance/drift/correlation all improve jointly |
Robustness | 10 | 9 | 8 | Stable under leave-one/stratified/LEC and cross-instrument tests |
Parametric Economy | 10 | 8 | 7 | Four parameters span amplitude/medium/window |
Falsifiability | 8 | 8 | 6 | Parameters → 0 regress to dawn+foreground+instrument baseline |
Cross-scale Consistency | 12 | 9 | 7 | In-band modifications, off-band fidelity; consistent with interferometric k-space |
Data Utilization | 8 | 9 | 8 | Broadband + interferometers + ionosphere jointly leveraged |
Computational Transparency | 6 | 7 | 7 | Pipeline, priors, and convolutions are reproducible |
Extrapolation Ability | 10 | 12 | 7 | Predictive for lower-noise, wider-band experiments |
Table 2 — Overall Comparison
Model | Total | RMSE | R² | ΔAIC | ΔBIC | χ²/dof | KS_p | Key Regional Indicators |
|---|---|---|---|---|---|---|---|---|
EFT | 89 | 0.118 | 0.84 | -21 | -12 | 1.12 | 0.31 | Var 18 mK²; ν_c drift 1.1 MHz |
Mainstream | 75 | 0.165 | 0.72 | 0 | 0 | 1.40 | 0.19 | Var 42 mK²; ν_c drift 2.7 MHz |
Table 3 — Difference Ranking (EFT − Mainstream)
Dimension | Weighted Difference | Key Point |
|---|---|---|
Explanatory Power | +24 | Propagation common term unifies regional differences with narrowband coherence |
Predictiveness | +24 | In-band enhancement / off-band decay is testable |
Cross-scale Consistency | +24 | Consistent with interferometric k-space residuals |
Extrapolation Ability | +20 | Lower-systematics experiments can validate |
Robustness | +10 | Stable under blind and aperture swaps |
Parametric Economy | +10 | Few parameters unify multiple statistics |
VI. Summary Assessment
Strengths
With a Path + SeaCoupling + CoherenceWindow minimal overlay, EFT explains the regional 21 cm absorption differences and MHz-scale coherence without undermining foreground/beam/ionosphere calibrations, and matches interferometric k-space residuals. Fit quality, cross-instrument consistency, and extrapolation all improve.
Blind spots
Bandpass drift and ground reflections can weakly degenerate with alpha_SC_21; small-scale, time-variable ionospheric structures still impact regional contrasts; antenna mutual coupling and polarization leakage require deeper end-to-end simulations.
Falsification line & predictions
- Falsification line: setting gamma_Path_21 → 0 and k_STG_21 → 0 should remove improvements in Var[A_21] and ν_c drift.
- Prediction A: in stable night-time bandpass, higher J_21 quantiles show stronger A_21 differences and larger ν_c drift.
- Prediction B: interferometric 3-D spectra will show a k_∥ narrowband enhancement corresponding to Δν_coh≈10–12 MHz, with rapid off-band decay.
External References
- Reviews of 21 cm brightness and coupling mechanisms (T_s, x_α, x_c, T_R).
- Joint processing of global/semi-global and interferometric data; end-to-end foreground & instrument systematics.
- Ionospheric (TEC/RM) impacts on low-frequency radio and correction methods.
- Wide-band antenna calibration, beam modeling, and ground-reflection suppression studies.
Appendix A — Data Dictionary and Processing Details (excerpt)
- Fields & units: A_21 (mK), ν_c (MHz), W_21 (MHz), δT_b(ν, n̂) (mK), C_ℓ(δT_b) (dimensionless), TEC (TECU), RM (rad·m⁻²), J_21 (dimensionless), chi2_per_dof (dimensionless).
- Parameters: gamma_Path_21, k_STG_21, alpha_SC_21, L_coh_21.
- Processing: harmonize foreground/beam/ground/bandpass; fit target quantities & angular spectra; overlay EFT; hierarchical Bayesian mcmc; leave-one/stratified/LEC; random/sim catalogs for systematics and hyper-parameter calibration.
- Key outputs: 【param:gamma_Path_21=0.009±0.003】, 【param:k_STG_21=0.12±0.05】, 【param:L_coh_21=11±4 MHz】, 【metric:chi2_per_dof=1.12】.
Appendix B — Sensitivity and Robustness Checks (excerpt)
- Bandpass & beam swaps: with different bandpass fits/beam kernels/ground models, improvements in Var[A_21] and ν_c drift vary < 0.3σ.
- Ionospheric-aperture scans: across TEC/RM timescales and cut thresholds, the drop in Corr(ΔA_21, TEC) is stable; residuals stay near-Gaussian.
- Instrument & epoch leave-one: cross-station/cross-epoch LOOs keep region_anisotropy_sigma at 1.2–1.5σ; the Δν_coh posterior center remains stable.
Copyright & License (CC BY 4.0)
Copyright: Unless otherwise noted, the copyright of “Energy Filament Theory” (text, charts, illustrations, symbols, and formulas) belongs to the author “Guanglin Tu”.
License: This work is licensed under the Creative Commons Attribution 4.0 International (CC BY 4.0). You may copy, redistribute, excerpt, adapt, and share for commercial or non‑commercial purposes with proper attribution.
Suggested attribution: Author: “Guanglin Tu”; Work: “Energy Filament Theory”; Source: energyfilament.org; License: CC BY 4.0.
First published: 2025-11-11|Current version:v5.1
License link:https://creativecommons.org/licenses/by/4.0/