{"id":"d229c73e-3ac5-404a-a088-913cde7cf42c","arxiv_id":"2501.03323","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Raising the interstellar radiation field and cosmic-ray ionisation rate in cloud simulations makes both the core mass function and the stellar-system mass function top-heavy.","lead":"Simulations of star-forming clouds with 10 to 1,000 times the usual radiation and cosmic-ray levels produce fewer, heavier cores but richer clusters of stellar-system particles, with top-heavy mass functions. The result suggests stars may form differently in the Milky Way's centre and in starburst galaxies than in the solar neighbourhood.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The top-heavy CMF/SMF trend may not generalize to starburst clouds: the grid fixes initial density at n=10^3 cm^-3, but at CMZ-like pressures the ISRF/CRIR-driven Jeans-mass increase is offset, so the fragmentation shift could vanish.","rationale":"The paper is an honest parameter study: the fiducial run reproduces the Salpeter slope, the trends are monotonic in gamma_SFR and appear in both turbulent seeds where shown, and the authors explicitly list limitations. I do not find an internal inconsistency. The load-bearing gap is external validity. The central claim asserts a change in fragmentation in starburst environments, but the simulations hold fixed the initial density (n = 10^3 cm^-3), geometry (4.1 pc uniform sphere), and an undriven, solar-neighbourhood-like turbulent field (Section 2.2). The characteristic fragmentation mass M_J depends on both temperature and density. In a pressure-confined cloud, n is set by the external pressure, and M_J varies as T^2 / sqrt(P_ext). Real CMZ and starburst clouds have external pressures roughly 100-10,000 times the solar neighbourhood, so the Jeans mass there need not increase when the ISRF and CRIR are raised; it can remain similar to or even below the fiducial value. The high-gamma clouds in these runs develop a dense, compressive envelope, but that is a response of an isolated cloud to internal heating rather than a model of external confinement. Consequently, the top-heavy CMF and SMF may be specific to the low-density, isolated initial conditions rather than generic to high-ISRF/CRIR environments. The reader's weakest assumption identifies the same point, and the proposed density/pressure test would settle whether the result transfers. The verdict should remain CONDITIONAL, as the internal result is plausible but the generalisation to starbursts is not yet established.","tokens_in":24615,"tokens_out":12387,"duration_ms":131874,"concrete_test":"Rerun the gamma1 and gamma1000 cases with initial densities of n0 = 10^4 and 10^5 cm^-3 (or equivalently with a background pressure P/k of roughly 10^6-10^8 K cm^-3), keeping the virial parameter, Mach number, and turbulent seeds fixed. If the high-gamma CMF and SMF no longer become top-heavy as the initial density increases, the conclusion does not transfer to pressure-confined starburst/CMZ clouds.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is phrased as a statement about starburst environments, but the simulations vary only the ISRF and CRIR while fixing initial density, geometry, and an undriven, solar-neighbourhood-like turbulent velocity field (Section 2.2, Table 1). The paper itself notes that these initial conditions are 'more typical of the solar neighbourhood than of clouds in the CMZ or starbursts', yet the abstract and conclusions apply the result to those environments. The load-bearing assumption is that the ISRF/CRIR increase alone controls the fragmentation scale. Quantitatively, the Jeans mass scales as M_J proportional to T^(3/2) n^(-1/2); for a cloud in pressure equilibrium with an external medium, n approximately P_ext/(kT), so M_J proportional to T^2 P_ext^(-1/2). CMZ and starburst clouds have external pressures orders of magnitude above the solar neighbourhood, so even with the elevated temperatures produced by a 1000x ISRF/CRIR, the absolute Jeans mass can be comparable to or smaller than in the fiducial run. The dense envelope that compresses the high-gamma clouds (Section 3.1) is a consequence of heating in an isolated, low-pressure setup, not of realistic external confinement. Thus the reported top-heavy CMF and SMF may be a property of heated low-density clouds rather than a universal feature of high-ISRF/CRIR environments. This is an external-validity gap, not an internal inconsistency, but it is central to the paper's stated goal of describing star formation in starbursts.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents adaptive-mesh-refinement hydrodynamical simulations of isolated, virialised molecular clouds in which the interstellar radiation field (ISRF) and cosmic-ray ionisation rate (CRIR) are increased together by factors of 1, 10, 100 and 1000 relative to solar-neighbourhood values. The simulations use a modified version of arepo with the sgchem astrochemistry network, sink particles with 180 AU interaction radii, and two turbulent seeds per parameter choice, although most analysis is shown for one seed. The authors find that higher ISRF/CRIR values heat the gas and dust, delay the transition to cold molecular gas, and increase the local Jeans mass. This produces fewer but more massive cores, a bottom-light and top-heavy core mass function, richer clusters of sink particles per core, enhanced but less competitive accretion, and a top-heavy system mass function relative to the fiducial run. The fiducial run reproduces a Salpeter-like high-mass tail, which the authors use as an internal consistency check. Runs with a column-density-dependent cosmic-ray attenuation are also presented and show little effect at solar-neighbourhood values but a notable effect on the high-ISRF/CRIR core mass function. The paper concludes that star formation in high-SFR environments differs from that in the solar neighbourhood.","tokens_in":24982,"tokens_out":5076,"duration_ms":53666,"significance":"If the reported trend is robust, the paper provides a useful step toward understanding whether the stellar initial mass function and core mass function are universal, and it connects to observations of top-heavy mass functions in the Galactic centre and starburst regions. Its strengths are the relatively comprehensive treatment of ISRF and CRIR variations, the internal check provided by the fiducial Salpeter-like tail, the inclusion of two turbulent seeds for at least some diagnostics, and the public availability of analysis code and simulation snapshots on request. The interpretation is careful to distinguish stellar systems from individual stars and to avoid overclaiming an IMF result. However, the generality of the central conclusion is limited by the solar-neighbourhood-like initial conditions, as discussed below, and by the limited treatment of stochastic variations and unresolved substructure.","major_comments":[{"comment":"The central claim that high ISRF/CRIR environments produce top-heavy CMFs and SMFs rests on simulations that vary only gamma_SFR while keeping the initial density, geometry, and turbulent velocity dispersion fixed at solar-neighbourhood-like values. The manuscript itself states in §2.2 that these initial conditions are \"more typical of the solar neighbourhood than of clouds in the CMZ or starbursts,\" yet the abstract and conclusions generalize the result to starbursts. For gas in pressure equilibrium with an external medium, the Jeans mass in Eq. (B1) scales as M_J ∝ T^(3/2) n^(-1/2) ∝ T^2 P_ext^(-1/2). CMZ and starburst clouds have external pressures orders of magnitude above the solar neighbourhood, so the higher temperatures produced by a 1000x ISRF/CRIR do not necessarily translate into larger absolute Jeans masses. The compression of the high-gamma clouds by their heated outer envelope, described in §3.1, is an internal effect of heating in an isolated, low-pressure setup rather than a realistic external confinement. This is an external-validity gap rather than an internal inconsistency, but it is load-bearing for the paper's stated goal. I recommend either restricting the conclusions to the regime actually simulated, or adding simulations with pressure-matched initial conditions (e.g., higher initial density or an external pressure term) to demonstrate how the top-heavy trend depends on P_ext at fixed gamma_SFR.","section":"§2.2, §3.1, Eq. (B1)"},{"comment":"The quantitative mass-function slopes and the power-law exponents in Table 2 appear to be based on simulations with one turbulent seed for the main figures. Section 2.2 notes that results are usually shown for one seed unless strongly affected by the seed, and Figure 6 shows that the sink formation histories differ substantially between the two seeds. A trend in alpha_Fit based on a single seed per gamma value cannot be assigned a robust uncertainty, and it is possible that the ordering of the slopes could change with a different seed. Please report both seeds for the SMF and CMF fits, or explicitly state which seed is used and provide a quantitative estimate of seed-to-seed scatter, for example by showing the fit exponents from both runs in Table 2.","section":"§2.2, Figures 4 and 8, Table 2"},{"comment":"The sink particles have an interaction radius of 180 AU and therefore represent stellar systems rather than individual stars, and protostellar discs are unresolved. The authors are appropriately cautious about not claiming an IMF result, but the interpretation in §5.1 and §7.2 that cores fragment into richer clusters of sinks depends on the assumption that fragmentation below the sink scale would not alter the multiplicity or the resulting system mass function. The manuscript itself acknowledges in §7.4 that some fragmentation between sink insertion and optically thick core formation may be missed. Given that the new result is a shift in the SMF and CMF, this resolution caveat should be elevated from a limitation to a tested assumption, for example by a resolution study or a sub-resolution model of disc fragmentation, or the claims about cluster richness should be softened.","section":"§2.1.1, §5.1, §7.4"}],"minor_comments":[{"comment":"The caption states that outliers are omitted but does not define the outlier criterion; please specify the interquartile range or percentile rule used.","section":"Figure 7 caption"},{"comment":"The sentence \"Fragmentation only slows when the clouds become isothermal\" is ambiguous, since the clouds are roughly isothermal at high densities; rewording to \"when the gas becomes isothermal again\" or similar would clarify the point.","section":"§3.2"},{"comment":"The initial H2 abundance for gamma10 is listed as 0.363, which is higher than the fiducial value; this is plausible but worth a sentence in the text explaining why the equilibrium abundance at n=10^3 cm^-3 is not monotonic in gamma_SFR.","section":"Table A1"},{"comment":"Equation (C1) is written with a piecewise definition, but the middle line appears to be missing an exponentiation operator or parentheses; please check the typesetting so the functional form is unambiguous.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid simulation study with a clear internal result, and the external-validity issue identified in the stress-test note is genuine but addressable. I would not reject the manuscript, because the fiducial Salpeter check and the systematic variation of ISRF/CRIR provide a useful basis even if the conclusions need to be narrowed. The main concern for the editor is scope: the abstract and conclusions currently overstate the applicability to starburst and CMZ clouds given the solar-neighbourhood-like initial conditions. A major revision that either reframes the claims or adds pressure-matched tests would make the paper suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe short version: this is a solid parameter study showing that raising the ISRF and CRIR together shifts both the core mass function and the stellar-system mass function top-heavy in idealized low-density clouds. The two-scale fragmentation picture is coherent and new in this combination. But the leap to real starburst environments is only partly earned, because the initial conditions are solar-neighbourhood-like and the high external pressure of CMZ/starburst clouds is not modelled.\n\nWhat is actually new: previous grids varied the ISRF (Guszejnov et al. 2022, 2023) without the CRIR; here both are scaled together, cosmic-ray attenuation runs are included, and the same simulations are analysed for dendrogram cores and sink-particle systems. The fiducial run reproduces a Salpeter-like tail, which is a good internal consistency check. The interpretation—fewer, more massive cores that then fragment more richly into stellar systems, with weakened competitive accretion—is supported by the accretion-rate diagnostics. The paper is also honest: it says plainly that the sinks represent systems, not stars; that discs are unresolved; that feedback and magnetic fields are omitted; and that the initial density, geometry and velocity dispersion are more typical of the solar neighbourhood than of the CMZ or starbursts. The analysis code is on GitHub, which helps.\n\nThe soft spot that matters is external validity. The clouds are 4.1 pc uniform spheres at n = 10^3 cm^-3 in an isolated box, with undriven turbulence. Real CMZ and starburst clouds are under external pressures orders of magnitude higher. For pressure-confined clouds, M_J ∝ T^2 P_ext^{-1/2}, so the elevated temperatures from a 1000x ISRF/CRIR can be offset by the pressure term. The dense envelope that compresses the γ1000 cloud is generated by heating in the low-pressure setup, not by a realistic surrounding medium. So the reported top-heavy mass functions might be a property of heated low-density clouds, not a universal feature of high-ISRF/CRIR environments. This is not an internal inconsistency; the simulations show what they show. But the abstract and conclusions overreach slightly when they apply the result to starbursts without qualification.\n\nMinor issues: only two turbulent seeds, with most figures showing one; the 180 AU sinks are a real limit for IMF-style claims (the authors phrase it as 'system mass function', which is fair); and the snapshot data is on request, though the analysis code is public.\n\nWho it is for: anyone working on IMF universality, CMZ star formation, or numerical star formation in extreme environments. It deserves a serious referee; the referee should push on the initial-condition and external-pressure point and ask for more seeds or a clearly scoped claim. I would take it for review, expecting revisions.\n\nBest,","headline":"A careful parameter study showing top-heavy CMF/SMF under high ISRF/CRIR in low-density isolated clouds, but the starburst generalization is only partly earned because external pressure is not modelled.","tokens_in":25548,"tokens_out":3297,"would_cite":true,"duration_ms":29373,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Starburst radiation yields fewer heavy cores and richer star clusters","keywords":["star formation","molecular clouds","starburst galaxies","interstellar radiation field","cosmic ray ionisation","initial mass function","core mass function","numerical simulations"],"falsifier":"A resolved core and protostellar census of a high-radiation starburst cloud, for example in the Central Molecular Zone, that found a standard power-law core mass function with abundant sub-solar cores and a steep system mass function would directly contradict the predicted top-heavy shift.","tokens_in":1842,"feed_emoji":"🌟","tokens_out":3179,"duration_ms":103240,"temperature":0.7,"pith_summary":"This paper argues that star formation in starburst-like environments is not a scaled-up version of star formation in the solar neighbourhood. Using hydrodynamical simulations with a coupled chemical network, the authors raise the interstellar radiation field and cosmic ray ionisation rate in steps of 10, 100, and 1000 times the local values and follow the collapse of a ten-thousand-solar-mass cloud. The extra heating raises the Jeans mass of the gas, so the cloud fragments into fewer, more massive cores; those cores then fragment into richer groups of stellar systems that accrete faster and less competitively. The result is a core mass function and a stellar system mass function that are both bottom-light and top-heavy, deviating from the standard initial mass function slope that the fiducial run reproduces. If correct, the environments where most cosmic star formation actually occurred would produce systematically different stellar populations from those seen in the Milky Way's neighbourhood.","feed_headline":"Starburst radiation yields fewer heavy cores and richer star clusters","feed_subtitle":"Hydrodynamic simulations show high-radiation environments shift core and system mass functions to top-heavy shapes.","key_machinery":"The mechanism is the thermodynamic response of the gas to elevated heating, expressed through the Jeans mass (the minimum mass a region needs to overcome thermal support and collapse) and the effective Mach number of the turbulence. A stronger interstellar radiation field and cosmic ray ionisation rate heat the gas at all densities, raising the Jeans mass and making turbulent shocks less efficient at generating density structure; delayed dust-gas coupling keeps the gas warm to higher densities, and cosmic ray heating displaces photoelectric heating as the dominant heat source in the most extreme runs. These changes set the mass scale on which the cloud fragments into cores, and the larger core masses in turn set up richer fragmentation into stellar systems during collapse. The sink particles, which represent stellar systems rather than individual stars, allow the simulation to report a stellar system mass function that can be compared with observed initial mass functions.","core_discovery":"The central claim is that raising the interstellar radiation field and cosmic ray ionisation rate together changes fragmentation in opposite senses on two scales. On the scale of cores and clumps, the warmer, higher-pressure gas resists shock compression and has a larger Jeans mass, so fewer cores form and the ones that do are heavier, shifting the core mass function to higher masses and suppressing sub-solar cores. On the scale of stellar systems, those more massive cores are more Jeans-unstable as they collapse and fragment into larger groups of sink particles, which grow rapidly through enhanced, less competitive accretion from a plentiful reservoir. The net effect is that both the core mass function and the system mass function become top-heavy, with the high-mass slope flattening from roughly the canonical power law at fiducial conditions to a distinctly shallower slope at the highest irradiation. The paper interprets this as a picture where high-$\\gamma$ clouds fragment less on the scale of cores and clumps but more on the scale of stellar systems.","pith_inferences":["If the system-level top-heaviness survives when individual stars are resolved, the stellar initial mass function itself may be non-universal, which would change estimates of the stellar and metal content contributed by the galaxies that dominate cosmic star formation.","The combined raising of ISRF and CRIR, rather than either parameter alone, appears to be what produces the strongest effect; a natural next step would be to run a similar grid with denser, more compact initial conditions to see whether the top-heavy trend persists under Galactic-centre-like densities.","An observable prediction that could be tested with current interferometers is that dense-core surveys toward the Galactic centre should find fewer, more massive clumps per unit mass than in the solar neighbourhood, together with a deficit of low-mass cores.","The reduced rate of overall sink formation in high-$\\gamma$ clouds suggests that star formation in such environments may be spread over a longer timescale or proceed more in bursts, which could affect interpretations of the star formation efficiency in starburst galaxies."],"forward_implications":["The peak of the core mass function shifts upward by roughly an order of magnitude between the fiducial and the most extreme runs, while the high-mass tail of the system mass function flattens from a power-law slope near $-1.25$ to near $-0.70$.","Sink formation is delayed and the overall sink formation rate decreases in high-$\\gamma$ clouds, because the gas is more stable against collapse on large scales even though individual cores produce richer clusters.","Cores in the extreme runs fragment into significantly richer embedded clusters, with median distances to the tenth nearest neighbour falling below one Jeans length, implying ten or more stellar systems per core.","Including cosmic ray attenuation barely changes the fiducial cloud but substantially restores low-mass core formation in the extreme cloud, showing that cosmic ray heating is a regulator of the core mass function."],"supporting_citations":[{"why":"Provides the moving-mesh hydrodynamics code on which all simulations are run.","marker":"Springel (2010)"},{"why":"Supplies the sink particle implementation and the updated astrochemistry modifications to the code.","marker":"Hunter et al. (2023)"},{"why":"Supplies the chemical network that the simulations' astrochemistry model is based on.","marker":"Gong et al. (2017)"},{"why":"Gives the observational evidence for a shallow, top-heavy core mass function in a Galactic bar starburst that the simulations aim to reproduce.","marker":"Motte et al. (2018)"},{"why":"The closest prior simulation study varying the ISRF, whose weaker IMF shift and differing star formation rate are compared with these results.","marker":"Guszejnov et al. (2022)"},{"why":"Earlier starburst-environment simulations with a polytropic equation of state that already found a top-heavy IMF.","marker":"Klessen et al. (2007)"},{"why":"Predicts a raised minimum star-forming mass under stronger background radiation, supporting the bottom-light system mass function found here.","marker":"Whitworth et al. (2024)"},{"why":"Justifies the high ISRF and CRIR values adopted for the Galactic centre runs.","marker":"Clark et al. (2013)"},{"why":"Defines the canonical stellar mass function power law that the fiducial run must match for the high-$\\gamma$ deviations to be meaningful.","marker":"Salpeter (1955)"}],"fun_headline_variants":["Starburst clouds fragment less into cores, more into clusters","High radiation flips fragmentation: fewer cores, heavier clusters","Starburst radiation makes fewer, heavier cores and richer clusters","Fragmentation splits: starburst zones yield massive cores, rich clusters","Fewer cores, richer clusters: high-radiation clouds shift mass functions"],"cache_read_input_tokens":27520,"weakest_assumption_plain":"The work assumes that a uniform, low-density, virialised spherical cloud with solar-neighbourhood-like turbulence stands in for the clouds of the Galactic centre and starburst galaxies, an assumption the authors explicitly flag as better matched to the solar neighbourhood than to those extreme environments.","fun_headline_variants_meta":{"raw":{"variants":["Starburst clouds fragment less into cores, more into clusters","High radiation flips fragmentation: fewer cores, heavier clusters","Starburst radiation makes fewer, heavier cores and richer clusters","Fragmentation splits: starburst zones yield massive cores, rich clusters","Fewer cores, richer clusters: high-radiation clouds shift mass functions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000673,"raw_usage":{"total_tokens":3083,"prompt_tokens":985,"completion_tokens":2098,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":601,"completion_tokens_details":{"reasoning_tokens":2008}},"tokens_in":601,"tokens_out":2098,"duration_ms":12734,"temperature":1.0,"reasoning_tokens":2008,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:52:28.125739+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A resolved core and protostellar census of a high-radiation starburst cloud, for example in the Central Molecular Zone, that found a standard power-law core mass function with abundant sub-solar cores and a steep system mass function would directly contradict the predicted top-heavy shift.","supporting_citations":[],"review_version":1}