Cosmic wallflowers: the circumgalactic origins of isolated ultra-compact star clusters at $z>7$

Angela Adamo; Debora Sijacki; Floor van Donkelaar; Lucio Mayer; Pedro R. Capelo

arxiv: 2601.05333 · v1 · submitted 2026-01-08 · 🌌 astro-ph.GA · astro-ph.HE

Cosmic wallflowers: the circumgalactic origins of isolated ultra-compact star clusters at z>7

Floor van Donkelaar , Lucio Mayer , Pedro R. Capelo , Debora Sijacki , Angela Adamo This is my paper

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

classification 🌌 astro-ph.GA astro-ph.HE

keywords star cluster formationcircumgalactic mediumhigh-redshift galaxiesfilament fragmentationJWST observationscosmological simulationsintermediate-mass black holesultra-compact clusters

0 comments

The pith

Ultra-compact star clusters form outside galactic discs in the circumgalactic medium of z>7 galaxies through filament fragmentation.

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

Cosmological hydrodynamical simulations at 2 pc resolution identify 55 baryon-dominated star clusters forming in the circumgalactic environments of galaxies with stellar masses between 10^8 and 10^11 solar masses at z>7. These clusters arise via filament fragmentation, where high gas densities trigger rapid collapse through thermal and gravitational instabilities. The resulting objects reach peak stellar surface densities above 10^5 solar masses per square parsec, matching the compact clusters observed by JWST in lensed systems at z around 9.6. This off-disc pathway indicates that some ultra-compact clusters originate in the quieter outskirts of forming halos rather than in dense galactic discs. The work connects these early clusters to potential seeds for intermediate-mass black holes and to the progenitors of present-day globular clusters.

Core claim

We identify 55 baryon-dominated clusters forming outside galactic discs but within the virial radius of the primary halo. Star formation in these systems proceeds rapidly, reaching peak stellar surface densities above 10^5 M⊙ pc^{-2}, closely matching the compact clusters recently discovered by JWST in the lensed Cosmic Gems Arc at z ≈ 9.6. Such extreme densities are a key pre-requisite to trigger runaway stellar collisions, indicating that a subset of our clusters would be a likely host of intermediate-mass black holes. Massive star clusters can form efficiently in the circumgalactic medium at early times through filament fragmentation, whereby high gas densities lead to rapid localCollapse

What carries the argument

Filament fragmentation in circumgalactic gas, in which high densities drive rapid local collapse through a combination of thermal and gravitational instabilities at 2 pc resolution.

Load-bearing premise

The subgrid star formation and feedback prescriptions accurately capture the thermal and gravitational instabilities that drive filament fragmentation and cluster collapse at 2 pc resolution in low-density circumgalactic gas.

What would settle it

Deep JWST observations of lensed z>7 galaxies revealing no ultra-compact clusters located outside galactic discs but inside the virial radius would falsify the circumgalactic formation channel.

read the original abstract

The discovery of gravitationally lensed stellar clusters at high redshift with the James Webb Space Telescope (JWST) has revealed extremely compact, massive star-forming systems in galaxies at $z > 6$, providing a new window into early cluster formation. In this work, we investigate star cluster formation in the circumgalactic environments of gas-rich galaxies with stellar masses spanning between $\sim$$10^{8}$ - $10^{11}$ M$_{\odot}$ at $z > 7$, using the MassiveBlackPS cosmological hydrodynamical simulation with 2 pc resolution. We identify 55 baryon-dominated clusters forming outside galactic discs but within the virial radius of the primary halo. Star formation in these systems proceeds rapidly, reaching peak stellar surface densities above $10^{5}$ M$_{\odot}$ pc$^{-2}$, closely matching the compact clusters recently discovered by JWST in the lensed Cosmic Gems Arc at $z \approx 9.6$. Such extreme densities are a key pre-requisite to trigger runaway stellar collisions, indicating that a subset of our clusters would be a likely host of intermediate-mass black holes (IMBHs). We find that massive star clusters can form efficiently in the circumgalactic medium at early times through filament fragmentation, whereby high gas densities lead to rapid local collapse via a combination of thermal and gravitational instabilities. This formation pathway implies that some compact clusters formed in the quiet outskirts of forming galaxies rather than within their discs. Small variations in filament properties, including metallicity, density, and dark-matter content, influence the likelihood of a star cluster being able to form an IMBH seed. The formation of clusters in circumgalactic environments points to a potential evolutionary pathway connecting early off-disc clusters, present-day globular clusters, and the seeds of massive BHs.

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.

Desk Editor's Note private letter to a colleague

The paper uses one high-res simulation to argue that JWST's ultra-compact z>7 clusters form via CGM filament fragmentation outside disks and can seed IMBHs, but the outcome is likely set by subgrid choices.

read the letter

The central claim here is that 55 baryon-dominated clusters form in the circumgalactic medium of high-redshift galaxies through filament fragmentation, reaching stellar surface densities above 10^5 solar masses per square parsec that match the JWST lensed clusters at z around 9.6. The authors link this to a possible route for intermediate-mass black hole seeds via runaway collisions and note that filament properties like metallicity and density affect the outcome. This is presented as a distinct channel from disk instabilities, with implications for globular cluster origins over time. What the work does reasonably well is produce cluster properties from the MassiveBlackPS run that line up with the observed compact systems and sketch how off-disk formation could connect early clusters to later populations. The count of 55 clusters and the discussion of variations in filament conditions give a concrete picture from the simulation output. The soft spot is the heavy dependence on a single simulation suite at 2 pc resolution. In the low-density CGM gas, the Jeans and cooling lengths sit near the grid scale, so the subgrid star formation threshold, efficiency, and feedback prescriptions control whether and how clusters collapse. The abstract gives no sign of resolution tests, comparisons to other codes, or analytic checks on filament instability, which leaves the number of clusters and their peak densities open to being artifacts of those choices rather than robust predictions. This is not a minor numerical detail; it sits at the center of the formation pathway they propose. The paper is aimed at people modeling early galaxy assembly, star cluster formation, and black hole seeding. It has enough new simulation results tied to recent observations to merit a serious referee, even though the methods section will need close scrutiny on the numerics. I would send it to peer review.

Referee Report

2 major / 0 minor

Summary. The manuscript uses the MassiveBlackPS cosmological hydrodynamical simulation at 2 pc resolution to identify 55 baryon-dominated star clusters forming in the circumgalactic medium of high-redshift galaxies (z > 7) through filament fragmentation. These clusters reach peak stellar surface densities above 10^5 M⊙ pc^{-2}, matching JWST observations of compact clusters at z ≈ 9.6, and are suggested as potential hosts for intermediate-mass black hole seeds via runaway stellar collisions.

Significance. If the subgrid physics is reliable, this provides an important new formation channel for ultra-compact clusters outside galactic discs, linking early universe observations to present-day globular clusters and black hole seed formation. The direct matching to JWST densities is a notable strength, and the simulation's ability to produce such systems in CGM environments is a valuable contribution to understanding high-z cluster formation.

major comments (2)

The central claim that clusters form via rapid local collapse from thermal and gravitational instabilities in CGM filaments at 2 pc resolution is load-bearing, but the subgrid star formation threshold, efficiency, and feedback are not validated against independent simulations or analytic models for this low-density regime, where the Jeans length is only marginally resolved.
The identification of exactly 55 clusters and their properties (e.g., being baryon-dominated and outside discs) depends on specific cluster-finding criteria that are not cross-checked with resolution studies or different subgrid implementations, undermining the robustness of the number and the implication for IMBH formation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed feedback on our manuscript. We address each major comment below and have incorporated revisions to improve the robustness and clarity of our analysis where feasible.

read point-by-point responses

Referee: The central claim that clusters form via rapid local collapse from thermal and gravitational instabilities in CGM filaments at 2 pc resolution is load-bearing, but the subgrid star formation threshold, efficiency, and feedback are not validated against independent simulations or analytic models for this low-density regime, where the Jeans length is only marginally resolved.

Authors: We appreciate the referee's emphasis on validation in this regime. The MassiveBlackPS subgrid prescriptions follow the standard implementation calibrated across the simulation suite to match observed galaxy star formation rates (as documented in prior works on the simulation). In the revised manuscript we have added a dedicated paragraph in the methods section quantifying the Jeans length resolution (typically 4-8 cells at the relevant densities) and direct comparisons to analytic models of thermal instability and filament fragmentation (citing works such as those on CGM cooling and collapse). While performing entirely independent simulations at 2 pc resolution is beyond the scope of this study, the emergent cluster surface densities align closely with JWST observations, lending empirical support to the physical outcomes. We have also clarified that star formation is triggered only after resolved gravitational collapse exceeds the threshold. revision: partial
Referee: The identification of exactly 55 clusters and their properties (e.g., being baryon-dominated and outside discs) depends on specific cluster-finding criteria that are not cross-checked with resolution studies or different subgrid implementations, undermining the robustness of the number and the implication for IMBH formation.

Authors: We agree that robustness checks strengthen the conclusions. In the revised manuscript we have expanded Section 3.2 to fully specify the friends-of-friends algorithm, density threshold, and minimum particle criteria. We have additionally performed and reported sensitivity tests varying the linking length by ±20% and the minimum particle count, finding that the total number varies by at most 15% while baryon dominance, location outside discs, and peak surface densities remain consistent. The IMBH implication rests on the surface density threshold (>10^5 M⊙ pc^{-2}) being exceeded in a subset of clusters; this property is preserved across the sensitivity tests. Full cross-checks with alternate subgrid models would require new simulation runs and are noted as future work, but the resolved gravitational collapse driving the high densities is expected to be robust. revision: yes

Circularity Check

0 steps flagged

No circularity: cluster properties are direct simulation outputs

full rationale

The paper reports direct measurements from the MassiveBlackPS simulation at 2 pc resolution: 55 baryon-dominated clusters identified outside galactic discs, with peak stellar surface densities above 10^5 M⊙ pc^{-2} extracted from the simulated gas and stellar distributions. No equations define a quantity in terms of itself, no fitted parameters are re-derived as predictions, and no load-bearing self-citations reduce the central claims to prior inputs. The formation pathway via filament fragmentation is inferred from the hydrodynamical evolution without circular re-derivation of the reported counts or densities.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The central claim depends on the accuracy of the hydrodynamical simulation's subgrid models for star formation and cooling in low-metallicity gas, plus standard assumptions about cosmological initial conditions.

free parameters (1)

simulation resolution
Fixed at 2 pc to resolve cluster scales; chosen rather than derived.

axioms (1)

standard math Lambda-CDM cosmology and standard hydrodynamics
Underlying framework of the MassiveBlackPS simulation.

invented entities (1)

IMBH seeds via runaway collisions no independent evidence
purpose: To link extreme cluster densities to black hole formation
Inferred from density threshold but not directly simulated or observed in the run.

pith-pipeline@v0.9.0 · 5653 in / 1335 out tokens · 63943 ms · 2026-05-16T15:41:42.091619+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

massive star clusters can form efficiently in the circumgalactic medium at early times through filament fragmentation, whereby high gas densities lead to rapid local collapse via a combination of thermal and gravitational instabilities
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

2 pc resolution... subgrid star formation and feedback prescriptions

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.