Statistical analysis of the mass-to-flux ratio in turbulent cores: effects of magnetic field reversals and dynamo amplification

We study the mass-to-flux ratio (M/\Phi) of clumps and cores in simulations of supersonic, magnetohydrodynamical turbulence for different initial magnetic field strengths. We investigate whether the (M/\Phi)-ratio of core and envelope, R = (M/\Phi){core}/(M/\Phi){envelope} can be used to distinguish between theories of ambipolar diffusion and turbulence-regulated star formation. We analyse R for different Lines-of-Sight (LoS) in various sub-cubes of our simulation box. We find that, 1) the average and median values of |R| for different times and initial magnetic field strengths are typically greater, but close to unity, 2) the average and median values of |R| saturate at average values of |R| ~ 1 for smaller magnetic fields, 3) values of |R| < 1 for small magnetic fields in the envelope are caused by field reversals when turbulence twists the field lines such that field components in different directions average out. Finally, we propose two mechanisms for generating values |R| ~< 1 for the weak and strong magnetic field limit in the context of a turbulent model. First, in the weak field limit, the small-scale turbulent dynamo leads to a significantly increased flux in the core and we find |R| ~< 1. Second, in the strong field limit, field reversals in the envelope also lead to values |R| ~< 1. These reversals are less likely to occur in the core region where the velocity field is more coherent and the internal velocity dispersion is typically subsonic.

💡 Research Summary

The paper investigates whether the mass‑to‑flux ratio (M/Φ) of dense cores relative to their surrounding envelopes, expressed as the dimensionless quantity R = (M/Φ)_core / (M/Φ)_envelope, can serve as a diagnostic to discriminate between two competing theories of star formation: ambipolar‑diffusion‑regulated collapse and turbulence‑regulated collapse. Using three‑dimensional, ideal magnetohydrodynamic (MHD) simulations performed with the FLASH v2.5 code, the authors generate a suite of supersonic turbulent clouds (rms Mach number ≈ 10) with five different initial plasma‑β values (β₀ = 0.01, 0.1, 1, 10, 100), corresponding to magnetic field strengths ranging from 44 µG down to 0.44 µG. The simulations are isothermal (cₛ = 2 km s⁻¹) and driven by a solenoidal Ornstein‑Uhlenbeck forcing that injects energy on the largest scales (k L/2π = 1–3). Each run is evolved on a fixed 256³ grid (with additional resolution tests at 128³ and 512³) until a statistically steady turbulent state is reached (t ≈ 2 T, where T = L/(2Mcₛ) is the eddy turnover time). The analysis focuses on three snapshots (t = 2.0 T, 2.4 T, 2.8 T) to assess temporal stability.

Core–envelope identification is performed in two complementary ways. First, a simple geometric prescription selects the densest cell in each clump, surrounds it with a spherical core of diameter five grid cells (≈ 0.08 pc) and defines a thin shell of two grid‑cell thickness as the envelope. Second, the friend‑of‑friend algorithm “clumpfind” extracts connected high‑density regions; within each clump the top two‑thirds of the density distribution are labeled as core, the bottom third as envelope. Both methods yield comparable statistical ensembles (≈ 100 clumps per snapshot in 3‑D, ≈ 40 in projected 2‑D), and the authors demonstrate that the main results are insensitive to the exact definition.

The mass‑to‑flux ratio R is computed using two distinct approaches. (i) The “average‑field” method calculates a mass‑weighted mean line‑of‑sight magnetic field B_LOS for the core and for the envelope, then forms R = (Σ/B_LOS)_core / (Σ/B_LOS)_envelope, where Σ is the column density integrated along the line of sight. (ii) The “pixel‑wise” method evaluates R for every possible pair of core and envelope pixels, takes the absolute value, and averages logarithmically over the entire set, thereby capturing the broad, roughly log‑normal distribution of R values. Both techniques produce distributions centered near |R| ≈ 1, but with systematic trends that depend on the initial magnetic field strength.

Key findings are as follows. For strong initial fields (β₀ = 0.01, 0.1), the median |R| is essentially unity (|R| ≈ 1.0 ± 0.2) and shows little variation with time or resolution. For intermediate fields (β₀ = 1) the distribution widens slightly, with a median |R| ≈ 1.1. For weak fields (β₀ = 10, 100), the median |R| falls below unity (≈ 0.8–0.9), yet remains close to one. The authors also quantify magnetic‑field reversals along the line of sight by counting the fraction of cells with opposite sign B_LOS in the envelope; this reversal parameter X ranges from –1 (all negative) to +1 (all positive). In strong‑field runs, X is near zero in the envelope, indicating frequent cancellations, whereas in the core X stays close to ±1, reflecting a more coherent field.

The physical interpretation hinges on two mechanisms. (1) Field reversals: In highly magnetised turbulence, the envelope experiences strong shear and vorticity, causing the line‑of‑sight component of the magnetic field to flip sign on scales comparable to the envelope thickness. When averaging, positive and negative contributions cancel, reducing the effective magnetic flux in the envelope and driving |R| toward unity despite the underlying strong field. The core, being denser and typically sub‑sonic, retains a coherent field, so its flux is not similarly diminished. (2) Small‑scale turbulent dynamo: In weak‑field regimes, the turbulent cascade efficiently amplifies magnetic energy on scales comparable to or smaller than the core size. The dynamo can increase the core magnetic flux by factors of 10–100 relative to the envelope, thereby lowering the core’s mass‑to‑flux ratio and producing |R| < 1. The simulations show that this dynamo action is active even without explicit non‑ideal MHD terms, consistent with previous studies (e.g., Sur et al. 2010; Federrath et al. 2011).

The authors compare their synthetic results with observational studies that measured R via Zeeman splitting (Crutcher et al. 2009; Lunttila et al. 2008). Those observations typically report |R| ≈ 1, which has been interpreted as evidence for turbulence‑regulated star formation. The present work demonstrates that the same observational signature can arise from a combination of field reversals (in strong‑field clouds) and dynamo‑enhanced flux (in weak‑field clouds). Consequently, a simple threshold of |R| > 1 versus |R| < 1 is insufficient to discriminate between ambipolar diffusion and turbulence as the dominant star‑formation regulator. Additional diagnostics—such as direct measurements of field topology, velocity dispersion profiles, or independent estimates of the turbulent power spectrum—are required.

In summary, the paper provides a thorough statistical analysis of the core‑to‑envelope mass‑to‑flux ratio in supersonic, driven MHD turbulence. It shows that, across a wide range of initial magnetic field strengths and resolutions, the average |R| remains close to unity because (i) magnetic‑field reversals in the envelope suppress the envelope flux in strong‑field cases, and (ii) small‑scale turbulent dynamo amplification raises the core flux in weak‑field cases. These findings caution against using R alone as a decisive test of ambipolar diffusion versus turbulence‑regulated star formation, and they highlight the need for more nuanced observational strategies that can disentangle the underlying magnetic‑field dynamics.