{"id":"e7b56619-3d5e-4b17-a75b-9f38e71d3ba1","arxiv_id":"2507.22470","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Dense cores in Orion A that form binary or multiple systems have higher density and Mach number than cores forming single stars, while magnetic field strength does not significantly differ.","lead":"Using archival radio and submillimeter observations, this paper compares dense star-forming cores in Orion A that host one protostar versus those hosting two or more. The cores that fragment into multiple stars are denser and more turbulent, and magnetic fields appear to play little role.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Core-averaged density and Mach number are measured with beams larger than the smallest cores; the stated 'minimal impact' is unverified and could systematically drive the single-vs-multiple differences.","rationale":"The reader's weakest_assumption correctly identifies the resolution-matching issue as the most fragile link in the argument. My reading of the paper confirms that the two parameters driving the central claim—density and Mach number—are the ones most exposed to beam-smearing effects. The JCMT, Nobeyama, and Herschel beams are each larger than the smallest cores, and the paper's dismissal of resolution effects in Section 3 is a one-sentence assertion, not a demonstrated test. The concrete test I propose would settle whether the apparent group differences survive when the analysis is restricted to well-resolved cores or when the pipeline is calibrated on synthetic sources. The other concerns raised by the reader (multiple-testing correction, marginal density p-value, weak magnetic-field null) are real but secondary: they affect the strength of the statistical claim, whereas the resolution issue threatens the physical interpretation of the measured quantities themselves. I agree with the reader's CONDITIONAL verdict because the concern is substantial enough to require additional analysis, but not so severe as to warrant rejection without further testing. The paper does include some supportive checks, such as excluding very close multiples (§3.1), quantifying disk flux contamination (§3.2), and testing outflow contamination via radial line-width profiles (§5.1, Figure 8), which indicate careful work. However, none of these checks addresses the beam-vs-core-size mismatch, so the central claim remains conditional pending the resolution test.","tokens_in":22989,"tokens_out":6219,"duration_ms":77392,"concrete_test":"Run a synthetic-aperture injection test: place Gaussian model cores with known density, temperature, and velocity dispersion at the positions of the 53 sample cores in the JCMT 850 μm, Herschel temperature, and Nobeyama N2H+ data, convolve each map with its actual beam and add the observed noise, then run the same astrodendro/Gaussian-fitting and hyperfine-fitting pipeline to recover core radii, masses, densities, and Mach numbers. Compare the recovered values to the injected values as a function of core radius. If the recovery for cores with R < 6000 au shows systematic biases (e.g., Mach number inflated by beam-averaged external gas or density biased by beam-convolved radii), the reported single-vs-multiple differences in density and Mach number cannot be considered robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (abstract; §5.1) rests on two measurements: core density (Eq. 1, §3.2) and Mach number (Eqs. 4–5, §3.5). Both are derived from data whose beams are comparable to or larger than the cores. The JCMT 850 μm beam is 14.6″ (~5800 au at 400 pc), the Nobeyama N2H+ beam is 23.4″ (~9400 au), and the Herschel temperature map has 37″ resolution (~14800 au), while the sample core radii are 3200–16000 au (§3.1). For the smallest cores, e.g., radius 3200 au (index 45), the N2H+ and Herschel beams average emission over regions ~3–5× the core radius. The manuscript asserts (§3) that averaging over the cores makes resolution effects 'minimal', but this is not demonstrated. If the surrounding gas in Orion A has higher velocity dispersion or different temperature than the core, the inferred Mach number and mass (via T_d in Eq. 1) will be contaminated. Moreover, §3.1 states that cores smaller than the beam are excluded, yet several listed radii (e.g., 3200, 3400, 3700 au) are smaller than the JCMT beam FWHM (~5800 au), indicating either that the reported radii are not beam-deconvolved or the exclusion was not applied. Since the density is M/(4/3πR^3), a beam-convolved R for unresolved cores would bias density. The paper does not report a resolution-matched subsample or a deconvolution test, so the apparent higher density and Mach number in fragmented cores could be a resolution artifact rather than a physical property of the cores.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper compares the physical properties of 38 dense cores hosting single protostars and 15 dense cores hosting binary or multiple systems in Orion A, using JCMT 850 micron continuum and polarimetry, Herschel dust temperatures, Nobeyama N2H+ line data, and VANDAM multiplicity measurements. The authors derive core mass, density, Jeans length and mass, Mach number, turbulent pressure, DCF magnetic field strength, and mass-to-flux ratio, and then compare the two groups with cumulative distributions and KS tests. The central claim is that fragmented cores have significantly higher density and Mach number, and larger ratios of core radius/mass to Jeans radius/mass, while the energy ratios of turbulence and magnetic field to gravity show no significant difference. The paper interprets this as evidence that density and supersonic turbulence drive core fragmentation in Orion A and that magnetic support is not a dominant factor.","tokens_in":23343,"tokens_out":7148,"duration_ms":85268,"significance":"If the central result holds, the paper provides a valuable observational constraint on the fragmentation of low-mass cores: it supports the turbulent core fragmentation scenario and challenges a strong role for magnetic fields at the ~0.1 pc scale in Orion A. The study extends the work of Luo et al. (2022) by adding magnetic field measurements and by using the high-resolution VANDAM multiplicity catalog. Strengths include a clearly described sample selection, full tables of core properties, Monte Carlo treatment of CDF uncertainties, and a transparent statistical comparison. The main novel contribution is the comparison of the mass-to-flux ratio and turbulent-to-gravitational energy ratio between fragmented and unfragmented cores. However, the headline conclusions depend on resolution-related assumptions and on the treatment of multiple statistical tests, both of which need to be demonstrated before the result can be considered robust.","major_comments":[{"comment":"The paper states that unresolved cores smaller than the JCMT beam are excluded, yet Table 1 lists many cores with radii well below the 14.6 arcsec FWHM beam (~5800 au at 400 pc), e.g., indices 45 (3200 au), 52 (3400 au), 53 (4400 au). It is not stated whether the reported radii are beam-deconvolved. If they are not, the fitted 2D Gaussian sizes should be at least the beam's own 2-sigma width (~4900 au), so either the radii are deconvolved by an unspecified procedure or the claimed exclusion was not applied. Since the core density in Section 3.2 is computed as M/(4/3 pi R^3), a biased R for the smallest cores could directly produce the apparent density difference between the single and multiple groups. In addition, the N2H+ beam (23.4 arcsec) and Herschel temperature beam (37 arcsec) average over regions 3-5 times the radius of the smallest cores, so the core-averaged Mach number and dust temperature may be contaminated by surrounding gas. The statement in Section 3 that the resolution impact is 'minimal' is asserted rather than demonstrated. Please provide a resolution-matched subsample analysis, a synthetic beam-convolution test, or a clear description of any deconvolution procedure.","section":"§3.1–3.2, Eq. (1)"},{"comment":"The KS tests in Table 3 involve 11 parameters, and no correction for multiple comparisons is applied. The density p-value is 0.04 with an upper uncertainty of 0.15, and the Mach number p-value is 0.006; with a Bonferroni threshold of 0.0045 neither would be individually significant. The abstract and §5.1 highlight these parameters as the key evidence, so the paper should either apply a multiple-testing correction or reframe the conclusion around the more robust Jeans-ratio test (p=0.002) and the mass test (p=0.005), noting that the density and Mach number results are suggestive trends. The Spearman correlations in §4.5 and Figure 7 also involve multiple tests and are reported without adjustment.","section":"§4, Table 3"},{"comment":"Magnetic field strengths and mass-to-flux ratios are measured for only 27 of the 53 cores, with different detection fractions for the two groups: 17 of 38 single systems (45%) and 10 of 15 multiple systems (67%). Because polarization detection requires sufficiently bright and polarized emission, the magnetic subsample may be biased toward the same dense, massive cores that drive the main result, so the null KS results for B_pos (p=0.14) and lambda (p=0.64) could be a selection artifact rather than a physical absence of magnetic effects. The authors should test whether restricting the density and Mach number comparisons to the 27 cores with magnetic detections reproduces their headline differences, and should discuss the impact of the unequal detection fractions, possibly by using upper limits for the nondetections.","section":"§3.6 and §4.4, Table 2"}],"minor_comments":[{"comment":"The sentence 'we note that dense core with sizes smaller than the JCMT beam size may not be fully resolved' contains a typo ('core' should be 'cores'), and the exclusion criterion should be stated precisely, e.g., whether 'smaller than the beam' refers to the FWHM or to the 2-sigma width of the fitted Gaussian.","section":"§3.1"},{"comment":"When two velocity components are fitted, the paper selects the component with higher optical depth; the potential bias introduced by this choice should be briefly justified or tested, since the higher-optical-depth component may preferentially trace the densest gas and thus affect the derived velocity dispersion and Mach number.","section":"§3.4"},{"comment":"The angular dispersion is measured after subtracting a smoothed field with a 3x3 pixel (36 arcsec) kernel that is said to correspond to the median core size; the sensitivity of the derived B_pos and mass-to-flux ratio to the kernel choice is not discussed and could be checked, especially for the smallest cores.","section":"§3.6"},{"comment":"The discussion of mass-to-flux ratios in the context of other surveys would benefit from stating explicitly that the DCF-derived B_pos is a plane-of-sky lower limit, so the conclusion that magnetic fields are relatively weak is partly projection-dependent; the text makes this point earlier but the summary and abstract do not.","section":"§5.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and addresses a timely question. The main risk is that the headline density and Mach number result may not be robust to the resolution and multiple-testing concerns described above, while the magnetic-field conclusion rests on a small and possibly biased subsample. I recommend asking for the specific robustness tests rather than rejecting, because the data analysis is otherwise careful and the question is well motivated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — the genuinely new piece is the magnetic-field comparison: 27 Orion A cores with JCMT POL-2 polarization, split by multiplicity, with DCF field strengths and mass-to-flux ratios. That analysis is new and comes out null. It is a useful data point for the 'B-fields suppress fragmentation' debate. The density and Mach-number part is less new — Luo et al. (2022) already reported higher density and Mach number in fragmented cores on a partially overlapping sample — but this paper uses a separate core sample and a cleaner protostellar census from VANDAM, so the confirmation has value.\n\nThe paper is honestly put together. Core identification is transparent, the N2H+ hyperfine fitting is careful, and they check that outflows aren't driving the line-width excess (Figure 8). The KS testing with Monte Carlo uncertainties is standard. I also see no circularity: the PT/uG null against a Mach-number difference is a real consequence of the measured quantities.\n\nThe soft spots are the usual ones. No multiple-comparison correction across 11 KS tests: mass (p=0.005) and Mach number (p=0.006) are individually significant but would not survive Bonferroni (0.0045), and density (p=0.04) clearly wouldn't. Given the parameters are correlated, a permutation test would be more appropriate than a simple p-value list. Second, the magnetic-field subsample is 27 of 53 cores with uneven fractions and large DCF uncertainties; the abstract states the null more strongly than the data support. Third, there is a concrete internal inconsistency: Section 3.1 says cores smaller than the beam are excluded, but Table 1 lists radii of 3200–5000 au against a 14.6\\\" (5800 au) beam. Either the reported radii are not beam-deconvolved or the exclusion wasn't applied. That needs a clear statement. The stress-test worry that this creates the central result is probably wrong — beam-smearing of compact cores would dilute density differences, not produce them, and the radii distributions of the two groups are similar — but the N2H+ and Herschel beams are broad enough that environment contamination of Mach number and temperature for the smallest cores is a real possibility. The manuscript's assertion that the impact is 'minimal' is not demonstrated; a resolution-matched subsample or deconvolution test would settle it.\n\nOverall this is a solid incremental paper with one new result and one replicated result. It deserves a serious referee and is publishable after moderate revision. I wouldn't cite it in my own work unless I were working directly on Orion multiplicity, but I'd send it out.","headline":"Worth a serious referee: new magnetic-field null in Orion A core multiplicity, but the density/Mach result partly replicates Luo et al. and has unresolved beam issues.","tokens_in":23901,"tokens_out":4119,"would_cite":false,"duration_ms":46125,"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 paper claims that in Orion A, dense cores hosting binary or multiple protostars are significantly denser and more supersonically turbulent than cores hosting single protostars, while turbulence-to-gravity and magnetic-field-to-gravity…","keywords":["star formation","dense cores","fragmentation","multiplicity","turbulence","magnetic fields","Orion A","Mach number"],"falsifier":"Re-measure the density and Mach number of the same 53 cores after convolving all data to a common beam that fully resolves the smallest cores, or observe the sample with a matched high-resolution line survey; if the density and Mach number differences between single and multiple cores disappear, the central claim fails.","tokens_in":1830,"feed_emoji":"⭐","tokens_out":1883,"duration_ms":54979,"temperature":0.7,"pith_summary":"The paper compares 15 dense cores that host binary or multiple protostars with 38 cores hosting a single protostar in the Orion A cloud, using archival JCMT 850 micron continuum, Herschel dust temperature, and Nobeyama N2H+ line data plus JCMT polarization maps. Its central claim is that the fragmented cores have significantly higher density and Mach number, while their ratios of turbulent pressure to gravitational energy and their mass-to-flux ratios are indistinguishable from single-star cores. If correct, this pins core fragmentation on high density and supersonic turbulence and dismisses magnetic support as a decisive factor in these cores. The result would sharpen the search for what sets stellar multiplicity in low-mass star-forming regions.","feed_headline":"Cores that birth multiple stars are denser and more turbulent","feed_subtitle":"Survey of 53 Orion A cores finds no link between magnetic support and fragmentation.","key_machinery":"The argument rests on a catalogue of 53 dense cores extracted from the JCMT 850 micron map with astrodendro, classified as single or binary/multiple using the VANDAM ALMA survey's protostar positions at 0.1 arcsecond resolution. Core-averaged density comes from 850 micron flux with Herschel dust temperatures, Mach number and non-thermal velocity dispersion from hyperfine fits to Nobeyama N2H+ spectra, and magnetic field strength from POL-2 polarization using the Davis-Chandrasekhar-Fermi method with unsharp-masked angular dispersions. The statistical engine is a bootstrap-resampled cumulative distribution comparison, with Kolmogorov-Smirnov tests applied to each physical parameter to decide which properties separate the two groups.","core_discovery":"The paper's central discovery, stated in Sections 5.1 and 6, is that dense cores forming binary or multiple systems in Orion A are statistically denser and more supersonically turbulent than cores forming a single protostar. Kolmogorov-Smirnov tests give p-values of 0.005 for core mass, 0.04 for density, 0.006 for Mach number, and 0.002 and 0.009 for the ratios of core radius to Jeans radius and core mass to Jeans mass, all favoring the binary/multiple group. In contrast, the p-values for core radius (0.21), aspect ratio (0.50), turbulent-to-gravitational energy ratio (0.27), magnetic field strength (0.14), and mass-to-flux ratio (0.64) show no significant difference. The number of protostars per core also correlates with density, Mach number, and the Jeans-based ratios. The authors interpret these results as evidence that high density and supersonic turbulence promote local collapse and fragmentation, while the magnetic field has limited influence on whether a core fragments in Orion A.","pith_inferences":["A direct test this paper leaves implicit is a matched-resolution reanalysis: convolving all data to a common beam that covers the smallest cores and re-measuring density and Mach number would show whether the single-versus-multiple differences are robust or partly a resolution artifact.","If the trend holds generally, multiplicity fraction should increase toward the densest, most turbulent parts of any star-forming cloud, offering a cloud-scale diagnostic for predicting where wide binaries and higher-order multiples are born.","The authors argue that outflow or infall energy is unlikely to inflate the line widths in multiple cores; an independent check using outflow tracers such as CO or SiO would verify whether any of the Mach number excess is feedback-related rather than intrinsic turbulence.","Extending the same comparison to other clouds with measured mass-to-flux ratios would test whether Orion A's null magnetic-field result generalizes or is peculiar to its magnetic environment."],"forward_implications":["If high density and supersonic turbulence drive fragmentation, surveys of core multiplicity can use density and Mach number as predictors of which cores will split into binary or multiple systems.","The absence of a mass-to-flux ratio trend suggests that magnetic support is not the dominant regulator of multiplicity in Orion A, so models of low-mass star formation should treat turbulence and gravity, not magnetic fields, as the main fragmentation switches in such cores.","The correlation between protostar number and Mach number supports turbulent core fragmentation scenarios in which supersonic fluctuations create multiple density peaks that exceed the local Jeans mass.","The null result on the turbulent-to-gravitational energy ratio indicates that turbulence acts as a fragmentation promoter rather than as a stabilizing support, because gravity dominates in most of these cores.","Combined with the earlier sample of Luo et al. (2022), the two studies claim that 91 cores across the Orion complex point to the same drivers of fragmentation: gas density and Mach number."],"supporting_citations":[{"why":"Supplies the VANDAM survey catalogue of protostar positions and multiplicities in Orion A, used to classify each dense core as single or binary/multiple.","marker":"Tobin et al. 2020"},{"why":"Provides the Nobeyama 45-m N2H+ (1-0) data cube from which line widths and non-thermal velocity dispersions are measured.","marker":"Tatematsu et al. 2008"},{"why":"Presents part of the JCMT POL-2 Orion A polarization data used to infer magnetic field orientations and strengths.","marker":"Pattle et al. 2017"},{"why":"Is the previous study of core fragmentation in the Orion complex that found higher density and Mach number in multiple systems, which the present paper extends with a complementary sample.","marker":"Luo et al. 2022"},{"why":"Supplies the dust opacity at 850 micron used to convert continuum flux to core mass.","marker":"Johnstone et al. 2017"},{"why":"Is the original method, as adapted via the Davis-Chandrasekhar-Fermi relation, used to estimate magnetic field strengths from polarization angular dispersion.","marker":"Chandrasekhar & Fermi 1953"},{"why":"Provides the numerical correction factor Q=0.5 adopted in the Davis-Chandrasekhar-Fermi estimate of magnetic field strength.","marker":"Ostriker et al. 2001"}],"fun_headline_variants":["Orion A cores spawning multiple stars are denser and more turbulent","Dense, supersonic cores fragment into multiple stars; magnetic fields don't","No magnetic link to fragmentation in Orion A's dense cores","Binary-forming cores in Orion A are denser and more turbulent"],"cache_read_input_tokens":25856,"weakest_assumption_plain":"The inferred core-averaged density and Mach number faithfully represent the pre-fragmentation physical conditions, even though the JCMT, Nobeyama, and Herschel beams are comparable to or larger than the smallest cores in the sample.","fun_headline_variants_meta":{"raw":{"variants":["Orion A cores spawning multiple stars are denser and more turbulent","Dense, supersonic cores fragment into multiple stars; magnetic fields don't","No magnetic link to fragmentation in Orion A's dense cores","Binary-forming cores in Orion A are denser and more turbulent"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000313,"raw_usage":{"total_tokens":1813,"prompt_tokens":1014,"completion_tokens":799,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":724}},"tokens_in":630,"tokens_out":799,"duration_ms":8896,"temperature":1.0,"reasoning_tokens":724,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:38:02.238482+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-measure the density and Mach number of the same 53 cores after convolving all data to a common beam that fully resolves the smallest cores, or observe the sample with a matched high-resolution line survey; if the density and Mach number differences between single and multiple cores disappear, the central claim fails.","supporting_citations":[{"cited_title":"2017, , 836, 132, 10.3847/1538-4357/aa5b95","cited_arxiv_id":null,"evidence_quote":"Supplies the dust opacity at 850 micron used to convert continuum flux to core mass."}],"review_version":1}