The relation between the turbulent Mach number and observed fractal dimensions of turbulent clouds
Pith reviewed 2026-05-25 10:34 UTC · model grok-4.3
The pith
The turbulent Mach number of molecular clouds can be estimated from the fractal dimension of their column density maps alone.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the turbulent Mach number M can be recovered from the fractal dimension D of the column density via the fitted relation log M(D) = ξ1 (erfc^{-1}[(D - D_min)/Ω] + ξ2), where D_min = 1.55 ± 0.13, Ω = 0.22 ± 0.07, ξ1 = 0.9 ± 0.1 and ξ2 = 0.2 ± 0.2. This mapping is constructed from six simulations and validated by recovering consistent Mach numbers for observed subregions without using velocity information.
What carries the argument
The empirical relation log M(D) = ξ1 (erfc^{-1}[(D - D_min)/Ω] + ξ2) that maps the fractal dimension of column density to the turbulent Mach number.
If this is right
- Mach number estimates become possible using only column density geometry.
- The relation supplies cloud kinematic information without line-of-sight velocity data.
- Star formation rate predictions can be made from observed density structure alone.
- The minimum fractal dimension of column density is bounded near 1.55 for the simulated range.
Where Pith is reading between the lines
- The same formula could be applied to column density maps from other telescopes or wavelengths to produce Mach number maps across entire clouds.
- Varying the simulation driving mechanism or adding magnetic fields might shift the fitted parameters and test the relation's robustness.
- If the inverse error function form persists across different tracers, it could link fractal geometry to other turbulence statistics.
Load-bearing premise
The fractal dimension computed from simulated column density projections has the same meaning and is measured in the same way as the fractal dimension extracted from real observational maps.
What would settle it
Independent measurements of velocity dispersion in the saxophone and quiet subregions that yield Mach numbers far from 10 and 2, respectively.
read the original abstract
Supersonic turbulence is a key player in controlling the structure and star formation potential of molecular clouds (MCs). The three-dimensional (3D) turbulent Mach number, $\mathcal{M}$, allows us to predict the rate of star formation. However, determining Mach numbers in observations is challenging because it requires accurate measurements of the velocity dispersion. Moreover, observations are limited to two-dimensional (2D) projections of the MCs and velocity information can usually only be obtained for the line-of-sight component. Here we present a new method that allows us to estimate $\mathcal{M}$ from the 2D column density, $\Sigma$, by analysing the fractal dimension, $\mathcal{D}$. We do this by computing $\mathcal{D}$ for six simulations, ranging between $1$ and $100$ in $\mathcal{M}$. From this data we are able to construct an empirical relation, $\log\mathcal{M}(\mathcal{D}) = \xi_1(\text{erfc}^{-1} [(\mathcal{D}-\mathcal{D}_{\text{min}})/\Omega] + \xi_2),$ where $\text{erfc}^{-1}$ is the inverse complimentary error function, $\mathcal{D}_{\text{min}} = 1.55 \pm 0.13$ is the minimum fractal dimension of $\Sigma$, $\Omega = 0.22 \pm 0.07$, $\xi_1 = 0.9 \pm 0.1$ and $\xi_2 = 0.2 \pm 0.2$. We test the accuracy of this new relation on column density maps from $Herschel$ observations of two quiescent subregions in the Polaris Flare MC, `saxophone' and `quiet'. We measure $\mathcal{M} \sim 10$ and $\mathcal{M} \sim 2$ for the subregions, respectively, which is similar to previous estimates based on measuring the velocity dispersion from molecular line data. These results show that this new empirical relation can provide useful estimates of the cloud kinematics, solely based upon the geometry from the column density of the cloud.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims that an empirical relation log M(D) = ξ1 (erfc^{-1}[(D - D_min)/Ω] + ξ2) can be constructed from fractal dimensions D measured on 2D column-density projections of six supersonic turbulence simulations (M ranging from 1 to 100). The four fitted parameters (D_min = 1.55 ± 0.13, Ω = 0.22 ± 0.07, ξ1 = 0.9 ± 0.1, ξ2 = 0.2 ± 0.2) are then used to infer M ≈ 10 and M ≈ 2 from Herschel column-density maps of the 'saxophone' and 'quiet' subregions of the Polaris Flare, values stated to be consistent with prior velocity-dispersion estimates.
Significance. If the relation is robust and the D measurements are demonstrably comparable, the work supplies a concrete functional form that converts observed column-density geometry into a Mach-number estimate without requiring line-of-sight velocity data. The explicit parametrization and the direct comparison to two observational fields constitute the main strengths; the small calibration set (six points) and the four free parameters are noted as limiting factors for the claimed generality.
major comments (3)
- [Abstract] Abstract and methods: the relation is calibrated on only six simulation points yet employs four free parameters (D_min, Ω, ξ1, ξ2). No description is given of the fitting procedure, the objective function, or how uncertainties were propagated, so it is unclear whether the erfc^{-1} functional form is uniquely required by the data or simply adopted.
- [Abstract] Abstract, application to Polaris Flare: the transfer of the four fitted parameters to the Herschel maps rests on the premise that D is computed with identical box-counting parameters, spatial filtering, thresholding, and noise treatment in both the simulation projections and the observations. No verification or cross-check of this equivalence is reported, yet any systematic offset in D maps directly into M via the inverse-erfc mapping.
- [Abstract] Abstract: the quoted uncertainties on the four parameters are supplied, but the manuscript does not state how many independent D measurements were extracted per simulation, whether the six points are independent across the M range, or whether the fit residuals justify the claimed precision of the relation.
minor comments (2)
- Notation: the symbol M is used for the Mach number while the abstract also employs script-M; consistent use of a single symbol throughout would improve readability.
- [Abstract] The abstract states that the relation 'can provide useful estimates' but does not quantify the expected uncertainty on the inferred M values when D is measured from real maps; adding a brief error-propagation estimate would strengthen the claim.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments highlight important omissions in the description of our fitting procedure, the equivalence of fractal dimension measurements between simulations and observations, and the statistical details of the calibration. We agree that these points require clarification and will revise the manuscript accordingly to strengthen the presentation. Our point-by-point responses follow.
read point-by-point responses
-
Referee: [Abstract] Abstract and methods: the relation is calibrated on only six simulation points yet employs four free parameters (D_min, Ω, ξ1, ξ2). No description is given of the fitting procedure, the objective function, or how uncertainties were propagated, so it is unclear whether the erfc^{-1} functional form is uniquely required by the data or simply adopted.
Authors: We agree that the manuscript does not describe the fitting procedure. The erfc^{-1} form was adopted empirically because it captures the observed saturation of D toward both low and high Mach numbers; it is not claimed to be the unique functional form. We will add a methods subsection specifying that parameters were obtained via least-squares minimization, with uncertainties from the covariance matrix of the fit. The small number of calibration points (six) and four parameters will be discussed explicitly as a limitation on generality. revision: yes
-
Referee: [Abstract] Abstract, application to Polaris Flare: the transfer of the four fitted parameters to the Herschel maps rests on the premise that D is computed with identical box-counting parameters, spatial filtering, thresholding, and noise treatment in both the simulation projections and the observations. No verification or cross-check of this equivalence is reported, yet any systematic offset in D maps directly into M via the inverse-erfc mapping.
Authors: The referee is correct that explicit verification is missing. The same box-counting implementation, grid resolution, and relative thresholding (above the mean column density) were applied to both the projected simulation maps and the Herschel data; simulations contain no observational noise, so no additional filtering was used. We will insert a paragraph in the methods section that tabulates the processing steps for simulations versus observations to demonstrate equivalence and will note any remaining differences. revision: yes
-
Referee: [Abstract] Abstract: the quoted uncertainties on the four parameters are supplied, but the manuscript does not state how many independent D measurements were extracted per simulation, whether the six points are independent across the M range, or whether the fit residuals justify the claimed precision of the relation.
Authors: One D value was extracted per simulation (the full 2D projection), producing six independent points spanning M = 1–100. Parameter uncertainties were obtained from the fit covariance; we will report the number of points, confirm independence, and include the fit residuals (or reduced chi-squared) to support the quoted precision. The limited sample size will be acknowledged as restricting the robustness of the relation. revision: yes
Circularity Check
No significant circularity; empirical fit applied to independent data
full rationale
The paper computes fractal dimensions D from 2D projections of six simulations spanning known Mach numbers M=1–100, then explicitly fits the four parameters of the erfc^{-1} functional form to those (D, M) pairs. The resulting relation is then applied to D values measured from separate Herschel column-density maps of the Polaris Flare. This is a standard empirical calibration followed by out-of-sample application; the observational estimates are not forced by construction from the same dataset. No self-citations, self-definitional steps, or renamings of known results appear in the derivation chain. The assumption that the box-counting algorithm yields comparable D values is a methodological premise, not a circular reduction of the claimed relation itself.
Axiom & Free-Parameter Ledger
free parameters (4)
- D_min =
1.55 ± 0.13
- Omega =
0.22 ± 0.07
- xi_1 =
0.9 ± 0.1
- xi_2 =
0.2 ± 0.2
axioms (1)
- domain assumption The fractal dimension of 2D column density maps correlates with the 3D Mach number across the simulated range in a manner captured by the chosen erfc functional form
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.