pith. sign in

arxiv: 2606.08494 · v1 · pith:XF6WCH5Bnew · submitted 2026-06-07 · 🌌 astro-ph.SR · astro-ph.GA

Virial-based extraction of structures in numerical simulations: The vibes tool

Pith reviewed 2026-06-27 18:06 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GA
keywords star formationcore extractionvirial theoremnumerical simulationscore mass functioninitial mass functionstructure identificationSTARFORGE
0
0 comments X

The pith

Virial theorem defines core boundaries in star formation simulations by tracking energy evolution during growth.

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

The paper presents a method to extract cores from 3D numerical simulations of star formation by applying the virial theorem iteratively around density peaks. Structures grow step by step and their boundaries are set where the total energy changes in a way that marks a physical limit rather than an arbitrary density value. This addresses the sensitivity of existing tools to user-chosen density thresholds, which affects how the core mass function is measured and linked to the stellar initial mass function. A reader would care because cores extracted this way are intended to match the definition of gas that will form one star or a close multiple system. The approach yields structures that remain consistent across different parameter choices in the tested simulations.

Core claim

The vibes tool builds structures iteratively from density peaks and applies the virial theorem at each step to monitor the structure's total energy; the boundary is chosen from the evolution of that energy, producing cores that are stable with respect to shape constraints, iteration step size, and peak selection while being less sensitive than density-threshold methods such as hop and dendrogram.

What carries the argument

Iterative virial-theorem evaluation on growing structures, which sets the boundary from the spatial evolution of the structure's total energy.

If this is right

  • Extracted cores remain coherent with each other and physically motivated across different simulation snapshots.
  • The method is far less sensitive to its working parameters than density-based algorithms.
  • Cores defined this way are expected to better represent the reservoirs that collapse into single or multiple stars.
  • The core mass function derived from such extractions should be more robust for studies connecting it to the stellar initial mass function.

Where Pith is reading between the lines

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

  • The same energy-based boundary criterion could be adapted to observational data cubes if line-width and column-density maps allow an approximate virial analysis.
  • Running the extraction on time-series snapshots might show how core boundaries migrate as turbulence and gravity evolve.
  • If the method proves reliable, it could be used to re-analyze existing simulation suites and test whether the resulting CMF shape changes systematically.

Load-bearing premise

The evolution of a structure's total energy as it grows spatially provides a reliable physical boundary that corresponds to the gas reservoir forming a single star or close multiple system.

What would settle it

In the same simulations, compare the material inside each extracted core at the identified boundary with the stars or multiples that actually form from that material at later times.

Figures

Figures reproduced from arXiv: 2606.08494 by Benjamin Thomasson, Daniel J. Price, Estelle Moraux, Fabien Louvet, Fr\'ed\'erique Motte, Isabelle Joncour, Marta Gonz\'alez-Garcia, Maxime Valeille-Manet, No\'e Brucy, Pierre Didelon, Simon Chevalier, Yann Bernard.

Figure 1
Figure 1. Figure 1: Structure energy components as a function of the equiva [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Illustration of a single cell iteration. Structures are built [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Illustration of the boundary selection process. The dashed [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Number of objects relative to the reference number of [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Number of peaks kept after sorting relative to the total [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Relative number of extracted objects with respect to the [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: Relative number of extracted objects with respect to the [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 10
Figure 10. Figure 10: Mass distribution of the structures extracted with [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Matrix giving the p values for Kolmogorov-Smirnov [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
Figure 13
Figure 13. Figure 13: Equivalent radius with respect to the mass of the struc [PITH_FULL_IMAGE:figures/full_fig_p012_13.png] view at source ↗
Figure 12
Figure 12. Figure 12: Peak density with respect to the number of cells per [PITH_FULL_IMAGE:figures/full_fig_p012_12.png] view at source ↗
read the original abstract

The processes that determine the stellar initial mass function (IMF) and its connection to the core mass function (CMF) are among the major open questions in star formation. The definition of a core remains unclear, yet the way they are extracted from simulations and observations critically shapes the CMF. Nowadays, cores are mostly detected through their density or intensity only. We aim to explore a new way to define cores in 3D numerical simulations based on a direct application of the virial theorem, and break free from some limitations induced by density-based methods. We intend to improve the accuracy and the physical meaning of the extracted cores. We developed vibes, an innovative method that makes full use of the virial theorem to extract overdensities in simulation snapshots. It works by building structures iteratively around density peaks, and applying the virial theorem to the structure at each iteration. Then, the structure boundary is set from the evolution of the its energy as it spatially grows. We used STARFORGE simulations to test the sensitivity of the extraction process to the main working parameters (constraints on the structure shape, iteration step, and peak selection criteria). This sensitivity is observed to be low. We compared our extraction with two density-based extraction algorithms, hop and dendrogram, that are observed to be very sensitive to their input density threshold parameter. Vibes returns structures that are coherent to each other and physically motivated, and it appears much more stable than existing 3D extraction tools. By defining the boundary of the cores on a physical criterion rather than on a user-defined set of density parameters, we expect such extracted cores to be closer to their forsaken definition: gas reservoirs that will form a single star or a close multiple system.

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

2 major / 2 minor

Summary. The paper introduces the vibes tool for extracting structures from 3D star-formation simulations. It builds overdensities iteratively around density peaks, applies the virial theorem at each step, and sets the boundary where the structure's total energy evolution indicates a physical transition. Tests on STARFORGE snapshots show low sensitivity to iteration step, peak selection, and shape constraints; comparisons indicate greater stability than hop or dendrogram methods, which depend strongly on density thresholds. The authors expect the resulting cores to better match the definition of gas reservoirs that form single stars or close multiples.

Significance. If the energy-evolution boundary reliably selects Lagrangian material that collapses to single stars, the method would offer a parameter-light, physically motivated alternative to density-threshold extractions and could tighten the CMF–IMF connection. The reported low parameter sensitivity and direct use of the virial theorem are concrete strengths; however, the manuscript provides no forward integration or multiplicity statistics to test the central physical correspondence.

major comments (2)
  1. [Abstract] Abstract: the expectation that energy-based boundaries yield cores 'closer to their forsaken definition' as single-star reservoirs rests on an untested premise; the manuscript reports only algorithmic robustness tests (parameter sweeps and hop/dendrogram comparisons) and contains no particle-tracing, forward collapse, or multiplicity comparison that would link the extracted volumes to actual star-formation outcomes.
  2. [Method] Method (iterative virial application): the boundary is chosen from the evolution of total energy as the structure grows, yet no quantitative criterion, error propagation, or sensitivity of this choice to the free parameters (iteration step, peak selection) is supplied; without this, it is unclear whether the reported low overall sensitivity extends to the physically decisive step.
minor comments (2)
  1. [Abstract] The abstract states that sensitivity 'is observed to be low' but does not report numerical measures (e.g., fractional variation in core mass or radius across the tested parameter ranges).
  2. [Method] Notation for the virial quantities (kinetic, potential, and surface terms) should be defined explicitly when first introduced to allow readers to reproduce the energy-evolution diagnostic.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments. We address each major comment below.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the expectation that energy-based boundaries yield cores 'closer to their forsaken definition' as single-star reservoirs rests on an untested premise; the manuscript reports only algorithmic robustness tests (parameter sweeps and hop/dendrogram comparisons) and contains no particle-tracing, forward collapse, or multiplicity comparison that would link the extracted volumes to actual star-formation outcomes.

    Authors: We agree that the manuscript contains no particle-tracing, forward integration, or multiplicity statistics. The abstract uses 'we expect' to signal a physically motivated hypothesis based on the virial theorem rather than a demonstrated result. We will revise the abstract to distinguish clearly between the demonstrated algorithmic stability and the anticipated physical correspondence to single-star reservoirs. revision: yes

  2. Referee: [Method] Method (iterative virial application): the boundary is chosen from the evolution of total energy as the structure grows, yet no quantitative criterion, error propagation, or sensitivity of this choice to the free parameters (iteration step, peak selection) is supplied; without this, it is unclear whether the reported low overall sensitivity extends to the physically decisive step.

    Authors: The reported sensitivity tests cover iteration step and peak selection, but we acknowledge that an explicit quantitative criterion for boundary selection from the energy evolution, together with associated sensitivity analysis, is not detailed. We will expand the methods section to define the boundary criterion precisely and add targeted tests of its sensitivity to the free parameters. revision: yes

Circularity Check

0 steps flagged

No circularity: direct application of virial theorem to external simulation data

full rationale

The derivation applies the standard virial theorem iteratively to structures grown from density peaks in STARFORGE snapshots. Boundary selection follows from computed total-energy evolution at each spatial growth step. This is a direct computation on input data without parameter fitting, self-referential definitions, or load-bearing self-citations. No equation reduces to its own input by construction, and the method is independent of the interpretive claim about correspondence to single-star collapse outcomes. The paper remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim depends on the standard virial theorem holding for growing structures in the simulations and on the interpretation that the resulting boundaries match the physical core definition. Two working parameters are noted but their sensitivity is reported as low.

free parameters (2)
  • iteration step
    Listed among the main working parameters whose sensitivity was tested.
  • peak selection criteria
    Listed among the main working parameters whose sensitivity was tested.
axioms (1)
  • domain assumption The virial theorem can be applied iteratively to growing overdensities to determine physically meaningful boundaries.
    This is the core mechanism described for setting structure boundaries.

pith-pipeline@v0.9.1-grok · 5897 in / 1307 out tokens · 26185 ms · 2026-06-27T18:06:18.905558+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

52 extracted references · 50 canonical work pages · 25 internal anchors

  1. [1]

    Powerlaw: a Python package for analysis of heavy-tailed distributions

    Powerlaw: a Python package for analysis of heavy-tailed distributions , volume =. PLoS ONE , author =. doi:10.1371/journal.pone.0085777 , shorttitle =. 1305.0215 [physics] , keywords =

  2. [2]

    , keywords =

    Six myths on the virial theorem for interstellar clouds. , keywords =. doi:10.1111/j.1365-2966.2006.10880.x , archivePrefix =. astro-ph/0606071 , primaryClass =

  3. [3]

    , keywords =

    On the structure of molecular clouds. , keywords =. doi:10.1111/j.1365-2966.2012.22130.x , archivePrefix =. 1209.3346 , primaryClass =

  4. [4]

    Gravity or turbulence? IV. Collapsing cores in out-of-virial disguise

    Gravity or turbulence? - IV. Collapsing cores in out-of-virial disguise. , keywords =. doi:10.1093/mnras/sty1515 , archivePrefix =. 1710.07384 , primaryClass =

  5. [5]

    FellWalker - a Clump Identification Algorithm

    FellWalker-A clump identification algorithm. Astronomy and Computing , keywords =. doi:10.1016/j.ascom.2014.11.004 , archivePrefix =. 1411.6267 , primaryClass =

  6. [6]

    Towards a more realistic sink particle algorithm for the RAMSES code

    Towards a more realistic sink particle algorithm for the RAMSES CODE. , keywords =. doi:10.1093/mnras/stu2005 , archivePrefix =. 1409.6528 , primaryClass =

  7. [7]

    Star formation: statistical measure of the correlation between the prestellar core mass function and the stellar initial mass function

    Star Formation: Statistical Measure of the Correlation between the Prestellar Core Mass Function and the Stellar Initial Mass Function. , keywords =. doi:10.1088/2041-8205/725/1/L79 , archivePrefix =. 1011.1185 , primaryClass =

  8. [8]

    , year = 1953, month = jul, volume =

    Problems of Gravitational Stability in the Presence of a Magnetic Field. , year = 1953, month = jul, volume =. doi:10.1086/145732 , adsurl =

  9. [9]

    Power-law distributions in empirical data

    Power-Law Distributions in Empirical Data. SIAM Review , keywords =. doi:10.1137/070710111 , archivePrefix =. 0706.1062 , primaryClass =

  10. [10]

    , keywords =

    Cloud properties across spatial scales in simulations of the interstellar medium. , keywords =. doi:10.1051/0004-6361/202348983 , archivePrefix =. 2403.00512 , primaryClass =

  11. [11]

    , keywords =

    The fragmentation of molecular clouds in starburst environments. , keywords =. doi:10.1093/mnras/stae2812 , archivePrefix =. 2501.03323 , primaryClass =

  12. [12]

    Journal of Computational Physics , year = 2002, month = jan, volume =

    Hyperbolic Divergence Cleaning for the MHD Equations. Journal of Computational Physics , year = 2002, month = jan, volume =. doi:10.1006/jcph.2001.6961 , adsurl =

  13. [13]

    The Virial Balance of Clumps and Cores in Molecular Clouds

    The Virial Balance of Clumps and Cores in Molecular Clouds. , keywords =. doi:10.1086/513708 , archivePrefix =. astro-ph/0607362 , primaryClass =

  14. [14]

    Core Mass Function: The Role of Gravity

    Core Mass Function: The Role of Gravity. , keywords =. doi:10.1086/588608 , archivePrefix =. 0801.2257 , primaryClass =

  15. [15]

    HOP: A New Group-Finding Algorithm for N-body Simulations

    HOP: A New Group-Finding Algorithm for N-Body Simulations. , keywords =. doi:10.1086/305535 , archivePrefix =. astro-ph/9712200 , primaryClass =

  16. [16]

    A High Order Godunov Scheme with Constrained Transport and Adaptive Mesh Refinement for Astrophysical MHD

    A high order Godunov scheme with constrained transport and adaptive mesh refinement for astrophysical magnetohydrodynamics. , keywords =. doi:10.1051/0004-6361:20065371 , archivePrefix =. astro-ph/0607230 , primaryClass =

  17. [17]

    , keywords =

    Smoothed particle hydrodynamics: theory and application to non-spherical stars. , keywords =. doi:10.1093/mnras/181.3.375 , adsurl =

  18. [18]

    Dense core formation in supersonic turbulent converging flows

    Dense Core Formation in Supersonic Turbulent Converging Flows. , keywords =. doi:10.1088/0004-637X/729/2/120 , archivePrefix =. 1101.2650 , primaryClass =

  19. [19]

    , keywords =

    STARFORGE: Towards a comprehensive numerical model of star cluster formation and feedback. , keywords =. doi:10.1093/mnras/stab1347 , archivePrefix =. 2010.11254 , primaryClass =

  20. [20]

    , keywords =

    Effects of the environment and feedback physics on the initial mass function of stars in the STARFORGE simulations. , keywords =. doi:10.1093/mnras/stac2060 , archivePrefix =. 2205.10413 , primaryClass =

  21. [21]

    , keywords =

    The Physical Origin of the Stellar Initial Mass Function. , keywords =. doi:10.1146/annurev-astro-052622-031748 , archivePrefix =. 2404.07301 , primaryClass =

  22. [22]

    GIZMO: A New Class of Accurate, Mesh-Free Hydrodynamic Simulation Methods

    A new class of accurate, mesh-free hydrodynamic simulation methods. , keywords =. doi:10.1093/mnras/stv195 , archivePrefix =. 1409.7395 , primaryClass =

  23. [23]

    A Constrained-Gradient Method to Control Divergence Errors in Numerical MHD

    A constrained-gradient method to control divergence errors in numerical MHD. , keywords =. doi:10.1093/mnras/stw1578 , archivePrefix =. 1509.07877 , primaryClass =

  24. [24]

    The Formation of Stellar Clusters: Mass Spectra from Turbulent Molecular Cloud Fragmentation

    The Formation of Stellar Clusters: Mass Spectra from Turbulent Molecular Cloud Fragmentation. , keywords =. doi:10.1086/321626 , archivePrefix =. astro-ph/0104127 , primaryClass =

  25. [25]

    A census of dense cores in the Aquila cloud complex: SPIRE/PACS observations from the Herschel Gould Belt survey

    A census of dense cores in the Aquila cloud complex: SPIRE/PACS observations from the Herschel Gould Belt survey. , keywords =. doi:10.1051/0004-6361/201525861 , archivePrefix =. 1507.05926 , primaryClass =

  26. [26]

    The Core Mass Function across Galactic Environments. II. Infrared Dark Cloud Clumps. , keywords =. doi:10.3847/1538-4357/aacb7c , archivePrefix =. 1806.02213 , primaryClass =

  27. [27]

    The W43-MM1 mini-starburst ridge, a test for star formation efficiency models

    The W43-MM1 mini-starburst ridge, a test for star formation efficiency models. , keywords =. doi:10.1051/0004-6361/201423603 , archivePrefix =. 1404.4843 , primaryClass =

  28. [28]

    , keywords =

    Strong dependence of the physical properties of cores on spatial resolution in observations and simulations. , keywords =. doi:10.1051/0004-6361/202040053 , archivePrefix =. 2107.05813 , primaryClass =

  29. [29]

    Core mass function in the high-mass star formation regime

    ALMA-IMF: XV. Core mass function in the high-mass star formation regime. , keywords =. doi:10.1051/0004-6361/202345986 , archivePrefix =. 2407.18719 , primaryClass =

  30. [30]

    , keywords =

    A numerical approach to the testing of the fission hypothesis. , keywords =. doi:10.1086/112164 , adsurl =

  31. [31]

    , keywords =

    On the Virial Theorem for Turbulent Molecular Clouds. , keywords =. doi:10.1086/171946 , adsurl =

  32. [32]

    , keywords =

    Multiscale, multiwavelength extraction of sources and filaments using separation of the structural components: getsf. , keywords =. doi:10.1051/0004-6361/202039913 , archivePrefix =. 2102.11565 , primaryClass =

  33. [33]

    Source extraction and photometry for the far-infrared and sub-millimeter continuum in the presence of complex backgrounds

    Source extraction and photometry for the far-infrared and sub-millimeter continuum in the presence of complex backgrounds. , keywords =. doi:10.1051/0004-6361/201014752 , archivePrefix =. 1011.3946 , primaryClass =

  34. [34]

    , keywords =

    A refined particle method for astrophysical problems. , keywords =

  35. [35]

    , keywords =

    The initial conditions of star formation in the rho Ophiuchi main cloud: wide-field millimeter continuum mapping. , keywords =

  36. [36]

    The unexpectedly large proportion of high-mass star-forming cores in a Galactic mini-starburst

    The unexpectedly large proportion of high-mass star-forming cores in a Galactic mini-starburst. Nature Astronomy , keywords =. doi:10.1038/s41550-018-0452-x , archivePrefix =. 1804.02392 , primaryClass =

  37. [37]

    High-Mass Star and Massive Cluster Formation in the Milky Way

    High-Mass Star and Massive Cluster Formation in the Milky Way. , keywords =. doi:10.1146/annurev-astro-091916-055235 , archivePrefix =. 1706.00118 , primaryClass =

  38. [38]

    The core and stellar mass functions in massive collapsing filaments

    Core and stellar mass functions in massive collapsing filaments. , keywords =. doi:10.1051/0004-6361/201834094 , archivePrefix =. 1902.05744 , primaryClass =

  39. [39]

    , keywords =

    The Life and Times of Star-forming Cores: An Analysis of Dense Gas in the STARFORGE Simulations. , keywords =. doi:10.3847/1538-4357/adb71d , archivePrefix =. 2502.15057 , primaryClass =

  40. [40]

    The Perils of Clumpfind: The Mass Spectrum of Sub-structures in Molecular Clouds

    The Perils of Clumpfind: The Mass Spectrum of Substructures in Molecular Clouds. , keywords =. doi:10.1088/0004-637X/699/2/L134 , archivePrefix =. 0906.0331 , primaryClass =

  41. [41]

    ALMA-IMF. III. Investigating the origin of stellar masses: top-heavy core mass function in the W43-MM2&MM3 mini-starburst. , keywords =. doi:10.1051/0004-6361/202142951 , archivePrefix =. 2203.03276 , primaryClass =

  42. [42]

    Journal of Computational Physics , year = 1999, month = sep, volume =

    A Solution-Adaptive Upwind Scheme for Ideal Magnetohydrodynamics. Journal of Computational Physics , year = 1999, month = sep, volume =. doi:10.1006/jcph.1999.6299 , adsurl =

  43. [43]

    The density variance -- Mach number relation in supersonic, isothermal turbulence

    The Density Variance-Mach Number Relation in Supersonic, Isothermal Turbulence. , keywords =. doi:10.1088/2041-8205/727/1/L21 , archivePrefix =. 1010.3754 , primaryClass =

  44. [44]

    Structural Analysis of Molecular Clouds: Dendrograms

    Structural Analysis of Molecular Clouds: Dendrograms. , keywords =. doi:10.1086/587685 , archivePrefix =. 0802.2944 , primaryClass =

  45. [45]

    The ALMA Survey of 70 m Dark High-mass Clumps in Early Stages (ASHES). I. Pilot Survey: Clump Fragmentation. , keywords =. doi:10.3847/1538-4357/ab45e9 , archivePrefix =. 1909.07985 , primaryClass =

  46. [46]

    doi:10.1111/j.1365-2966.2009.15868.x , eprint =

    The simultaneous formation of massive stars and stellar clusters. , keywords =. doi:10.1111/j.1365-2966.2009.15621.x , archivePrefix =. 0908.3910 , primaryClass =

  47. [47]

    , keywords =

    The highly variable time evolution of star-forming cores identified with dendrograms. , keywords =. doi:10.1093/mnras/staa2253 , archivePrefix =. 2004.01263 , primaryClass =

  48. [48]

    Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES

    Cosmological hydrodynamics with adaptive mesh refinement. A new high resolution code called RAMSES. , keywords =. doi:10.1051/0004-6361:20011817 , archivePrefix =. astro-ph/0111367 , primaryClass =

  49. [49]

    Towards a comprehensive scenario

    Global hierarchical collapse in molecular clouds. Towards a comprehensive scenario. , keywords =. doi:10.1093/mnras/stz2736 , archivePrefix =. 1903.11247 , primaryClass =

  50. [50]

    , keywords =

    Determining Structure in Molecular Clouds. , keywords =. doi:10.1086/174279 , adsurl =

  51. [51]

    ALMA-IMF. XX. Core Fragmentation in the W51 High-mass Star-forming Region. , keywords =. doi:10.3847/1538-4357/ae0619 , archivePrefix =. 2509.06749 , primaryClass =

  52. [52]

    , keywords =

    Dense Cores in IRDC G14.225-0.506 Revealed by ALMA Observations. , keywords =. doi:10.3847/1538-4357/adc398 , archivePrefix =. 2502.16413 , primaryClass =