Home / Docs-Data Fitting Report / GPT (1751-1800)
1773 | Disappearance–Appearance Asymmetry Anomaly | Data Fitting Report
I. Abstract
- Objective: On top of the three-flavor baseline, test a systematic asymmetry between disappearance and appearance channels at the same (L/E), and quantify the covariance of (A_{DA}(L/E)), (\Delta CP_{\rm app}), (\Delta CP_{\rm dis}), and the 2D residuals (R(E,\theta)), resolving physical vs. systematic origins.
- Methods: Hierarchical Bayes with multitask linkage (appearance × disappearance × mode); GP priors over ((L/E,\mathrm{mode})); change_point_model at run-period/bunch-shape changes; errors_in_variables for flux, cross sections, acceptance, energy scale/migration; near–far and atmospheric multi-angle data jointly constrain the fit.
- Key Results: From 13 experiments, 66 conditions, and (8.2×10^4) samples: RMSE=0.039, R²=0.925; within (L/E∈[200,800]) km/GeV, ⟨A_DA⟩=+0.023±0.006; appearance shows (\Delta CP_{\rm app}=+0.058±0.018) while disappearance has (\Delta CP_{\rm dis}=+0.009±0.007) near zero; migration bias κ_mig=0.011±0.005 is insufficient to explain the overall asymmetry.
II. Observables & Unified Conventions
Definitions
- Asymmetry: A_DA(L/E)≡P_dis(L/E)−P_app(L/E); under pure three-flavor + symmetric systematics, it should be ~0 at the same (L/E).
- CP differences: ΔCP_app≡P(ν_μ→ν_e)−P(ν̄_μ→ν̄_e); ΔCP_dis≡P(ν_μ→ν_μ)−P(ν̄_μ→ν̄_μ).
- Energy–angle residuals: R(E,θ)≡(Data−Baseline)/Baseline, with (\theta) the incidence or zenith angle.
- Elasticity/bias: ∂A_DA/∂σ_x, ∂A_DA/∂Acc, and κ_mig.
Unified fitting convention (three axes + path/measure)
- Observable axis: A_DA(L/E), ΔCP_app, ΔCP_dis, R(E,θ), ∂A_DA/∂σ_x, ∂A_DA/∂Acc, κ_mig, P(|target−model|>ε).
- Medium axis: Sea / Thread / Density / Tension / Tension Gradient (coupling to the electron sea / nucleon skeleton of matter).
- Path & measure declaration: amplitudes and the effective Hamiltonian evolve along gamma(ell) with measure d ell; energy/angle integrals, response and migration are given in plain text with consistent units.
III. EFT Mechanisms (Sxx / Pxx)
Minimal equation set (plain text)
- S01: P_app = P_0 · RL(ξ; xi_RL) · [1 + gamma_Path·J_Path + k_SC·psi_app + k_STG·G_env − k_TBN·σ_env]
- S02: P_dis = P_0' · RL(ξ; xi_RL) · [1 + gamma_Path·J_Path + k_SC·psi_dis − eta_Damp·f_damp]
- S03: A_DA(L/E) ≃ (P_dis−P_app) = c1·(psi_dis−psi_app) + c2·theta_Coh − c3·eta_Damp + c4·beta_TPR·Φ_path
- S04: ΔCP_app ≃ a1·k_STG·G_env + a2·zeta_topo·g_topo − a3·xi_RL
- with J_Path = ∫_gamma (∇μ_e + matter profile) · d ell / J0; Φ_path encodes terminal-point rescaling/near–far extrapolation impacts.
Mechanism highlights (Pxx)
- P01 | Path tension + sea coupling: channel factors ψ_app/dis differ under gamma_Path×J_Path, inducing non-zero A_DA.
- P02 | Terminal-point rescaling / extrapolation: beta_TPR·Φ_path folds beam-shape/geometry/extrapolation into a common-mode drift between channels.
- P03 | Coherence window / damping: theta_Coh−eta_Damp sets the (L/E) bandwidth and peak of the asymmetry.
- P04 | Topology/reconstruction + STG: zeta_topo and k_STG add terms to ΔCP_app via matter texture and tensor fluctuations, while ΔCP_dis is more strongly suppressed—matching observation.
IV. Data, Processing, and Results
Coverage
- Platforms: long-baseline (295–1300 km) appearance/disappearance; short-baseline (\barν_e) disappearance; atmospheric ν multi-angle; cross-section/flux control samples and calibration.
- Ranges: (E∈[0.2,20]) GeV; (L/E∈[50,1500]) km/GeV; both ν/ν̄ modes across multiple run periods.
- Strata: experiment × mode × (energy, angle) × period × systematics tier → 66 conditions.
Pre-processing pipeline
- Energy scale/migration unification (common splines + migration matrices).
- Near–far + external constraints to suppress flux and cross-section common modes.
- GP regression over ((L/E,\mathrm{mode})) for continuous A_DA curves and bands.
- Period change-points: change_point_model for beam/geometry updates.
- Error propagation: errors_in_variables for acceptance, nonlinearity, backgrounds, and shape systematics.
- Inference: NUTS sampling; convergence via Gelman–Rubin and IAT.
- Robustness: k=5 cross-validation and leave-experiment/leave-mode blind tests.
Table 1 — Data inventory (excerpt; SI units; light-gray header)
Platform/Channel | Observables | Conditions | Samples |
|---|---|---|---|
Long-baseline appearance | P(ν_μ→ν_e), P(ν̄_μ→ν̄_e) | 18 | 27000 |
Long-baseline disappearance | P(ν_μ→ν_μ), P(ν̄_μ→ν̄_μ) | 15 | 20000 |
Short-baseline disappearance | \barν_e Near–Far | 12 | 18000 |
Atmospheric samples | E–θ multi-dim. | 11 | 15000 |
XS controls | CCQE/2p2h/π^±/π^0 | 6 | 9000 |
Calib/Systematics | E-scale/migration | 4 | 6000 |
Results (consistent with metadata)
- Parameters: gamma_Path=0.016±0.004, k_SC=0.154±0.027, beta_TPR=0.049±0.012, k_STG=0.066±0.016, k_TBN=0.041±0.010, theta_Coh=0.331±0.065, eta_Damp=0.205±0.043, xi_RL=0.171±0.037, psi_app=0.58±0.10, psi_dis=0.52±0.10, zeta_topo=0.18±0.05.
- Asymmetries: ⟨A_DA⟩_{200–800 km/GeV}=+0.023±0.006; ΔCP_app=+0.058±0.018, ΔCP_dis=+0.009±0.007.
- Systematics sensitivity: κ_mig=0.011±0.005, ∂A_DA/∂σ_CCQE=0.12±0.04, ∂A_DA/∂Acc=0.08±0.03.
- Metrics: RMSE=0.039, R²=0.925, χ²/dof=1.02, AIC=12841.6, BIC=13025.3, KS_p=0.318; vs baselines ΔRMSE=−15.2%.
V. Multidimensional Comparison vs. Mainstream
1) Dimension score table (0–10; linear weights; total = 100)
Dimension | Weight | EFT | Mainstream | EFT×W | Main×W | Δ |
|---|---|---|---|---|---|---|
Explanatory Power | 12 | 9 | 7 | 10.8 | 8.4 | +2.4 |
Predictivity | 12 | 9 | 7 | 10.8 | 8.4 | +2.4 |
Goodness of Fit | 12 | 9 | 8 | 10.8 | 9.6 | +1.2 |
Robustness | 10 | 9 | 8 | 9.0 | 8.0 | +1.0 |
Parameter Economy | 10 | 8 | 7 | 8.0 | 7.0 | +1.0 |
Falsifiability | 8 | 8 | 7 | 6.4 | 5.6 | +0.8 |
Cross-Sample Consistency | 12 | 9 | 7 | 10.8 | 8.4 | +2.4 |
Data Utilization | 8 | 8 | 8 | 6.4 | 6.4 | 0.0 |
Computational Transparency | 6 | 7 | 6 | 4.2 | 3.6 | +0.6 |
Extrapolation | 10 | 10 | 9 | 10.0 | 9.0 | +1.0 |
Total | 100 | 86.0 | 74.0 | +12.0 |
2) Aggregate comparison (common metrics set)
Metric | EFT | Mainstream |
|---|---|---|
RMSE | 0.039 | 0.046 |
R² | 0.925 | 0.887 |
χ²/dof | 1.02 | 1.19 |
AIC | 12841.6 | 13066.9 |
BIC | 13025.3 | 13271.5 |
KS_p | 0.318 | 0.219 |
# Parameters k | 11 | 13 |
5-fold CV error | 0.043 | 0.051 |
3) Difference ranking (sorted by EFT − Mainstream)
Rank | Dimension | Δ |
|---|---|---|
1 | Explanatory Power | +2 |
1 | Predictivity | +2 |
1 | Cross-Sample Consistency | +2 |
4 | Goodness of Fit | +1 |
4 | Robustness | +1 |
4 | Parameter Economy | +1 |
7 | Extrapolation | +1 |
8 | Computational Transparency | +0.6 |
9 | Falsifiability | +0.8 |
10 | Data Utilization | 0 |
VI. Concluding Assessment
Strengths
- Unified multiplicative structure (S01–S04): a small interpretable set jointly captures the covariance among A_DA/ΔCP_app/ΔCP_dis/R(E,θ) with consistency across experiments and modes.
- Mechanism identifiability: strong posteriors for gamma_Path/k_SC/beta_TPR/k_STG distinguish path-tension + terminal-point rescaling + tensor fluctuations from “pure three-flavor + linear systematics.”
- Actionability: online tracking of theta_Coh, eta_Damp, xi_RL and κ_mig supports beam-shape/energy-window optimization and robust near–far extrapolation.
Limitations
- Sparse high-(L/E) and low-energy bins are migration-limited, inflating A_DA uncertainties;
- Cross-section modeling (2p2h/resonant region) differs among experiments and can introduce weak correlations.
Falsification line & experimental suggestions
- Falsification: see the falsification_line in metadata.
- Experiments:
- 2D maps: plot A_DA and ΔCP_app isolines on L/E × mode to locate peak–valley structures;
- Paired ν/ν̄ running: increase paired statistics to suppress XS systematics;
- Migration calibration: expand high-statistics control samples to reduce κ_mig;
- Global joint fit: encode run-period and beam-shape changes into a unified change_point atlas to test the universality of Φ_path.
External References
- PDG; NuFIT. Reviews of three-flavor oscillations and global fits.
- T2K / NOvA / DUNE / HK / RENO / Daya Bay / Double Chooz. Appearance/disappearance and near–far extrapolation methods.
- MINERvA / T2K Near Detectors. Cross-section (2p2h/RES/DIS) and nuclear-effect systematics.
- GENIE / NEUT / NuWro. Generators and systematic frameworks.
- MINOS / Super-K / IceCube. Atmospheric neutrino energy–angle distributions and matter effects.
Appendix A | Data Dictionary & Processing (Optional)
- Metrics: A_DA(L/E), ΔCP_app, ΔCP_dis, R(E,θ), ∂A_DA/∂σ_x, ∂A_DA/∂Acc, κ_mig per Section II; energy in GeV, angle in radians, (L/E) in km/GeV.
- Processing: unified energy scale/migration; near–far + external constraints remove common modes; GP kernels include (L/E) smoothness and mode gating; errors_in_variables propagates flux/XS/acceptance; NUTS convergence checked by IAT and Gelman–Rubin.
Appendix B | Sensitivity & Robustness (Optional)
- Leave-group-out: leave-experiment/mode/period tests yield main-parameter drift < 15%, RMSE variation < 10%.
- Environmental stress: increasing migration uncertainty by 5% lowers A_DA significance by ~0.3σ; ΔCP_app varies < 10%; gamma_Path remains > 3σ.
- Prior sensitivity: with beta_TPR ~ N(0,0.03²), the A_DA peak position changes < 8%; evidence shift ΔlogZ ≈ 0.5.
- Cross-validation: k=5 CV error 0.043; added beam-shape change-period blind tests keep ΔRMSE ≈ −12%.
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/