Ambipolar diffusion and the molecular abundances in prestellar cores
Pith reviewed 2026-05-25 02:30 UTC · model grok-4.3
The pith
Line intensity profiles of HCN and CH3OH differ enough between models to trace whether prestellar cores are magnetically super- or sub-critical.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Sub-critical models produce greater central molecular depletion than super-critical ones because the collapse lasts longer. Observable molecule-to-hydrogen column density ratios remain broadly similar. N2H+ and HCO+ show qualitative profile differences on 0.01 pc scales but require independent knowledge of the hydrogen column and yield similar line intensities after radiative transfer. HCN and CH3OH produce line intensity profiles that differ more strongly between the two regimes.
What carries the argument
Post-processing of non-ideal MHD simulations with a time-dependent gas-grain chemical code and subsequent radiative transfer to generate observable line intensity profiles.
If this is right
- Sub-critical cores exhibit longer collapse times and therefore stronger central depletion for most species.
- Molecule-to-hydrogen column density ratios are generally too similar to discriminate super- from sub-critical models.
- N2H+ and HCO+ display qualitative profile differences on 0.01 pc scales that would require separate hydrogen column measurements to interpret.
- Predicted line intensities for N2H+ and HCO+ remain comparable between models, while those for HCN and CH3OH differ more strongly.
Where Pith is reading between the lines
- Targeted mapping of HCN and CH3OH in nearby cores with current interferometers could test whether magnetic support sets the collapse timescale.
- If the differences persist under varied initial conditions, these molecules could become standard diagnostics for the role of ambipolar diffusion in core evolution.
Load-bearing premise
The chemical reaction rates, initial elemental abundances, and simplified treatment of ambipolar diffusion in the MHD runs produce molecular distributions and line intensities that match reality without dominant systematic errors.
What would settle it
High-resolution observations of HCN or CH3OH line intensity profiles toward a prestellar core that fail to show the predicted model-dependent differences on 0.01 pc scales.
read the original abstract
We investigate differences in the molecular abundances between magnetically super- and sub-critical prestellar cores, performing three-dimensional non-ideal magnetohydrodynamical (MHD) simulations with varying densities and magnetic field strengths, and post-processing the results with a time-dependent gas-grain chemical code. Most molecular species show significantly more central depletion in subcritical models, due to the longer duration of collapse. However, the directly observable quantities - the molecule to hydrogen column density ratios - are generally too similar for observational data to discriminate between models. The profiles of N$_2$H$^+$ and HCO$^+$ show qualitative differences between supercritical and subcritical models on scales of $0.01 \, {\rm pc}$, which may allow the two cases to be distinguished. However, this requires knowledge of the hydrogen column density, which is not directly measureable, and predicted line intensity profiles from radiative transfer modelling are similar for these molecules. Other commonly observed species, such as HCN and CH$_3$OH, have line intensity profiles which differ more strongly between models, and so are more promising as tracers of the mechanism of cloud collapse.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper performs 3D non-ideal MHD simulations of prestellar cores with varying density and magnetic field strength (super- vs. sub-critical), followed by post-processing with a time-dependent gas-grain chemical network and radiative transfer. It reports greater central depletion in subcritical models for most species due to longer collapse timescales, but finds that molecule-to-hydrogen column-density ratios are generally too similar to discriminate observationally. N2H+ and HCO+ show qualitative profile differences on 0.01 pc scales, while HCN and CH3OH exhibit stronger differences in predicted line intensity profiles, positioning the latter as promising tracers of collapse mechanism.
Significance. If the reported line-intensity differences prove robust, the work would supply concrete, observationally testable predictions for distinguishing magnetically regulated vs. non-regulated collapse, a longstanding question in star formation. The forward-modeling approach (MHD + chemistry + RT) is a methodological strength that generates falsifiable outputs rather than post-hoc fits.
major comments (2)
- [Abstract and §4] Abstract and §4: the central claim that HCN and CH3OH line intensity profiles 'differ more strongly' and are 'more promising as tracers' is based on a single chemical-network realization and single RT setup; no sensitivity runs are reported that vary key rate coefficients (depletion/desorption), initial C/O ratio, or cosmic-ray ionization rate, all of which are known to shift abundances by factors of several on the 0.01 pc scales under discussion.
- [§2.3] §2.3: the gas-grain network is fixed without any ensemble or uncertainty quantification; therefore it is impossible to determine whether the qualitative differences in HCN/CH3OH intensities survive plausible variations in the chemical parameters that dominate the post-processing pipeline.
minor comments (1)
- Notation for column-density ratios vs. line intensities could be clarified in the figure captions to avoid reader confusion between the two quantities.
Simulated Author's Rebuttal
We thank the referee for their positive evaluation of the work's significance and for the constructive comments on the chemical modeling. We address the two major comments below, which both concern the lack of sensitivity tests. While we agree that such tests would be desirable, the reported differences arise principally from the contrasting dynamical timescales between super- and sub-critical models; we therefore maintain that the qualitative conclusions remain useful even without an ensemble of networks.
read point-by-point responses
-
Referee: [Abstract and §4] Abstract and §4: the central claim that HCN and CH3OH line intensity profiles 'differ more strongly' and are 'more promising as tracers' is based on a single chemical-network realization and single RT setup; no sensitivity runs are reported that vary key rate coefficients (depletion/desorption), initial C/O ratio, or cosmic-ray ionization rate, all of which are known to shift abundances by factors of several on the 0.01 pc scales under discussion.
Authors: We acknowledge that the results rest on a single chemical network and radiative-transfer setup. The dominant physical distinction between the models is the longer free-fall time available for depletion and freeze-out in the sub-critical case. Because this is a first-order effect of integration time rather than a fine-tuned rate coefficient, we expect the relative contrast in central depletion (and therefore in the emergent line profiles of HCN and CH3OH) to persist under moderate variations of the parameters listed by the referee. Nevertheless, we agree that an explicit discussion of this limitation is warranted and will add a short paragraph in §4 (and a corresponding sentence in the abstract) noting that absolute abundances are uncertain but that the timescale-driven contrast is expected to be robust. No new simulations are performed, so the revision is limited to this clarification. revision: partial
-
Referee: [§2.3] §2.3: the gas-grain network is fixed without any ensemble or uncertainty quantification; therefore it is impossible to determine whether the qualitative differences in HCN/CH3OH intensities survive plausible variations in the chemical parameters that dominate the post-processing pipeline.
Authors: We agree that the network is fixed and that no ensemble or uncertainty quantification is provided. As explained in the response to the first comment, the qualitative difference in line-intensity profiles is driven by the order-of-magnitude difference in collapse duration rather than by the precise numerical values of individual rate coefficients. We will therefore expand the methods section (§2.3) to state explicitly that the adopted network is the standard gas-grain network of the cited reference and that future work should explore parameter variations. This addition constitutes a partial revision; we do not claim to have performed the sensitivity study requested. revision: partial
Circularity Check
No circularity; forward simulation outputs independent of inputs
full rationale
The paper executes non-ideal MHD simulations with ambipolar diffusion from stated initial densities and magnetic field strengths, followed by time-dependent gas-grain chemistry and radiative transfer post-processing. Reported line intensity profile differences (e.g., for HCN and CH3OH) are direct numerical outputs of this pipeline rather than quantities fitted to or defined in terms of the target observables. No self-definitional equations, fitted-input predictions, or load-bearing self-citations appear in the abstract or described chain; the derivation remains self-contained against external benchmarks.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.