pith. sign in

arxiv: 2604.12570 · v2 · submitted 2026-04-14 · 🌌 astro-ph.GA · astro-ph.IM

Deriving volume density profiles of filaments from observed surface densities

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

Pith reviewed 2026-05-15 07:01 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.IM
keywords filamentsmolecular cloudssurface density profilesvolume density profilesPlummer functionfitting methodsstar formation
0
0 comments X

The pith

Observed surface density profiles of filaments cannot be translated to volume densities using the simple Plummer assumptions that have been standard.

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

The paper introduces a fitting method for filament surface density profiles that explicitly incorporates finite cylindrical geometry rather than assuming infinite cylinders. Traditional Plummer fits assume the volume density slope is exactly one unit steeper than the observed surface slope and that the two widths are nearly identical, but these relations break down for shallow profiles. The new approach derives empirical relations that link the observed parameters to the true volume density slope and width, showing large deviations especially when profiles are shallow or the filament is compact or extended. Accurate recovery demands the observed width to be at least ten times the beam size; below that threshold slopes are badly overestimated and the filament appears unresolved. Because current far-infrared data usually fall short of this resolvedness, previously reported filament properties may need revision.

Core claim

Fitting observed surface density profiles Σ(r) with slope γ and width H by the Plummer function under the assumptions β = γ + 1 and h ≈ H produces systematic errors that grow for shallow slopes; a new method that accounts for finite cylindrical geometry supplies self-consistent empirical relations between the observed parameters and the volume density parameters β and h, with the slope difference δ = β − γ falling below unity for shallow compact cases and the width ratio h/H differing by more than an order of magnitude for extended shallow filaments.

What carries the argument

Finite-cylindrical-geometry fitting method that derives empirical relations between the parameters of the projected surface density profile Σ(r) and the intrinsic volume density profile ρ(r).

If this is right

  • The slope difference δ ≡ β − γ falls below unity for shallow (β ≲ 2) and compact profiles.
  • Widths h and H can differ by over an order of magnitude for extended filaments with shallow slopes.
  • Accurate parameter recovery requires resolvedness R ≡ H/O ≳ 10; lower values cause severe slope overestimation.
  • The traditional Plummer function yields systematically overestimated slopes.
  • Current far-infrared observations typically lack the resolution needed for reliable deconvolution.

Where Pith is reading between the lines

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

  • Many existing catalogs of filament slopes and widths derived from low-resolution data would shift if reprocessed with the geometry-aware relations.
  • The circular requirement of knowing true parameters to deconvolve observations can only be broken by pushing angular resolution well above the current typical values.
  • The same finite-cylinder correction could be tested on other elongated structures such as jets or galactic arms where projection effects are similarly important.

Load-bearing premise

The model profiles chosen to build the empirical relations accurately represent the true three-dimensional geometry and density structure of real molecular-cloud filaments.

What would settle it

Apply the new method to synthetic surface-density maps generated from known volume-density cylinders with β < 2 and compare the recovered β and h against the input values; systematic mismatch would falsify the claimed relations.

Figures

Figures reproduced from arXiv: 2604.12570 by Alexander Men'shchikov, Guo-Yin Zhang.

Figure 1
Figure 1. Figure 1: Analytical profiles of volume and surface densities (solid and dashed curves, respectively) for structures with a circular boundary at R = 103 au and selected values of ξ ≡ R/h (defined in Sect. 2.1) and β. Both profile types are computed from the analytical Eqs. (3) and (4) and are shown together to directly illustrate the key point that the slopes and widths of ρ(r) and Σ(r) differ systematically dependi… view at source ↗
Figure 2
Figure 2. Figure 2: Analytical dependencies of surface and volume density parameters (Eqs. (5)–(8)) on β for selected values of ξ. The dotted curve displays the relation γ = β − 1 adopted in previous studies, illustrating that this approximation becomes invalid for β ≲ 2. Our experience has shown that γ is the most sensitive fitting parameter; therefore, we devised a more elaborate scheme for es￾timating its initial value γ0.… view at source ↗
Figure 3
Figure 3. Figure 3: Accuracy of the analytical profiles and their parameters, as functions of βT and ξT. The left panels show the maximum and median relative deviations of the analytical surface density profiles from the numerical models. Relative errors in the main parameters h, β, w, and γ (middle and right), derived from Eqs. (5), (7), and (8), are defined with respect to the true values hT and βT used to compute the numer… view at source ↗
Figure 4
Figure 4. Figure 4: Accuracies of all parameters derived from fitting the numerical model profiles. The accuracy ratios with respect to the true model values are shown as functions of the true βT for fitted surface density profiles without noise (left) and with 10% noise added (right). Filled circles represent the parameters indicated at the edge of the right panels; their vertical spread reflects variations due to different … view at source ↗
Figure 5
Figure 5. Figure 5: Accuracy of the widths and slopes derived from fitting surface density profiles of a model filament convolved with increasingly larger beams. The accuracy ratios hh−1 T , ββ−1 T , and ww−1 T are represented by filled circles connected by solid lines as functions of the profile resolvedness R for different values of βT and ξT. Dashed lines show accuracies for modified profiles truncated at the 0.1ΣC level a… view at source ↗
Figure 6
Figure 6. Figure 6: Application of the fitting method to three filaments in the Cal￾ifornia molecular cloud. The surface density profiles ΣO(r) of selected filaments and the corresponding fitted profiles are shown with dots and solid curves, respectively. The derived β and ξ of the fits are specified, and the half-maximum widths H and h of the Σ(r) and ρ(r) profiles, respectively, are indicated next to the curves at the radia… view at source ↗
Figure 7
Figure 7. Figure 7: Simulation of the traditional approach of fitting median surface density profiles of entire filaments with the Plummer function. Dotted blue, green, and brown curves show 8 representative original model pro￾files with R = 0.5 pc, ξT = {1, 2}, and βT = {1, 2, 3, 6}, as well as the same profiles scaled to 0.7R and 1.5R. Overlapping gray curves show 400 randomly scaled original profiles in the range [0.7, 1.5… view at source ↗
Figure 8
Figure 8. Figure 8: Accuracies of traditional fitting with the Plummer function of Eq. (1) for numerical model profiles (cf [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Accuracies of traditional fitting with the Plummer function of Eq. (1) for surface density profiles of a model filament convolved with increasingly larger beams (cf [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Application of the traditional fitting method of Eq. (1) to the observed filament profiles shown in [PITH_FULL_IMAGE:figures/full_fig_p013_10.png] view at source ↗
read the original abstract

Accurate characterization of filamentary structures in star-forming clouds is essential for understanding star formation. Traditional methods fit observed surface density profiles $\Sigma(r)$ with slope $\gamma$ and width $H$ using the Plummer function, assuming $\beta=\gamma+1$ and $h\approx H$ for the volume density slope and width. These assumptions break down for shallow profiles, with the slope and width relations deviating progressively more for compact and extended filaments, respectively. We present a new fitting method that explicitly accounts for finite cylindrical geometry and establishes self-consistent empirical relationships between the parameters of $\Sigma(r)$ and those of the volume density profile $\rho(r)$ with slope $\beta$ and width $h$. The method was validated on model profiles and applied to selected filaments in the California molecular cloud. The slope difference $\delta\equiv\beta-\gamma$ falls below unity for shallow ($\beta\lesssim 2$) and compact profiles; $h$ and $H$ can differ by over an order of magnitude for extended filaments with shallow slopes. Accurate parameter recovery requires high resolvedness $R\equiv H/O\gtrsim 10$ (where $O$ is the beam width); at lower resolvedness, slopes are severely overestimated and filaments remain unresolved even when $H\gg O$. The traditional Plummer function yields systematically overestimated slopes. Accurate deconvolution requires a priori knowledge of the true parameters, creating a fundamental circular problem whose only robust solution is obtaining sufficiently high angular resolution. Current far-infrared observations typically lack sufficient resolution, and some previously reported filament properties may require reinterpretation.

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 / 2 minor

Summary. The paper claims to introduce a new fitting method that explicitly accounts for finite cylindrical geometry when deriving volume density profile parameters (slope β and width h) from observed surface density profiles (slope γ and width H). It establishes self-consistent empirical relationships between these parameters, validates the approach on families of model profiles, applies it to selected filaments in the California molecular cloud, and argues that traditional Plummer-function assumptions (β=γ+1, h≈H) lead to overestimated slopes, that accurate recovery requires resolvedness R≡H/O ≳10, and that deconvolution suffers from a fundamental circularity problem resolvable only by high angular resolution.

Significance. If the empirical Σ–ρ mappings prove robust beyond the idealized models, the work would meaningfully improve quantitative characterization of filament density structures, which are central to star-formation theory in molecular clouds. The explicit treatment of finite-cylinder projection effects and the resolution requirement are useful additions. The paper correctly flags the circularity in traditional deconvolution and the need for high-resolution data, both of which have direct observational implications.

major comments (3)
  1. [Abstract] Abstract and validation description: the method is stated to have been validated on model profiles, yet no quantitative recovery metrics (RMS errors, bias, or scatter in recovered β and h versus input values) are supplied. This absence is load-bearing because the central claim is that the new empirical relations enable accurate parameter recovery.
  2. [Method] Derivation of empirical relations: the self-consistent mappings between (γ,H) and (β,h) are obtained by fitting the same class of idealized finite-cylinder models that are later used for validation. This circularity, already flagged in the manuscript for traditional deconvolution, is not resolved by the new procedure and directly affects the reliability of the reported δ≡β−γ and h/H offsets.
  3. [Application to California Cloud] Application section: the reinterpretation of California Cloud filament properties rests on the assumption that real filaments match the constant-β, constant-h, perfectly cylindrical models used to calibrate the corrections. No test is shown for the impact of longitudinal density variations or non-circular cross-sections, which are known to exist in observed filaments and would propagate into biased β and h.
minor comments (2)
  1. The definition of resolvedness R≡H/O appears only in the abstract; it should be introduced with an explicit equation in the main text before the resolution-requirement discussion.
  2. Notation for the Plummer-function parameters (γ,H) versus the intrinsic (β,h) should be tabulated once for quick reference, as the many comparisons between them are central to the argument.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight important aspects of validation, model assumptions, and applicability that we address below. We have revised the manuscript to incorporate quantitative metrics and expanded discussion of limitations.

read point-by-point responses

  1. Referee: [Abstract] Abstract and validation description: the method is stated to have been validated on model profiles, yet no quantitative recovery metrics (RMS errors, bias, or scatter in recovered β and h versus input values) are supplied. This absence is load-bearing because the central claim is that the new empirical relations enable accurate parameter recovery.

    Authors: We agree that explicit quantitative metrics strengthen the validation. In the revised manuscript we have added a new table (Table 2) reporting RMS errors, mean biases, and standard deviations in recovered β and h as functions of input β, h, and resolvedness R. These metrics show that the empirical mappings recover input values with bias <0.05 and scatter <0.15 for R ≳ 10, while confirming the systematic overestimation under the traditional Plummer assumptions. The abstract has been updated to reference these results. revision: yes

  2. Referee: [Method] Derivation of empirical relations: the self-consistent mappings between (γ,H) and (β,h) are obtained by fitting the same class of idealized finite-cylinder models that are later used for validation. This circularity, already flagged in the manuscript for traditional deconvolution, is not resolved by the new procedure and directly affects the reliability of the reported δ≡β−γ and h/H offsets.

    Authors: The empirical relations are constructed by forward-projecting a grid of finite-cylinder models (varying β and h) to obtain synthetic Σ profiles, measuring γ and H from those profiles, and then fitting the inverse mappings. Validation then applies the fitted mappings to recover the original β and h from the same synthetic Σ profiles. While this uses the model family for both steps, it is not circular in the same sense as traditional deconvolution: the latter assumes a functional form to invert the projection, whereas here the mapping is an empirical calibration whose accuracy can be tested on the models themselves. We have added text in Section 3 clarifying this distinction and noting that the relations remain model-dependent; users should treat them as approximate corrections when real filaments depart from the assumed geometry. revision: partial

  3. Referee: [Application to California Cloud] Application section: the reinterpretation of California Cloud filament properties rests on the assumption that real filaments match the constant-β, constant-h, perfectly cylindrical models used to calibrate the corrections. No test is shown for the impact of longitudinal density variations or non-circular cross-sections, which are known to exist in observed filaments and would propagate into biased β and h.

    Authors: We acknowledge that the California Cloud application is an illustrative demonstration rather than a comprehensive re-analysis. Real filaments exhibit longitudinal variations and non-circular cross-sections that are not captured by our idealized models. In the revised manuscript we have expanded Section 5 to discuss these limitations explicitly, citing observational evidence for such complexities and noting that they would likely increase scatter in recovered parameters. We have also added a forward-looking statement that more realistic 3-D simulations will be needed to quantify the resulting biases. The original application is retained as a proof-of-concept showing how the new relations alter inferred slopes relative to Plummer fits. revision: yes

Circularity Check

1 steps flagged

Empirical Σ–ρ mappings calibrated and validated exclusively on the same family of idealized finite-cylinder models

specific steps
  1. fitted input called prediction [Abstract]
    "establishes self-consistent empirical relationships between the parameters of Σ(r) and those of the volume density profile ρ(r) with slope β and width h. The method was validated on model profiles and applied to selected filaments in the California molecular cloud."

    The relationships are obtained by fitting the very same class of model profiles (finite cylinders with fixed β, h) that are later used for validation; parameter recovery on those models is therefore forced by the calibration procedure rather than independently verified.

full rationale

The paper derives its central empirical relationships by generating and fitting families of model surface-density profiles that assume perfect finite cylinders with constant β and h, then recovers the input (β, h) values from the fitted (γ, H) parameters. Validation is performed on the identical model set, confirming the mapping by construction within the assumed geometry. While the paper correctly diagnoses circularity in traditional Plummer deconvolution, its own correction inherits the same model-class dependence; success on synthetic data does not independently test applicability to real filaments with longitudinal variations or non-circular cross-sections. This produces partial circularity in the load-bearing empirical step, though the underlying radiative-transfer modeling itself remains non-circular.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the assumption that filaments are well-approximated by finite cylinders and that synthetic profiles generated from power-law volume densities capture the essential projection effects; no new physical entities are introduced.

free parameters (1)
  • empirical mapping coefficients
    Coefficients in the self-consistent relations between (γ, H) and (β, h) are determined from model fits and therefore constitute fitted parameters.
axioms (1)
  • domain assumption Filaments can be modeled as finite-length cylinders with power-law volume density profiles
    Invoked to generate the model surface density profiles used for validation and to derive the empirical relations.

pith-pipeline@v0.9.0 · 5588 in / 1348 out tokens · 43636 ms · 2026-05-15T07:01:35.645458+00:00 · methodology

Review history (2 revisions) →

discussion (0)

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

Forward citations

Cited by 3 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.

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

    astro-ph.GA 2026-04 unverdicted novelty 5.0

    Filament widths in molecular clouds scale with spatial scale following power laws H ∝ Y^0.50 and h ∝ Y^0.37, with volume density slopes of 2.1-2.4 and contrasts 17-52 versus 1.1-2.7 for surface densities.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages · cited by 2 Pith papers

  1. [1]

    true” widths of power-law surface density profiles of filaments (or sources) when applying the widely used Gaussian (“naive

    André, P., Di Francesco, J., Ward-Thompson, D., et al. 2014, in Protostars and Planets VI, ed. H. Beuther, R. S. Klessen, C. P. Dullemond, & Th. Henning, Space Science Series, 27–51 André, P., Men’shchikov, A., Bontemps, S., et al. 2010, A&A, 518, L102+ André, P., Palmeirim, P., & Arzoumanian, D. 2022, A&A, 667, L1 Arzoumanian, D., André, P., Didelon, P.,...

    work page 2014

  2. [2]

    in- finite

    The inaccuracies of Eq. (A.2) are mostly within a few percent, except near the lower bounds of resolved- ness (see below), where errors become larger. For Gaussian pro- files (γ→ ∞), the formula is accurate within 12%. The concept of resolvedness becomes more complex for “in- finite” power-law structures (Men’shchikov 2023), becauseR has its lower bound o...

    work page 2023