{"id":"f1810eeb-145d-4bd3-9c17-647eec6df793","arxiv_id":"2501.07626","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A high-resolution MHD re-simulation shows a young protostellar disc formed from cloud-scale initial conditions is highly turbulent (sonic Mach ~2), near magnetically equipartitioned (beta ~ 1), and accretes episodically.","lead":"This paper simulates the birth of a protostellar disc from a turbulent, magnetized molecular cloud core, resolving the disc down to sub-AU scales over the first 10,000 years after star formation. It finds the young disc is highly turbulent, accretes in bursts, and develops a magnetic field strong enough to rival the thermal pressure.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed equipartition (β~1, M_A~2) is not supported by the paper's own Eq. 15: with M_A~2 and M_s~2, β~2; the abstract overstates the result, and the ideal-MHD limitation in Sec. 4.4 reinforces the need for a corrected claim.","rationale":"The paper is a well-executed zoom-in MHD simulation with appropriate refinement and resolution studies; I do not question the qualitative picture of a turbulent, thick, episodically accreting disc. The load-bearing problem is that one of the headline quantitative claims—equipartition of thermal and magnetic pressure with β~1 and M_A~2—is contradicted by the paper's own defining relation. Eq. (15) is β = 2 M_A^2 / M_s^2; substituting the quoted M_A~2 and M_s~2 gives β~2, and the measured M_A range (1.5–3.5) makes β vary substantially across the disc. The abstract and conclusion present β~1 as a robust result, which is not supported by the data shown. The ideal-MHD caveat is not a separate fatal flaw—it is an additional reason to distrust the equipartition number, and the authors themselves acknowledge in Sec. 4.4 that β could be higher by factors of a few. Because the central claim about magnetic field influence and its dynamical consequences rests on this number, the paper needs a corrected, qualified statement before acceptance. The reader's CONDITIONAL verdict is appropriate. I differ slightly from the reader's framing of the weakest assumption: the internal inconsistency of Eq. (15) is more immediately falsifiable, although non-ideal MHD is the larger physical uncertainty.","tokens_in":28183,"tokens_out":7792,"duration_ms":76738,"concrete_test":"Recompute the radial β profile using Eq. (15) from the published M_s and M_A data in Figs. 11 and 14. If the representative or disc-averaged β is ~2 rather than ~1, or varies by more than a factor of ~3 across the disc, the abstract and conclusions should be revised to say 'β of order a few' rather than 'equipartition, β~1'. A single quantitative statement of the radial range of β would settle whether the overclaim is real.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing issue is the headline magnetic-state claim: 'equipartition of the thermal and magnetic pressure' with Alfvén Mach number ~2 and β~1. This is not what the paper's own equations and measurements give. Using Eq. (15), β = 2 M_A^2 / M_s^2. With the reported sonic and Alfvén Mach numbers both ~2, β ~ 2, not 1. The β~1 value appears only near the radius where M_A ~ 1.4 while M_s ~ 2; the radial profile in Fig. 14 shows M_A ranging from ~1.5 to ~3.5, which would give β from ~1 to ~6 across the disc if M_s ~ 2. Thus the abstract and conclusion point (vii) overstate an equipartition. This is not a semantic quibble: the 'substantial magnetic influence' and fragmentation-inhibition statements are built on this number. The ideal-MHD assumption (Eq. 1, no non-ideal terms) makes it worse: in the disc's high-density, low-ionization regime, ambipolar diffusion and Ohmic dissipation decouple the field from the gas, and the authors themselves note in Sec. 4.4 that the field may be over-amplified and the scale height overestimated; Sec. 3.4.3 concedes β could be higher by factors of a few in non-ideal MHD. So the central magnetic-equipartition claim is both internally inconsistent with the paper's own beta formula and not robust to the acknowledged physics omission.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a high-resolution adaptive mesh refinement (AMR) ideal MHD simulation of the collapse of a dense core extracted from a (2 pc)^3 molecular cloud simulation, re-simulated in a (0.1 pc)^3 domain at up to 0.63 AU resolution for 10 kyr after sink (protostar) formation. The authors characterise the resulting disc: it grows to a radius of roughly 50 AU (diameter ~100 AU), has a mass of 0.12 M_sun around a 0.15 M_sun protostar, accretes episodically on ~100 yr timescales, is geometrically thick and strongly turbulent with sonic Mach number ~2, develops a density profile that steepens toward a power-law index ~1, and shows magnetic field amplification with turbulent Alfvén Mach number ~2 and plasma beta claimed to be ~1. The paper also reports intermittent sub-Alfvénic regions above and below the mid-plane that produce magnetic 'bubbles' rather than a coherent jet.","tokens_in":28366,"tokens_out":3620,"duration_ms":38606,"significance":"If the findings hold, the paper provides a useful case study of protostellar disc formation from realistic, cloud-scale initial conditions, capturing sub-AU structure that is difficult to achieve in cloud-scale simulations. The strengths include the careful use of inherited initial conditions, the inclusion of box-size and resolution convergence studies in Appendices A and C, quantitative power-law fits to density and surface density profiles, and a generally detailed description of the numerical setup. The central magnetic-state claim, however, is not internally consistent with the paper's own Eq. (15), and the acknowledged ideal-MHD assumption further weakens it. The paper is therefore a valuable numerical experiment whose headline conclusion about equipartition needs substantial correction and qualification before publication.","major_comments":[{"comment":"The claim that plasma beta is ~1, i.e., thermal and magnetic pressure are in equipartition, is not supported by the paper's own numbers. Equation (15) gives beta = 2 M_A^2 / M_s^2; with the quoted M_s ~ 2 and M_A ~ 2, this gives beta ~ 2, not 1. The radial profile in Fig. 14 shows M_A varying from about 1.5 to 3.5 across the disc, which for M_s ~ 2 gives beta ranging from roughly 1 to 6, with beta ~ 1 only near r ~ 40 AU. The abstract's wording 'equipartition ... i.e., leading to an Alfvén Mach number of ~2' is also internally inconsistent, since beta = 1 and M_s = 2 would imply M_A ~ 1.4. This is not a semantic issue because conclusion (vii) and the abstract use the equipartition statement to argue for significant magnetic influence on disc dynamics and fragmentation. The claim should be corrected to beta ~ 2 (or the appropriate radial range) and the presentation adjusted accordingly.","section":"Section 3.4.3, Eq. (15); Abstract; Section 5 (vii)"},{"comment":"The magnetic field amplification, Alfvén Mach number, and beta are computed under the ideal MHD approximation. The authors themselves state in Sec. 4.4 that non-ideal MHD effects may be highly relevant in the dense disc and that the disc scale height may be overestimated, and in Sec. 3.4.3 they note that beta could be higher by factors of a few with non-ideal MHD. Because the paper's headline claim about equipartition and magnetic influence depends directly on the amplified field under ideal MHD, the conclusion should be presented as an ideal-MHD result that is an upper limit on magnetic field strength, not as a robust physical finding. A brief quantitative test or at least a strongly worded uncertainty statement is needed.","section":"Section 4.4; Section 3.4.3"}],"minor_comments":[{"comment":"The sentence 'oursimulationsisconsistent with the expected MMSN' contains a typo; it should read 'our simulation is consistent'.","section":"Section 5 (iii)"},{"comment":"The caption appears to have a copy-paste error in the colorbar label: 'Alfvén Mach number over 5 AU along LoS □ g cm□3' should be a dimensionless quantity, not a density unit.","section":"Figure 15 caption"},{"comment":"The reference 'Federrathet al. 2010b' is misspelled; it should be 'Federrath et al. 2010b'.","section":"Section 2.2.2"},{"comment":"The abstract states the disc grows to a diameter of approximately 100 AU, while Fig. 3 shows r_disc ~ 50 AU at 10 kyr; this is consistent, but the abstract could clarify that this is the 80%-mass radius definition, to avoid confusion with observational radius definitions that may differ.","section":"Abstract and Section 3.2.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent single-case numerical study with clear methods and useful convergence tests. The main concern is that the headline magnetic-equipartition claim is internally inconsistent with the paper's own Eq. (15) and is further weakened by the acknowledged ideal-MHD limitation. These issues are correctable within the scope of the manuscript, so major revision rather than rejection seems appropriate. The editor may also wish to consider whether the single-realisation nature of the study is sufficiently emphasised for the journal's readership."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this is a solid re-simulation study with real value, but the headline magnetic claim does not survive contact with the paper's own equation. The abstract says the disc reaches equipartition (beta ~ 1) with an Alfven Mach number of ~2, yet Eq. 15 gives beta = 2 M_A^2 / M^2, and with M_A ~ 2 and M ~ 2 that is beta ~ 2, not 1. The radial profile in Fig. 14 shows M_A ranging from ~1.5 to ~3.5, which would put beta somewhere between roughly 1 and 6 across the disc. The authors need to correct this before publication.\n\nWhat is actually new: a single protostellar disc at 0.63 AU effective resolution, re-simulated from the 2 pc cloud-scale GTMJR simulation (Appel et al. 2023) and followed for 10 kyr after sink formation. The density-profile steepening, episodic accretion on ~100 yr timescales, sonic Mach ~2, and the wound-up azimuthal field are all measured in a self-consistent environment rather than hand-crafted initial conditions. The box-size and resolution convergence studies are real, not decorative. The paper is careful about numerical definitions and the analysis choices are post-hoc but not circular. The qualitative picture is plausible and consistent with earlier re-simulation work (Kuffmeier et al. 2017; Bate 2018; Lebreuilly et al.).\n\nSoft spots, in proportion: (1) The equipartition overclaim is the main issue. It is not semantic, because the conclusions lean on \"substantial magnetic influence\" and fragmentation inhibition. (2) Ideal MHD only. The authors flag in Sec. 4.4 that non-ideal effects may matter in the dense disc and that beta might shift by factors of a few, but they do not quantify it for this run. That is an honest limitation, not a fatal flaw, but it means the magnetic claims should be lower-key. (3) Reproducibility: 10 TB of data on request, no public analysis scripts. Understandable, but weak.\n\nMinor: \"equipartition of thermal and magnetic pressure\" is ambiguous in the abstract. What they compute is the turbulent beta; they should say beta is of order a few, with a range.\n\nWho this is for: people working on protostellar disc formation, episodic accretion, and disc magnetism. It deserves a serious referee. I would send it back with a request to fix the beta claim, add error bars on the Mach numbers and beta, and qualify the ideal-MHD dependence in abstract and conclusions. Conditional accept.\n\nBest,\n[you]","headline":"A capable re-simulation study whose magnetic equipartition headline is undercut by the paper's own Eq. 15 — worth revising, not rejecting.","tokens_in":29099,"tokens_out":2494,"would_cite":true,"duration_ms":24139,"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":"This paper claims that a protostellar disc formed from cloud-scale initial conditions is born turbulent (Mach ~ 2) and bursty on ~100 yr timescales, with magnetic pressure reaching near equipartition with thermal pressure (beta ~ 1).","keywords":["protostellar discs","magnetohydrodynamics","star formation","turbulence","accretion variability","adaptive mesh refinement","cloud-scale initial conditions","magnetic field amplification"],"falsifier":"Re-run the same 0.1 pc re-simulation with non-ideal MHD (ambipolar diffusion plus Ohmic resistivity, with or without the Hall term) at the same maximum resolution and compare the turbulent Alfvén Mach number, plasma $\\beta$, and disc scale height at 10 kyr. If $\\beta$ rises well above 1, $\\mathcal{M}_A$ falls well below 2, or $H/r$ drops significantly, the claimed equipartition and thick, turbulent disc structure are artifacts of ideal MHD rather than properties of real protostellar discs. Resolved magnetic-field or turbulence measurements of a Class 0 disc showing magnetic energy far below thermal or turbulent energy would also contradict the claim.","tokens_in":27801,"feed_emoji":"🌌","tokens_out":13977,"duration_ms":117512,"temperature":0.7,"pith_summary":"The paper tries to show what a protostellar disc looks like when it forms from the turbulent, magnetised gas left over by a molecular-cloud simulation, rather than from an idealised spherical cloud. Following the collapse of one dense core for 10,000 years after the protostar appears, it claims the disc grows to roughly 100 AU across, accretes mass in episodic bursts on $\\sim100$-year timescales, and is strongly turbulent with a sonic Mach number near 2. It further claims that accretion events amplify the magnetic field until thermal and magnetic pressure reach rough equipartition (plasma $\\beta \\sim 1$, Alfvén Mach number $\\sim 2$), with intermittent sub-Alfvénic outflow regions. A sympathetic reader would care because these numbers set the physical conditions in which planets begin to form and tie observed accretion variability to the turbulence of the parent cloud rather than to internal disc physics alone.","feed_headline":"Infant discs accrete in bursts with near-equipartition magnetic fields","feed_subtitle":"Parent-cloud turbulence drives 100-year accretion bursts and magnetic field growth to equipartition.","key_machinery":"The load-bearing machinery is the re-simulation technique: extract a pre-collapse dense core from an existing (2 pc)^3 molecular-cloud simulation, inheriting its turbulence, magnetic field, and density structure, and re-run it in a 0.1 pc box with adaptive mesh refinement that keeps the Jeans length resolved by 30 to 60 cells, modelling the protostar as a sink particle of radius 1.6 AU. The simulations solve the ideal MHD equations closed by a piecewise polytropic equation of state. The physical argument is carried by two dimensionless ratios computed from the simulated fields, the sonic Mach number $\\mathcal{M} = \\sigma_v/c_s$ and the turbulent Alfvén Mach number $\\mathcal{M}_A = \\sigma_v/v_A$, combined through $\\beta = 2 c_s^2/v_A^2 = 2 \\mathcal{M}_A^2/\\mathcal{M}^2$ to conclude that turbulent magnetic pressure is comparable to thermal pressure in the disc.","core_discovery":"Starting from a 0.1 pc box carved out of a (2 pc)^3 turbulent, magnetised molecular-cloud simulation, the authors follow the collapse of one 1.34 solar-mass core and the first 10 kyr of the star+disc system at 0.63 AU maximum resolution. They find that the disc grows to a radius of about 50 AU (diameter ~100 AU) and a mass of 0.12 solar masses around a 0.15 solar-mass protostar, with the disc-to-star mass ratio near unity. Accretion onto both the disc and the star is episodic: brief bursts at roughly $10^{-5}$ solar masses per year on $\\sim100$ yr timescales alternate with lulls, and the disc radius fluctuates on the same timescale. The disc is geometrically thick and turbulent, with velocity dispersions comparable to the local sound speed and a sonic Mach number near 2; its density profile steepens from nearly flat at formation to a power law with exponent $\\sim 1$ at 10 kyr, while the surface density remains shallower than the minimum-mass solar nebula. The magnetic field, wound up by rotation, is amplified so that the turbulent Alfvén Mach number is $\\sim 2$ and the plasma $\\beta$ is $\\sim 1$, implying thermal and magnetic pressures are in rough equipartition; regions above and below the midplane intermittently become sub-Alfvénic and launch magnetic bubbles rather than a coherent jet.","pith_inferences":["Beyond the paper's claims, if a Mach $\\sim 2$, $\\beta \\sim 1$ state is generic for young discs, dust settling and planetesimal-formation models that assume laminar, weakly magnetised discs start from the wrong initial conditions; dust-evolution simulations run on top of this turbulent state would be a direct test.","The $\\sim100$ yr episodic burst timescale implies short monitoring campaigns could easily misread stochastic turbulent accretion as periodic variability; comparing burst waiting-time statistics from the simulation with long-baseline protostellar light curves would test the mechanism.","The absence of a coherent jet in a disc that still shows intermittent sub-Alfvénic outflows suggests jets may turn on only after the disc settles, so a search for jets preferentially in older, less turbulent Class I sources would be a consistent observational test."],"forward_implications":["At 10 kyr the disc holds 0.12 $M_\\odot$ around a 0.15 $M_\\odot$ protostar, so young discs can carry as much mass as their stars, setting a massive early reservoir for planet formation.","Accretion onto both star and disc is episodic, with bursts reaching roughly $10^{-5}\\,M_\\odot\\,\\mathrm{yr}^{-1}$ separated by lulls and radius fluctuations on ~100 yr timescales, tying young-protostar variability to the turbulence of the parent core.","The disc is geometrically thick ($H/r \\sim 0.1$–$0.5$) and flared, with scale heights of a few to ~10 AU at tens of AU, so early discs are not thin Keplerian structures.","Magnetic field amplification by collapse and rotation brings the turbulent Alfvén Mach number to ~2 and plasma beta to ~1, making magnetic pressure dynamically comparable to thermal pressure and capable of suppressing fragmentation despite strong spiral features.","The surface-density profile steepens from a nearly flat distribution at formation to a power law with exponent $\\sim 1$ at 10 kyr, still shallower than the minimum-mass solar nebula but moving toward it as mass accumulates inward."],"supporting_citations":[{"why":"Supplies the (2 pc)^3 molecular-cloud simulation from which the progenitor core and its turbulent, magnetised initial conditions are extracted.","marker":"Appel et al. (2023)"},{"why":"Defines the sink-particle creation and accretion technique used to model the protostar.","marker":"Federrath et al. (2010b)"},{"why":"Adds the sink-particle gravity implementation and jet/outflow feedback framework used in the base simulation.","marker":"Federrath et al. (2014)"},{"why":"Earlier similar re-simulation study used for comparison of the accretion-rate evolution and disc densities.","marker":"Kuffmeier et al. (2017)"},{"why":"Provides the disc-radius definition and statistical disc properties that the disc mass and size results are compared with.","marker":"Bate (2018)"},{"why":"Establishes the Jeans-length resolution criterion underlying the refinement strategy.","marker":"Truelove et al. (1997)"},{"why":"Justifies resolving the Jeans length by 30-60 cells to capture turbulent dynamics and magnetic amplification.","marker":"Federrath et al. (2011)"},{"why":"Quantifies the expected upward shift in plasma beta when non-ideal MHD is included, marking the main caveat on the equipartition claim.","marker":"Masson et al. (2016)"}],"fun_headline_variants":["Turbulent young discs accrete in bursts with magnetic fields near equipartition","Episodic accretion and magnetic equipartition in turbulent protostellar discs","Cloud turbulence drives 100-year accretion bursts and near-equipartition fields","Simulated protostellar discs show bursty accretion and magnetic equipartition"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything about the magnetic state of the disc — the amplification to $\\beta \\sim 1$, the Alfvén Mach number $\\mathcal{M}_A \\sim 2$, and the sub-Alfvénic bubbles — depends on the ideal MHD approximation, which omits ambipolar diffusion, Ohmic dissipation, and the Hall effect; in the dense, poorly ionised gas of a real disc those effects can alter field strength and geometry substantially, and the paper itself notes the disc scale height may be overestimated as a result.","fun_headline_variants_meta":{"raw":{"variants":["Turbulent young discs accrete in bursts with magnetic fields near equipartition","Episodic accretion and magnetic equipartition in turbulent protostellar discs","Cloud turbulence drives 100-year accretion bursts and near-equipartition fields","Simulated protostellar discs show bursty accretion and magnetic equipartition"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000764,"raw_usage":{"total_tokens":3485,"prompt_tokens":1133,"completion_tokens":2352,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":749,"completion_tokens_details":{"reasoning_tokens":2268}},"tokens_in":749,"tokens_out":2352,"duration_ms":16848,"temperature":1.0,"reasoning_tokens":2268,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:38:47.871873+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the same 0.1 pc re-simulation with non-ideal MHD (ambipolar diffusion plus Ohmic resistivity, with or without the Hall term) at the same maximum resolution and compare the turbulent Alfvén Mach number, plasma $\\beta$, and disc scale height at 10 kyr. If $\\beta$ rises well above 1, $\\mathcal{M}_A$ falls well below 2, or $H/r$ drops significantly, the claimed equipartition and thick, turbulent disc structure are artifacts of ideal MHD rather than properties of real protostellar discs. Resolved magnetic-field or turbulence measurements of a Class 0 disc showing magnetic energy far below thermal or turbulent energy would also contradict the claim.","supporting_citations":[],"review_version":1}