REVIEW 2 major objections 4 minor 2 references
Early Planet Formation in Embedded Disks (eDisk) XXII: Keplerian disk, disk structures and jets/outflows in the Class 0 protostar IRAS 04166+2706
T0 review · 2 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read This paper claims that the Class 0 protostar IRAS 04166+2706 hosts a compact, rotationally supported disk orbiting a low-mass star, and that its jets fired their latest knots within the last 20–25 years.
desk verdict A careful eDisk paper whose new 22 au disk image and innermost jet knots are solid, but the central stellar mass rests on a one-sided PV analysis that needs robustness testing before the 'Keplerian disk' title earns its keep. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing measurement is the edge/ridge fit of position–velocity (PV) diagrams for $^{13}$CO(2–1) and C$^{18}$O(2–1) along the dust disk's major axis, using the form $V_{\rm rot} = v_b (r/r_b)^{-p}$. A fitted exponent near $p = 0.5$ identifies Keplerian rotation, and the stellar mass follows from $M_\star = v_b^2 r_b / (G \sin^2 i)$ with inclination $i = 47°$ fixed by the deconvolved dust shape under a vertically thin circular disk assumption. The ridge and edge fits bracket the mass between 0.15 and 0.39 $M_\odot$ because ridge fitting is known to underestimate and edge fitting to overestimate the true stellar mass; the authors adopt these as lower and upper limits.
What would settle it
Take the same PV cuts at higher velocity resolution with a midplane tracer (e.g., C$^{18}$O 3–2 or CS) and test whether the power-law index $p$ remains consistent with 0.5 over a wider radius range, and whether the inferred $M_\star$ matches a radiative-transfer model of the dust asymmetry; in parallel, re-image the continuum at higher signal-to-noise to check whether the ~14 au gap is a true axisymmetric ring or an azimuthal-average artifact.
Extended reading notes
Core claim
The paper reports ALMA 1.3 mm observations resolving a dust disk of radius ~22 au around IRAS 04166+2706. Position–velocity diagrams of $^{13}$CO(2–1) and C$^{18}$O(2–1) along the disk major axis, fitted with a power-law rotation curve $V = v_b(r/r_b)^{-p}$, give exponents $p = 0.58$–$0.67$ consistent within uncertainties with Keplerian rotation ($p = 0.5$). The same fits, combined with the inclination $i = 47°$ from the deprojected dust shape, place the central star at $0.15\,M_\odot < M_\star < 0.39\,M_\odot$. H$_{2}$CO and CH$_3$OH emission near the star shows velocity gradients aligned with the disk major axis and is consistent with the same Keplerian pattern. The dust continuum shows a
Load-bearing premise
The 0.15–0.39 $M_\odot$ range presumes that the $^{13}$CO and C$^{18}$O emission in the fitted PV windows rotates purely about the disk axis in a vertically thin, circular, uncontaminated Keplerian disk tilted at the 47° inclination derived from the dust; a flared disk, a warped emitting layer, or leaked infall/outflow emission would shift the mass range.
Editorial extensions
If this is right
- If true, rotationally supported, disk-mediated accretion begins already at the Class 0 stage, so the clock for planet formation can start when the central star is only a few tenths of a solar mass.
- The disk's mass budget (0.008–0.015 $M_\odot$; roughly 25–50 $M_\oplus$ of dust) is enough to grow a 4–7 $M_\oplus$ core at canonical efficiency, with envelope accretion still resupplying the disk.
- The ~73-year episodic ejection cycle and the ~20–25-year-old innermost knots imply pulsed accretion and angular-momentum removal during the deeply embedded phase.
- Toomre $Q$ values of ~1.4–5.6 put the disk near gravitational marginal stability, so modest future accretion could trigger spiral structure or fragmentation.
- The candidate ~14 au ring/gap, if confirmed, would be one of the smallest-radius substructures claimed around a Class 0 protostar, but the paper itself flags it as unconfirmed.
Reading between the lines
- Editorial: because the mass scaling goes as $\sin^{-2} i$, the flared-disk geometry invoked to explain the minor-axis brightness asymmetry would steepen the true inclination and push $M_\star$ above the quoted upper limit; a self-consistent radiative-transfer fit of the asymmetry could test this.
- Editorial: the noted misalignment between the $^{13}$CO velocity gradient and the dust major axis suggests the CO emission may trace a warped or elevated layer rather than the disk midplane, so the fitted rotation curve could blend velocities from different radii even if the midplane is perfectly Keplerian.
- Editorial: the ~73-year jet period and the candidate gap at ~14 au are near commensurate with the orbital period expected for a ~0.2 $M_\odot$ star at that radius, suggesting a testable link between episodic accretion and growing substructure, though the paper does not claim it.
- Editorial: the innermost knots B0/R0 at ~100–300 au can serve as a clock: multi-epoch proper-motion observations would measure the jet speed and launch time independently of the assumed inclination, checking the 20–25-year age.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper reports ALMA eDisk observations of Class 0 protostar IRAS 04166+2706. From 1.3 mm continuum at 0.05″ resolution the authors find a compact (~22 au) disk-like structure with PA=122° and inclination i=47°. Using SLAM edge/ridge fits to 13CO(2-1) and C18O(2-1) PV diagrams cut along the dust major axis, they fit power-law rotation curves with indices p=0.49–0.67 and derive a stellar mass range 0.15–0.39 M⊙, interpreting the rotation as possibly Keplerian. H2CO and CH3OH PV diagrams are overlaid with Keplerian curves as a consistency check. The paper also reports a minor-axis continuum brightness asymmetry (interpreted as a flared, optically thick disk), possible unconfirmed ring/gap features, and knotty, wiggling 12CO/SiO jets with inner knots R0/B0 implying mass ejection within the last ~20–25 years.
Significance. If confirmed, the paper provides an important datum for the eDisk program: a low-mass Class 0 protostar with a compact rotationally supported disk, a disk-scale estimate of M⋆, and a well-resolved view of the youngest jet knots. The manuscript is carefully hedged in several places: ring/gap features are explicitly noted as unconfirmed (Section 4.2), the Keplerian identification is qualified as 'possibly', and the mass range is presented as an edge/ridge systematic bracket. The analysis uses a public tool (SLAM) and presents all moment maps and PV diagrams, which is a strength. However, the quantitative mass range in the abstract rests on kinematic assumptions that are not yet demonstrated robust against infall/outflow contamination and cut geometry; this is the main obstacle to publication.
major comments (2)
- [§3.2, §4.3, Table 3] The PV fits in §4.3 are made along the dust major axis (PA=122°), yet §3.2 reports that the 13CO and C18O velocity gradient 'appears to be misaligned' from this axis and interprets the kinematics as a mixture of rotation and infall. If the actual kinematic major axis differs from the chosen cut by θ, the observed projected velocity is reduced by a factor cosθ; since M⋆ ∝ v²r/sin²i (Eq. 2 and the Keplerian formula in §4.3), the mass estimate is biased by cos²θ. The ~30° misalignment visible in Figure 5 corresponds to a ~25% bias, comparable to the spread between edge and ridge values in Table 3. In addition, the adopted i=47° assumes a vertically thin circular disk (§3.1) while §5.1 argues for a flared disk on the basis of the brightness asymmetry; M⋆ scales as sin⁻²i. Please quantify the sensitivity of p and M⋆ to cut position angle and to inclination, or otherwise show that the quoted r
- [§4.3, Figure 9] The fit windows are selected so that only one redshifted quadrant (v≥8.5 km/s) and one blueshifted quadrant (v≤4.5 km/s) enter the fit; however, the full PV diagrams show emission in all four quadrants, which the authors attribute to infall/outflow. The paper does not test whether an infalling-rotating envelope model, without a Keplerian disk, can reproduce these one-sided features at the selected velocities. Because the derived power-law indices are not tightly constrained (e.g., C18O ridge p=0.49±0.22), the claim that the rotation is 'closely consistent with Keplerian' and the subsequent stellar mass estimate are not uniquely supported. I ask for an explicit infalling-rotating-envelope fit to the full PV data (or an equivalent test) to show that the one-sided features require a Keplerian component. Without this, the mass range in the abstract should be presented as conditional.
minor comments (4)
- [Table 1 footnote] 'robuts' should be 'robust'.
- [Eq. (3)] The projection formula is unclear. Directly, 11″ at 156 pc corresponds to 1716 AU; dividing by 3×40 km/s would give ~68 yr. The printed expression with cos(i) and /sin(i) is dimensionally ambiguous and appears to give a different number. Please state the deprojection explicitly.
- [Section 6, bullet 3] The summary says C18O 'likely traces the Keplerian motion', whereas the abstract says 'possibly a Keplerian disk.' Please use consistent hedging.
- [Figure 8] The outer ring is only ~1.6 au wide, i.e., a fraction of the beam. The text appropriately cautions that it is unconfirmed; consider moving it to the Discussion to avoid it being read as a detection.
Circularity Check
No significant circularity: the stellar mass is obtained from a standard PV-diagram power-law fit, and the H2CO/CH3OH Keplerian overlay is an independent consistency check rather than a fitted prediction.
full rationale
The paper's central quantitative result is the protostellar mass range 0.15–0.39 M_sun, derived in Section 4.3 by fitting the 13CO and C18O PV diagrams with a power-law Vrot = vb (r/rb)^-p and then interpreting the resulting p ~ 0.5 as consistent with Keplerian rotation. This is a standard fitting-and-interpretation procedure: the data are the observed PV positions, the free parameters are vb, rb, p, and the mass is computed afterward from v = sqrt(GM/r) under the stated inclination assumption. The power-law index is not fixed to 0.5 in advance, so the Keplerian conclusion is not guaranteed by construction. The subsequent overlay of Keplerian curves for 0.15 and 0.39 M_sun on the independent H2CO and CH3OH PV diagrams is presented as a separate consistency check, not as a fit to those lines; no fitted parameter from H2CO/CH3OH is fed back into the mass estimate. The inclination i = 47° is derived geometrically from the continuum axis ratio under an explicitly stated thin-circular-disk assumption; the later flared-disk discussion in Section 5.1 is a caveat about that assumption, not a circular re-use of the mass result. The self-citations that appear (eDisk overview Ohashi et al. 2023, SLAM Aso & Sai 2024, ridge/edge bias Aso et al. 2015, and RADMC-3D model Takakuwa et al. 2024) are methodological, code-descriptive, or interpretive. The RADMC-3D model was not fitted to IRAS 04166+2706, and the paper explicitly notes that alternative explanations for the minor-axis asymmetry 'cannot be ruled out' and that a 'more quantitative comparison' is needed. No uniqueness theorem from the same authors is invoked to forbid alternative interpretations, and no equation reduces by construction to its own inputs. The analysis is therefore not circular; concerns about cut orientation, infall contamination, or the flared-disk geometry are robustness and model-selection issues, not derivational circularity.
Assumptions & free parameters
free parameters (5)
- Dust temperature T_d =
20 K (34 K alternative)
- Dust mass opacity κ_1.3mm =
2.3 cm²/g (1.0 cm²/g alternative)
- Gas-to-dust ratio =
100
- Distance =
156 pc
- Inclination i =
47°
assumptions (6)
- domain assumption Dust emission is optically thin and isothermal (Eq. 1) for the disk mass estimate.
- domain assumption Gas in the fitted PV regions follows Keplerian rotation, v = sqrt(GM⋆/r), at inclination 47°.
- domain assumption The 2D Gaussian deconvolution of the continuum gives the disk geometry (radius, PA, inclination) of a single thin circular disk.
- ad hoc to paper The selected PV velocity ranges isolate rotation from infall/outflow.
- domain assumption Dust temperature follows T_d = 43(L_bol/L_⊙)^0.25 for the alternative mass estimate.
- domain assumption The flared, optically thick disk model of Takakuwa et al. 2024 explains the minor-axis brightness asymmetry.
Cite this review
Pith. "Pith review of Early Planet Formation in Embedded Disks (eDisk) XXII: Keplerian disk, disk structures and jets/outflows in the Class 0 protostar IRAS 04166+2706." pith.science (2026). https://pith.science/paper/YALBQTCI
@misc{pith2026250807212,
author = {Pith},
title = {Pith review of: Early Planet Formation in Embedded Disks (eDisk) XXII: Keplerian disk, disk structures and jets/outflows in the Class 0 protostar IRAS 04166+2706},
year = {2026},
howpublished = {\url{https://pith.science/paper/YALBQTCI}},
note = {Machine review of arXiv:2508.07212}
}
abstract
We present ALMA observations of the Class 0 protostar IRAS 04166+2706, obtained as part of the ALMA large program Early Planet Formation in Embedded Disks (eDisk). These observations were made in the 1.3 mm dust continuum and molecular lines at angular resolutions of $\sim 0.05''$ ($\sim 8$ au) and $\sim 0.16''$ ($\sim25$ au), respectively. The continuum emission shows a disk-like structure with a radius of $\sim22$ au. Kinematical analysis of $^{13}$CO(2-1), C$^{18}$O(2-1), H$_2$CO (3$_{0,3}$-2$_{0,2}$), CH$_3$OH (4$_2$-3$_1$) emission demonstrates that these molecular lines trace the infalling-rotating envelope and possibly a Keplerian disk, enabling us to estimate the protostar mass to be $0.15 \rm{M_\odot} < \rm{M_\star} < 0.39 M_\odot$. The dusty disk is found to exhibit a brightness asymmetry along its minor axis in the continuum emission, probably caused by a flared distribution of the dust and the high optical depth of the dust emission. In addition, the CO(2-1) and SiO(5-4) emissions show knotty and wiggling motions in the jets. Our high angular resolution observations revealed the most recent mass ejection events, which have occurred within the last $\sim 25$ years.
Figures
Figures from the paper (12 more)
Reference graph
Works this paper leans on
-
[1]
Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 369, 229, doi: 10.1111/j.1365-2966.2006.10294.x Alves, F. O., Cleeves, L. I., Girart, J. M., et al. 2020, ApJL, 904, L6, doi: 10.3847/2041-8213/abc550 Andrews, S. M., & Williams, J. P. 2005, ApJ, 631, 1134, doi: 10.1086/432712 Aso, Y., & Sai, J. I. C. 2024, Publications of the Korean Astronomi...
arXiv 2006
-
[27]
2015, ApJ, 812, 27, doi: 10.1088/0004-637X/812/1/27 Beckwith, S
http://pkas.kas.org/journal/article.php?code=91201 Aso, Y., Ohashi, N., Saigo, K., et al. 2015, ApJ, 812, 27, doi: 10.1088/0004-637X/812/1/27 Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 1990, AJ, 99, 924, doi: 10.1086/115385 Booth, R. A., & Clarke, C. J. 2016, MNRAS, 458, 2676, doi: 10.1093/mnras/stw488 Chen, H., Myers, P. C., Ladd, E....
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.