Extending velocity-centroid statistics in position–position–velocity space beyond the Gaussian, weak-absorption regime: intermittency, density–velocity correlations, non-linear self-absorption, and non-degenerate inversion for MHD mode content and magnetization

The analytic theory of velocity centroid statistics in position–position–velocity (PPV) space derived by Kandel, Lazarian & Pogosyan (2017) predicts when centroids recover the underlying three-dimensional velocity spectrum, how self-absorption and absorption lines modify the scalings, and how the anisotropy multipoles of the centroid structure function decompose into Alfvén, slow, and fast magnetohydrodynamic modes — but it is valid only under Gaussian statistics, negligible density–velocity correlations ($\sigma_\rho/\rho_0\lesssim1$), a perturbative (linear) treatment of self-absorption, and marginal coupling between the mode cascades. The theory's own stated limitations, left open by the subsequent literature, define the problem: quantify and correct the effects of intermittent, non-Gaussian statistics at high sonic Mach number; incorporate the density–velocity cross terms that dominate in supersonic turbulence; treat the non-linear strong-self-absorption regime in which the linear expansion fails; and, on that basis, construct a demonstrably non-degenerate inversion of observed centroid statistics into mode fractions, Alfvén Mach number, and line-of-sight magnetic-field angle. Later work has numerically tested the mode-decomposition and opacity-robustness predictions (Hernández-Padilla et al. 2020, 2024), supplied partial numerical workarounds for supersonic density contamination (Xu & Hu 2021; Yuen, Ho & Lazarian 2021), and resolved the inclination–magnetization degeneracy for channel-intensity anisotropy (Hu, Lazarian & Xu 2021), but the analytic extensions and the centroid-based mode-fraction inversion remain open.

Background

Turbulence in the ionized and molecular interstellar medium (ISM) is studied through spectroscopic line observations. A mapping of an emission line such as H I 21 cm or a rotational CO transition delivers a data cube of intensity indexed by two sky coordinates and the line-of-sight (LOS) velocity inferred from the Doppler shift — position–position–velocity (PPV) space. Because Doppler shifts probe only the LOS velocity component and radiative transfer (self-absorption by intervening gas along the same sightline, plus thermal broadening) reshapes the line profile, recovering the true three-dimensional (3D) density and velocity statistics from a PPV cube is a non-trivial inference problem.

A classical statistic for this purpose is the velocity centroid: the intensity-weighted mean LOS velocity at each sky position. Its unnormalized form (UVC, the first spectral moment without division by the integrated intensity) is analytically tractable, and numerical studies showed the normalization does not matter for recovering velocity statistics; the modified velocity centroid (MVC) additionally subtracts a term proportional to the column-density structure function to suppress density contamination. The basic measure is the second-order structure function — the variance of the difference of centroid values at two sky positions separated by a lag vector — whose dependence on the lag magnitude recovers the velocity spectrum: for optically thin, subsonic turbulence the UVC structure function scales as $R^{1+\nu}$, where $\nu$ is the exponent of the velocity structure function (Kolmogorov $\nu=2/3$). The analytic backbone for all such statistics is the PPV fluctuation theory of Lazarian & Pogosyan: velocity channel analysis (VCA; statistics of intensity fluctuations in velocity slices of finite width), its extension to absorbing media, and the velocity coordinate spectrum (VCS; statistics along the velocity axis). Kandel, Lazarian & Pogosyan (2016) extended the VCA to anisotropy studies, and Kandel, Lazarian & Pogosyan (2017) — the source of this problem — applied the same machinery to velocity centroids.

The 2017 centroid theory provides: regime maps for when centroids recover 3D velocity spectra (reliable mainly for subsonic turbulence, with an absorption-induced 'universal' regime in which spectral information is lost, mirroring the VCA); a perturbative treatment of self-absorption valid only when $\alpha^2[d_s(R,v)-d_s(0,v)]\ll1$, where $\alpha$ is an absorption parameter and $d_s$ the PPV density dispersion (the practical criterion is the source's eq. (23)); restricted velocity centroids (RVCs) built from line wings for absorption-line studies; and a decomposition of centroid anisotropy into contributions from the Alfvén, slow, and fast modes of compressible magnetohydrodynamic (MHD) turbulence. Mode decomposition rests on the Goldreich–Sridhar picture of anisotropic Alfvénic cascades and on the procedure of Cho & Lazarian (2002, 2003), which is accurate when coupling between the cascades is moderate. The anisotropy is encoded in the angular multipole moments $D_n(R)$ of the centroid structure function (moments of its dependence on the position angle of the lag), with a one-to-one correspondence between real-space and Fourier-space multipoles that enables interferometric measurement. The mode-dependent multipoles were proposed as diagnostics of the medium's magnetization — the Alfvén Mach number $M_{\rm A}$ (ratio of turbulent to Alfvénic velocity) — and of the angle between the mean magnetic field and the LOS.

The theory is explicitly built on assumptions: Gaussian velocity statistics (and a Gaussian approximation for the PPV density when treating absorption); uncorrelated density and velocity, justified only when the relative density dispersion $\sigma_\rho/\rho_0\lesssim1$; the linear self-absorption expansion; and marginal inter-mode coupling. Numerical work before and after — Esquivel & Lazarian (2005) and Esquivel, Lazarian, Horibe, Cho, Ossenkopf & Stutzki (2007) — established that density–velocity correlations and non-Gaussian density statistics distort centroid statistics once $\sigma_\rho/\rho_0\gtrsim1$, typical of supersonic turbulence where the sonic Mach number $M_{\rm s}$ (flow velocity relative to the sound speed) exceeds unity and the density spectrum is shallow, and where intermittency (non-Gaussian, scale-dependent fluctuations concentrated in rare structures) becomes prominent. What does not exist is the corresponding quantitative theory: how intermittent statistics modify the centroid scalings, multipoles, and regime boundaries; how to include the density–velocity cross terms; what happens in the non-linear strong-absorption regime; and how to invert observed centroid statistics into mode fractions, $M_{\rm A}$, and LOS angle without degeneracy.

Problem Statement

Extend the analytic theory of velocity centroid statistics in PPV space beyond the validity assumptions under which it was derived, and make the extended statistics invertible for the physical parameters of the turbulence. The components below share this single objective.

(1) Intermittency and non-Gaussianity: quantify how non-Gaussian, intermittent velocity and density statistics — prominent at high sonic Mach number — modify the centroid (and modified-centroid) structure-function scalings, the anisotropy multipoles $D_n(R)$, and the boundaries between the regimes of the 2017 theory (e.g. between the usual $R^{1+\nu}$ scaling, the absorption-induced universal regime, and the non-linear regime), either by analytic extension of the PPV fluctuation theory or by corrections calibrated on and validated against mode-decomposed MHD simulations.

(2) Density–velocity correlations: incorporate the density–velocity cross terms into the centroid formulae for $\sigma_\rho/\rho_0\gtrsim1$ and shallow density spectra, specifying under what conditions velocity spectra and anisotropies remain recoverable from centroids and what corrections or data-selection strategies (beyond the existing low-column-density sampling workaround) are required, with their domain of validity.

(3) Non-linear strong self-absorption: treat the regime in which the linear-expansion condition behind eq. (23) of the source fails — derive, or numerically establish with controlled simulations, the behavior of the centroid structure function and its multipoles there, in particular whether power-law asymptotes survive and with what exponents and regime boundaries.

(4) Non-degenerate inversion: on the basis of (1)–(3), construct and validate a robust inversion of observed centroid statistics (structure functions, multipole ratios, possibly in synergy with velocity-channel and velocity-coordinate data) into the Alfvén/slow/fast mode fractions, the Alfvén Mach number $M_{\rm A}$, and the LOS angle of the magnetic field — demonstrably breaking the confirmed degeneracies (Alfvén versus slow modes at small viewing angles; magnetization versus viewing angle) — and demonstrate it on synthetic observations with known ground truth.

All four components are required for this single objective: components (1)–(3) supply the extension of the theory beyond its stated assumptions, and component (4) renders those extended statistics invertible, so a submission that resolves only a proper subset of the components — for example an inversion resting on the unextended theory, or any single one of (1)–(3) without the others — is partial progress rather than resolution of the objective. Where a component permits either an analytic derivation or corrections calibrated on and validated against simulations, either route suffices to satisfy that component. The fourth validity assumption of the theory — marginal coupling between the mode cascades — is retained rather than relaxed: relaxing it is not one of the required components, and the mode-fraction inversion of component (4) need only be validated where the Alfvén/slow/fast decomposition it rests on is accurate.

Instrumental optimization (tracer choice, channel width, resolution needed to access a reliable inertial range) enters as the operational context of component (4), not as an independent deliverable.

References

1. Kandel, D., Lazarian, A., & Pogosyan, D. (2017). Study of velocity centroids based on the theory of fluctuations in position–position–velocity space. Monthly Notices of the Royal Astronomical Society, 464, 3617–3635. DOI: 10.1093/mnras/stw2512. https://doi.org/10.1093/mnras/stw2512 (arXiv:1607.04316, https://arxiv.org/abs/1607.04316)
2. Lazarian, A., & Pogosyan, D. (2000). Velocity Modification of H I Power Spectrum. The Astrophysical Journal, 537, 720–748. DOI: 10.1086/309040. https://doi.org/10.1086/309040
3. Lazarian, A., & Pogosyan, D. (2004). Velocity Modification of the Power Spectrum from an Absorbing Medium. The Astrophysical Journal, 616, 943–965. DOI: 10.1086/422462. https://doi.org/10.1086/422462
4. Lazarian, A., & Pogosyan, D. (2006). Studying Turbulence Using Doppler-broadened Lines: Velocity Coordinate Spectrum. The Astrophysical Journal, 652, 1348–1365. DOI: 10.1086/508012. https://doi.org/10.1086/508012
5. Lazarian, A., & Esquivel, A. (2003). Statistics of Velocity from Spectral Data: Modified Velocity Centroids. The Astrophysical Journal, 592, L37–L40. DOI: 10.1086/377427. https://doi.org/10.1086/377427
6. Esquivel, A., & Lazarian, A. (2005). Velocity Centroids as Tracers of the Turbulent Velocity Statistics. The Astrophysical Journal, 631, 320–350. DOI: 10.1086/432458. https://doi.org/10.1086/432458
7. Esquivel, A., Lazarian, A., Horibe, S., Cho, J., Ossenkopf, V., & Stutzki, J. (2007). Statistics of velocity centroids: effects of density-velocity correlations and non-Gaussianity. Monthly Notices of the Royal Astronomical Society, 381, 1733–1744. DOI: 10.1111/j.1365-2966.2007.12359.x. https://doi.org/10.1111/j.1365-2966.2007.12359.x
8. Goldreich, P., & Sridhar, S. (1995). Toward a theory of interstellar turbulence. 2: Strong Alfvénic turbulence. The Astrophysical Journal, 438, 763. DOI: 10.1086/175121. https://doi.org/10.1086/175121
9. Cho, J., & Lazarian, A. (2002). Compressible sub-Alfvénic MHD turbulence in low-beta plasmas. Physical Review Letters, 88, 245001. DOI: 10.1103/PhysRevLett.88.245001. https://doi.org/10.1103/PhysRevLett.88.245001 (arXiv:astro-ph/0205282, https://arxiv.org/abs/astro-ph/0205282)
10. Cho, J., & Lazarian, A. (2003). Compressible magnetohydrodynamic turbulence: mode coupling, scaling relations, anisotropy, viscosity-damped regime and astrophysical implications. Monthly Notices of the Royal Astronomical Society, 345, 325–339. DOI: 10.1046/j.1365-8711.2003.06941.x. https://doi.org/10.1046/j.1365-8711.2003.06941.x (arXiv:astro-ph/0301062, https://arxiv.org/abs/astro-ph/0301062)
11. Kandel, D., Lazarian, A., & Pogosyan, D. (2016). Extending velocity channel analysis for studying turbulence anisotropies. Monthly Notices of the Royal Astronomical Society, 461, 1227–1259. DOI: 10.1093/mnras/stw1296. https://doi.org/10.1093/mnras/stw1296 (arXiv:1604.06102, https://arxiv.org/abs/1604.06102)
12. Hernández-Padilla, D., Esquivel, A., Lazarian, A., Pogosyan, D., Kandel, D., & Cho, J. (2020). Anisotropy of Velocity Centroids and the Signature of Different MHD Modes in the Turbulent ISM. The Astrophysical Journal, 901, 11. DOI: 10.3847/1538-4357/abad9e. https://doi.org/10.3847/1538-4357/abad9e (arXiv:2008.05393, https://arxiv.org/abs/2008.05393)
13. Xu, S., & Hu, Y. (2021). Measuring Magnetization with Rotation Measures and Velocity Centroids in Supersonic MHD Turbulence. The Astrophysical Journal, 910, 88. DOI: 10.3847/1538-4357/abe403. https://doi.org/10.3847/1538-4357/abe403 (arXiv:2102.05647, https://arxiv.org/abs/2102.05647)
14. Hu, Y., Lazarian, A., & Xu, S. (2021). Anisotropic Turbulence in Position–Position–Velocity Space: Probing Three-dimensional Magnetic Fields. The Astrophysical Journal, 915, 67. DOI: 10.3847/1538-4357/ac00ab. https://doi.org/10.3847/1538-4357/ac00ab (arXiv:2104.02842, https://arxiv.org/abs/2104.02842)
15. Yuen, K. H., Ho, K. W., & Lazarian, A. (2021). Technique for Separating Velocity and Density Contributions in Spectroscopic Data and Its Application to Studying Turbulence and Magnetic Fields. The Astrophysical Journal, 910, 161. DOI: 10.3847/1538-4357/abe4d4. https://doi.org/10.3847/1538-4357/abe4d4 (arXiv:2012.15776, https://arxiv.org/abs/2012.15776)
16. Hernández-Padilla, D., Esquivel, A., Lazarian, A., Velázquez, P. F., & Cho, J. (2024). Effects of Radiative Transfer on the Observed Anisotropy in Magnetohydrodynamic Turbulent Molecular Simulations. The Astrophysical Journal, 972, 93. DOI: 10.3847/1538-4357/ad57c8. https://doi.org/10.3847/1538-4357/ad57c8 (arXiv:2406.13926, https://arxiv.org/abs/2406.13926)
17. Lazarian, A., Yuen, K. H., & Pogosyan, D. (2024). Gradient Technique Theory: Tracing Magnetic Field and Obtaining Magnetic Field Strength. The Astrophysical Journal, 974, 237. DOI: 10.3847/1538-4357/ad6d62. https://doi.org/10.3847/1538-4357/ad6d62 (arXiv:2406.01545, https://arxiv.org/abs/2406.01545)
18. Hu, Y. (2025). Characterizing 3D Magnetic Fields and Turbulence in H I Clouds. arXiv:2505.07422. https://arxiv.org/abs/2505.07422