PaperlayerView exact revision on arXiv

Source-linked AI summary

Structural Analysis of Molecular Clouds: Dendrograms

E. W. Rosolowsky, J. E. Pineda, J. Kauffmann, A. A. Goodman

arXiv:0802.2944v1astro-ph

TL;DR

Molecular-line analyses need a way to represent nested cloud structure and measure its properties across scales while accounting for noise and position-position-velocity ambiguities. The paper uses dendrograms to encode changing isosurface topology, applies multiscale molecular-gas measurements, and tests the method on molecular-cloud data. It finds self-gravitating structure across L1448 scales, recovers σv ∝ R^0.58, and identifies GMCs in blended Orion-Monoceros emission.

  • Problem

    Interpreting intensity-defined regions in position-position-velocity data cubes as physical objects is difficult because line-of-sight superposition and broad line profiles can make the mapping ambiguous.

  • Method

    The paper represents nested isosurfaces with dendrograms, measures molecular-gas properties on their branches, and accounts for noise in the resulting hierarchy.

  • Results

    Self-gravitating structures appear on all spatial scales in L1448, the size-linewidth relation is σv ∝ R^0.58, and dendrograms identify three GMCs in blended Orion-Monoceros data.

  • Takeaways & Limitations

    Dendrograms provide a multiscale structural reduction that links molecular-line emission to cloud properties and identifies physically defined GMCs in blended data.

  • Takeaways & Limitations

    Noise makes merge levels uncertain by approximately 2σrms, while position-position-velocity bijections can fail under line-of-sight superposition and broad line profiles.

Abstract

from arXiv · show

We demonstrate the utility of dendrograms at representing the essential features of the hierarchical structure of the isosurfaces for molecular line data cubes. The dendrogram of a data cube is an abstraction of the changing topology of the isosurfaces as a function of contour level. The ability to track hierarchical structure over a range of scales makes this analysis philosophically different from local segmentation algorithms like CLUMPFIND. Points in the dendrogram structure correspond to specific volumes in data cubes defined by their bounding isosurfaces. We further refine the technique by measuring the properties associated with each isosurface in the analysis allowing for a multiscale calculation of molecular gas properties. Using COMPLETE 13CO(1-0) data from the L1448 region in Perseus and mock observations of a simulated data cube, we identify regions that have a significant contribution by self-gravity to their energetics on a range of scales. We find evidence for self-gravitation on all spatial scales in L1448 though not in all regions. In the simulated observations, nearly all of the emission is found in objects that would be self-gravitating if gravity were included in the simulation. We reconstruct the size-line width relationship within the data cube using the dendrogram-derived properties and find it follows the standard relation: s_v ~ R^0.58. Finally, we show that constructing the dendrogram of CO J=1-0 emission from the Orion-Monoceros region allows for the identification of giant molecular clouds in a blended molecular line data set using only a physically motivated definition (self-gravitating clouds with masses 5x10^4 Msun.

1. INTRODUCTION

Molecular clouds have hierarchical structure that is closely tied to star formation, but existing analyses use different scales and struggle to connect observed molecular emission to physical structure. This paper introduces dendrograms as a multiscale representation of nested isosurfaces in three-dimensional molecular-line data cubes.

  • Motivation: Molecular-cloud structure helps determine the locations, numbers, and masses of newly formed stars.Its role is connected to the initial mass function and local star-formation rate.
  • Molecular-cloud hierarchy: Dense, small-scale features are nested within larger envelopes of lower-density gas, producing a non-trivial hierarchy across scales.There are more small-scale, dense structures than large-scale, sparse structures at a given scale.
  • Observational challenges: Connecting molecular clouds to surrounding atomic gas is difficult because 21-cm observations suffer foreground and background confusion and poorer spatial resolution.Specific geometries, self-absorption, or photodissociation-region modeling can enable some cases to be studied.
  • Contribution: Dendrograms represent nested isosurfaces in position-position-velocity data cubes and support simultaneous measurement of gas properties across physical scales.They also reduce large data sets to essential structural features and simple models.
  • Relation to prior work: The approach extends structure-tree ideas from Houlahan and Scalo by applying them to three-dimensional isosurfaces and emphasizing their properties.The paper distinguishes this use from dendrograms in statistical analysis.
  • Paper scope: The paper develops dendrograms, addresses noise and physical-domain complications, and applies the technique to self-gravity and blended GMC identification.The applications use L1448 and Orion-Monoceros molecular-line data.

2. THE ANALYSIS OF MOLECULAR LINE DATA

Molecular-line data analysis commonly either summarizes an entire data set statistically or segments it into presumed physically relevant structures. Segmentation is useful for substructure studies but is sensitive to observational resolution and sensitivity.

  • Analytic approaches: Molecular-line analyses generally follow either whole-data statistical description or segmentation into physically relevant structures.The two paths study global emission statistics or property distributions of identified objects.
  • Segmentation: Segmentation methods identify connected emission regions above thresholds and are used to construct GMC catalogs and clump populations.CLUMPFIND is a prominent watershed-segmentation example.
  • Observational dependence: Segmentation results are controlled by the sensitivity and resolution of the data set.This dependence is a stated limitation of threshold-based identification approaches.

3. THE DENDROGRAM TECHNIQUE

Dendrograms encode how connected structures split and merge as contour level changes, representing isosurfaces and their hierarchy in multidimensional data. The method selects significant local maxima while suppressing noise, but its tree structure is reliable only above noise-dependent amplitude changes.

  • Dendrogram representation: Each dendrogram point represents a three-dimensional isosurface bounding a distinct volume in the data cube at a specified contour level.Distinct surfaces are identified by the local maxima they contain.
  • Data representation: The dendrogram abstracts hierarchical structure while retaining spatial relationships that can be used to label regions on maps.The plotted leaf ordering itself does not encode spatial position.
  • Noise treatment: Noise suppression rejects local maxima likely caused by noise before constructing the initial dendrogram leaves.Real-data noise can create spurious maxima that do not correspond to physical emission structure.
  • Noise treatment: Local maxima are decimated using merge-level volumes, temperature contrasts, and minimum-volume criteria to retain distinct structures.The method uses Nmin and ΔTmax criteria, together with spatial and velocity windows Dmax and ΔVmax.
  • Dendrogram representation: Thresholding an emission profile at different contour levels changes the number of connected objects, and the dendrogram records the levels where structures separate or merge.In the schematic, I1 yields one object, I2 yields two, and Icrit separates those regimes.
  • Parameter dependence: Relaxing noise-suppression conditions produces more independent leaves while preserving the dendrogram’s basic structure and retaining the original isosurfaces as a subset.This comparison is shown for L1448 13CO emission.

4. MEASURING CLOUD PROPERTIES IN DENDROGRAMS

The paper measures cloud properties along dendrogram isosurfaces, while recognizing that mapping position-position-velocity emission to physical regions is interpretation-dependent. It compares bijection, clipping, and extrapolation treatments, then adopts bijection for substructure and uses the virial parameter as an approximate dynamical diagnostic.

  • Dendrogram properties: Each dendrogram point corresponds to a unique isosurface whose bounded emission is used to calculate physical properties.Properties vary continuously along branches and change abruptly where branches merge.
  • Property measurements: Cloud size, line width, luminosity, and mass are estimated from intensity-weighted moments and integrated emission within each isosurface.The method uses brightness-temperature weights, second spatial moments, and zeroth-moment flux measurements.
  • Interpretation: The physical meaning of an isosurface is ambiguous because PPV-space regions may not correspond directly to individual physical objects.This motivates comparing multiple interpretations of the emission and their resulting cloud properties.
  • Bijection paradigm: The bijection treatment assumes each data-cube pixel maps one-to-one onto a physical volume, but line-of-sight superposition and broad line profiles can violate that assumption.These effects can place multiple objects in one pixel or make one volume appear at multiple velocities.
  • Virial analysis: The virial parameter is treated as an approximate, mainly relative indicator of self-gravity, with regions having α ≲ 2 associated with significant gravitational potential.Magnetic fields, surface pressures, and bulk motions are omitted from this simple estimate.
  • Paradigm choice: The authors adopt bijection for substructure, while extrapolation produces larger radii and can add inferred emission beyond what the data cube contains.The extrapolated radius is most different at high contour levels and approaches approximately 1.5 pc at low thresholds.

5. THE HIERARCHICAL SUBSTRUCTURE OF A MOLECULAR

The paper applies dendrogram analysis to observed L1448 data and a matched turbulent simulation to characterize hierarchical molecular-cloud structure.

  • Observational data: The L1448 data cover a 3.1 pc × 3.1 pc region centered on the star-forming region, with 13CO emission spanning a 10 km s−1 interval.The observations have 46′′ angular resolution, 23′′ pixels, and uniform noise rms of 0.3 K.
  • Simulation data: The simulation uses synthetic 13CO emission from a 6 pc MHD turbulence simulation with mean density n = 10^3 cm−3 and mean Mach number M = 6.The simulated emission was generated by Monte Carlo radiative transfer from a snapshot of the simulated density and velocity fields.
  • Simulation data: The simulated data were convolved, resampled, and noise-matched to the observational data for direct comparison.The trial simulation region was selected from the most compact identifiable emission feature.

5.3. The 13CO-to-H2 Conversion Factor

The paper calibrates a linear 13CO-to-H2 conversion factor by comparing 13CO integrated intensity with near-infrared extinction-derived H2 column density.

  • Calibration: The conversion calibration matches the 13CO integrated-intensity map to the 48′′ extinction map derived from NICER near-infrared observations.The extinction map saturates above AV ∼22 mag.
  • Calibration: The adopted conversion is a single inverse-variance-weighted mean of N(H2)/W(13CO) across the data.The simple ratio is chosen because the conversion is applied to individual velocity channels.
  • Caveat: The linear conversion is adequate over the full range but systematically underestimates column density in high-brightness regions where 13CO saturates.Virial-parameter estimates are therefore likely overestimates in those regions.
  • Adopted factor: For the mass calculations, the paper adopts X2 = 4.0, comparable to the more sophisticated total-intensity result X2 = 2.1.The sophisticated analysis is not applicable to individual channels.

5.4. The Dynamical State of L1448

The L1448 dendrogram reveals self-gravitating structure across multiple spatial scales, including individual leaves and larger connected regions, while channel maps expose a low-velocity component hidden in integrated intensity.

  • Data visualization: The integrated-intensity map shows the L1448 emission, while selected channel maps reveal a low-velocity feature at vLSR ∼0.5 km s−1.The feature is not discernible in the integrated-intensity map.
  • Dynamical state: The L1448 dendrogram contains leaves with evidence for self-gravitation on small scales and larger structures that remain self-gravitating across the region.Several leaves associated with individual local maxima show self-gravitation, while the left-hand branch is largely self-gravitating until merger at contour levels ≲1.5 K.
  • Dynamical state: Using α ≤2 as the threshold for significant self-gravity identifies objects on a variety of scales, including nearly all of the left-hand branch and three distinct sub-branches.Four self-gravitating leaves are located and spatially delineated in the data cube.
  • Dendrogram analysis: Dendrogram properties are calculated for vertical branches representing isosurfaces, whereas horizontal merger branches carry no reported physical-property data.The dendrogram branches are color-coded by virial parameter, with values suppressed when formal errors exceed 50%.

5.5. The Dynamical State of the Turbulent Simulation

Applying the same dendrogram analysis to a turbulent simulation shows substantially more self-gravitating structure than in L1448, despite similar numbers of leaves and comparable antenna-temperature ranges.

  • Analysis setup: The simulation adopts a 13CO-to-H2 conversion factor of X2 = 10.9 derived by comparing synthetic 13CO data with simulated column density.The calibration uses the same procedure applied to the observational data.
  • Dendrogram comparison: The simulation and observations contain 39 and 26 leaves, respectively, with similar antenna-temperature spans.Most simulated mergers occur at higher intensity levels than observed mergers.
  • Dynamical state: Far more of the simulated data cube corresponds to self-gravitating objects than the observed L1448 data.The simulation dendrogram contains nearly all structure in self-gravitating objects, with only a few unbound leaves.
  • Dynamical state: The simulation contains ∼4 times as much molecular mass as L1448, while that additional mass is spread over a similar volume.The paper links this mass difference to the contrasting dynamical states of the two data cubes.

5.6. Interpretation of Dendrogram Properties

The interpretation of dendrogram-derived virial parameters depends strongly on how observed structures are mapped to physical properties. The clipping paradigm finds no self-gravitating structure in L1448, whereas extrapolation finds more than the bijection paradigm.

  • The clipping paradigm finds no self-gravitating structure anywhere in L1448.The authors regard this result as overly conservative given the region’s estimated virial parameter and star-forming small-scale clumps.
  • The extrapolation paradigm identifies more self-gravitating structure than the bijection paradigm.This follows because extrapolation applies a larger luminosity correction than corrections to radius and line width.

5.7. The Scale of Self-Gravity

Dendrograms measure self-gravity across spatial scales rather than selecting structures through local segmentation. In L1448, the self-gravitating fraction is small at small scales and increases toward larger scales.

  • Scale-dependent self-gravity: Dendrogram analysis tracks self-gravitating structure across multiple spatial scales, unlike local segmentation methods such as CLUMPFIND.The analysis measures virial parameters for all isosurfaces as a function of size scale.
  • Scale-dependent self-gravity: A small fraction of L1448 structure at small scales is self-gravitating, and the fraction grows at larger scales.The fraction is based on luminosity contained within isosurfaces with virial parameter α < 2.
  • Scale-dependent self-gravity: The analysis computes the luminosity fraction in size bins of ΔR = 0.2 dex for isosurfaces corresponding to self-gravitating objects.Each isosurface contributes through its luminosity, radius, and virial parameter.

5.8. The Size-Line Width Relationship in L1448

Dendrogram-derived isosurface properties probe the size-line width relation on intermediate scales within molecular clouds. For L1448, the fitted relation is σNT = (0.62 ± 0.04)R^0.58±0.04, while small-scale measurements become unreliable because of noise and instrumental effects.

  • Relationship and method: Dendrogram isosurfaces provide a multiscale probe of the size-line width relationship inside molecular clouds.The method measures each isosurface’s spatial and velocity extent and can represent independent dendrogram branches.
  • Relationship and method: σNT = (0.62 ± 0.04)R^0.58±0.04 describes the fitted L1448 size-line width relationship.The fit uses the gray-circled characteristic size and line width for independent dendrogram branches.
  • Limitations: Small-scale data become unreliable as thermal noise and instrumental convolution effects make isosurface properties poorly defined.The plotted sizes and line widths were corrected for beam convolution and thermal contributions.
  • Relationship and method: The dendrogram contains 26 leaves and 51 independent branches, yielding one representative point for each significantly distinct isosurface set.The branch count follows 2N − 1 for binary mergers.

6. IDENTIFYING GIANT MOLECULAR CLOUDS

The dendrogram identifies giant molecular clouds in blended Orion-Monoceros CO data as the largest-scale self-gravitating structures. Applying a mass threshold separates three GMCs from the remaining emission.

  • Physical definition: GMCs are defined as the largest-scale self-gravitating structures, identified through virial parameters near unity.The approach targets dynamically stable structures rather than relying only on segmentation.
  • Orion-Monoceros application: The Orion-Monoceros complex lies within one Tmb = 0.4 K isosurface, which the dendrogram decomposes into constituent clouds.The analysis uses 12CO data and an adopted CO-to-H2 conversion factor.
  • Orion-Monoceros application: Three GMCs are identified automatically as distinct self-gravitating regions with masses M > 5 × 10^4 M⊙.These three regions segregate naturally from the rest of the emission.

7. SUMMARY

Dendrograms provide a minimally model-dependent, multiscale representation of hierarchical molecular-line structure while supporting measurements of physical properties across nested isosurfaces. Applications to L1448, simulations, and Orion-Monoceros demonstrate analyses of self-gravity, size–line width scaling, and physically motivated GMC identification, although observed-to-physical-domain ambiguities remain.

  • Technique: Each dendrogram point corresponds to an isosurface, enabling molecular-gas properties, virial parameters, and size–line width relations to be measured across multiple scales.This provides a uniform way to examine energetics and characteristic structure from small to large scales.
  • Technique: Unlike CLUMPFIND-style segmentation, dendrograms preserve hierarchical emission structure and can support physically motivated object segmentation.The approach tracks nested structures over a range of scales, while user-selected simplification still governs the actual output.
  • Limitations: Relating observed structures to the physical domain remains ambiguous, and the paper states that no satisfactory universally applicable method resolves these ambiguities.The authors retain caveats about interpretation while arguing that they do not undermine the technique’s applicability.
  • Applications: L1448 contains self-gravitating structures on all spatial scales, although not in every region; most small-scale emission is not self-gravitating, while larger scales are often gravity-influenced.The analysis used COMPLETE 13CO(1 →0) observations of L1448.
  • Applications: Within a single molecular cloud, the recovered size–line width relation is σv ∝R^0.58.The relation is obtained from characteristic sizes and line widths of constituent isosurfaces.
  • Applications: GMCs are defined as clouds with M > 5 × 10^4 that are self-gravitating but not bound to their surrounding medium, identifying three Orion-Monoceros GMCs.The dendrogram identifies these objects in blended line data without including dynamically unrelated low-mass material.
  • Technique: Dendrograms reduce three-dimensional hierarchical data to a two-dimensional representation that preserves essential emission topology with minimal model dependence.The reduction follows the intrinsic structure of nested isosurfaces rather than imposing a fixed segmentation model.
Loading 0802.2944v1…