REVIEW 3 major objections 5 minor 1 cited by
Impact of the Tayler magnetic instability on the surface abundance of boron in massive stars
T0 review · 3 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Purely hydrodynamic rotating models overestimate boron depletion in fast B-type stars, and adding asteroseismically calibrated magnetic Tayler transport fixes the discrepancy.
desk verdict A useful, honest test of Tayler transport against boron data, but the main conclusion stands only if you already buy GENEC's advective treatment of meridional circulation. 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 central object is the pair of transport coefficients in the rotating-star equations: the shear diffusion coefficient $D_{\rm shear}$ and the effective chemical diffusion coefficient $D_{\rm eff}$ from meridional advection combined with horizontal turbulence. In non-magnetic models, meridional circulation builds radial differential rotation, making $D_{\rm shear}$ dominate chemical transport in the outer layers and driving rapid boron destruction. Adding the magnetic Tayler instability, a magnetohydrodynamic instability that transports angular momentum in stably stratified radiative zones, with the asteroseismically calibrated strength ($n=1$, $C_T=216$) flattens the rotation profile, suppresses $D_{\rm shear}$, and leaves $D_{\rm eff}$ as the main mixing agent. Because $D_{\rm eff}$ grows only as the four-thirds power of the meridional velocity, the faster meridional flow in magnetic models only partially compensates, so the net boron depletion is reduced.
What would settle it
Asteroseismic measurement of a moderate-to-fast rotating B-type star that shows strong radial differential rotation while its surface boron is only mildly depleted would break the claimed link between differential rotation, shear mixing, and boron destruction.
Extended reading notes
Core claim
Using stellar evolution models of $9$, $12$, and $15\,M_\odot$ stars with and without magnetic transport, the paper finds that purely hydrodynamic rotating models overestimate boron depletion for stars with high rotation rates, in disagreement with observed surface abundances. The excessive mixing is traced to the shear instability in the outer radiative layers, which is fed by the strong radial differential rotation that meridional circulation creates in the advective treatment. When the magnetic Tayler instability is added using the asteroseismically calibrated prescription ($n=1$, $C_T=216$), angular-momentum transport keeps the rotation profile nearly flat, shear mixing drops, and the predicted surface boron abundances match the evolutionary states and projected velocities of moderately and fast-rotating B-type stars without adjusting any free parameter. At low projected velocities ($v\sin i\lesssim50$ km/s), the magnetic models deplete too little boron, which the paper interprets as a sign that current Tayler prescriptions may overestimate angular-momentum transport in slow rotators.
Load-bearing premise
The entire boron result follows from the asteroseismically calibrated Tayler-instability prescription flattening the rotation profile of 9-15 solar-mass main-sequence B stars; if that prescription is wrong in this mass and temperature range, the predicted boron abundances and the criticism of hydrodynamic models would change.
Editorial extensions
If this is right
- For initial rotation rates $v_{\rm ini}/v_{\rm crit}\gtrsim0.3$, purely hydrodynamic models predict more boron depletion than observed, and the discrepancy grows as the star evolves on the main sequence.
- Magnetic models built on the asteroseismically calibrated Tayler transport reproduce surface boron, surface gravity, and projected rotation velocity together for moderate and fast rotators, with no free parameter adjusted.
- The magnetic models keep fast rotators at higher equatorial velocities than non-magnetic models, matching the fastest projected velocities in the observed sample.
- At $v\sin i\lesssim50$ km/s the non-magnetic models fit the boron data better, suggesting the Tayler prescription over-transports angular momentum in slow rotators.
- The surface-boron conclusions are nearly independent of which Tayler-Spruit dynamo prescription is used, because the outer rotation profiles are similar for the original and revised prescriptions.
Reading between the lines
- Extending the comparison to beryllium and lithium, which are destroyed at lower temperatures than boron, would map the transport profile at several depths and provide a sharper test of the magnetic models.
- If slow rotators really need weaker angular-momentum transport, the asteroseismic calibration of the Tayler instability could be made rotation-rate dependent; the low-$v\sin i$ boron data would then constrain that dependence.
- Because advective and diffusive treatments of meridional circulation respond oppositely to magnetic transport, conclusions about magnetic mixing drawn from diffusive stellar-evolution codes should be re-examined before being generalised to advective codes.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper computes Geneva stellar evolution models of 9, 12, and 15 solar-mass stars with an advective treatment of meridional circulation, with and without transport by the magnetic Tayler instability, and compares the predicted surface boron abundances with B-star observations from Jin et al. (2024) and Proffitt et al. (2024). The central claim is that purely hydrodynamic rotating models over-deplete boron at moderate and high rotation rates because they develop strong differential rotation and efficient shear mixing, while models with an asteroseismically calibrated Tayler-instability transport (Eggenberger et al. 2022, n=1, CT=216) flatten the rotation profile, reduce shear mixing, and reproduce the observed boron abundances of moderately and fast-rotating B-type stars without tuning parameters to the boron data. An appendix (A.2) shows that with the purely diffusive MESA scheme used by the observational comparison papers, non-magnetic hydrodynamic models have flat rotation profiles and do not over-deplete boron, and magnetic transport barely changes the boron abundance. The paper also examines robustness across different Tayler-instability prescriptions and post-main-sequence evolution.
Significance. If the central claim holds, the paper is significant because it offers a single framework that reconciles asteroseismic constraints on internal angular momentum transport with boron surface-abundance constraints on mixing in main-sequence B stars. The paper has genuine strengths: it uses recently measured boron abundances of fast rotators, it explicitly includes an appendix (A.2) that demonstrates the scheme dependence of the non-magnetic result, it checks robustness across the Spruit (2002), Fuller et al. (2019), and Eggenberger et al. (2022) prescriptions (Appendix B.2), and it extends the comparison to post-main-sequence stars (Appendix C). The significance as a model-independent observational conclusion is, however, not yet established because the key claim depends on the choice of the advective treatment of meridional circulation, and the paper does not provide an independent validation of that treatment for B-type stars.
major comments (3)
- [§3, §4, Appendix A.2] The central claim is not model-independent as stated. In the diffusive MESA scheme with fc=0.017 and fmu=0.1, the same choices used by the observational comparison papers, non-magnetic 15 solar-mass models have nearly flat rotation profiles (Fig. A.4) and do not over-deplete boron (Fig. A.6); magnetic transport barely changes the boron abundance. The statements in Sections 3 and 4 that 'models with only hydrodynamic transport processes overestimate the amount of boron depletion' and that boron abundances 'indicate that a more efficient AM transport is needed' are therefore valid only within the advective framework adopted in GENEC. The paper argues that the advective scheme is more physical, but it provides no independent test of that scheme for main-sequence B stars. Since Proffitt et al. (2024) and Jin et al. (2024) use the diffusive scheme, the disagreement is between two modelling choices rather than a robust model-observation discrepancy. The authors should either provide a direct test of the advective scheme's rotation and mixing properties for these stars or explicitly and prominently restrict the conclusions to the advective framework, with a corresponding revision of the abstract and conclusion.
- [§3, Figs. 1 and 2] The claimed agreement of the magnetic models and the incompatibility of the non-magnetic models are assessed only visually. No quantitative goodness-of-fit measure, treatment of observational upper limits, or propagation of stellar parameter uncertainties is presented. Since the letter's main positive claim is that the magnetic models are 'in good agreement' with the boron data, a simple quantitative statistic (for example, a likelihood that accounts for upper limits and log g uncertainties) or an explicit statement of the limitations of the visual comparison is needed to support the strength of the conclusion.
- [§4 and Appendix B.2] The sentence 'we find that this conclusion does not depend on the prescription adopted for the exact modelling of the TS dynamo' is stronger than what Appendix B.2 demonstrates. Figure B.3 shows that the Spruit (2002) prescription gives results similar to the calibrated one in the external layers, but the quantitative differences at the end of the main sequence are not evaluated, and the statement rests on visual similarity. A quantitative comparison of the boron predictions across the three prescriptions would make the robustness claim precise.
minor comments (5)
- [§2] The parameters n and CT of the Eggenberger et al. (2022) prescription are introduced without a definition; since the robustness of the conclusion across prescriptions is emphasized, a one-sentence definition of these parameters would be helpful.
- [Appendix B.1] In the expression Dh = Ar(rΩV(2V−αU))^(1/3), the symbol V is used without definition; it appears only in this appendix, and U was defined earlier as the vertical component of the meridional circulation velocity.
- [Fig. 2] The repeated '9M', '12M', '15M' labels in the left panel are confusing and overlap the tracks; placing the mass labels on individual tracks or in the caption would improve readability.
- [§4] 'We have showed' should be 'We have shown'. In addition, 'without the need to adjust any free parameter' should be qualified as 'without adjusting parameters to the boron data', because CT=216 is itself a parameter calibrated to asteroseismic data.
- [References] The reference to Proffitt et al. (2024) in the bibliography lacks journal, volume, and page information (or a DOI); it should be completed.
Circularity Check
No significant circularity: boron surface abundances are predicted outputs compared with external observations, not fitted inputs.
full rationale
The paper's derivation chain is not circular. The adopted Tayler-instability transport prescription (Eggenberger et al. 2022, n=1, CT=216) was calibrated on asteroseismic measurements of internal rotation, not on boron abundances, so the subsequent comparison with observed boron abundances is an external test rather than a fit. The central result that magnetic models reproduce boron abundances is also shown in Appendix B.2 to be nearly insensitive to the choice of TS-dynamo prescription (Spruit 2002; Fuller et al. 2019; Eggenberger et al. 2022), so the conclusion does not rest solely on a self-citation. Appendix A.2 explicitly shows that non-magnetic MESA models with a diffusive treatment of meridional circulation predict flat rotation profiles and only moderate boron depletion, stating that 'the fundamental difference in the treatment of meridional circulation between the advective scheme used in the present work with GENEC models and the purely diffusive scheme used, for instance, in MESA can lead to opposite results for the transport of AM and chemicals.' This is a genuine model-dependence caveat rather than a circular reduction: the paper acknowledges this opposite result and argues on physical grounds for the advective scheme. No equation or parameter is defined in terms of the target boron abundances, and no boron-derived quantity is renamed as a prediction. The self-citation to Eggenberger et al. (2022) is therefore not load-bearing in the sense prohibited by the review rules.
Assumptions & free parameters
free parameters (2)
- CT (Tayler-Spruit angular momentum transport calibration) =
216
- fc and fmu in MESA comparison models =
fc = 0.017, fmu = 0.1
assumptions (6)
- domain assumption Shellular rotation: angular velocity is approximately constant on isobars (Zahn 1992).
- domain assumption Advective treatment of meridional circulation is physically appropriate, while the purely diffusive approximation misapplies meridional velocities.
- domain assumption The adopted expressions for D_shear (Talon & Zahn 1997) and D_h (Maeder 2003) describe the turbulent transport.
- domain assumption The Eggenberger et al. (2022) Tayler instability prescription (n=1, CT=216), calibrated on asteroseismic data, is valid for main-sequence B-type stars of 9-15 solar masses.
- domain assumption Solar metallicity is representative of the majority of the observed sample.
- domain assumption Boron is destroyed by proton captures at temperatures below about 6 million kelvin, so surface boron depletion traces mixing from the outer radiative zone.
Cite this review
Pith. "Pith review of Impact of the Tayler magnetic instability on the surface abundance of boron in massive stars." pith.science (2026). https://pith.science/paper/4X5ODN3K
@misc{pith2026250704267,
author = {Pith},
title = {Pith review of: Impact of the Tayler magnetic instability on the surface abundance of boron in massive stars},
year = {2026},
howpublished = {\url{https://pith.science/paper/4X5ODN3K}},
note = {Machine review of arXiv:2507.04267}
}
read the original abstract
Context: The surface abundances of massive stars show evidence of internal mixing, while asteroseismic data suggest that efficient angular momentum (AM) transport occurs in stellar interiors. It is of interest to find a consistent physical framework that is able to account for both of these effects simultaneously. Aims: We investigate the impact of the Tayler instability on the surface abundance of boron in massive B-type stars as predicted by rotating stellar models accounting for the advective nature of meridional currents. Methods: We used the Geneva stellar evolution code to compute models of 9, 12, and 15 Msun stars at different rotational velocities and with and without magnetic fields. We compared the surface boron abundances predicted by these models with those of observed B-type stars. Results: We find that models with only hydrodynamic transport processes overestimate the amount of boron depletion for stars with high rotation rates, in disagreement with observational constraints. We show that this excessively high mixing efficiency is a consequence of the high degree of differential rotation predicted by purely hydrodynamic models. We thus conclude that surface abundances of boron indicate that a more efficient AM transport is needed in stellar radiative zones. We then studied the impact of the Tayler instability as a possible physical explanation to this issue. Models including this instability are found to be in good agreement with constraints on the surface boron abundances, the evolutionary state, and the projected rotational velocity of moderately and fast-rotating B-type stars. Finally, we note that at low rotational velocities, models with magnetic fields do not predict sufficient depletion to be consistent with the observations. This could suggest that the current prescriptions for the Tayler instability may overestimate the AM transport in slow-rotating B-type stars.
Figures
Forward citations
Cited by 1 Pith paper
-
The IACOB project XIX. Revisiting massive-star evolution with empirical TAMS constraints: updated models, overshoot calibration, and the population of blue supergiants
Massive-star models require mass-dependent core overshoot (α_ov ≈ 0.18–0.45) to match the empirical TAMS, but still fail to explain the velocity dependence of the TAMS and the observed blue supergiant population.
Reference graph
Works this paper leans on
-
[1]
2014, The Astrophysical Journal, 788, 93
Cantiello, M., Mankovich, C., Bildsten, L., Christensen-Dalsgaard, J., & Paxton, B. 2014, The Astrophysical Journal, 788, 93
work page 2014
-
[2]
Ceillier , T., Eggenberger , P., Garc \' a , R. A., & Mathis , S. 2013, , 555, A54
work page 2013
- [3]
-
[4]
W., Eggenberger , P., & Hirschi , R
den Hartogh , J. W., Eggenberger , P., & Hirschi , R. 2019, , 622, A187
work page 2019
-
[5]
Eggenberger , P., Buldgen , G., & Salmon , S. J. A. J. 2019, , 626, L1
work page 2019
-
[6]
2008, , 316, 43
Eggenberger , P., Meynet , G., Maeder , A., et al. 2008, , 316, 43
2008
-
[7]
2012, Astronomy & Astrophysics, 544, L4
Eggenberger, P., Montalbán, J., & Miglio, A. 2012, Astronomy & Astrophysics, 544, L4
work page 2012
-
[8]
Eggenberger, P., Moyano, F. D., & den Hartogh, J. W. 2022, Astronomy & Astrophysics, 664, L16
work page 2022
Show all 29 references
-
[9]
2010, , 522, A39
Frischknecht, U., Hirschi, R., Meynet, G., et al. 2010, , 522, A39
2010
-
[10]
L., & Jermyn, A
Fuller, J., Piro, A. L., & Jermyn, A. S. 2019, Monthly Notices of the Royal Astronomical Society, 485, 3661–3680
2019
-
[11]
J., & Proffitt, C
Jin, H., Langer, N., Lennon, D. J., & Proffitt, C. R. 2024, , 690, A135
2024
-
[12]
2009, Proceedings of the International Astronomical Union, 5, 411–420
Langer, N., Brott, I., Cantiello, M., et al. 2009, Proceedings of the International Astronomical Union, 5, 411–420
2009
-
[13]
2003, , 399, 263–269
Maeder, A. 2003, , 399, 263–269
2003
-
[14]
2009, Physics, Formation and Evolution of Rotating Stars ( Springer )
Maeder , A. 2009, Physics, Formation and Evolution of Rotating Stars ( Springer )
2009
-
[15]
& Meynet, G
Maeder, A. & Meynet, G. 2005, Astronomy & Astrophysics, 440, 1041–1049
2005
-
[16]
& Zahn , J.-P
Maeder , A. & Zahn , J.-P. 1998, , 334, 1000
1998
-
[17]
P., Goupil , M
Marques , J. P., Goupil , M. J., Lebreton , Y., et al. 2013, , 549, A74
2013
-
[18]
2011, , 192, 3
Paxton , B., Bildsten , L., Dotter , A., et al. 2011, , 192, 3
2011
-
[19]
2013, , 208, 4
Paxton , B., Cantiello , M., Arras , P., et al. 2013, , 208, 4
2013
-
[20]
2015, , 220, 15
Paxton , B., Marchant , P., Schwab , J., et al. 2015, , 220, 15
2015
-
[21]
B., et al
Paxton , B., Schwab , J., Bauer , E. B., et al. 2018, , 234, 34
2018
-
[22]
2019, , 243, 10
Paxton , B., Smolec , R., Schwab , J., et al. 2019, , 243, 10
2019
-
[23]
R., Jin, H., Daflon, S., et al
Proffitt, C. R., Jin, H., Daflon, S., et al. 2024, Boron Abundances in Early B Dwarfs of the Galactic Open Cluster NGC 3293
2024
-
[24]
R., Lennon , D
Proffitt , C. R., Lennon , D. J., Langer , N., & Brott , I. 2016, , 824, 3
2016
-
[25]
F., Meynet, G., & Maeder, A
Przybilla, N., Firnstein, M., Nieva, M. F., Meynet, G., & Maeder, A. 2010, Astronomy and Astrophysics, 517, A38
2010
-
[26]
Salmon , S. J. A. J., Moyano , F. D., Eggenberger , P., Haemmerl \'e , L., & Buldgen , G. 2022, , 664, L1
2022
-
[27]
Spruit , H. C. 2002, , 381, 923
2002
-
[28]
& Zahn , J
Talon , S. & Zahn , J. P. 1997, , 317, 749
1997
-
[29]
Zahn , J. P. 1992, , 265, 115
1992
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.