pith. machine review for the scientific record. sign in

arxiv: 2511.10486 · v2 · submitted 2025-11-13 · 🌌 astro-ph.GA

Revealing the Connection Between the Filamentary Hierarchy and Star Cluster Formation in a Simulated NGC 628 Galaxy

Pith reviewed 2026-05-17 22:33 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords filamentsstar clustersgravitational fragmentationmass functionpower-law indexgalaxy simulationNGC 628interstellar medium
0
0 comments X

The pith

Young star cluster masses in a simulated galaxy inherit their distribution directly from the mass spectrum of the gas filaments that fragment to form them.

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

The paper runs a magnetohydrodynamic simulation of a galaxy resembling NGC 628 and uses a filament-finding algorithm to measure the lengths and masses of the gas filaments that form inside it. It then applies a clustering algorithm to locate the star clusters that appear along those filaments between 268 and 278 Myr and tracks how the clusters change over the next 60 Myr. The mass distribution of the filaments follows a power law with index -1.35; the newly formed clusters show exactly the same index. This match indicates that the initial mass spectrum of the clusters is set by how the filaments break apart under their own gravity. As the clusters expand and lose mass they become unbound and the power-law index of their mass function steepens to -1.55, reproducing values seen in both other simulations and real observations.

Core claim

The mass function of young star clusters originates from the mass function of their parent filaments through gravitational fragmentation, demonstrated by identical power-law indices of -1.35 at the moment of cluster formation; continued evolution then drives the cluster mass function to a steeper index of -1.55 after 60 Myr as clusters expand and lose mass.

What carries the argument

Gravitational fragmentation of filaments, tracked by comparing the mass probability density functions of filaments and of the clusters they produce, which yields matching power-law slopes at birth.

If this is right

  • The initial mass spectrum of star clusters is fixed at the moment their host filaments fragment, before any later dynamical evolution occurs.
  • Cluster mass functions become steeper over time because clusters lose mass and grow in radius as they become unbound.
  • The observed power-law indices around -1.5 in real galaxies can be explained as the evolved state of an initially flatter distribution inherited from filaments.
  • Hierarchical filamentary structure in the interstellar medium directly imprints on the demographics of the first generation of star clusters.

Where Pith is reading between the lines

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

  • If filament mass functions vary across different galactic environments, the birth mass functions of clusters should vary in the same way.
  • Simulations that do not resolve individual filaments may need to impose the observed filament mass spectrum as an initial condition to recover realistic cluster populations.
  • Future maps of filament masses in nearby galaxies could be used to predict the mass distribution of the youngest clusters before dynamical evolution alters it.

Load-bearing premise

The filament-finding tool and clustering algorithm identify real physical structures rather than introducing selection biases that could force the two mass distributions to match by accident.

What would settle it

A high-resolution simulation or JWST observation of NGC 628 in which the filament mass power-law index differs from the initial cluster mass power-law index by more than the reported uncertainties.

Figures

Figures reproduced from arXiv: 2511.10486 by Rachel Pillsworth, Ralph E. Pudritz, Tamara Koletic.

Figure 1
Figure 1. Figure 1: Left: Image of NGC 628 taken by JWST (ESA/Webb, NASA & CSA, J. Lee and the PHANGS-JWST Team. Acknowledgement: J. Schmidt). Right: A plot of the gas density of our NGC 628-like simulation at 327 Myr zoomed-in to the center. The white lines are filaments identified using FilFinder which trace the densest gas concentrated along the spiral arms. NGC 628 simulated filament data is available on GitHub at https:/… view at source ↗
Figure 2
Figure 2. Figure 2: The gas density of our NGC 628-like sim￾ulation at 327 Myr. The white lines are filaments identified using FilFinder which trace the densest gas concentrated in the spiral arms. The clusters are of star particle formed between 307-327 Myr, semi-evenly spaced along the densest part of the arms. NGC 628 simulated filament data is available on GitHub at https://github.com/Tkoletic/NGC628-Filaments. skeletoniz… view at source ↗
Figure 3
Figure 3. Figure 3: Log-log plot of filament length (pc) vs mass (M⊙). Purple points are filaments from our NGC 628 simulation and grey are observed filaments compiled in A. Hacar et al. (2022) and reproduced with permission from the correspond￾ing PIs of A. Hacar et al. (2022); E. Schisano et al. (2019). NGC 628 simulated filament data is available on GitHub at https://github.com/Tkoletic/NGC628-Filaments. The sub￾set of Hac… view at source ↗
Figure 5
Figure 5. Figure 5: Mass distribution of filaments. Filaments with lengths less than 25 pc were cut to avoid structures near the resolution limit of our simulation, therefore these are also absent from the mass distribution. The black line shows the power-law fit to the distribution and the dotted and dash– dotted lines indicate the median and mean filament masses. NGC 628 simulated filament data is available on GitHub at htt… view at source ↗
Figure 4
Figure 4. Figure 4: Length distribution of filaments. Lengths less than 25 pc were cut to avoid structures near the resolu￾tion limit of our simulation. The black line shows the pow￾er-law fit to the distribution and the dotted and dash-dot￾ted lines indicating the median and mean filament lengths. NGC 628 simulated filament data is available on GitHub at https://github.com/Tkoletic/NGC628-Filaments. We employ the same proces… view at source ↗
Figure 6
Figure 6. Figure 6: The first frame depicts star particle formed between 268 − 278 Myr, subsequent frames show this same population after the galaxy has evolved for 20, 30, 40, 50, and 60 Myr. clusters are formed out of dense clouds, it follows that if the clouds fragment out of filaments and adopt their structural properties, clusters do as well. As the galaxy evolves over 50 Myr and the star particles start to drift from th… view at source ↗
Figure 7
Figure 7. Figure 7: Six plots of the mass distribution of filaments (green) and clusters (teal) identified in 3D. The filament mass distribution is of all filaments at a time of 327 Myr. The cluster mass distribution is of the star particle population formed between 268-278 Myr. Each frame evolves the star particle population by 10 Myr, ending at 327 Myr. The black lines show the power-law fit to the filament and cluster mass… view at source ↗
Figure 8
Figure 8. Figure 8: Six plots depicting the probability distribution function of all cluster radii at the 6 different star particle evolution times. Custom cuts were made to each time bin to remove the population of merging clusters identified by HDBSCAN as a single structure. The purple dash-dotted line is the mean cluster radius with values located in table 1 [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Radius growth rate (purple) and mass loss rate (teal) of clusters as the galaxy evolves [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Average eccentricities of clusters measured along the XY, YZ and ZX axes versus time. two-arm spiral of the galaxy, with very few inter-arm filamentary networks. Furthermore, since we have controlled the supernova rate and the column density of our NGC 628-like disk to be the same as in the Milky Way analogue, we can verify that the differences in the filamentary ISM’s structure between the two are exclus… view at source ↗
Figure 11
Figure 11. Figure 11: 6 select clusters 5 Myr after the star particles were formed. The first column shows the face-on view of identified clusters where the large pink triangle represents the location of the selected cluster. The middle column shows the star particles belonging to the cluster: each point = one star particle. We identify the cluster center with a small blue star. The last column depicts the mass density profile… view at source ↗
Figure 12
Figure 12. Figure 12: Cluster 0 time evolution over 20 Myr in steps of 5 Myr. The first row shows the location (pink triangle) of the cluster in the galaxy. The second row shows the scatter plot of individual star particles in the cluster. The last row shows the mass density profile where the dotted line represents the half-mass radius and the dash-dotted line is the standard deviation radius. The solid blue line is the Plumme… view at source ↗
Figure 13
Figure 13. Figure 13: shows six histograms in steps of 10 Myr of the distribution in the lengths of clusters measured along the x, y, and z axes. We see no shift along a single axis meaning that clusters increase in radius approximately isotropically [PITH_FULL_IMAGE:figures/full_fig_p018_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Cluster 1. Description of plots is the same as figure 12 [PITH_FULL_IMAGE:figures/full_fig_p019_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Cluster 2. Description of plots is the same as figure 12 [PITH_FULL_IMAGE:figures/full_fig_p019_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Cluster 3. Description of plots is the same as figure 12 [PITH_FULL_IMAGE:figures/full_fig_p020_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Cluster 4. Description of plots is the same as figure 12 [PITH_FULL_IMAGE:figures/full_fig_p020_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Cluster 5. Description of plots is the same as figure 12 [PITH_FULL_IMAGE:figures/full_fig_p022_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Six plots of the mass distribution of filaments (green) and clusters (teal) identified in 2D. The filament mass distribution is of all filaments at a time of 327 Myr. The cluster mass distribution is of the star particle population formed between 268-278 Myr and identified by applying HDBSCAN to the projected particle positions. Each frame evolves the star particle population by 10 Myr, ending at 327 Myr.… view at source ↗
read the original abstract

There is abundant observational evidence for the hierarchical, interconnected nature of filaments in the interstellar medium (ISM) extending from galactic down to sub-parsec scales. New JWST images of NGC 628 in particular, show clusters forming along the two spiral arms of this galaxy. In this paper we investigate filament and cluster properties in an NGC 628-like multi-scale high-resolution magnetohydrodynamic simulation. We use a filament finding tool to identify filaments and derive the probability density functions (PDFs) for the filament lengths and masses. Using a clustering algorithm we identify star clusters formed between 268 to 278 Myr and follow this population as the galaxy evolves for 60 Myr, calculating their mass PDFs, average radius growth rate, and average mass loss rate. We find a power-law index of alpha_m = -1.35 for the filament masses. Calculating the power-law index from our cluster mass PDF, we find a value of alpha_{c,m} = -1.35 when the clusters first form, exactly our filament mass power-law index. This shows that properties of young clusters arise from the gravitational fragmentation of their host filaments. We track the post-formation evolution of the clusters as they become unbound, increase in radius and decrease in mass yielding an ever steeper mass power-law index. After 60 Myr, the mass power-law index is alpha_{c,m} = -1.55, matching other simulations and observations.

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

Summary. The manuscript reports results from a multi-scale high-resolution MHD simulation of an NGC 628-like galaxy. A filament-finding tool is used to identify filaments and derive PDFs for filament lengths and masses, yielding a mass power-law index α_m = -1.35. A clustering algorithm identifies star clusters formed between 268 and 278 Myr; their mass PDF shows the identical index α_{c,m} = -1.35 at formation. The clusters are then evolved for 60 Myr, during which they become unbound, grow in radius, lose mass, and the mass power-law index steepens to α_{c,m} = -1.55. The authors interpret the initial index match as direct evidence that young cluster properties originate from gravitational fragmentation of their host filaments.

Significance. If the reported index match is shown to be physical rather than methodological, the result would link galactic filament fragmentation directly to the initial mass distribution of young star clusters, providing a numerical basis for the hierarchical nature of star formation observed in galaxies such as NGC 628. The tracking of post-formation evolution (radius growth and mass loss leading to a steeper index) adds a useful temporal dimension that aligns with other simulations and observations. The work is therefore potentially significant for models of clustered star formation, but its impact depends on demonstrating that the identification pipelines do not introduce correlated biases.

major comments (2)
  1. [Abstract] Abstract: The central claim is that the exact match α_m = -1.35 (filaments) = α_{c,m} = -1.35 (newly formed clusters at 268-278 Myr) demonstrates that cluster properties arise from gravitational fragmentation of host filaments. However, the abstract supplies no details on the filament-finding algorithm, the clustering method, the fitting procedure for the power laws, error bars, or possible selection biases. This omission leaves the equality difficult to evaluate and is load-bearing for the interpretation.
  2. [Abstract] Abstract and methods description: The interpretation assumes the filament finder (density ridges or similar) and cluster finder (proximity grouping of star particles) recover physically independent structures. No cross-validation is reported, such as spatial overlap statistics, results from alternative finders, or null tests on randomized particle distributions. Without these, the matching indices could arise from shared implicit thresholds or spatial scales rather than from the underlying fragmentation physics.
minor comments (1)
  1. [Abstract] The time interval 268-278 Myr is described as the formation window; clarifying whether this is a single snapshot or an integrated period over which clusters are identified would improve reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript. The comments highlight important aspects of clarity and validation that we address below. We have revised the manuscript to incorporate additional details and tests as outlined in our point-by-point responses.

read point-by-point responses
  1. Referee: [Abstract] Abstract: The central claim is that the exact match α_m = -1.35 (filaments) = α_{c,m} = -1.35 (newly formed clusters at 268-278 Myr) demonstrates that cluster properties arise from gravitational fragmentation of host filaments. However, the abstract supplies no details on the filament-finding algorithm, the clustering method, the fitting procedure for the power laws, error bars, or possible selection biases. This omission leaves the equality difficult to evaluate and is load-bearing for the interpretation.

    Authors: We agree that the abstract would benefit from greater methodological transparency to support evaluation of the central claim. In the revised manuscript we have expanded the abstract to include brief descriptions of the filament-finding algorithm (density ridge detection on the gas distribution), the clustering method (proximity-based grouping of star particles), the power-law fitting procedure (maximum-likelihood estimation on the cumulative distribution function), and the approach to uncertainties (bootstrap resampling). We also direct readers to the methods section for discussion of selection biases. These changes make the reported index match more readily evaluable while preserving the abstract's brevity. revision: yes

  2. Referee: [Abstract] Abstract and methods description: The interpretation assumes the filament finder (density ridges or similar) and cluster finder (proximity grouping of star particles) recover physically independent structures. No cross-validation is reported, such as spatial overlap statistics, results from alternative finders, or null tests on randomized particle distributions. Without these, the matching indices could arise from shared implicit thresholds or spatial scales rather than from the underlying fragmentation physics.

    Authors: This concern about possible methodological correlation is well taken. Although the filament identification operates on the gas density field and the cluster identification on the stellar particle distribution, we acknowledge that the original submission did not include explicit cross-validation. To address this, the revised manuscript now includes spatial overlap statistics between the identified filaments and young clusters, results from an alternative filament-finding algorithm, and null tests performed on randomized star-particle distributions. These additions demonstrate that the power-law index match is not reproduced under randomization and is therefore unlikely to stem from shared thresholds or scales. revision: yes

Circularity Check

0 steps flagged

No circularity: mass power-law indices are independently measured from distinct identification pipelines on simulation output.

full rationale

The paper extracts alpha_m = -1.35 directly from the filament mass PDF produced by its filament-finding tool and alpha_{c,m} = -1.35 from the cluster mass PDF produced by its separate clustering algorithm applied to star particles at formation time. These are reported as parallel measurements whose numerical agreement is then interpreted as evidence for gravitational fragmentation. No equation defines one index in terms of the other, no parameter is fitted to enforce the match, and the abstract invokes no self-citation or prior ansatz to justify the values. The derivation chain therefore consists of two independent post-processing steps on the same simulation snapshot and remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumption that the numerical tools extract unbiased physical filaments and clusters from the MHD data and that the chosen time windows capture representative formation and evolution phases.

free parameters (2)
  • Cluster formation time window = 268 to 278 Myr
    The interval 268-278 Myr is selected to identify newly formed clusters.
  • Post-formation evolution duration = 60 Myr
    60 Myr is used to track cluster radius growth and mass loss.
axioms (2)
  • domain assumption The filament finding tool accurately identifies gravitationally relevant structures in the simulated ISM.
    Invoked when deriving PDFs for filament lengths and masses.
  • domain assumption The clustering algorithm correctly groups stars into bound clusters and tracks their unbound evolution.
    Used to compute cluster mass PDFs, average radius growth rate, and mass loss rate.

pith-pipeline@v0.9.0 · 5567 in / 1473 out tokens · 40265 ms · 2026-05-17T22:33:21.030565+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

59 extracted references · 59 canonical work pages · 1 internal anchor

  1. [1]

    2021, The Astrophysical Journal, 916, 83, doi: 10.3847/1538-4357/ac07a1

    Abe, D., Inoue, T., Inutsuka, S.-i., & Matsumoto, T. 2021, The Astrophysical Journal, 916, 83, doi: 10.3847/1538-4357/ac07a1

  2. [2]

    2014, PLOS ONE, 9, 1, doi: 10.1371/journal.pone.0085777 Andr´ e, P

    Alstott, J., Bullmore, E., & Plenz, D. 2014, PLOS ONE, 9, 1, doi: 10.1371/journal.pone.0085777 Andr´ e, P. 2017, Comptes Rendus. G´ eoscience, 349, 187, doi: 10.1016/j.crte.2017.07.002 Andr´ e, P., Di Francesco, J., Ward-Thompson, D., et al. 2014, From Filamentary Networks to Dense Cores in Molecular Clouds: Toward a New Paradigm for Star Formation (Unive...

  3. [3]

    C., Arnaboldi, M., et al

    Aniyan, S., Freeman, K. C., Arnaboldi, M., et al. 2018, Monthly Notices of the Royal Astronomical Society, 476, 1909, doi: 10.1093/mnras/sty310

  4. [4]

    Arzoumanian, D., Arakawa, S., Kobayashi, M. I. N., et al. 2023, The Astrophysical Journal Letters, 947, L29, doi: 10.3847/2041-8213/acc849

  5. [5]

    2019, A&A, 621, A42, doi: 10.1051/0004-6361/201832725 Astropy Collaboration, Robitaille, T

    Arzoumanian, D., Andr´ e, Ph., K¨ onyves, V., et al. 2019, A&A, 621, A42, doi: 10.1051/0004-6361/201832725 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10.3847/1538-3881/aabc4f Astropy Coll...

  6. [6]

    Campello, R. J. G. B., Moulavi, D., & Sander, J. 2013, in Advances in Knowledge Discovery and Data Mining, ed. J. Pei, V. S. Tseng, L. Cao, H. Motoda, & G. Xu (Berlin, Heidelberg: Springer Berlin Heidelberg), 160–172

  7. [7]

    Chen, Y., Li, H., & Gnedin, O. Y. 2025, The Astrophysical Journal Letters, 980, L18, doi: 10.3847/2041-8213/adaf93

  8. [8]

    Chevance, M., Kruijssen, J. M. D., Hygate, A. P. S., et al. 2019, Monthly Notices of the Royal Astronomical Society, 493, 2872, doi: 10.1093/mnras/stz3525

  9. [9]

    M., Chastenet, J., et al

    Chiang, I.-D., Sandstrom, K. M., Chastenet, J., et al. 2021, The Astrophysical Journal, 907, 29, doi: 10.3847/1538-4357/abceb6 ESA/Webb, CSA, N. ., Lee, J., & Team, P.-J. 2022, Webb Inspects the Heart of the Phantom Galaxy,

  10. [10]

    2016, Monthly Notices of the Royal Astronomical Society, 457, 375, doi: 10.1093/mnras/stv2880

    Federrath, C. 2016, Monthly Notices of the Royal Astronomical Society, 457, 375, doi: 10.1093/mnras/stv2880

  11. [11]

    S., & Shaviv, N

    Fiege, J. D., & Pudritz, R. E. 2000, Monthly Notices of the Royal Astronomical Society, 311, 85, doi: 10.1046/j.1365-8711.2000.03066.x Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1, doi: 10.1051/0004-6361/202243940 Galv´ an-Madrid, R., Zhang, Q., Keto, E., et al. 2010, The Astrophysical Journal, 725, 17, doi: 10.1088/0004-6...

  12. [12]

    2022, Initial Conditions for Star Formation: A Physical Description of the Filamentary ISM, https://arxiv.org/abs/2203.09562

    Hacar, A., Clark, S., Heitsch, F., et al. 2022, Initial Conditions for Star Formation: A Physical Description of the Filamentary ISM, https://arxiv.org/abs/2203.09562

  13. [13]

    2024, Astronomy &; Astrophysics, 694, A69, doi: 10.1051/0004-6361/202450779

    Hacar, A., Konietzka, R., Seifried, D., et al. 2024, Astronomy &; Astrophysics, 694, A69, doi: 10.1051/0004-6361/202450779

  14. [14]

    R., Millman, K

    Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357, doi: 10.1038/s41586-020-2649-2

  15. [15]

    2023, Monthly Notices of the Royal Astronomical Society, 521, 5152–5159, doi: 10.1093/mnras/stad852

    Hoemann, E., Heigl, S., & Burkert, A. 2023, Monthly Notices of the Royal Astronomical Society, 521, 5152–5159, doi: 10.1093/mnras/stad852

  16. [16]

    S., Pudritz, R

    Howard, C. S., Pudritz, R. E., Sills, A., & Harris, W. E. 2019, Monthly Notices of the Royal Astronomical Society, 486, 1146–1155, doi: 10.1093/mnras/stz924

  17. [17]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, doi: 10.1109/MCSE.2007.55

  18. [18]

    Inutsuka, S.-i., & Miyama, S. M. 1997, The Astrophysical Journal, 480, 681, doi: 10.1086/303982

  19. [19]

    A., Kim, H., et al

    Kahre, L., Walterbos, R. A., Kim, H., et al. 2018, The Astrophysical Journal, 855, 133, doi: 10.3847/1538-4357/aab101

  20. [20]

    S., & Sills, A

    Karam, J., Fujii, M. S., & Sills, A. 2025, Dynamics of Star Cluster Formation: The Effects of Ongoing Star Formation and Stellar Feedback, https://arxiv.org/abs/2503.21716

  21. [21]

    2008, Proceedings of SPIE - The International Society for Optical Engineering, doi: 10.1117/12.788372

    Kashyap, V., Lee, H., Siemiginowska, A., et al. 2008, Proceedings of SPIE - The International Society for Optical Engineering, doi: 10.1117/12.788372

  22. [22]

    Kim, J., Chevance, M., Kruijssen, J. M. D., et al. 2022, Monthly Notices of the Royal Astronomical Society, 516, 3006–3028, doi: 10.1093/mnras/stac2339

  23. [23]

    2016, ApJ, 833, 202, doi: 10.3847/1538-4357/833/2/202

    Kim, J.-h., Agertz, O., Teyssier, R., et al. 2016, The Astrophysical Journal, 833, 202, doi: 10.3847/1538-4357/833/2/202

  24. [24]

    C., Bourke, T

    Kirk, H., Myers, P. C., Bourke, T. L., et al. 2013, The Astrophysical Journal, 766, 115, doi: 10.1088/0004-637x/766/2/115

  25. [25]

    W., & Rosolowsky, E

    Koch, E. W., & Rosolowsky, E. W. 2015, Monthly Notices of the Royal Astronomical Society, 452, 3435, doi: 10.1093/mnras/stv1521

  26. [26]

    Inutsuka, S. I. 2020, Astronomy &; Astrophysics, 642, A87, doi: 10.1051/0004-6361/202038232 17

  27. [27]

    2024, A&A, 681, A45, doi: 10.1051/0004-6361/202347268

    Lacchin, E., Mastrobuono-Battisti, A., Calura, F., et al. 2024, A&A, 681, A45, doi: 10.1051/0004-6361/202347268

  28. [28]

    2025, The Astrophysical Journal, 989, 22, doi: 10.3847/1538-4357/ade9b5

    Laverde-Villarreal, E., Sills, A., Cournoyer-Cloutier, C., & Arias Callejas, V. 2025, The Astrophysical Journal, 989, 22, doi: 10.3847/1538-4357/ade9b5

  29. [29]

    C., Sandstrom, K

    Lee, J. C., Sandstrom, K. M., Leroy, A. K., et al. 2023, ApJL, 944, L17, doi: 10.3847/2041-8213/acaaae

  30. [30]

    K., Hughes, A., Liu, D., et al

    Leroy, A. K., Hughes, A., Liu, D., et al. 2021, ApJS, 255, 19, doi: 10.3847/1538-4365/abec80

  31. [31]

    Malicious User Experience Design Research for Cybersecurity

    McKee, C. F., & Ostriker, E. C. 2007, Annual Review of Astronomy and Astrophysics, 45, 565–687, doi: 10.1146/annurev.astro.45.051806.110602

  32. [32]

    Myers, P. C. 2009, ApJ, 700, 1609, doi: 10.1088/0004-637X/700/2/1609

  33. [33]

    2024, NGC 628,

    NASA, ESA, CSA, et al. 2024, NGC 628,

  34. [34]

    1964, ApJ, 140, 1056, doi: 10.1086/148005 pandas development team, T

    Ostriker, J. 1964, ApJ, 140, 1056, doi: 10.1086/148005 pandas development team, T. 2025, pandas-dev/pandas: Pandas, v2.3.3 Zenodo, doi: 10.5281/zenodo.17229934

  35. [35]

    A., Andr´ e, P., et al

    Peretto, N., Fuller, G. A., Andr´ e, P., et al. 2014, Astronomy &; Astrophysics, 561, A83, doi: 10.1051/0004-6361/201322172

  36. [36]

    2013, Astronomy &; Astrophysics, 559, A38, doi: 10.1051/0004-6361/201322134

    Pfalzner, S., & Kaczmarek, T. 2013, Astronomy &; Astrophysics, 559, A38, doi: 10.1051/0004-6361/201322134

  37. [37]

    E., & Koch, E

    Pillsworth, R., Roscoe, E., Pudritz, R. E., & Koch, E. W. 2025a, Filamentary Hierarchies and Superbubbles I: Characterizing filament properties across a simulated spiral galaxy, https://arxiv.org/abs/2504.01099

  38. [38]

    E., & Koch, E

    Pillsworth, R., Roscoe, E., Pudritz, R. E., & Koch, E. W. 2025b, The Astrophysical Journal, 989, 233, doi: 10.3847/1538-4357/adf62c

  39. [39]

    Plummer, H. C. 1911, MNRAS, 71, 460, doi: 10.1093/mnras/71.5.460

  40. [40]

    S., Goodman, A

    Rice, T. S., Goodman, A. A., Bergin, E. A., Beaumont, C., & Dame, T. M. 2016, The Astrophysical Journal, 822, 52, doi: 10.3847/0004-637X/822/1/52

  41. [41]

    A., & Pudritz, R

    Robinson, H., Wadsley, J., Sellwood, J. A., & Pudritz, R. E. 2025, Star Formation and Magnetic Field Amplification due to Galactic Spirals, https://arxiv.org/abs/2506.16515

  42. [42]

    2019, Monthly Notices of the Royal Astronomical Society, 492, 5420, doi: 10.1093/mnras/stz3466

    Schisano, E., Molinari, S., Elia, D., et al. 2019, Monthly Notices of the Royal Astronomical Society, 492, 5420, doi: 10.1093/mnras/stz3466

  43. [43]

    Schneider, S., & Elmegreen, B. G. 1979, ApJS, 41, 87, doi: 10.1086/190609

  44. [44]

    2023, A&A, 672, A133, doi: 10.1051/0004-6361/202140857

    Shimajiri, Y., Andr´ e, Ph., Peretto, N., et al. 2023, A&A, 672, A133, doi: 10.1051/0004-6361/202140857

  45. [45]

    J., Shetty, R., Stutz, A

    Smith, R. J., Shetty, R., Stutz, A. M., & Klessen, R. S. 2012, The Astrophysical Journal, 750, 64, doi: 10.1088/0004-637x/750/1/64

  46. [46]

    Tang, J., Grasha, K., & Krumholz, M. R. 2024, Monthly Notices of the Royal Astronomical Society, 532, 4583, doi: 10.1093/mnras/stae1799

  47. [47]

    M., Kavanagh, P

    Temim, T., Laming, J. M., Kavanagh, P. J., et al. 2024, The Astrophysical Journal Letters, 968, L18, doi: 10.3847/2041-8213/ad50d1

  48. [48]

    2002, A&A, 385, 337, doi: 10.1051/0004-6361:20011817

    Teyssier, R. 2002, Astronomy &; Astrophysics, 385, 337, doi: 10.1051/0004-6361:20011817 The pandas development team. 2020, pandas-dev/pandas: Pandas, latest Zenodo, doi: 10.5281/zenodo.3509134

  49. [49]

    A., Lee, J

    Thilker, D. A., Lee, J. C., Deger, S., et al. 2023, The Astrophysical Journal Letters, 944, L13, doi: 10.3847/2041-8213/acaeac van der Giessen, S. A., Matsumoto, K., Relano, M., et al. 2024, Radial properties of dust in galaxies: Comparison between observations and isolated galaxy simulations, https://arxiv.org/abs/2410.21430

  50. [50]

    E., et al

    Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261, doi: 10.1038/s41592-019-0686-2

  51. [51]

    2015, Monthly Notices of the Royal Astronomical Society, 450, 4043, doi: 10.1093/mnras/stv735

    Wang, K., Testi, L., Ginsburg, A., et al. 2015, Monthly Notices of the Royal Astronomical Society, 450, 4043, doi: 10.1093/mnras/stv735

  52. [52]

    Wells, M. R. A., Beuther, H., Molinari, S., et al. 2024, Astronomy &; Astrophysics, 690, A185, doi: 10.1051/0004-6361/202449794

  53. [53]

    C., Chandar, R., Rodr´ ıguez, M

    Whitmore, B. C., Chandar, R., Rodr´ ıguez, M. J., et al. 2023, The Astrophysical Journal Letters, 944, L14, doi: 10.3847/2041-8213/acae94

  54. [54]

    J., Srinivasan, S., Pudritz, R

    Whitworth, D. J., Srinivasan, S., Pudritz, R. E., et al. 2025, On the relation between magnetic field strength and gas density in the interstellar medium: A multiscale analysis, https://arxiv.org/abs/2407.18293

  55. [55]

    G., Sun, J., Barnes, A

    Williams, T. G., Sun, J., Barnes, A. T., et al. 2022, The Astrophysical Journal Letters, 941, L27, doi: 10.3847/2041-8213/aca674

  56. [56]

    G., Lee, J

    Williams, T. G., Lee, J. C., Larson, K. L., et al. 2024, PHANGS-JWST: Data Processing Pipeline and First Full Public Data Release, https://arxiv.org/abs/2401.15142

  57. [57]

    2024, Filamentary Hierarchies and Superbubbles: Galactic Multiscale MHD Simulations of GMC to Star Cluster Formation, https://arxiv.org/abs/2405.18474

    Wadsley, J. 2024, Filamentary Hierarchies and Superbubbles: Galactic Multiscale MHD Simulations of GMC to Star Cluster Formation, https://arxiv.org/abs/2405.18474

  58. [58]

    2024, A&A, 691, A204, doi: 10.1051/0004-6361/202452064

    Zhou, Jian-Wen, Kroupa, Pavel, & Wu, Wenjie. 2024, A&A, 691, A204, doi: 10.1051/0004-6361/202452064

  59. [59]

    Zucker, C., Battersby, C., & Goodman, A. 2018, The Astrophysical Journal, 864, 153, doi: 10.3847/1538-4357/aacc66 18 APPENDIX A.EVOLUTION OF CLUSTER RADII BY AXIS Figure 13 shows six histograms in steps of 10 Myr of the distribution in the lengths of clusters measured along the x, y, and z axes. We see no shift along a single axis meaning that clusters in...