A Unified Magnetohydrodynamic Scaling Relation for the Multiphase Interstellar Medium
- 1Nanjing University
- 2Yunnan University
- 3Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University)
Abstract
The interplay of magnetic fields, turbulence, and gravity governs the structural evolution of the interstellar medium (ISM) and the initial conditions of star formation, yet observational gaps have long enforced a “broken power-law” description of the magnetic field–density relation. Here, we assemble a unified dataset spanning 10 orders of magnitude in density (10 ⁻²⁶ –10 ⁻¹⁶ g cm ⁻³ ) by combining pulsar and Zeeman observations. The unified data are consistent with a single, continuous magnetic-field evolution organised by the Alfvén Mach number M A = E K / E B . Within this interpretation, the low-density gas is magnetically dominated ( M A < 1 ), whereas the high-density gas becomes kinetically dominated ( M A > 1 ) as gravity increasingly contributes to the kinetic-energy budget, with magnetic tension continuing to influence the collapse geometry. Within the Gradual Transition interpretation, the empirical break density traces the vicinity of the trans-Alfvénic equipartition point, M A = 1 . This Gradual Transition model makes three predictions tested here. Its exponent and background field are fixed in advance by turbulent physics and recovered by the fits ( β ≈ 0.15–0.21 against a predicted 0.147; B c ≈ 2.0 μ G). A dense-gas fit, extrapolated blindly across four decades, passes through the diffuse pulsar data. And, under the adopted scale mappings, the implied magnetic-energy spectrum approaches a k −5/3 -like scaling on large scales and departs from this extrapolation on small scales, where gravitational compression amplifies the field. The broken power law can therefore be viewed as a piecewise approximation to the continuous magnetic equation of state, with its fitted transition density potentially retaining a physical connection to the onset of gravity-driven motions.
Related
- Unified Scaling Relation for Coronal Temperature Anomaly and Solar Wind Acceleration
- The atomic multiphase interstellar medium of galaxies in the COLIBRE simulations
- The multiphase interstellar medium as a common origin for magnetic misalignment and TB parity violation
- The σ_s^2- relation in the multi-phase interstellar medium. Exploring the density distribution with the Cloud Factory simulations
- On the relation between magnetic field strength and gas density in the interstellar medium. II. Density uncertainties and diffuse gas constraints
BibTeX
@article{zhao2026unified,
title = {{A Unified Magnetohydrodynamic Scaling Relation for the Multiphase Interstellar Medium}},
author = {Mengke Zhao and Guang-Xing Li and Keping Qiu},
journal = {The Astrophysical Journal},
year = {2026},
month = oct,
doi = {10.3847/1538-4357/aea203},
url = {https://iopscience.iop.org/article/10.3847/1538-4357/aea203},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - A Unified Magnetohydrodynamic Scaling Relation for the Multiphase Interstellar Medium
AU - Mengke Zhao
AU - Guang-Xing Li
AU - Keping Qiu
PY - 2026
DA - 2026/10/06
JO - The Astrophysical Journal
DO - 10.3847/1538-4357/aea203
UR - https://iopscience.iop.org/article/10.3847/1538-4357/aea203
AB - The interplay of magnetic fields, turbulence, and gravity governs the structural evolution of the interstellar medium (ISM) and the initial conditions of star formation, yet observational gaps have long enforced a “broken power-law” description of the magnetic field–density relation. Here, we assemble a unified dataset spanning 10 orders of magnitude in density (10 ⁻²⁶ –10 ⁻¹⁶ g cm ⁻³ ) by combining pulsar and Zeeman observations. The unified data are consistent with a single, continuous magnetic-field evolution organised by the Alfvén Mach number M A = E K / E B . Within this interpretation, the low-density gas is magnetically dominated ( M A < 1 ), whereas the high-density gas becomes kinetically dominated ( M A > 1 ) as gravity increasingly contributes to the kinetic-energy budget, with magnetic tension continuing to influence the collapse geometry. Within the Gradual Transition interpretation, the empirical break density traces the vicinity of the trans-Alfvénic equipartition point, M A = 1 . This Gradual Transition model makes three predictions tested here. Its exponent and background field are fixed in advance by turbulent physics and recovered by the fits ( β ≈ 0.15–0.21 against a predicted 0.147; B c ≈ 2.0 μ G). A dense-gas fit, extrapolated blindly across four decades, passes through the diffuse pulsar data. And, under the adopted scale mappings, the implied magnetic-energy spectrum approaches a k −5/3 -like scaling on large scales and departs from this extrapolation on small scales, where gravitational compression amplifies the field. The broken power law can therefore be viewed as a piecewise approximation to the continuous magnetic equation of state, with its fitted transition density potentially retaining a physical connection to the onset of gravity-driven motions.
SN - 0004-637X
SN - 1538-4357
ER -