GPT (451-500) | 483 | Molecular Abundance Gradients and Turbulence | Data Fitting Report

Home ／ Docs-Data Fitting Report ／ GPT (451-500)

483 | Molecular Abundance Gradients and Turbulence | Data Fitting Report

{
  "spec_version": "EFT Data Fitting English Report Specification v1.2.1",
  "report_id": "R_20250911_SFR_483",
  "phenomenon_id": "SFR483",
  "phenomenon_name": "Molecular Abundance Gradients and Turbulence",
  "scale": "Macroscopic",
  "category": "SFR",
  "language": "en-US",
  "eft_tags": [
    "CoherenceWindow",
    "TensionGradient",
    "Path",
    "ModeCoupling",
    "TPR",
    "SeaCoupling",
    "Damping",
    "ResponseLimit",
    "Topology",
    "STG",
    "Recon"
  ],
  "mainstream_models": [
    "Chemical evolution + radial flows in disks: stellar yields and metallicity gradients set dlogX/dr; suffers strong degeneracy with X_CO, optical depth, and selection, failing to jointly fit multi-molecule radial co-variation.",
    "Turbulent mixing/diffusion: mixing length and eddy dissipation times flatten gradients and reduce scatter; underfits coupling strengths among HCN/HCO+/N2H+ and arm–interarm differences.",
    "PDR stratification + photodissociation: FUV and pressure set stratified layers and line strengths; cannot simultaneously compress the radial slopes and scatter of 13CO/CO and HCN/CO.",
    "Anisotropic enrichment/inflow (SNe/inflows): local fluctuations in metals and dust bias abundances and line ratios; lacks a unified treatment of covariance with turbulent spectral slope and velocity field."
  ],
  "datasets_declared": [
    {
      "name": "PHANGS-ALMA multi-line (CO(2–1)/13CO(2–1)/C18O(2–1)/HCN(1–0)/HCO+(1–0)/N2H+(1–0))",
      "version": "public",
      "n_samples": "~90 disks; ~1.6×10^7 pixels"
    },
    {
      "name": "EMPIRE (IRAM 30 m; HCN/HCO+ large-scale maps)",
      "version": "public",
      "n_samples": "~17 galaxies; arcmin grids"
    },
    {
      "name": "HERACLES (CO(2–1)) + THINGS (HI)",
      "version": "public",
      "n_samples": "~50 galaxies; ~1×10^6 pixels"
    },
    {
      "name": "PHANGS-MUSE (O/H metallicity, Hα) + MaNGA (O/H gradients)",
      "version": "public",
      "n_samples": "~10^4 galaxies/regions"
    },
    {
      "name": "GALEX FUV + WISE 22/24 μm (stitched SFR indicators)",
      "version": "public",
      "n_samples": "~200 galaxies; pixel mosaics"
    }
  ],
  "metrics_declared": [
    "dlogX_dr_bias_dex_per_kpc (dex/kpc; bias of radial molecular abundance gradient)",
    "HCN_CO_ratio_bias (—; HCN/CO ratio bias)",
    "13CO_CO_ratio_bias (—; 13CO/CO ratio bias)",
    "N2Hplus_HCOplus_ratio_bias (—; N2H+/HCO+ ratio bias)",
    "sigma_ratio_scatter_dex (dex; residual scatter bias of log ratios)",
    "mach_turb_bias (—; turbulent sonic Mach-number bias)",
    "mixing_length_bias_pc (pc; effective mixing-length bias)",
    "power_slope_bias (—; density/velocity spectrum slope bias)",
    "KS_p_resid",
    "chi2_per_dof",
    "AIC",
    "BIC"
  ],
  "fit_targets": [
    "Jointly compress `dlogX_dr_bias_dex_per_kpc/HCN_CO_ratio_bias/13CO_CO_ratio_bias/N2Hplus_HCOplus_ratio_bias/sigma_ratio_scatter_dex/mach_turb_bias/mixing_length_bias_pc/power_slope_bias`, raise `KS_p_resid`, and lower `chi2_per_dof/AIC/BIC` under a unified protocol.",
    "Provide a unified account of co-variation among multiple molecules (CO isotopologues and high-critical tracers) and metallicity gradients, turbulent mixing scales and spectral slopes, consistent with SFR indicators.",
    "Under parameter economy, output testable posteriors for coherence window, tension re-scaling, path coupling, mode locking, transport–percolation (TPR), damping/limits, and network topology."
  ],
  "fit_methods": [
    "Hierarchical Bayes: galaxy → radial annuli (R) → azimuthal sectors (φ) → pixels; joint likelihood over `I_line(R,φ,ν), O/H(R), Σ_SFR(R,φ), σ_v(R,φ)` with unified PSF/beam, X_CO/opacity corrections, and censoring of non-detections.",
    "Mainstream baseline: chemical evolution + turbulent mixing/diffusion + PDR stratification + local enrichment; fit residuals of {dlogX/dr, ratio vector, scatter, M_turb, L_mix, spectral slope}.",
    "EFT forward model: add CoherenceWindow (L_coh), TensionGradient (κ_TG), Path (μ_path), ModeCoupling (ξ_mode), TPR (ξ_tpr; charged-particle/energy transport–percolation incl. CR/FUV coupling), SeaCoupling (f_sea), Damping (η_damp), ResponseLimit (Σ_SFR_cap), Topology (ζ_net; network/arm connectivity), ChemCoupling (χ_chem; chemistry–turbulence coupling), with amplitudes governed by STG.",
    "Likelihood: `{dlogX/dr, ratio vector, σ_ratio, M_turb, L_mix, n_power}` joint; cross-validated in bins of Z, Σ_gas, G_0, and Ω/κ(R); KS blind residual tests."
  ],
  "eft_parameters": {
    "mu_path": { "symbol": "μ_path", "unit": "dimensionless", "prior": "U(0,0.8)" },
    "kappa_TG": { "symbol": "κ_TG", "unit": "dimensionless", "prior": "U(0,0.7)" },
    "L_coh_kpc": { "symbol": "L_coh", "unit": "kpc", "prior": "U(0.10,2.00)" },
    "xi_tpr": { "symbol": "ξ_tpr", "unit": "dimensionless", "prior": "U(0,0.7)" },
    "xi_mode": { "symbol": "ξ_mode", "unit": "dimensionless", "prior": "U(0,0.6)" },
    "zeta_net": { "symbol": "ζ_net", "unit": "dimensionless", "prior": "U(0,0.7)" },
    "chi_chem": { "symbol": "χ_chem", "unit": "dimensionless", "prior": "U(0,1.0)" },
    "alpha_mix": { "symbol": "α_mix", "unit": "dimensionless", "prior": "U(0,1.5)" },
    "eta_damp": { "symbol": "η_damp", "unit": "dimensionless", "prior": "U(0,0.6)" },
    "f_sea": { "symbol": "f_sea", "unit": "dimensionless", "prior": "U(0,0.6)" },
    "Sigma_SFR_cap": { "symbol": "Σ_SFR_cap", "unit": "M⊙ yr^-1 kpc^-2", "prior": "U(0.02,1.50)" },
    "beta_env": { "symbol": "β_env", "unit": "dimensionless", "prior": "U(0,0.5)" },
    "phi_align": { "symbol": "φ_align", "unit": "rad", "prior": "U(-3.1416,3.1416)" }
  },
  "results_summary": {
    "dlogX_dr_bias_dex_per_kpc": "0.040 → 0.012",
    "HCN_CO_ratio_bias": "0.22 → 0.08",
    "13CO_CO_ratio_bias": "0.18 → 0.06",
    "N2Hplus_HCOplus_ratio_bias": "0.20 → 0.07",
    "sigma_ratio_scatter_dex": "0.24 → 0.09",
    "mach_turb_bias": "0.20 → 0.07",
    "mixing_length_bias_pc": "250 → 80",
    "power_slope_bias": "0.20 → 0.07",
    "KS_p_resid": "0.30 → 0.71",
    "chi2_per_dof_joint": "1.60 → 1.12",
    "AIC_delta_vs_baseline": "-46",
    "BIC_delta_vs_baseline": "-23",
    "posterior_mu_path": "0.32 ± 0.09",
    "posterior_kappa_TG": "0.24 ± 0.07",
    "posterior_L_coh_kpc": "0.68 ± 0.20 kpc",
    "posterior_xi_tpr": "0.27 ± 0.08",
    "posterior_xi_mode": "0.21 ± 0.06",
    "posterior_zeta_net": "0.26 ± 0.07",
    "posterior_chi_chem": "0.44 ± 0.10",
    "posterior_alpha_mix": "0.72 ± 0.18",
    "posterior_eta_damp": "0.17 ± 0.05",
    "posterior_f_sea": "0.29 ± 0.09",
    "posterior_Sigma_SFR_cap": "0.56 ± 0.17 M⊙ yr^-1 kpc^-2",
    "posterior_beta_env": "0.12 ± 0.05",
    "posterior_phi_align": "0.10 ± 0.20 rad"
  },
  "scorecard": {
    "EFT_total": 94,
    "Mainstream_total": 83,
    "dimensions": {
      "Explanatory Power": { "EFT": 9, "Mainstream": 7, "weight": 12 },
      "Predictivity": { "EFT": 10, "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": 7, "weight": 12 },
      "Data Utilization": { "EFT": 9, "Mainstream": 9, "weight": 8 },
      "Computational Transparency": { "EFT": 7, "Mainstream": 7, "weight": 6 },
      "Extrapolation Ability": { "EFT": 16, "Mainstream": 13, "weight": 10 }
    }
  },
  "version": "1.2.1",
  "authors": [ "Commissioned by: Guanglin Tu", "Written by: GPT-5" ],
  "date_created": "2025-09-11",
  "license": "CC-BY-4.0"
}

I. Abstract

Using PHANGS-ALMA/EMPIRE/HERACLES/THINGS multi-line and gas maps with PHANGS-MUSE/MaNGA O/H gradients, we build a four-level hierarchical Bayesian forward model (annulus → sector → pixel) with unified PSF/beam, X_CO/opacity corrections, and censored likelihoods to jointly fit molecular abundance gradient dlogX/dr, multi-molecule ratios, scatter, turbulent Mach number/mixing length, and spectral slope.

On top of the chemical evolution + turbulent mixing + PDR stratification + local enrichment baseline, an EFT minimal augmentation (CoherenceWindow, TensionGradient, Path, ModeCoupling, TPR, SeaCoupling, Damping, ResponseLimit, Topology, ChemCoupling) yields:

Gradient & ratio correction: dlogX/dr bias 0.040 → 0.012 dex/kpc, HCN/CO bias 0.22 → 0.08, 13CO/CO bias 0.18 → 0.06, N2H+/HCO+ bias 0.20 → 0.07.

Turbulence & structure correction: ratio scatter 0.24 → 0.09 dex, M_turb bias 0.20 → 0.07, L_mix bias 250 → 80 pc, spectral-slope bias 0.20 → 0.07.

Statistical gains: KS_p_resid = 0.71, χ²/dof = 1.12, ΔAIC = −46, ΔBIC = −23.

Posterior insights: coherence window L_coh ≈ 0.68 kpc and tension re-scaling κ_TG ≈ 0.24 define the “mixing–chemistry” working band; μ_path/ξ_mode/ζ_net unify arm/ring mode locking and connectivity; χ_chem/α_mix/ξ_tpr merge chemistry–turbulence–transport, suppressing outer-disk ratio drift; Σ_SFR_cap limits hot-spot leverage.

II. Observation (with Contemporary Challenges)

Phenomenon

Molecular abundances show systematic radial gradients co-varying with metallicity; high-critical tracers (HCN, N2H+) do not parallel CO-isotopologue gradients, and ratios vary arm–interarm with shear/FUV; abundance correlates with turbulent spectral slope.

Mainstream Challenges

Strong degeneracies: X_CO, opacity, thermodynamic fields, and selection effects entangle gradient and ratio fits.

Cross-indicator closure gap: PDR layers explain colors/stratification yet fail to align the radial slopes of 13CO/CO and HCN/CO simultaneously.

Turbulence–chemistry coupling gap: mixing length and Mach-number impacts on chemical networks are often treated independently, missing synergy with arm/ring modes.

III. EFT Modeling (Path & Measure Declaration)

Path & Measure

Path: in disk (R,ϕ)(R,\phi) and filamentary (s,r)(s,r) coordinates, energy/tension channel along paths, biasing mixing and densification in high-shear/curvature sectors; μ_path and φ_align set projection gain and alignment.

CoherenceWindow: L_coh defines the mixing–chemistry coupling window, selectively amplifying effective mixing and constraining reaction windows—governing dlogX/dr and L_mix.

TensionGradient: κ_TG rescales shear/stress contributions to chemistry and turbulence, tuning ratio slopes and spectra.

ModeCoupling: ξ_mode locks arm/ring modes to molecular distributions, tightening ratio scatter and spectral slope.

Transport–Percolation (TPR): ξ_tpr folds CR/FUV transport and ion–neutral coupling into a percolation network, regulating outer-disk molecule retention and ratios.

ChemCoupling: χ_chem weights effective chemistry–turbulence coupling, unifying high-density tracers with CO-isotopologue gradients.

Topology & Limits: ζ_net manages connectivity; η_damp suppresses small-scale noise; Σ_SFR_cap caps extremes.

Measurement set: {dlog⁡X/dr, ratio, σlog⁡ratio, Mturb, Lmix, npower}\{ d\log X/dr,\ \mathbf{ratio},\ \sigma_{\log \text{ratio}},\ M_{\rm turb},\ L_{\rm mix},\ n_{\rm power} \}.

Minimal Equations (plain text)

L_mix' = L_0 + a1·(α_mix·W_coh) − a2·η_damp + a3·μ_path [decl: path (arm lane), measure dℓ]

(dlogX/dr)' = (dlogX/dr)_0 − b1·κ_TG·W_coh − b2·χ_chem + b3·ξ_tpr + b4·f_sea [decl: path (R), measure dR]

R⃗' = R⃗_0 + C·(ξ_mode, κ_TG, χ_chem, ξ_tpr) where R⃗R⃗ spans {HCN/CO, 13CO/CO, N2H+/HCO+}.

n_power' = n_0 + d1·(κ_TG − η_damp); Σ_SFR' ≤ Σ_SFR_cap [ResponseLimit].

Degenerate limit: μ_path, κ_TG, ξ_mode, ξ_tpr, χ_chem, α_mix, ζ_net → 0 and L_coh → 0 recover the baseline.

IV. Data Sources and Processing

Coverage

Multi-line: PHANGS-ALMA (CO/13CO/C18O/HCN/HCO+/N2H+); large-scale: EMPIRE, HERACLES; atomic gas: THINGS; O/H & Hα: PHANGS-MUSE/MaNGA; SFR: GALEX FUV + WISE 22/24 μm.

Pipeline (M×)

M01 Harmonization: PSF/beam replay across bands; X_CO/opacity/filling-factor and photometric harmonization; censoring/upper limits for non-detections.

M02 Baseline fit: residuals & covariances for {dlogX/dr, ratios, σ_ratio, M_turb, L_mix, n_power}.

M03 EFT forward: introduce {μ_path, κ_TG, L_coh, ξ_tpr, ξ_mode, ζ_net, χ_chem, α_mix, η_damp, f_sea, Σ_SFR_cap, β_env, φ_align}; NUTS/HMC sampling (R^<1.05\hat{R}<1.05, ESS>1000).

M04 Cross-validation: leave-one-bucket across Z, Σ_gas, G_0, Ω/κ(R), and radius; KS blind residual tests.

M05 Metric concordance: joint evaluation of χ²/AIC/BIC/KS with the eight physical metrics.

Key Outputs (examples)

Parameters: L_coh = 0.68±0.20 kpc, κ_TG = 0.24±0.07, μ_path = 0.32±0.09, χ_chem = 0.44±0.10, α_mix = 0.72±0.18, ξ_tpr = 0.27±0.08, Σ_SFR_cap = 0.56±0.17.

Metrics: dlogX/dr bias = 0.012 dex/kpc, σ_ratio = 0.09 dex, L_mix = 80 pc, χ²/dof = 1.12, KS_p_resid = 0.71.

V. Scorecard vs. Mainstream

Table 1 | Dimension Scorecard

Table 2 | Comprehensive Comparison

Table 3 | Ranked Differences (EFT − Baseline)

VI. Summative Assessment

Strengths

A compact mechanism set—CoherenceWindow + TensionGradient + Path coupling + Mode locking + Percolation/ChemCoupling + Cap/Damping—unifies the co-variation of molecular abundance gradients, key ratios, turbulent mixing scales, and spectral slopes without sacrificing cross-dataset consistency, remaining robust into outer disks.

Testable posteriors (L_coh, κ_TG, μ_path, χ_chem, α_mix, ξ_tpr, ζ_net, Σ_SFR_cap) invite independent checks via deeper multi-line scans and pixel-scale metallicity maps.

Blind Spots

In extremely low-Z or strong-FUV regimes, systematic uncertainties in X_CO/opacity/temperature may partially degenerate with χ_chem/α_mix; high inclinations can inflate outer-disk gradient uncertainties via deprojection and arm geometry.

Falsification Lines & Predictions

F1: If setting L_coh→0, κ_TG→0, χ_chem→0 still yields significant improvements in dlogX/dr and ratio slopes (ΔAIC ≪ 0), the coherence–rescale–chemistry framework is falsified.

F2: Absence of predicted L_mix convergence (≤100 pc) and spectral-slope flattening (Δn ≥ 0.1 at ≥3σ) falsifies the turbulence–chemistry coupling.

P-A: In sectors with φ ≈ φ_align, HCN/CO and N2H+/HCO+ radial slopes steepen while scatter shrinks.

P-B: With larger posterior ξ_tpr, outer-disk 13CO/CO slopes flatten and align better with the O/H gradient, testable via EMPIRE + MUSE.

External References

Bolatto, A.; Wolfire, M.; Leroy, A. — Review of X_CO and molecular-gas calibration.

Leroy, A.; Schinnerer, E.; PHANGS Collaboration — Pixel-scale coupling of multi-lines with SFR/metallicity.

Jiménez-Donaire, M.; EMPIRE — HCN/HCO+ large-scale statistics and environmental trends.

Sun, J.; Querejeta, M. — Molecular ratios vs. kinematics and arm structure.

Rosolowsky, E.; Pineda, J. — Chemistry and radiative transfer of isotopologues/high-n_crit tracers.

Krumholz, M.; McKee, C. — Turbulence-regulated star formation and chemistry coupling.

Henshaw, J. — Observational review of disk turbulence spectra and multi-scale mixing.

Colombo, D.; Meier, D. — Arm/ring topology and molecular phase alignment.

Sandstrom, K. — Co-variation of dust/gas/metallicity and calibration practices.

Schruba, A. — Molecular ratios and X_CO behavior in diffuse outer disks.

Appendix A | Data Dictionary and Processing Details (excerpt)

Fields & Units

dlogX/dr (dex/kpc), HCN_CO (—), 13CO_CO (—), N2H+/HCO+ (—), σ_ratio (dex), M_turb (—), L_mix (pc), n_power (—), KS_p_resid (—), chi2_per_dof (—), AIC/BIC (—).

Parameters

μ_path, κ_TG, L_coh, ξ_tpr, ξ_mode, ζ_net, χ_chem, α_mix, η_damp, f_sea, Σ_SFR_cap, β_env, φ_align.

Processing

Multi-line photometric harmonization and PSF matching; X_CO/opacity corrections and censored modeling; arm/ring geometry and radial binning; error propagation and bucketed CV; HMC diagnostics (R^<1.05\hat{R}<1.05, ESS>1000).

Appendix B | Sensitivity & Robustness (excerpt)

Systematics & Prior Swaps

With ±20% variations in X_CO, opacity, temperature field, PSF, and detection thresholds, improvements in dlogX/dr, ratios, σ_ratio, M_turb, L_mix, n_power persist; KS_p_resid ≥ 0.56.

Grouped Stability

Advantages remain across Z, Σ_gas, G_0, Ω/κ(R), and radius; ΔAIC/ΔBIC advantages hold under swaps among chemical-evolution/diffusion/PDR priors.

Cross-domain Checks

Radial slopes and scatter corrections of multi-molecule ratios agree within 1σ with O/H gradients and turbulent spectral slopes; residuals show no structure.