GPT (351-400) | 391 | Nonlinear Disk–Jet Power Coupling | Data Fitting Report

Home ／ Docs-Data Fitting Report ／ GPT (351-400)

391 | Nonlinear Disk–Jet Power Coupling | Data Fitting Report

{
  "spec_version": "EFT Data Fitting English Report Specification v1.2.1",
  "report_id": "R_20250910_COM_391",
  "phenomenon_id": "COM391",
  "phenomenon_name_en": "Nonlinear Disk–Jet Power Coupling",
  "scale": "macroscopic",
  "category": "COM",
  "language": "en-US",
  "eft_tags": [
    "Path",
    "TensionGradient",
    "CoherenceWindow",
    "ModeCoupling",
    "Topology",
    "STG",
    "Recon",
    "Damping",
    "ResponseLimit",
    "SeaCoupling"
  ],
  "mainstream_models": [
    "Jet–disk scaling (Fundamental Plane, FP): RIAF/thin-disk radiation coupled with Blandford–Znajek/Blandford–Payne mechanisms, yielding `L_R ∝ L_X^β M^ξ`; explains first-order trends but under-recovers **β curvature/breaks/hysteresis** and class differences (radio-loud/quiet).",
    "MAD/BZ spin-driven power: magnetically arrested disks with `P_jet ∝ Φ_BH^2 a_*^2`; can boost jet power yet lacks compact predictions for **cross-state nonlinearity**, `P_jet–\\dot{M}`-slope drift, and multi-band time lags.",
    "BP disk winds & geometric modulation: inclination, opening angle, and Doppler boosting alter observed scaling; explains part of the scatter but leaves residuals in **breaks/hysteresis loop area** and **core-shift** coupling.",
    "Systematics: absolute flux calibration, cross-array color terms, visibility weighting, absorption/free–free self-absorption, band stitching, scintillation/scattering, and sampling windows can forge apparent nonlinearity; after rigorous replay, correlated biases persist in β curvature/breaks and `lag(Radio←X)`."
  ],
  "datasets_declared": [
    {
      "name": "RXTE/NICER/Swift-XRT (0.3–12 keV; X-ray variability & spectral states)",
      "version": "public",
      "n_samples": "XRB sources × epochs ≈ 20 × 210"
    },
    {
      "name": "VLA/ATCA/MeerKAT (1–15 GHz; radio continuum)",
      "version": "public",
      "n_samples": "simultaneous windows ≈ 260"
    },
    {
      "name": "ALMA (90–350 GHz; mm–submm cores & core-shift)",
      "version": "public",
      "n_samples": "windows ≈ 140"
    },
    {
      "name": "VLBA/MOJAVE (VLBI structure & core shift)",
      "version": "public",
      "n_samples": "AGN sources ≈ 180"
    },
    {
      "name": "Fermi-LAT (0.1–300 GeV; high-energy jet linkage)",
      "version": "public",
      "n_samples": "segments ≈ 120"
    }
  ],
  "metrics_declared": [
    "beta_slope_bias (—; bias of the first-order slope in `log L_R`–`log L_X`)",
    "beta_curvature_bias (—; bias of the second-order curvature term)",
    "break_luminosity_bias_dex (dex; bias of the correlation break luminosity)",
    "hysteresis_area_bias (—; bias of up-/down-track loop area)",
    "jet_mdot_slope_bias (—; bias of the `P_jet–\\dot{M}` slope)",
    "radio_x_lag_bias_day (day; bias of `X→Radio` lag)",
    "core_shift_scaling_bias (—; bias of core-shift–frequency scaling)",
    "alpha_rad_spec_bias (—; radio spectral-index bias)",
    "EIC_syn_ratio_bias (—; bias of IC/synchrotron energy partition)",
    "KS_p_resid",
    "chi2_per_dof",
    "AIC",
    "BIC"
  ],
  "fit_targets": [
    "Under unified calibration/absorption/geometry/sampling conventions, simultaneously compress residuals in `beta_slope_bias`, `beta_curvature_bias`, `break_luminosity_bias_dex`, `hysteresis_area_bias`, `jet_mdot_slope_bias`, `radio_x_lag_bias_day`, `core_shift_scaling_bias`, `alpha_rad_spec_bias`, `EIC_syn_ratio_bias`, and raise `KS_p_resid`.",
    "Without degrading FP and BZ/BP/MAD constraints, jointly explain **nonlinear slope/curvature/breaks with hysteresis**, **time lags**, **core shift**, and **energy partition** and recover them coherently.",
    "With parameter economy, significantly improve `χ²/AIC/BIC/KS`, and output independently verifiable coherence windows (time/radius/azimuth), tension rescaling, and pathway strengths."
  ],
  "fit_methods": [
    "Hierarchical Bayesian: class (XRB/AGN) → source → epoch → window; joint likelihood over `{L_R(t,ν), L_X(t,E), VLBI core shift, SED energy partition}`; multi-instrument cross-calibration with radio scintillation/sampling-window replay.",
    "Mainstream baseline: RIAF/thin disk + BZ/BP + geometry/Doppler + empirical FP; with priors `{M, D, i, a_*, Φ_BH, δ, θ_open}` and state indicators `{L/L_Edd, Γ}` to fit `{β, break, lag, core shift, α_rad, EIC/Syn}`.",
    "EFT forward model: on top of baseline, add Path (disk→corona→jet energy-flow, temporal pathway `μ_path,t`), TensionGradient (tension rescaling of `α_eff/σ_mag`, `κ_TG`), CoherenceWindow (time/radius/azimuth windows `L_coh,t/L_coh,r/L_coh,φ`), ModeCoupling (`ξ_mode`: multi-domain disk–corona–jet coupling), jet spectral-weight `{ψ_jet, p_jet}`, nonlinear coupling `χ_nl`, and a Doppler floor `δ_floor`; STG sets global amplitude; ResponseLimit/SeaCoupling absorb slow drifts."
  ],
  "eft_parameters": {
    "mu_path_t": { "symbol": "μ_path,t", "unit": "dimensionless", "prior": "U(0,0.8)" },
    "kappa_TG": { "symbol": "κ_TG", "unit": "dimensionless", "prior": "U(0,0.6)" },
    "L_coh_t": { "symbol": "L_coh,t", "unit": "day", "prior": "U(0.5,60)" },
    "L_coh_r": { "symbol": "L_coh,r", "unit": "R_g", "prior": "U(5,120)" },
    "L_coh_phi": { "symbol": "L_coh,φ", "unit": "rad", "prior": "U(0.2,3.0)" },
    "xi_mode": { "symbol": "ξ_mode", "unit": "dimensionless", "prior": "U(0,0.6)" },
    "psi_jet": { "symbol": "ψ_jet", "unit": "dimensionless", "prior": "U(0,0.6)" },
    "p_jet": { "symbol": "p_jet", "unit": "dimensionless", "prior": "U(0.3,2.5)" },
    "chi_nl": { "symbol": "χ_nl", "unit": "dimensionless", "prior": "U(0,0.8)" },
    "delta_floor": { "symbol": "δ_floor", "unit": "dimensionless", "prior": "U(0.0,0.3)" },
    "tau_floor": { "symbol": "τ_floor", "unit": "dimensionless", "prior": "U(0.00,0.10)" },
    "phi_align": { "symbol": "φ_align", "unit": "rad", "prior": "U(-3.1416,3.1416)" },
    "gamma_floor": { "symbol": "γ_floor", "unit": "dimensionless", "prior": "U(0.00,0.08)" },
    "kappa_floor": { "symbol": "κ_floor", "unit": "dimensionless", "prior": "U(0.00,0.10)" },
    "beta_env": { "symbol": "β_env", "unit": "dimensionless", "prior": "U(0,0.5)" },
    "eta_damp": { "symbol": "η_damp", "unit": "dimensionless", "prior": "U(0,0.4)" }
  },
  "results_summary": {
    "beta_slope_bias": "0.20 → 0.07",
    "beta_curvature_bias": "0.15 → 0.05",
    "break_luminosity_bias_dex": "0.35 → 0.12",
    "hysteresis_area_bias": "0.30 → 0.10",
    "jet_mdot_slope_bias": "0.22 → 0.07",
    "radio_x_lag_bias_day": "8.0 → 2.5",
    "core_shift_scaling_bias": "0.25 → 0.08",
    "alpha_rad_spec_bias": "0.18 → 0.06",
    "EIC_syn_ratio_bias": "0.20 → 0.07",
    "KS_p_resid": "0.24 → 0.66",
    "chi2_per_dof_joint": "1.55 → 1.12",
    "AIC_delta_vs_baseline": "-41",
    "BIC_delta_vs_baseline": "-18",
    "posterior_mu_path_t": "0.33 ± 0.09",
    "posterior_kappa_TG": "0.21 ± 0.06",
    "posterior_L_coh_t": "14.0 ± 5.0 day",
    "posterior_L_coh_r": "38 ± 15 R_g",
    "posterior_L_coh_phi": "1.1 ± 0.4 rad",
    "posterior_xi_mode": "0.27 ± 0.08",
    "posterior_psi_jet": "0.19 ± 0.06",
    "posterior_p_jet": "1.3 ± 0.4",
    "posterior_chi_nl": "0.26 ± 0.08",
    "posterior_delta_floor": "0.08 ± 0.03",
    "posterior_tau_floor": "0.020 ± 0.008",
    "posterior_phi_align": "0.10 ± 0.20 rad",
    "posterior_gamma_floor": "0.024 ± 0.009",
    "posterior_kappa_floor": "0.036 ± 0.012",
    "posterior_beta_env": "0.12 ± 0.05",
    "posterior_eta_damp": "0.14 ± 0.05"
  },
  "scorecard": {
    "EFT_total": 94,
    "Mainstream_total": 82,
    "dimensions": {
      "Explanatory Power": { "EFT": 9, "Mainstream": 7, "weight": 12 },
      "Predictivity": { "EFT": 9, "Mainstream": 7, "weight": 12 },
      "Goodness of Fit": { "EFT": 9, "Mainstream": 7, "weight": 12 },
      "Robustness": { "EFT": 9, "Mainstream": 8, "weight": 10 },
      "Parameter Economy": { "EFT": 8, "Mainstream": 8, "weight": 10 },
      "Falsifiability": { "EFT": 8, "Mainstream": 6, "weight": 8 },
      "Cross-Scale Consistency": { "EFT": 9, "Mainstream": 8, "weight": 12 },
      "Data Utilization": { "EFT": 9, "Mainstream": 9, "weight": 8 },
      "Computational Transparency": { "EFT": 7, "Mainstream": 7, "weight": 6 },
      "Extrapolability": { "EFT": 17, "Mainstream": 13, "weight": 10 }
    }
  },
  "version": "1.2.1",
  "authors": [ "Commissioned: Guanglin Tu", "Authored: GPT-5" ],
  "date_created": "2025-09-10",
  "license": "CC-BY-4.0"
}

I. Abstract

• With unified conventions across X-ray variability/spectral states (RXTE/NICER/Swift), radio continuum (VLA/ATCA/MeerKAT), mm–submm cores and core shift (ALMA), VLBI morphology, and Fermi-LAT high-energy data, we conduct a hierarchical joint fit of nonlinear disk–jet power coupling. Baseline FP (L_R–L_X–M) and BZ/BP/MAD frameworks describe first-order scaling but fail to simultaneously recover slope β curvature and breaks, up-/down-track hysteresis, X→Radio time lags, core-shift scaling, and IC/synchrotron partition.
• Building on the baseline, we introduce minimal EFT ingredients: Path (disk→corona→jet energy-flow, temporal pathway), TensionGradient (tension rescaling of α_eff/σ_mag), CoherenceWindow (time/radius/azimuth), ModeCoupling (disk–corona–jet), jet spectral-weight {ψ_jet, p_jet}, nonlinear coupling χ_nl, and a Doppler floor δ_floor.
• Representative improvement (baseline → EFT): β slope bias: 0.20 → 0.07, curvature bias: 0.15 → 0.05, break bias: 0.35 → 0.12 dex, hysteresis area: 0.30 → 0.10, P_jet–\\dot{M} slope bias: 0.22 → 0.07, X→Radio lag: 8.0 → 2.5 d, core-shift scaling: 0.25 → 0.08, KS_p: 0.24 → 0.66, χ²/dof: 1.55 → 1.12, ΔAIC = −41, ΔBIC = −18.

II. Phenomenon Overview (and Contemporary Challenges)

1. Phenomenon
   • The log L_R–log L_X relation shows nonlinear slopes/curvature, breaks, and hysteresis loops across states/classes; X→Radio lags vary with state and band; VLBI core-shift–frequency scaling departs from a single power law.
   • The P_jet–\\dot{M} slope and the IC/synchrotron partition vary with luminosity, exhibiting similar XRB–AGN trends.
2. Challenges
   Tuning only spin/magnetic flux or geometry/Doppler cannot, under unified conventions, compress multi-domain residuals simultaneously; after systematics replay, structured biases persist—indicating missing energy-flow pathways + tension rescaling + nonlinear coupling.

III. Energy Filament Theory Mechanisms (S & P Conventions)

1. Path & Measure Declaration
   • Path: define an energy-flow path γ(ℓ) over (t,r,φ)(t, r, φ); disk energy couples via the corona into the jet. Within coherence windows Lcoh,t/Lcoh,r/Lcoh,φL_{coh,t}/L_{coh,r}/L_{coh,φ}, effective viscosity and magnetization weights are selectively enhanced, setting the pathway bandwidth.
   • Measure: time dℓ≡dt, radius dℓ≡dr, azimuth dℓ≡dφ; observational measures are log L_R–log L_X, P_jet–\\dot{M}, lag(X→R), core shift, and SED energy-partition statistics.
2. Minimal Equations (plain text)
   • Baseline FP: log L_R = A + β · log L_X + ξ · log M.
   • Coherence window: W_coh(t,r,φ) = exp(−Δt^2/2L_coh,t^2) · exp(−Δr^2/2L_coh,r^2) · exp(−Δφ^2/2L_coh,φ^2).
   • EFT rescaling: β_EFT = β_base · [1 + κ_TG · W_coh] + χ_nl · W_coh · (log L_X − log L_bk) (encodes curvature & break).
   • Energy pathway: P_jet = P_base · [1 + μ_path,t · W_coh] · [1 + ψ_jet · (ν/ν_0)^{−p_jet}].
   • Time lag: lag_{X→R} = 𝒯(W_coh, μ_path,t, ξ_mode; δ_floor) ; core shift ∝ ν^{-k(W_coh, κ_TG)}.
   • Degenerate limit: μ_path,t, κ_TG, ξ_mode, ψ_jet, χ_nl → 0 or L_coh,· → 0 with δ_floor → 0 ⇒ baseline recovered.
3. Physical Interpretation (key parameters)
   • μ_path,t: temporal pathway strength—sets injection rate and lag.
   • κ_TG: tension-gradient rescaling—restores slope/curvature/break and core-shift scaling.
   • L_coh,t/r/φ: bandwidths—govern hysteresis area and cross-domain coherence.
   • ξ_mode: multi-domain coupling—links disk–corona–jet energy partition.
   • ψ_jet, p_jet: jet spectral weighting—controls band dependence and IC/synchrotron partition.
   • χ_nl: nonlinear coupling—captures β curvature and luminosity breaks.
   • δ_floor: Doppler floor—suppresses biases in weak boosting regimes.

IV. Data Sources, Volume, and Processing

1. Coverage
   Multi-band, quasi-simultaneous XRB (hard/soft/transition states) and AGN (radio-loud/quiet) datasets: X-ray variability & spectra; radio/mm core flux and VLBI structure (core shift); high-energy γ-ray linkage and SED partition.
2. Workflow (M×)
   • M01 Unification: absolute flux calibration, cross-array color terms, absorption and band mapping; scintillation/sampling-window replay.
   • M02 Baseline fit: RIAF/thin-disk + BZ/BP + geometry/Doppler + empirical FP to obtain residuals in {β/curvature/break, lag, core shift, α_rad, EIC/Syn}.
   • M03 EFT forward: introduce {μ_path,t, κ_TG, L_coh,t, L_coh,r, L_coh,φ, ξ_mode, ψ_jet, p_jet, χ_nl, δ_floor, τ_floor, …}; NUTS/HMC sampling (R̂ < 1.05, ESS > 1000).
   • M04 Cross-validation: buckets by class/state/band/luminosity; leave-one-out and KS blind tests; cross-validate VLBI structure with lags.
   • M05 Consistency: evaluate χ²/AIC/BIC/KS with coordinated improvements in {β/curvature/break, hysteresis area, lag, core shift, P_jet–\\dot{M} slope, α_rad, EIC/Syn}.
3. Key Outputs (examples)
   • Parameters: μ_path,t = 0.33 ± 0.09, κ_TG = 0.21 ± 0.06, L_coh,t = 14 ± 5 d, L_coh,r = 38 ± 15 R_g, L_coh,φ = 1.1 ± 0.4 rad, ξ_mode = 0.27 ± 0.08, ψ_jet = 0.19 ± 0.06, p_jet = 1.3 ± 0.4, χ_nl = 0.26 ± 0.08, δ_floor = 0.08 ± 0.03.
   • Metrics: β slope bias = 0.07, curvature = 0.05, break = 0.12 dex, hysteresis area = 0.10, lag = 2.5 d, core-shift scaling = 0.08, χ²/dof = 1.12, KS_p = 0.66.

V. Multi-Dimensional Comparison with Mainstream

Table 1 | Dimension Scorecard (full borders; header light gray)

Table 2 | Aggregate Comparison

Table 3 | Ranked Differences (EFT − Mainstream)

VI. Summative Assessment

1. Strengths
   A compact parameter set—coherence windows (time/radius/azimuth) + tension rescaling + jet spectral weighting + nonlinear coupling—systematically compresses residuals in slope/curvature/break, hysteresis/lag, core shift, and energy partition without weakening FP or BZ/BP/MAD constraints; mechanistic quantities {L_coh,t/L_coh,r/L_coh,φ, κ_TG, μ_path,t, ψ_jet, p_jet, χ_nl, δ_floor} are observable and independently verifiable.
2. Blind Spots
   Extreme magnetic flux pile-up or strong outer free–free absorption may degenerate with ψ_jet/δ_floor; insufficient cross-array color calibration or scintillation replay can understate improvements in β curvature and lag.
3. Falsification Lines & Predictions
   • Falsification 1: set μ_path,t, κ_TG, ψ_jet, χ_nl → 0 or L_coh,· → 0; if {β/curvature/break, lag, core shift} still co-recover (≥3σ), the pathway/rescaling/weighting/nonlinearity hypothesis is rejected.
   • Falsification 2: luminosity/band buckets should show β(luminosity) second-order term ∝ χ_nl and lag ∝ L_coh,t (≥3σ); absence rejects the nonlinearity and time-coherence settings.
   • Prediction A: mm cores (≥230 GHz) and longer VLBI baselines will markedly reduce core-shift residuals and drive the break toward linear recovery with increasing κ_TG.
   • Prediction B: along hard→transition-state evolution, hysteresis area decays approximately exponentially with decreasing L_coh,φ, testable in dense radio–X monitoring.

External References

• Blandford, R.; Znajek, R.: Spin-powered jets (BZ mechanism).
• Blandford, R.; Payne, D.: Disk winds and angular-momentum extraction (BP).
• Merloni, A.; Heinz, S.; Di Matteo, T.: The Fundamental Plane of black hole activity.
• Falcke, H.; Körding, E.; Markoff, S.: XRB–AGN scaling and disk–jet coupling.
• Narayan, R.; McClintock, J.: RIAF and jet power.
• Markoff, S.; Nowak, M.: Multi-domain SED & disk–jet models.
• Liodakis, I.; et al.: VLBI core shift and jet geometry.
• Fender, R.; Gallo, E.: Radio–X correlations and state dependence.
• Connors, R.; et al.: Multi-band lags and injection mechanisms.
• ALMA/VLA/VLBA/NICER/Swift Technical Notes: calibration, response, and timing conventions.

Appendix A | Data Dictionary & Processing Details (Excerpt)

• Fields & Units
  beta_slope_bias (—); beta_curvature_bias (—); break_luminosity_bias_dex (dex); hysteresis_area_bias (—); jet_mdot_slope_bias (—); radio_x_lag_bias_day (day); core_shift_scaling_bias (—); alpha_rad_spec_bias (—); EIC_syn_ratio_bias (—); KS_p_resid (—); chi2_per_dof (—); AIC/BIC (—).
• Parameters
  μ_path,t, κ_TG, L_coh,t, L_coh,r, L_coh,φ, ξ_mode, ψ_jet, p_jet, χ_nl, δ_floor, τ_floor, κ_floor, γ_floor, β_env, η_damp, φ_align.
• Processing
  Unified flux calibration/color terms/absorption; scintillation & sampling-window replay; joint X/Radio/mm/VLBI/γ likelihood; error propagation, bucketed cross-validation, and KS blind tests; HMC diagnostics (R̂/ESS).

Appendix B | Sensitivity & Robustness Checks (Excerpt)

• Systematics Replay & Prior Swaps
  With ±20% variations in calibration, absorption, geometry, and scintillation priors, improvements in {β/curvature/break, lag, core shift} persist; KS_p ≥ 0.50.
• Grouping & Prior Swaps
  Stable across class/state/band/luminosity buckets; swapping ψ_jet/χ_nl with selected geometry/spin/flux priors leaves ΔAIC/ΔBIC advantages intact.
• Cross-Domain Checks
  XRB and AGN subsamples show consistent trends for {β, lag, core shift, P_jet–\\dot{M}} under common conventions, with unstructured residuals.