arxiv: 2512.10298 · v3 · submitted 2025-12-11 · 🌌 astro-ph.GA

Recognition: 2 theorem links

· Lean Theorem

Tracking Protostellar Variability in Massive Protoclusters with ALMA: I. Insights from QUARKS and MaMMOtH

Yuhan Yang , Tie Liu , Sheng-Yuan Liu , Doug Johnstone , Gregory Herczeg , Wenyu Jiao , Yu-Nung Su , Xiaofeng Mai

show 17 more authors

Fengwei Xu Dominique Meyer Siju Zhang Eduard Vorobyov Suinan Zhang Qiuyi Luo Guido Garay Xi Chen Yunfan Jiao Qi-lao Gu Yan-kun Zhang Ken'ichi Tatematsu Kee-Tae Kim Andrey Sobolev Sergey Parfenov Dmitry A. Ladeyschikov Leonardo Bronfman

Authors on Pith no claims yet

Pith reviewed 2026-05-16 23:45 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords protostellar variabilityALMA continuummassive protoclustersmillimeter observationshypercompact HII regionsvariable fraction

0 comments

The pith

ALMA observations of 22 massive protoclusters detect millimeter variability in five of 383 condensations, including a 68% intensity increase in one source over one year.

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

The paper uses repeated ALMA Band 6 observations of 22 massive protoclusters, spaced from hours to more than two years apart, to search for changes in the millimeter continuum emission from embedded protostars. A custom pipeline aligns the images and calibrates fluxes to high precision, after which the astrodendro algorithm extracts 383 condensations whose peak intensities are tracked across epochs. Standard deviation maps and difference images flag five variable sources, setting a lower limit of 1.3 percent on the variable fraction in this environment. One object, I13111-6228, which contains a hypercompact H II region, brightens by 68 percent in peak intensity over a single year. The work extends variability searches to high-mass protoclusters and shows that such changes can be measured systematically with interferometric data.

Core claim

Multi-epoch ALMA Band 6 continuum imaging at 0.3 arcsec resolution of 22 massive protoclusters yields 383 tracked condensations, of which five exhibit statistically significant variability; the standout source I13111-6228 shows a 68 percent rise in continuum peak intensity over one year.

What carries the argument

Custom data-reduction pipeline for image alignment and relative flux calibration, combined with the astrodendro algorithm to identify and monitor condensations across epochs.

If this is right

Millimeter variability occurs in high-mass protoclusters but at a lower rate than reported for low-mass regions.
Hypercompact H II regions can produce large, rapid continuum changes on yearly timescales.
Interferometric multi-epoch surveys can deliver statistical samples of protostellar variability beyond nearby clouds.
The detected changes are consistent with dynamical accretion events that also operate in more massive systems.

Where Pith is reading between the lines

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

Longer time baselines or additional epochs could reveal whether the variability is periodic, episodic, or tied to specific accretion bursts.
Cross-matching the variable sources with infrared or radio data would test whether the millimeter changes correlate with other accretion tracers.
The 1.3 percent lower limit implies that larger surveys of massive clusters are needed to build a robust sample of variable objects.

Load-bearing premise

The measured intensity changes are produced by intrinsic protostellar processes rather than residual calibration or alignment errors that survive the custom processing.

What would settle it

A fresh ALMA Band 6 observation of I13111-6228 at matching resolution and frequency that shows the peak intensity has returned to the pre-increase value would falsify the reported 68 percent change.

Figures

Figures reproduced from arXiv: 2512.10298 by Andrey Sobolev, Dmitry A. Ladeyschikov, Dominique Meyer, Doug Johnstone, Eduard Vorobyov, Fengwei Xu, Gregory Herczeg, Guido Garay, Kee-Tae Kim, Ken'ichi Tatematsu, Leonardo Bronfman, Qi-lao Gu, Qiuyi Luo, Sergey Parfenov, Sheng-Yuan Liu, Siju Zhang, Suinan Zhang, Tie Liu, Wenyu Jiao, Xiaofeng Mai, Xi Chen, Yan-kun Zhang, Yuhan Yang, Yunfan Jiao, Yu-Nung Su.

**Figure 1.** Figure 1: Processing workflow. Schematic overview of the four main processing steps adopted in this study. The process begins with data selection and calibration, including dataset cross-matching and ALMA pipeline calibration to prepare the data for further steps (Sect. 2). In the second step, imaging and self-calibration, we perform spectral Line flagging and self-calibration (Sect. 3.1), followed by smoothing all … view at source ↗

**Figure 2.** Figure 2: ALMA Band 6 continuum images of I15520–5234 obtained from four epochs spanning over two years (with the observing dates indicated in the figure). The cyan dashed circle in each panel marks the 0.5 primary beam FWHM region of the MaMMOtH survey (∼13.24′′), which is also used as the reference area for image alignment. Upper panels (a, b): MaMMOtH 1.3 mm continuum emission. Lower panels (c, d): QUARKS 1.3 mm … view at source ↗

**Figure 3.** Figure 3: Top: Cartoon illustrating the extraction of structures and their fluxes. Red and green contours represent condensations independently identified in the MaMMOtH and QUARKS surveys using the astrodendro algorithm. The union of these two sets defines a common mask, shown in yellow, which is consistently applied across all observing epochs. Bottom: Application to I15520–5234. (a) 1.3 mm continuum image before … view at source ↗

**Figure 4.** Figure 4: Relative flux calibration overview. Top: distribution of flux calibration factors. Bottom: corresponding uncertainties. each epoch by this mean. We choose a threshold of 0.1 for the standard deviation of the normalized intensities as a compromise between retaining more calibrator candidates and ensuring the reliability of the calibration. Sources with standard deviations below this threshold were classi… view at source ↗

**Figure 5.** Figure 5: Normalized standard deviation of peak intensity (SD/SDfid) as a function of the mean peak intensity for condensations across 22 protoclusters. Horizontal dashed lines indicate levels at 1, 3, and 5 times the fiducial expectation (SD/SDfid = 1, 3, 5). Sources exceeding the 5 times fiducial level are highlighted in red and flagged as candidate variables. standard deviation (SD). The uncertainty in measurin… view at source ↗

**Figure 6.** Figure 6: Peak intensity ratio versus S/N at the reference epoch for I13111–6228. Each point represents an individual condensation extracted from the I13111–6228. The shaded regions denote ±1, ±3, and ±5σratio. Points with S/N greater than 30 are shown in blue, while those with S/N less than 30 are shown in green. The red star marks a deviation beyond 5σratio. are colored green, reflecting their larger noise-domina… view at source ↗

**Figure 7.** Figure 7: 1.3 mm continuum images of I13111–6228 observed with ALMA at two epochs and their difference map. (a) Image from Epoch 1 observed on 2023 May 02. (b) Image from Epoch 2 observed on 2024 June 02. (c) Difference map between panels (a) and (b), produced after aligning the two images for visual comparison. The cyan circle marks a radius of 0.75′′ centered on the residual peak in the difference map. All three p… view at source ↗

**Figure 8.** Figure 8: Integrated intensity map of H30α for I13111–6228. The black contours are at levels of [3, 4, 5]×σrms, where σrms = 0.13 Jy beam−1 km s−1 . The synthesized beam is shown in the lower-left corner of panel, and a 0.01 pc scale bar is indicated in the lower-right corner. Condensations that exhibit H30α emission were classified as Hii regions, while those without H30α but exhibiting CH3CN emission were classi… view at source ↗

**Figure 9.** Figure 9: presents a scatter plot of the fractional amplitude as a function of the mean peak intensity for all sources with valid multi-epoch millimeter measurements. The background histogram displays the distribution of mean peak intensities across the full sample, providing context for the population’s brightness. Each gray point denotes an individual condensation, while the red circles highlight the detected v… view at source ↗

**Figure 10.** Figure 10: “crowding effect” test for I13111–6228. Top panels: ALMA 1.3 mm images from 2023 May and 2024 June convolved to 5.0′′, and their difference map (right). Bottom panels: same data convolved to 10.0′′. The beam size is shown in the lower-left panel of the first column of each row. White contours in the third column indicate the union mask used for flux extraction, and the cyan circle marks the variable sourc… view at source ↗

read the original abstract

Millimeter/submillimeter variability is often attributed to dynamical disk-mediated accretion, yet detection is limited to low-mass protostars in nearby clouds. Recent observations have also revealed significant (sub)millimeter variability in high-mass protostars, but the confirmed cases are scarce and lack systematic monitoring. In this work, we analyzed multi-epoch Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 (1.3 mm) continuum observations of 22 massive protoclusters, with epoch separations ranging from a few hours to more than two years, while achieving a consistent angular resolution of approximately 0.3 arcsec. These data allow us to track variability of protostars across a broader mass range and in an environment markedly different from nearby clouds. Using a custom processing pipeline for data reduction, image alignment, and relative flux calibration, we achieve high-precision flux measurements and, for the first time, investigate millimeter variability in massive protoclusters based on interferometric data in a statistical manner. Applying the astrodendro algorithm, we identified 383 condensations and tracked their variations in peak intensities. Standard deviation analysis and difference maps reveal five variable sources, corresponding to a lower limit of 1.3% on the variable fraction. Among these, I13111-6228 stands out as it hosts a hypercompact H II region that exhibits a 68% increase in continuum peak intensity over one year, with an uncertainty of 2%.

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

This paper gives the first statistical lower limit on mm variability in massive protoclusters from ALMA data but the custom pipeline lacks validation against systematics.

read the letter

The main thing here is that they took multi-epoch ALMA Band 6 observations of 22 massive protoclusters, ran astrodendro on the images to get 383 condensations, and flagged five as variable via standard deviation and difference maps. That yields a 1.3% lower limit on the variable fraction, with one hypercompact HII region showing a 68% continuum rise over a year at 2% quoted uncertainty. It is the first interferometric statistical look at this in high-mass environments rather than the usual low-mass clouds. They deserve credit for pulling together the QUARKS and MaMMOtH datasets with consistent 0.3 arcsec resolution and doing a uniform search across a mass range that matters for Galactic star formation. The numbers are concrete and the question is worth asking. The soft spot is exactly what the stress-test note flags: the custom pipeline for reduction, alignment, and relative calibration is described only at a high level, with no error budget, no injected-signal tests, and no control-sample check that non-variables stay flat at the few-percent level. If residual phase or amplitude errors after alignment reach a few percent, difference maps can generate false positives at the reported rate, and the 2% uncertainty on the 68% change would not capture that. Without those checks the standout result is hard to trust at face value. This is for people working on high-mass accretion and variability statistics. A reader who wants a first look at how common mm changes are in protoclusters will get something from the sample size, but anyone planning to cite the fraction or the 68% number will need the full methods to hold up. I would send it to peer review because the dataset is real and the claim is testable, even though revisions on the error analysis are likely needed.

Referee Report

2 major / 2 minor

Summary. The paper analyzes multi-epoch ALMA Band 6 (1.3 mm) continuum data for 22 massive protoclusters at ~0.3 arcsec resolution. Using a custom pipeline for reduction, alignment, and relative calibration, the authors apply astrodendro to extract 383 condensations and identify five variable sources via standard-deviation cuts and difference maps, yielding a 1.3% lower limit on the variable fraction. The standout case is source I13111-6228, whose hypercompact H II region shows a 68% rise in peak intensity over one year (quoted uncertainty 2%).

Significance. If the variability detections hold after rigorous validation, the work supplies the first statistical millimeter-variability census in massive protoclusters, extending earlier low-mass results to a higher-mass regime and different environment. The large sample size (383 sources) and the dramatic 68% change in I13111-6228 are concrete observational anchors that could constrain accretion-disk models for high-mass star formation.

major comments (2)

[Methods (custom processing pipeline)] Methods / custom pipeline description: no error budget, no injected-signal simulations, and no control-sample statistics are presented to demonstrate that residual phase/amplitude errors after alignment and relative calibration remain below the quoted 2% level; this directly affects the reliability of the 68% intensity change reported for I13111-6228 and the overall 1.3% variable fraction.
[Results (standard deviation analysis and difference maps)] Results / variability detection: the standard-deviation analysis and difference-map criteria that flag the five variable sources lack a quantitative false-positive rate or Monte-Carlo validation against the noise properties of the aligned images; without this, it is unclear whether the reported lower limit of 1.3% is robust.

minor comments (2)

[Abstract] Abstract: the 2% uncertainty on the 68% change should explicitly state whether it is statistical only or includes systematic contributions from the pipeline.
[Figure 1 / text] Figure captions and text: epoch time baselines are described only qualitatively; a compact table listing the exact time separations for each cluster would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments, which help improve the rigor of our analysis. We address each major comment below and will revise the manuscript to incorporate the requested validations.

read point-by-point responses

Referee: Methods / custom pipeline description: no error budget, no injected-signal simulations, and no control-sample statistics are presented to demonstrate that residual phase/amplitude errors after alignment and relative calibration remain below the quoted 2% level; this directly affects the reliability of the 68% intensity change reported for I13111-6228 and the overall 1.3% variable fraction.

Authors: We agree that the current manuscript does not include a full error budget or quantitative validation of residual calibration errors. In the revised version we will add an explicit error budget for the alignment and relative-calibration steps, perform injected-signal Monte-Carlo simulations to measure residual phase/amplitude errors, and analyze a control sample of non-variable condensations to demonstrate that the achieved precision is at or below the quoted 2% level. These additions will directly support the reliability of the 68% intensity change in I13111-6228 and the overall variable fraction. revision: yes
Referee: Results / variability detection: the standard-deviation analysis and difference-map criteria that flag the five variable sources lack a quantitative false-positive rate or Monte-Carlo validation against the noise properties of the aligned images; without this, it is unclear whether the reported lower limit of 1.3% is robust.

Authors: We acknowledge that a quantitative false-positive assessment is required. We will add Monte-Carlo simulations that generate noise realizations matching the observed image properties, apply the identical standard-deviation and difference-map criteria to non-variable sources, and report the resulting false-positive rate. This will allow us to evaluate and, if necessary, adjust the 1.3% lower limit in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No circularity: purely observational flux measurements

full rationale

The paper reports direct ALMA continuum observations of 383 condensations across multiple epochs. Variability is identified by applying standard deviation analysis and difference maps to measured peak intensities after a custom data-reduction pipeline. No equations, predictions, or derivations are present that reduce to fitted parameters, self-definitions, or self-citations. The central results (5 variable sources, 1.3% lower limit, 68% change in I13111-6228) are empirical counts and ratios computed from the observed data values themselves. The pipeline is described as a processing step, not a theoretical model whose outputs are then re-derived as inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Observational study; no new theoretical parameters, axioms, or entities are introduced. All quantities derive from standard ALMA data reduction and source extraction assumptions already present in the cited literature.

pith-pipeline@v0.9.0 · 5684 in / 1150 out tokens · 63146 ms · 2026-05-16T23:45:17.324990+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.

Applying the astrodendro algorithm, we identified 383 condensations

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

3 extracted references · 3 canonical work pages

[1]

, " * write output.state after.block = add.period write newline

ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

work page
[2]

write newline

" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

work page
[3]

l N؟ /4'5տȯ w_ ?՟ òJ0MN M<?d>_ ?i왳?#NQH뺕#n ӵ AÌ(nkYC<fx>D׎A- ^*Xk ͥ͟ki j G\ /f ͦYĪCXԈlr3

thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

work page doi:10.3847/1538-4365/adf847 2017