Pith. sign in

REVIEW 3 major objections 6 minor 26 references

Radial velocities of narrow emission line components in the spectra of T Tauri stars

T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

Pith's one-line read The paper argues that the 2-7 km/s red shifts of narrow helium emission lines in four T Tauri stars are not caused by infalling accreted gas, because such inflow would phase-shift the radial velocity curves of lines with different shifts…

desk verdict A clean phase-shift test that weakens the Doppler-infall reading of helium line shifts, but the categorical conclusion rests on a radial-infall assumption and a partly circular Stark fit. read the letter →

arxiv 2412.06362 v1 pith:XQMG4IWD submitted 2024-12-09 astro-ph.SR

classification astro-ph.SR
keywords TTauristarsaccretiondiscsshockwavesradialvelocitiesheliumemissionlinesStarkeffectrotationalmodulation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Narrow helium emission lines in classical T Tauri stars are known to be red-shifted by several km/s relative to the star, and this has usually been read as gas still settling after passing through the accretion shock. The paper tests that interpretation by measuring the radial velocity curves of helium and metal lines in four stars (BP Tau, DK Tau, EX Lup, and TW Hya) across many spectra. If the shifts were Doppler shifts of radially inflowing gas, lines with different mean shifts should have radial velocity curves that lag or lead each other by a calculable amount; the measured curves instead move in phase, within about 10 degrees. The paper concludes that the helium-line shifts are not gas motion: the neutral helium shifts can be fully explained by the Stark effect and large optical depth of the lines at the density of the accretion column base, while the HeII shift remains without a definite explanation. The result changes what the line shifts can be used to measure.

What carries the argument

The central object is the phase-shift test embodied in Eqs. (11)--(12) of the paper. For a spot on a rotating star, a line formed in gas at rest relative to the surface has a purely sinusoidal radial-velocity curve; adding a radial inflow term $B\cos\varphi + C$ puts the two sinusoids in quadrature, so the observed mean shift $C$ and the phase offset $\Delta\varphi$ of the combined curve are tied together through $\cos\Delta\varphi = v\sin i / \sqrt{(v\sin i)^2 + (C\tan i/\cos\theta)^2}$. Larger shifts therefore force larger phase offsets, and comparing the predicted lower limits with the observed, near-zero phase offsets is the decisive test. On the interpretation side, the companion mechanism is the internal wavelength structure of the neutral helium lines: HeI 5876 is a blend of fine-structure components, so its measured center moves with optical depth, and the density-sensitive Stark shift moves the line center at the same time, together reproducing the observed shifts at $N \approx 10^{15}$ cm$^{-3}$ and $\tau \gg 1$.

What would settle it

A decisive test would be to find a classical T Tauri star in which the narrow HeI 6678 line shows a mean red shift of about 2 km/s relative to the stellar frame and, over several rotation cycles, its radial velocity curve lags or leads the metal-line curve by the amount Eq. (12) predicts for the star's $v\sin i$ and spot latitude; such a phase offset would contradict the paper's central claim.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central claim is that the persistent 2-7 km/s red shifts of the narrow HeI 5876, HeI 6678, and HeII 4686 emission components in classical T Tauri stars are not produced by motion of gas relative to the stellar surface. The supporting argument is a phase test: for radial inflow, the velocity curve of a line is the sum of a rotationally modulated sinusoid and a constant inflow term, so the mean velocity shift $C$ and the phase offset $\Delta\varphi$ of the sinusoid are linked through Eq. (12); larger mean shifts require larger phase offsets. For the four stars observed, the phase offsets between helium and metal lines are consistent with zero at the $\sim 10^\circ$ level, whereas the lower limits implied by the observed shifts are 15--44 degrees. The authors therefore conclude that the shifts are not Doppler shifts. For neutral helium, they show that the observed shifts can be reproduced by a combination of the Stark effect at electron density near $10^{15}$ cm$^{-3}$ and the optical-thickness-dependent blending of the HeI 5876 fine-structure components; for HeII 4686, no non-Doppler explanation is identified, but the zero phase shift still argues against gas motion.

Load-bearing premise

The load-bearing premise is that the infalling gas moves radially toward the star, because the paper's own appendix shows that an oblique inflow could hide the phase shifts that the radial-inflow formula predicts, and the authors set oblique inflow aside as inconsistent with the standard accretion picture.

Editorial extensions

If this is right

  • If the phase test is right, the mean red shifts of the helium lines cannot be used as measures of accretion-column infall or settling velocities; the paper's Table 6 places upper limits near 0.4--2.6 km/s instead.
  • The HeI 6678 shift becomes a density diagnostic for the line-forming region: the observed ~2 km/s shift implies an electron density around $10^{15}$ cm$^{-3}$, consistent with hotspot models.
  • The HeI 5876 shift, combined with the HeI 6678 density, constrains the optical depth of the line to $\tau \gg 1$, meaning the narrow neutral helium lines form in optically thick gas at the base of the accretion column.
  • The phase test can be applied to any star with sufficiently accurate multi-line velocity curves as a general check of whether emission-line shifts are kinematic or non-kinematic in origin.
  • For BP Tau, the in-phase variability of emission and absorption lines points to photometric distortions such as Rossiter--McLaughlin-type eclipses by the inner disk rather than to a rotating hotspot.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: if the zero-phase-shift result holds over longer baselines, the common practice of estimating accretion-zone latitude from helium-line velocity amplitudes while treating the mean shift as a free radial velocity should be revisited; latitude estimates from earlier studies may be biased if the offset is not purely geometric.
  • Editorial inference: the same phase-shift test could be applied to redshifted absorption components in Balmer or other helium profiles, which are also attributed to infall; those components should show phase behavior distinct from narrow emission if they arise in the pre-shock flow.
  • Editorial inference: since the paper leaves HeII 4686 unexplained, a testable extension is to compare high-signal-to-noise HeII 4686 profiles across many rotational phases to look for a non-modulated red-shifted emission component or blue absorption whose subtraction would account for the shift.
  • Editorial inference: because the appendix shows that oblique infall can cancel the phase shift, a decisive way to close the loophole is to couple the phase test with independent magnetospheric geometry constraints, such as spectropolarimetric magnetic field maps, to check whether the radial-infall assumption is actually justified in these stars.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper measures radial velocities of narrow emission-line components (HeI 5876, HeI 6678, HeII 4686, and metal lines) in four classical T Tauri stars (BP Tau, DK Tau, EX Lup, TW Hya) using archival ESPaDOnS and FEROS spectra. The authors confirm previously reported redshifted helium-line velocities of 2–7 km/s relative to the stellar rest frame. They argue that if these shifts were Doppler signatures of radially infalling gas, radial velocity curves of lines with different mean shifts would be phase-shifted relative to each other (Eqs. 11–12). They report observed phase shifts consistent with zero (Table 5) against predicted lower limits of 15–44 deg (Table 4), and conclude that the helium-line shifts are not caused by gas inflow. For neutral helium, they propose that the shifts arise from the Stark effect at electron density near 10^15 cm^-3 plus large optical depth for HeI 5876. The HeII 4686 shift is left unexplained.

Significance. If the phase-shift test is valid, the paper offers a novel observational argument against the standard interpretation of redshifted helium lines as residual infall in T Tauri hotspots. The analysis has notable strengths: the use of a non-accreting template, careful treatment of line blending and veiling, a clean analytic prediction for radial inflow, and an honest admission that HeII 4686 remains unexplained. The phase-shift test is in principle falsifiable and is applied to archival data that are not tailored to the result. However, the central conclusion depends on the radial-infall assumption, and the phase-shift constraints are weak for EX Lup and TW Hya. The Stark/optical-depth explanation for the HeI lines is a plausible consistency check rather than a unique or parameter-free determination.

major comments (3)
  1. [Appendix and HeII interpretation] The phase-shift null test does not exclude gas inflow in general. Equation (13) of the appendix gives the velocity curve as v = A' sin φ + B' cos φ + C with B' = (vr_in sin θ - vθ_in cos θ) sin i. For a non-radial inflow with vθ_in = vr_in tan θ, B' = 0 while C = vr_in cos i / cos θ is nonzero, so a line can have a nonzero mean velocity shift and zero phase shift relative to a stationary line. The HeII section's statement that a radial or poloidal component 'in this scenario, phase shifts would occur' is therefore contradicted by the paper's own appendix. Since the observed mean shifts C are 2–7 km/s, the data are compatible with a family of oblique-infall solutions. The abstract's claim that the shifts are 'not associated with the inflow of accreted gas' is accordingly supported only under the radial-infall assumption, which the authors adopt as a prior rather than derive from observations. The conclusion should be restricted to radial infall, or an observational constraint on vθ_in should be provided.
  2. [Observed phase shifts, Tables 4 and 5] The significance of the null result is overstated for two of the four stars. For EX Lup HeI 6678, the observed phase shift 12.9°±7.8° is consistent with the Table 4 lower limit of 17° at about 1σ; for TW Hya HeI 5876, 9.8°±18.0° against 15° is entirely consistent, and HeI 6678 5.7°±10.2° against 15° is consistent at <1σ. Only DK Tau (all lines) and EX Lup HeII 4686 (-3.3°±16.8° versus 44°, about 2.8σ) provide a meaningful exclusion. The text should report the per-line discrepancy in units of σ and should temper the abstract's categorical statement that the observed phase shifts 'do not correspond to the observed line velocity shifts.'
  3. [Interpretation of the observed line shifts, HeI] The HeI Stark/optical-depth explanation is presented as if it were a confirmation, but it is a two-parameter fit to the same data being explained. The electron density N ≈ 10^15 cm^-3 is derived directly from the observed HeI 6678 shift using the linear Stark calibration of Dimitrijevic and Sahal-Brechot (1990); the HeI 5876 shift is then matched by adjusting τ at that density. The Fig. 8 curves are 'calculated based on the observed densities derived from the HeI 6678' line, so the agreement for HeI 5876 is not an independent test. The paper should explicitly describe this as a consistency check, state the uncertainties on the inferred (N, τ), and compare with the model values at the base of the accretion column rather than claiming the shifts can be 'completely explained' without qualification.
minor comments (6)
  1. [Throughout] The unit 'kms−1' should be written with a space as 'km s−1' (e.g., in the Introduction and Section 2).
  2. [Title/first page] The running header 'Astronomy Letters, 2100, vol 1, № 5, p. 1–1' appears to be a template placeholder and should be corrected to the actual year, volume, and page numbers.
  3. [Observations] The SNR range '11 to 75' is given before the statement that spectra with SNR < 20 were excluded; please clarify that the range refers to the full archive before the cut.
  4. [Observed phase shifts] The statement that the uncertainties are '∼10°' is not consistent with Table 5, where TW Hya HeI 5876 has an uncertainty of 18°; please quote per-line uncertainties.
  5. [Observed phase shifts, Table 5] The method for computing the phase-shift uncertainties in Table 5 is not described; please state whether they come from the least-squares covariance, the period uncertainty, or a bootstrap.
  6. [Measurement of radial velocities, Eq. (5)] Equation (5) defines the gf-weighted central wavelength for the optically thin limit; it would be helpful to explicitly note that for the optically thick and Stark-shifted cases discussed later, this definition is not the exact zero-velocity reference, although the paper's analysis accounts for this.

Circularity Check

0 steps flagged · score 2.0 of 10

No load-bearing circularity: the phase-shift null test is independent; the HeI Stark/optical-depth explanation is a consistency check with a free parameter, not a forced derivation.

full rationale

The central phase-shift test is self-contained and does not reduce to its inputs. Expected phase shifts are computed from Eq. (12) using the measured mean line shifts C, the stellar v sin i, and inclination i, and then compared with independently fitted phase lags in Table 5. The observed null phase shifts are not algebraically forced by C, so this is a genuine kinematic test under the stated radial-inflow assumption. The appendix Eq. (13) shows that oblique infall can yield zero phase shift with nonzero mean shift, so the abstract's categorical wording is stronger than the test supports; this is a modeling assumption/limitation rather than circularity. The HeI interpretation derives N ~ 10^15 cm^-3 from the observed HeI 6678 shift via the Stark calibration and then uses that N with a free optical depth tau to match HeI 5876; this is a consistency check with a free parameter, not an independent prediction, and the match is not strictly forced because tau is unconstrained. The expected accretion-column density is supported by a self-citation to Dodin 2018, but that citation is not the sole support for the central phase-shift claim. Therefore no load-bearing circularity is present; at most a minor non-load-bearing self-citation justifies a score of 2 rather than 0.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The central phase-shift argument rests on standard kinematics plus the domain assumption of radial infall; the HeI Stark explanation additionally relies on density inferred from one of the observed shifts and a free optical depth. No new physical entities are introduced.

free parameters (5)
  • Stellar parameters (v sin i, inclination i, period P) = See Table 2
    Adopted from literature or fitted to absorption-line RV curves; the predicted phase-shift lower limits in Table 4 depend on these values.
  • Density inferred from HeI 6678 Stark shift = ~1e15 cm^-3
    Derived from the observed HeI 6678 A shift using the linear Stark calibration of Dimitrijevic and Sahal-Brechot (1990); used to compute the HeI 5876 shift curves in Fig. 8.
  • Optical depth tau of HeI 5876 = tau >> 1
    Not measured; the observed HeI 5876 shift is matched at large optical depth in Fig. 8, making the Stark/optical-depth explanation a consistency check rather than a prediction.
  • Gaussian broadenings for template = sigma = 6.7 km/s (DK Tau), 3.4 km/s (BP Tau); TAP45 v sin i = 7 km/s
    Chosen to match the TAP45 template to target stars; affects the accuracy of absorption-line subtraction and hence measured emission-line velocities.
  • Periods for phase curves = DK Tau 9.3902 d, EX Lup 7.417 d, TW Hya 3.568 d, BP Tau 8.6681 d
    Selected from periodograms or fits; the phase shifts of RV curves are computed at these periods, and the paper notes alternate periods for DK Tau give the same null result.
assumptions (5)
  • domain assumption Infalling gas moves nearly radially onto the stellar surface
    Equations (10)-(12) and the phase-shift predictions in Fig. 3 assume radial gas motion; the appendix generalizes to arbitrary angles but the paper sets oblique infall aside.
  • domain assumption Metal emission lines are stationary tracers of the hotspot
    The paper uses the average of six metal lines as the zero-shift reference for computing phase shifts of helium lines.
  • domain assumption Gaussian decomposition of narrow and broad line components is adequate
    Equation (6) assumes each component is Gaussian and that the broad components of HeI 5876 and 6678 share velocity and width; this constrains the measured narrow-line velocities.
  • domain assumption Template subtraction does not bias line centers
    The paper subtracts TAP45 or pySME synthetic spectra with veiling to remove absorption blends; later it concedes that residual template mismatch could mimic shifts, especially for HeII.
  • standard math Shock jump conservation equations for an ideal gas
    Used in the introduction to argue that a 10 km/s post-shock velocity implies temperatures where neutral helium cannot exist.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Radial velocities of narrow emission line components in the spectra of T Tauri stars." pith.science (2026). https://pith.science/paper/XQMG4IWD

@misc{pith2026241206362,
  author       = {Pith},
  title        = {Pith review of: Radial velocities of narrow emission line components in the spectra of T Tauri stars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XQMG4IWD}},
  note         = {Machine review of arXiv:2412.06362}
}
read the original abstract

We studied rotational modulation of the radial velocities of narrow emission lines in four classical T Tauri stars. We found that the previously declared shift of the mean velocity of neutral and ionized helium lines relative to the mean radial velocity of the star is not associated with the inflow of accreted gas into the hotspot, since the radial velocity curves for lines with different velocity shifts should exhibit phase shifts relative to each other, while the observed phase shifts are absent within their uncertainties and do not correspond to the observed line velocity shifts. This means that the line shifts are not caused by the actual gas motion. For neutral helium lines, the shifts can be explained by the large optical thickness of the lines and the Stark effect at plasma parameters expected at the base of the accretion column of T Tauri stars.

Figures

Figures reproduced from arXiv: 2412.06362 by the authors.

Figure 1
Figure 1. Removing the stellar absorption lines for [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Schematic view of a star with a spot. where v represents the linear rotational velocity of the star at the equator, φ denotes the longitude, θ is the angle between the rotation axis and the direction to the spot (latitude), and i is the angle between the rotation axis and the line of sight (see [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Possible values of phase shifts of the He II [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: Phase curves of radial velocities of helium [PITH_FULL_IMAGE:figures/full_fig_p008_4.png]
Figure 6
Figure 6. Figure 6: The phase curves of the radial velocities [PITH_FULL_IMAGE:figures/full_fig_p009_6.png]
Figure 5
Figure 5. Figure 5: The phase curves of the radial velocities of [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]
Figure 7
Figure 7. Figure 7: The radial velocity curve of He II 4686 A˚ (black curve) and He I 5876 A˚ (blue dashed curve) for the star BP Tau at the epoch MJD0 = 56667.41252. is comparable to the value of v sin i of the stars. If these shifts are caused by gas accreting onto the star, then signif…
Figure 8
Figure 8. Figure 8: The shift in km s−1 of the He I 5876 A˚ line as a function of the logarithm of the optical depth log τ. Red lines correspond to a temperature of 20 000 K, while blue lines correspond to 10 000 K. The curves are calculated based on the observed densities derived from th…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    Bouvier, S

    J. Bouvier, S. Cabrit, M. Fernandez, E. L. Martin, and J. M. Matthews, Astron Astrophys.,272, 176 (1993)

  2. [2]

    Bouvier, A

    J. Bouvier, A. Chelli, S. Allain, et al., Astron Astrophys., 349, 619 (1999)

  3. [3]

    Colombo, S

    S. Colombo, S. Orlando, G. Peres, C. Argiroffi, and F. Reale, Astron Astrophys.,594, A93 (2016)

  4. [4]

    M. S. Dimitrijevic and S. Sahal-Brechot, Astron. Astrophys. Suppl. Ser.,82, 519 (1990)

  5. [5]

    Dodin, MNRAS,475(4), 4367 (2018)

    A. Dodin, MNRAS,475(4), 4367 (2018)

  6. [6]

    A. V. Dodin and E. A. Suslina, MNRAS,503(4), 5704 (2021)

  7. [7]

    J.F.Donati,S.G.Gregory,S.H.P.Alencar,etal., MNRAS, 417(1), 472 (2011)

  8. [8]

    J. F. Donati, M. M. Jardine, S. G. Gregory, et al., MNRAS, 386(3), 1234 (2008)

Show all 26 references
  1. [9]

    Hartmann, G

    L. Hartmann, G. Herczeg, and N. Calvet, ARAA, 54, 135 (2016)

  2. [10]

    G. J. Herczeg, Y. Chen, J.-F. Donati, et al., Astrophys. J., 956(2), 102 (2023)

  3. [11]

    G. J. Herczeg and L. A. Hillenbrand, Astrophys. J., 786(2), 97 (2014)

  4. [12]

    K´ osp´ al, M

    ´A. K´ osp´ al, M. Mohler-Fischer, A. Sicilia-Aguilar, et al., Astron Astrophys.,561, A61 (2014)

  5. [13]

    F. Long, G. J. Herczeg, D. Harsono, et al., Astrophys. J., 882(1), 49 (2019)

  6. [14]

    McGinnis, J

    P. McGinnis, J. Bouvier, and F. Gallet, MNRAS, 497(2), 2142 (2020)

  7. [15]

    Nelissen, P

    M. Nelissen, P. McGinnis, C. P. Folsom, et al., Astron Astrophys., 670, A165 (2023)

  8. [16]

    J. R. Percy, S. Grynko, R. Seneviratne, and W. Herbst, PASP,122(893), 753 (2010)

  9. [17]

    P. P. Petrov, G. F. Gahm, J. F. Gameiro, et al., Astron Astrophys., 369, 993 (2001)

  10. [18]

    Piskunov and J

    N. Piskunov and J. A. Valenti, Astron Astrophys., 597, A16 (2017)

  11. [19]

    Pouilly, M

    K. Pouilly, M. Audard, A. Lavail, and´A. K´ osp´ al, arXiv e-prints, arXiv:2409.03322 (2024)

  12. [20]

    Pouilly, J

    K. Pouilly, J. Bouvier, E. Alecian, et al., Astron Astrophys., 656, A50 (2021)

  13. [21]

    Ryabchikova, N

    T. Ryabchikova, N. Piskunov, R. L. Kurucz, et al., Physica Scripta, 90(5), 054005 (2015)

  14. [22]

    Sicilia-Aguilar, M

    A. Sicilia-Aguilar, M. Fang, V. Roccatagliata, et al., Astron Astrophys.,580, A82 (2015)

  15. [23]

    Singh, J

    K. Singh, J. P. Ninan, M. M. Romanova, et al., arXiv e-prints, arXiv:2404.05420 (2024)

  16. [24]

    V. Y. Terebigh, Analiz vremennykh ryadov v astrofizike (in russian)(M.:Nauka. Gl.red.fiz.-mat. lit., 1992)

  17. [25]

    J. A. Valenti and N. Piskunov, Astron. Astrophys. Suppl. Ser., 118, 595 (1996)

  18. [26]

    Wehrhahn, N

    A. Wehrhahn, N. Piskunov, and T. Ryabchikova, Astron Astrophys., 671, A171 (2023) APPENDIX. THE RADIAL VELOCITY OF THE LINE IN THE CASE OF GAS INFLOW AT AN ARBITRARY ANGLE TO THE SPOT. When matter falls onto the surface of a star at an arbitrary angle, the velocity vector of t...

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.