A Unified Magnetohydrodynamic Scaling Relation for the Multiphase Interstellar Medium
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 ten orders of magnitude in density (10−26--10−16 g cm−3) by combining pulsar and Zeeman observations. The unified data are consistent with a continuous magnetic-field evolution organised by the Alfvén Mach number ℳA=√(EK/EB). Within this interpretation, the low-density gas is magnetically dominated (ℳA<1), whereas the high-density gas becomes kinetically dominated (ℳ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, ℳ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$; Bc≈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
@misc{zhao2026unified,
title = {{A Unified Magnetohydrodynamic Scaling Relation for the Multiphase Interstellar Medium}},
author = {Mengke Zhao and Guang-Xing Li and Keping Qiu},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.06900},
eprint = {2609.06900},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.06900},
note = {CC BY 4.0},
}
RIS
TY - UNPB
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/09/07
PB - arXiv
DO - 10.48550/arxiv.2609.06900
UR - https://arxiv.org/abs/2609.06900
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 ten orders of magnitude in density (10^−26--10^−16 g cm^−3) by combining pulsar and Zeeman observations. The unified data are consistent with a continuous magnetic-field evolution organised by the Alfvén Mach number ℳ_A=√(E_K/E_B). Within this interpretation, the low-density gas is magnetically dominated (ℳ_A<1), whereas the high-density gas becomes kinetically dominated (ℳ_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, ℳ_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.
ER -