{"id":"1fa0d70b-8849-431b-8d52-4568b2471390","arxiv_id":"2607.23274","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"ALMA resolves the hot core MM1 in IRAS 17233-3606 into 11 hot molecular fragments whose mean separation (1.8e3 au) is about half the thermal Jeans length (3.3e3 au), implying fragmentation by thermal instability and subsequent gravitational contraction.","lead":"Using ALMA 1.3 mm data, the authors resolve the hot core MM1 of the massive protocluster IRAS 17233-3606 into 11 hot molecular fragments and sort them into four evolutionary phases. The fragment separations are about half the thermal Jeans length, which they interpret as thermal-Jeans fragmentation followed by global gravitational contraction and accretion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central Jeans-fragmentation argument is internally inconsistent: homologous contraction from a thermal-Jeans spacing raises d/λ_J, so the observed d/λ_J≈0.55 cannot support both thermal fragmentation and global contraction.","rationale":"The observational dataset is valuable: 11 HMFs are credibly extracted, with hot CH3CN/CH3OH temperatures, high column densities, and detailed outflow/maser/HII associations. The conditional verdict is appropriate. However, the most load-bearing weakness is not simply the projection or density choice; it is the internal logic of Section 4.1. If a clump fragments on the thermal Jeans scale and then contracts, the separation and the Jeans length evolve with different powers of density, so d/λ_J increases, not decreases. Thus the observed sub-Jeans mean spacing, combined with a claimed global contraction, is in tension with the claim that the fragmentation was initially thermal-Jeans. The virial parameter arithmetic discrepancy (0.84 vs ~1.4) adds a second concrete correctness issue in the same dynamical argument. These issues do not invalidate the observational detections or the evolutionary-phase classification, but they require either a corrected derivation or a softened interpretation, so the paper should remain conditional rather than be accepted as is.","tokens_in":21891,"tokens_out":14537,"duration_ms":133715,"concrete_test":"Run a Monte Carlo/homologous-contraction check. Initialize a uniform isothermal sphere at T=100 K, n=1.7e7 cm^-3, and place 11 cores with separations drawn from a thermal-Jeans-spaced distribution (mean ≈ λ_J = 3.3e3 au). Contract the sphere isotropically in steps (ρ → fρ, R → f^(−1/3)R, f = 1.1, ..., 10) and compute d_MST/λ_J at each step. If the ratio is ≥1 for all f > 1, the observed d_MST/λ_J ≈ 0.55 cannot arise from contraction of a thermal-Jeans configuration, and the Sec. 4.1 interpretation should be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest quantitative pillar of the evolutionary scenario (Sec. 4.1) is d_MST ≈ 1.8e3 au < λ_th^J ≈ 3.3e3 au, used to infer both 'fragmentation dominated by thermal Jeans instability' and subsequent global gravitational contraction. This inference has a logical problem independent of projection or the adopted n(H2). In a uniform isothermal medium, modes with wavelength below λ_J are pressure-stabilized; thermal Jeans fragmentation produces separations of order λ_J, not 0.55 λ_J. Moreover, if the clump contracts, density increases, so λ_J ∝ ρ^(−1/2) shrinks faster than the separation d ∝ ρ^(−1/3); the ratio d/λ_J rises monotonically under homologous contraction. Starting from a thermal-Jeans spacing (d≈λ_J), contraction cannot drive d/λ_J below 1; the observed 0.55 therefore points either to initially sub-Jeans separations (i.e., not thermal-Jeans fragmentation) or to density/temperature evolution not captured by the current single value n=1.7e7 cm^-3, T=100 K. The paper offers no mechanism to reconcile these. Separately, the reported α_vir=0.84 does not follow from the stated inputs: α_vir=5σ_tot^2 R_eff/(G M_gas) ≈ 1.4 with σ_tot=2.0 km/s, R_eff=0.025 pc, M_gas=81.3 M_sun. This makes the dynamical part of the central claim doubly shaky.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents ALMA Band 3 (ATOMS) and Band 6 (QUARKS) observations of the massive protocluster IRAS 17233-3606. The authors identify 11 hot molecular fragments (HMFs) in the MM1 hot core, derive temperatures, masses, and column densities via LTE XCLASS fitting with explicit treatment of blended transitions, and assign evolutionary phases based on masers, outflows, HII regions, and CH3CN/CH3OH abundance ratios. The central interpretation is that the mean MST separation of the HMFs (~1.8e3 au), being about half the thermal Jeans length (~3.3e3 au), together with Q = 0.77 and alpha_vir = 0.84, indicates fragmentation initially driven by thermal Jeans instability followed by global gravitational contraction and active accretion. The paper also discusses feedback from a B2 ZAMS star and the UC HII region MM2.","tokens_in":22308,"tokens_out":6933,"duration_ms":60814,"significance":"The observational material is valuable, and the XCLASS fitting procedure is methodologically careful. The authors also appropriately hedge the evolutionary sequence, admitting that the abundance-ratio clock is not monotonic and that the phase assignment is tentative. If the central fragmentation scenario were correct, the paper would provide a useful case study of high-mass star formation on ~1000 au scales. However, the quantitative dynamical argument contains a logical inconsistency and a numerical error: the claimed direction of d/lambda_J evolution under contraction is contradicted by standard homologous contraction, and the reported alpha_vir value does not follow from the stated inputs. These issues weaken the paper's principal conclusion and require substantial revision.","major_comments":[{"comment":"The inference that d_MST < lambda_J^th implies thermal Jeans fragmentation followed by contraction is internally inconsistent. In a uniform isothermal medium, only perturbations with wavelength > lambda_J grow, so thermal Jeans fragmentation produces separations of order lambda_J, not ~0.5 lambda_J. Moreover, under homologous contraction d ∝ R and lambda_J ∝ rho^(-1/2) ∝ R^(3/2), giving d/lambda_J ∝ R^(-1/2); this ratio increases as the region contracts. Starting from d ≈ lambda_J, contraction cannot produce d/lambda_J ≈ 0.55. The observed ratio instead suggests either initially sub-Jeans separations (e.g., other fragmentation mechanisms) or that the adopted current n(H2)=1.7e7 cm^-3 and T=100 K do not represent the pre-fragmentation conditions. Please provide a quantitative evolutionary model or revise the conclusion.","section":"Sec. 4.1"},{"comment":"Using the stated inputs (sigma_tot = 2.0 km/s, R_eff = 0.025 pc, M_gas = 81.3 M_sun), the formula alpha_vir = 5 sigma_tot^2 R_eff / (G M_gas) yields alpha_vir ≈ 1.4, not 0.84. The reported value appears to correspond to 3 sigma_tot^2 R_eff / (G M_gas). Please correct the numerical value or the formula. The corrected value is still below the usual critical threshold of ~2, but the quantitative claim and the abstract need revision.","section":"Sec. 4.1, alpha_vir"},{"comment":"The statement that 'even accounting for the projection effect, the deprojected mean separation remains smaller than lambda_J^th' is not supported by any deprojection model. The MST separation is a projected quantity, and the factor relating projected to 3D separations depends on geometry. Please provide the assumed deprojection, or remove the claim and treat d_MST as a lower limit with a stated uncertainty.","section":"Sec. 4.1, deprojection"}],"minor_comments":[{"comment":"Grammatical error: 'none of them was observed' should be 'none of them were observed'.","section":"Sec. 3.3"},{"comment":"Typo: 'forulated' should be 'formulated'.","section":"Figure 5 caption"},{"comment":"Grammatical error: 'The result provide' should be 'The results provide'.","section":"Sec. 4.1"},{"comment":"The value L_av ≈ 4000 au used in the Q parameter is introduced without a precise definition. Please specify how the mean separation length between all HMFs is computed and cite the relevant method.","section":"Sec. 4.1, Q parameter"},{"comment":"The adopted volume density n(H2) ≈ 1.7e7 cm^-3 is taken from Chen et al. (2025) without discussing its uncertainty or derivation. Since the Jeans length comparison is sensitive to this value, please provide the uncertainty and justify its application to MM1.","section":"Sec. 4.1, adopted density"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is based on genuine ALMA observations and the line-fitting methodology is sound. The main issues are interpretive: the d/lambda_J argument is logically flawed as written, and the alpha_vir value is an arithmetic error that should have been caught. I have no concerns about data provenance or citation behavior. If the authors can correct the numerical error and either provide a quantitative evolutionary model or substantially qualify the fragmentation claim, the paper may be publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real contribution here is observational: first 1.3 mm resolution of MM1 into 11 hot molecular fragments, with XCLASS rotational temperatures (100-310 K), column densities, and a phase I-IV classification tied to outflows, masers, and H II regions. That part is careful and credible, a genuine step beyond the earlier lower-resolution work. The paper also does several things well: it states how blended CH3OH transitions were excluded, checks the continuum optical depth, and explicitly hedges the abundance-ratio clock and the outflow-driving-source assignments. Honest reporting. The soft spots are in the dynamical interpretation, and they are not minor. The central claim, thermal Jeans fragmentation followed by global contraction, is anchored on d_MST about 1.8e3 au < lambda_th about 3.3e3 au, plus Q = 0.77 and alpha_vir = 0.84. The Jeans argument has an internal consistency problem the authors don't address: thermal Jeans fragmentation separates fragments by roughly lambda_J, and homologous contraction raises d/lambda_J (d shrinks as rho^-1/3 while lambda_J shrinks as rho^-1/2). Starting near lambda_J, contraction cannot drive d/lambda_J below 1. Observed d/lambda_J about 0.55 therefore points either to initially sub-Jeans spacing (i.e., not thermal fragmentation) or to large temperature evolution not captured by the single current T and rho. The paper offers no mechanism. Also, the reported alpha_vir = 0.84 does not come out of the formula as written: 5 sigma^2 R_eff / (G M_gas) with sigma=2.0 km/s, R_eff=0.025 pc, M=81.3 M_sun gives about 1.4. Still below 2, so the global-contraction conclusion survives qualitatively, but the number needs correcting or the inputs need explaining. Projection deprojection is also hand-waved, though that is the lesser issue. Overall: the data deserve a serious referee. The source study will be useful to HMSF observers regardless of the interpretive fix. But the fragmentation section should be rewritten, not just patched, and the virial discrepancy resolved before acceptance. I'd send it for peer review with a request for major revision on those two points.","headline":"Solid ALMA core-resolution paper whose Jeans-fragmentation narrative doesn't survive contact with its own numbers.","tokens_in":626,"tokens_out":845,"would_cite":true,"duration_ms":40922,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The 11 hot molecular fragments inside the massive core MM1 are spaced at about half the thermal Jeans length, indicating that thermal instability set their initial fragmentation, followed by global gravitational contraction and active accre","keywords":["hot molecular cores","fragmentation","Jeans length","minimum spanning tree","massive star formation","protocluster","IRAS 17233-3606","thermal instability"],"falsifier":"Measure the 3D separations of the 11 HMFs using kinematic distances from molecular-line gradients or proper motions, and check if the mean deprojected separation exceeds ~3.3e3 au. Alternatively, derive the initial gas density and temperature before fragmentation (e.g., from the dust emission of the parent core on larger scales) and recompute the thermal Jeans length; if that length is smaller than the observed mean separation, thermal instability alone cannot explain the fragmentation scale.","tokens_in":21818,"feed_emoji":"🔭","tokens_out":6939,"duration_ms":61333,"temperature":0.7,"pith_summary":"The paper resolves the massive hot core MM1 in the protocluster IRAS 17233-3606 into 11 hot molecular fragments and measures their temperatures, densities, motions, and evolutionary states. Its central claim is that the typical spacing of these fragments, about 1,800 astronomical units, is nearly half the thermal Jeans length of the core, implying that the fragmentation was initially set by thermal instability rather than by turbulence. The accompanying Q parameter (0.77) and virial parameter (0.84) indicate a subclustered spatial distribution and ongoing global gravitational contraction, so the fragments are now being drawn together as they accrete. The fragments span four evolutionary phases, from quiescent to HII-region-associated, showing that star formation in this protocluster is asynchronous. The result connects an observable spacing to a specific physical fragmentation mechanism in massive star formation.","feed_headline":"Thermal instability, not turbulence, spaced a hot core's 11 fragments","feed_subtitle":"The group sits at half the thermal Jeans length and is now contracting, pointing to how massive protoclusters begin.","key_machinery":"Minimum-spanning-tree (MST) mean separation, the thermal Jeans length, the Q parameter, and the virial parameter. The MST mean separation quantifies the typical projected distance between the 11 HMFs; comparing it with the thermal Jeans length tests whether the fragment spacing matches the scale of thermal gravitational instability. The Q parameter distinguishes a centrally condensed cluster (Q > 0.8) from a subclustered one, and the virial parameter measures whether the region is gravitationally bound and contracting.","core_discovery":"Using millimeter observations at ~0.3 arcsecond resolution, the paper identifies 11 hot molecular fragments (HMFs) within the previously known hot core MM1, with rotational temperatures of 100-310 K and H2 column densities above 10^23 cm^-2, all capable of forming massive stars. The mean minimum-spanning-tree separation of the HMFs is ~1.8e3 au, about half the thermal Jeans length (~3.3e3 au) calculated from the core's mean density of n(H2)=1.7e7 cm^-3 and a rotational temperature of ~100 K. Because the observed spacing is below the thermal Jeans length, and because the Q parameter is 0.77 (subclustered, near the 0.8 threshold) and the virial parameter is 0.84 (gravitationally bound), the pa","pith_inferences":["If the thermal fragmentation scenario holds, the fragment masses should follow a thermal Jeans mass distribution; a direct test would be to compare the measured HMF masses with the local Jeans masses calculated at each fragment's position.","The conclusion depends on treating the projected MST separation as the physical spacing; high-resolution molecular-line data (e.g., CH3CN velocity gradients) could yield kinematic distances and deprojected separations to determine whether the true spacing is also sub-Jeans.","The phase I-IV gradient suggests that mixed-age core populations may be common in massive hot cores; if so, chemical clocks must be calibrated against spatial and dynamical indicators rather than assumed single-burst evolution.","MM2-2, a cold core coincident with an ultracompact HII region, predicts that shielded 'island' cores can survive feedback; such objects could be searched for in other UC HII regions to test the shielding hypothesis."],"forward_implications":["MM1's fragments are forming stars asynchronously: 11 HMFs coexist in a 0.1 pc region and span evolutionary phases I-IV, so massive protoclusters need not form in a synchronized burst.","The spacing result implies that the initial fragmentation of this massive hot core was thermal, meaning models that rely primarily on turbulent fragmentation should be revisited for cores at this stage.","With alpha_vir = 0.84, the region is globally contracting; the fragments should continue to converge and accrete, which may lead to core coalescence or competitive accretion.","The f[CH3CN/CH3OH] ratio generally increases from phase I to IV (within uncertainties), suggesting hot-core chemistry is established before an HII region appears.","The SiO arc at the interface with the B2 ZAMS star shows stellar feedback is actively shaping the molecular gas of MM1, while MM2-2 demonstrates that dense cores can remain cold and unperturbed inside an evolved HII region."],"fun_headline_variants":["Hot core's 11 fragments sit at half Jeans length, hinting thermal start","Thermal instability, not turbulence, spaced a massive core's 11 fragments","Half-Jeans spacing of 11 hot fragments: thermal fragmentation then collapse","ALMA resolves 11 massive-star seeds in one hot core, spaced by thermal physics"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The argument treats the projected on-sky separation of the 11 fragments as a faithful measure of their physical spacing and compares it with a thermal Jeans length computed from the current mean density and temperature; if the true 3D separations are significantly larger than the projected values, or if the pre-fragmentation density and temperature differed from the adopted values, the observed spacing would no longer implicate thermal instability.","fun_headline_variants_meta":{"raw":{"variants":["Hot core's 11 fragments sit at half Jeans length, hinting thermal start","Thermal instability, not turbulence, spaced a massive core's 11 fragments","Half-Jeans spacing of 11 hot fragments: thermal fragmentation then collapse","ALMA resolves 11 massive-star seeds in one hot core, spaced by thermal physics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001927,"raw_usage":{"total_tokens":7478,"prompt_tokens":941,"completion_tokens":6537,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":685,"completion_tokens_details":{"reasoning_tokens":6459}},"tokens_in":685,"tokens_out":6537,"duration_ms":38132,"temperature":1.0,"reasoning_tokens":6459,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T23:51:51.389546+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the 3D separations of the 11 HMFs using kinematic distances from molecular-line gradients or proper motions, and check if the mean deprojected separation exceeds ~3.3e3 au. Alternatively, derive the initial gas density and temperature before fragmentation (e.g., from the dust emission of the parent core on larger scales) and recompute the thermal Jeans length; if that length is smaller than the observed mean separation, thermal instability alone cannot explain the fragmentation scale.","supporting_citations":[],"review_version":1}