pith. sign in

arxiv: 2604.11485 · v3 · submitted 2026-04-13 · 🌌 astro-ph.GA · astro-ph.SR

Scale-dependent surface and volume density properties of filaments in molecular clouds

Guo-Yin Zhang , Alexander Men'shchikov , Jin-Zeng Li This is my paper

Pith reviewed 2026-05-10 15:41 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.SR
keywords filamentsmolecular cloudssurface density profilesvolume density profilesscale dependencestar formationdensity contrasts
0
0 comments X

The pith

Filament widths in molecular clouds grow with spatial scale following power laws, rather than staying fixed at 0.1 pc.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies a multiscale extraction method to map filaments across seven nearby molecular clouds and derives both surface density and volume density profiles for a large sample. It reports that half-maximum widths increase systematically with spatial scale, following power laws of roughly scale to the 0.5 for surface density and scale to the 0.37 for volume density. This scale dependence directly contradicts the common assumption of a universal filament width and shows volume density slopes shallower than those expected for an isothermal cylinder in hydrostatic equilibrium. The analysis also finds much higher density contrasts in three dimensions than in projection and notes that measured widths depend on resolution and distance.

Core claim

The half-maximum widths of surface density profiles scale proportionally to Y^0.50 and volume density profiles to Y^0.37 across scales from 0.01 to 1 pc. Median volume density slopes lie between 2.1 and 2.4, below the value of 4 predicted for hydrostatic equilibrium. Volume density contrasts reach 17-52 while surface density contrasts remain 1.1-2.7, and shallow-profile widths in volume density fall far below those in surface density.

What carries the argument

The getsf multiscale extraction method, which isolates filaments at successive spatial scales Y to obtain surface density profiles Σ(r) and the corresponding volume density profiles ρ(r).

If this is right

  • Filaments lack a single characteristic width and instead show widths that grow with the spatial scale of the structure being measured.
  • Volume density slopes shallower than 4 indicate that filaments are not generally in simple hydrostatic equilibrium as isothermal cylinders.
  • Surface density widths can overestimate the true physical extent of filaments by one to two orders of magnitude when profiles are shallow.
  • Comparative studies of filaments across clouds must correct for angular resolution and distance to avoid systematic bias in widths and slopes.
  • Filaments are far more prominent as three-dimensional structures than their projected surface density appearance suggests.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Star formation theories that rely on a fixed filament width for fragmentation and core formation may need revision to incorporate scale-dependent widths.
  • Improved three-dimensional density reconstruction techniques could be tested against the reported contrast differences between surface and volume density.
  • The power-law scaling may link filament properties to the larger-scale turbulence or gravitational collapse that sets the cloud structure at each Y.

Load-bearing premise

The multiscale extraction separates genuine scale-dependent filament properties without introducing artifacts, and volume density profiles are recovered from surface density profiles without major unaccounted projection or geometric effects.

What would settle it

Re-observe the same set of filaments in one cloud at two or more different distances or angular resolutions and check whether the measured widths follow the reported power-law relations with scale.

Figures

Figures reproduced from arXiv: 2604.11485 by Alexander Men'shchikov, Guo-Yin Zhang, Jin-Zeng Li.

Figure 1
Figure 1. Figure 1: Zoomed-in surface density image of the Taurus molecular cloud, overlaid with scale-dependent skeletons traced by getsf on scales of 14–216′′. Widths of the colored skeletons are proportional to the scale (specified in the color bar) on which the filaments are detected. Darker skeleton colors indicate filament segments with acceptably good radial profiles (Sect. 3.4). Two dashed lines indicate the filament … view at source ↗
Figure 2
Figure 2. Figure 2: Scale-dependent distributions of surface and volume density contrasts CΣk (blue) and Cρk (red) of filaments in the molecular clouds, obtained from profiles Σ(r) and ρ(r), respectively (Sect. 3.5). Dashed lines indicate median values C˜ Σk and C˜ ρk , given in the panels as upper and lower values, respectively. Rightmost panels show distributions accumulated over all scales. The resulting filament profiles … view at source ↗
Figure 3
Figure 3. Figure 3: Scale-dependent distributions of the half-maximum widths Hk (blue) and hk (red) of filaments in the molecular clouds, derived for the surface and volume density profiles Σ(r) and ρ(r), respectively (Sect. 3.5). Also shown are the distributions of the deconvolved widths H˘ k (green), estimated from Eq. (4) for only those filaments that can be deconvolved with an accuracy better than 20%. Dot-dashed lines in… view at source ↗
Figure 4
Figure 4. Figure 4: Scale-dependent distributions of surface and volume density slopes γk (blue) and βk (red) of filaments in the molecular clouds, obtained from profiles Σ(r) and ρ(r), respectively (Sect. 3.5). Dashed lines indicate median values ˜γk and β˜ k , given in the panels as upper and lower values, respectively. Only profiles passing the quality criteria (Sect. 3.4) are included. Rightmost panels show distributions … view at source ↗
Figure 5
Figure 5. Figure 5: Scale-dependent distributions of linear densities ΛΣk (blue) and Λρk (red) of filaments in the molecular clouds, integrated from profiles Σ(r) and ρ(r), respectively (Sect. 3.5). The two linear density measures are practically identical (within 3%). Dashed lines in the seven upper rows indicate median values Λ˜ Σk and Λ˜ ρk , listed as the upper and lower values, respectively. Green dashed lines indicate t… view at source ↗
Figure 6
Figure 6. Figure 6: Median half-maximum widths Hk , profile slopes γk , and linear densities Λk of scale-dependent filaments in the molecular clouds. Two-letter name abbreviations are shown above the upper axes at the cloud distances. To avoid overlaps, results for Taurus and Ophiuchus are shifted by ±5 pc, respectively. The plotted quantities are computed from profiles meeting the quality criteria for filament segments on sp… view at source ↗
Figure 7
Figure 7. Figure 7: Fractions of filament segments with supercritical linear densities as a function of the spatial scale. The three thresholds Λc/2, Λc , and 2Λc are indicated in the panels. Colored lines represent the seven molecular clouds studied in this work. Fractions are computed for filament segments that meet the quality criteria (Sect. 3.5) and have well-determined γ values. 5.2. Scale dependence of filament propert… view at source ↗
Figure 8
Figure 8. Figure 8: Convergence tests for the molecular clouds. Plots show measured half-maximum widths Hk , profile slopes γk , and linear densities Λk of filaments as a function of cloud distance D (lower axes) or angular resolution O (upper axes). Colored lines connect median values at resolutions of 13.5, 27, 38, 54, 76, 108, and 216′′, with thick dashed lines representing medians over all scales (14–216′′). Gray-scale ba… view at source ↗
read the original abstract

We present a systematic analysis of scale-dependent properties of filamentary structures in seven nearby molecular clouds $-$ Taurus, Ophiuchus, Perseus, Orion A, California, IC 5146, Vela C $-$ using the multiscale extraction method $getsf$. Alongside the usual surface density profiles $\Sigma(r)$, we derived volume density profiles $\rho(r)$ for a large sample of filaments, providing new observational constraints on their three-dimensional structure. The half-maximum widths $H$ and $h$ of the surface and volume density profiles, respectively, systematically increase with the spatial scale, following power laws $\tilde{H} \propto Y^{0.50}$ and $\tilde{h} \propto Y^{0.37}$, with distributions spanning $\sim 0.01 - 1 $ pc across all scales, challenging the notion of a universal filament width of $\sim 0.1$ pc. The median volume density slopes $\tilde{\beta} \approx 2.1 - 2.4$ are systematically lower than the value $\beta = 4$ expected for an isothermal cylinder in hydrostatic equilibrium. For shallow profiles with $\beta \lesssim 1$, the volume density width $h$ falls below the surface density width $H$ by one to two orders of magnitude, demonstrating that surface density widths overestimate the true physical extent of filaments with shallow profiles. Volume density contrasts are substantially higher than surface density contrasts ($\tilde{C}_{\rho} \approx 17 - 52$ versus $\tilde{C}_{\Sigma} \approx 1.1 - 2.7$), confirming that filaments are substantially more prominent in three dimensions than their projected appearance suggests. Measured filament widths and slopes systematically depend on the angular resolution and distance, highlighting the importance of accounting for resolution bias in comparative filament studies.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 3 minor

Summary. The paper analyzes filamentary structures across seven nearby molecular clouds (Taurus, Ophiuchus, Perseus, Orion A, California, IC 5146, Vela C) using the getsf multiscale extraction method. It derives both surface density profiles Σ(r) and volume density profiles ρ(r), reporting that half-maximum widths systematically increase with extraction scale Y according to power laws H̃ ∝ Y^{0.50} and h̃ ∝ Y^{0.37} (spanning ~0.01-1 pc), that median volume density slopes are β̃ ≈ 2.1-2.4 (lower than the β=4 expected for isothermal hydrostatic cylinders), that volume density contrasts are much higher than surface density contrasts, and that widths and slopes depend on angular resolution and distance.

Significance. If the reported scale dependence and shallow β values are intrinsic rather than methodological artifacts, the results would provide important observational constraints on the three-dimensional structure of filaments, challenging assumptions of a universal ~0.1 pc width and simple hydrostatic equilibrium models. The multi-cloud sample and derivation of volume densities from surface densities add value for star formation studies, though the significance hinges on demonstrating that getsf does not induce the observed correlations.

major comments (3)
  1. [§4] §4 (results on scale dependence): The headline power-law relations H̃ ∝ Y^{0.50} and h̃ ∝ Y^{0.37} rest on filaments extracted by getsf at scale Y; because the algorithm decomposes the map into a hierarchy of scales and assigns structures accordingly, any correlation between the assigned Y and the fitted half-maximum widths could arise from the decomposition procedure itself rather than the underlying density field. A quantitative test (e.g., injecting fixed-width filaments into synthetic maps and re-extracting) is required to establish that the trend is physical.
  2. [§3.3] §3.3 (volume density derivation): The conversion from observed Σ(r) to ρ(r) assumes cylindrical geometry and a Plummer-like form with free parameter β. For the reported shallow median slopes β̃ ≈ 2.1-2.4, even modest inclination or line-of-sight projection effects can shift the recovered β by ~0.5, directly affecting the claimed systematic deviation from the isothermal-cylinder value β=4; the manuscript does not quantify these uncertainties or provide independent verification of the cylindrical assumption for shallow profiles.
  3. [§5] §5 (discussion of resolution/distance dependence): While the text notes that measured widths and slopes depend on angular resolution and distance, the quantification of this bias (e.g., via resolution-matched subsamples or distance-binned statistics) is insufficient to assess whether the reported power-law indices remain robust when these effects are controlled, undermining cross-cloud comparisons.
minor comments (3)
  1. [Abstract] The notation H̃, h̃, β̃, C̃_ρ, and C̃_Σ is used without an explicit definition in the abstract and early sections; a clear statement that the tilde denotes median or characteristic values would improve readability.
  2. [Figures] Figures showing width versus scale (e.g., the panels presenting the power-law fits) would be clearer if they overlaid the individual filament measurements with uncertainties rather than only the binned medians and fits.
  3. [§3] The sample sizes per cloud and per scale bin, together with the precise fitting procedure for the power-law indices (including any weighting or outlier rejection), should be stated explicitly in the methods or results tables.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed review. We address each major comment below and have revised the manuscript to incorporate additional tests and quantifications where the concerns are valid.

read point-by-point responses

  1. Referee: [§4] The headline power-law relations H̃ ∝ Y^{0.50} and h̃ ∝ Y^{0.37} rest on filaments extracted by getsf at scale Y; because the algorithm decomposes the map into a hierarchy of scales and assigns structures accordingly, any correlation between the assigned Y and the fitted half-maximum widths could arise from the decomposition procedure itself rather than the underlying density field. A quantitative test (e.g., injecting fixed-width filaments into synthetic maps and re-extracting) is required to establish that the trend is physical.

    Authors: We agree that validation against algorithmic artifacts is essential given the multiscale decomposition in getsf. Although the trends appear consistently across seven clouds, we have added a new subsection in §4 presenting synthetic tests: fixed-width filaments (0.05 pc and 0.2 pc) were injected into realistic noise and background maps, then re-extracted with getsf. The recovered widths exhibit no artificial power-law dependence on Y, supporting that the observed relations are intrinsic. Corresponding figures and discussion have been included in the revision. revision: yes

  2. Referee: [§3.3] The conversion from observed Σ(r) to ρ(r) assumes cylindrical geometry and a Plummer-like form with free parameter β. For the reported shallow median slopes β̃ ≈ 2.1-2.4, even modest inclination or line-of-sight projection effects can shift the recovered β by ~0.5, directly affecting the claimed systematic deviation from the isothermal-cylinder value β=4; the manuscript does not quantify these uncertainties or provide independent verification of the cylindrical assumption for shallow profiles.

    Authors: The sensitivity of shallow β to inclination and projection is a valid concern. We have revised §3.3 to include Monte Carlo simulations of inclined Plummer cylinders (β=2.2) showing that for inclinations ≤45°, the recovered β increases by at most 0.4, preserving median values well below 3. The cylindrical assumption follows standard practice in filament studies and is consistent with the high aspect ratios in our sample; we have added explicit caveats noting that full independent 3D verification lies beyond the current dataset. revision: yes

  3. Referee: [§5] While the text notes that measured widths and slopes depend on angular resolution and distance, the quantification of this bias (e.g., via resolution-matched subsamples or distance-binned statistics) is insufficient to assess whether the reported power-law indices remain robust when these effects are controlled, undermining cross-cloud comparisons.

    Authors: We concur that stronger controls are required. The revised §5 now includes resolution-matched subsamples (restricted to comparable angular resolutions) and distance-binned statistics. These confirm that the power-law indices for H and h remain stable within ~10% after controlling for resolution and distance, although absolute widths show the expected resolution scaling. Updated text and supplementary figures document these checks. revision: yes

Circularity Check

0 steps flagged

No significant circularity in empirical filament analysis

full rationale

The paper applies the getsf multiscale extraction method to observational surface density maps of molecular clouds, extracts filaments at discrete spatial scales Y, measures their half-maximum widths H and h from fitted profiles, converts to volume density profiles ρ(r) under cylindrical assumptions, and reports empirical power-law fits and median slopes. These steps consist of direct measurements and standard statistical summaries on the data; no equation or procedure reduces a claimed result to a quantity defined by the paper's own fitted parameters, self-citations, or ansatz. The reported relations (H̃ ∝ Y^0.50, h̃ ∝ Y^0.37, β̃ ≈ 2.1–2.4) are therefore independent outputs of the analysis rather than tautological restatements of inputs.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

Central claims rest on the reliability of the getsf extraction algorithm and the validity of inferring 3D volume densities from 2D surface-density maps; two power-law exponents are reported as measured fits to the data.

free parameters (2)
  • surface-density width power-law index = 0.50
    Exponent 0.50 obtained by fitting observed half-maximum widths across spatial scales in the sample.
  • volume-density width power-law index = 0.37
    Exponent 0.37 obtained by fitting observed volume half-maximum widths across spatial scales.
axioms (2)
  • domain assumption The getsf multiscale extraction method accurately identifies filaments at each spatial scale without introducing artificial scale dependence.
    Invoked for all width and slope measurements in the seven clouds.
  • domain assumption Volume density profiles can be reliably recovered from projected surface density profiles under standard geometric assumptions.
    Required to derive ρ(r) and the reported β slopes and contrasts.

pith-pipeline@v0.9.0 · 5643 in / 1646 out tokens · 39425 ms · 2026-05-10T15:41:39.438547+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. A Salpeter-like filament linear density function across nearby molecular clouds

    astro-ph.GA 2026-04 unverdicted novelty 6.0

    Filament linear density functions integrated over spatial scales in multiple molecular clouds follow a power law with slope ~1.3, matching the Salpeter IMF.

  2. A Salpeter-like filament linear density function across nearby molecular clouds

    astro-ph.GA 2026-04 unverdicted novelty 5.0

    Across seven molecular clouds the integrated filament linear density function follows a power law with slope 1.30-1.34, matching the Salpeter IMF slope of 1.35.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    Alves, J., Lombardi, M., & Lada, C. J. 2007, A&A, 462, L17 André, P., Di Francesco, J., Ward-Thompson, D., et al. 2014, Protostars and Plan- ets VI, 27 André, P., Men’shchikov, A., Bontemps, S., et al. 2010, A&A, 518, L102 André, P. J., Palmeirim, P., & Arzoumanian, D. 2022, A&A, 667, L1 Arzoumanian, D., André, P., Didelon, P., et al. 2011, A&A, 529, L6 A...

    work page internal anchor Pith review Pith/arXiv arXiv 2007

  2. [2]

    C.1.Surface density profiles along two parallel cuts across the Taurus filament (see Fig

    filament 14 27 38 54 76 108 153 216 0 150 300 450 600 1021 1022 filament + background 0 200 400 600 0 0.1 0.2 0.3 0.4 0 0.1 0.2 0.3 0.4 Distance along cut (arcsec) Distance along cut (pc) Fig. C.1.Surface density profiles along two parallel cuts across the Taurus filament (see Fig. 1).Top panels:Observed filament profiles (solid curves) and scale-dependen...

    work page 2026