Pith. sign in

REVIEW 4 major objections 5 minor 13 references

Interactions between tidal flows and magnetic fields in stellar/planetary convective envelopes

T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Magnetism switches a rotating convective shell between a linear and a nonlinear tidal-dissipation regime: strong dipolar fields suppress zonal flows, weak fields let them persist and reshape the field.

desk verdict First nonlinear MHD simulations of tidally forced waves with a dipole field show two plausible regimes, but the quantitative threshold is provisional until the forcing amplitude and a full paper appear. read the letter →

arxiv 2411.16534 v1 pith:5BYVJPP3 submitted 2024-11-25 astro-ph.SR astro-ph.EP

Aurélie Astoul , Adrian J. Barker This is my paper · ORCID
classification astro-ph.SRastro-ph.EP
keywords tidalinteractionsmagnetohydrodynamicsconvectiveenvelopeszonalflowsLehnertnumbertorsionalAlfvénwavesdissipationlow-massstars
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

The paper tries to establish that magnetism changes how efficiently a rotating star or planet dissipates tidal energy. In nonlinear simulations of tidally forced waves in a spherical convective shell, a weak initial dipolar magnetic field leaves the tidally generated zonal flows intact, so dissipation matches the nonlinear hydrodynamic rate and the flows wind the field into a toroidal component. A strong enough dipole (Lehnert number above about $10^{-3}$ at magnetic Prandtl number one) destroys the zonal flow through Maxwell stresses and torsional Alfvén waves, restoring dissipation to the linear hydrodynamic prediction. The paper argues that both regimes are plausible in low-mass stars because their turbulent magnetic Prandtl number is near unity.

What carries the argument

The controlling parameter is the Lehnert number $Le$, the ratio of the Alfvén speed to the rotation speed, which sets the strength of the initial dipolar field. The controlling mechanism is the competition between Reynolds stresses, which build the axisymmetric zonal flow from nonlinear tidal-wave self-interactions, and Maxwell stresses, which act to destroy it. Torsional Alfvén waves carry this stress balance at intermediate field strengths and produce oscillations in the zonal flow and in the dissipation. In the weak-field branch the $\Omega$-effect, the stretching of poloidal field by the zonal shear, creates toroidal field, and the quadrupolar tidal flow converts part of it back to poloidal field. The simulations use Ekman number $Ek=10^{-5}$, magnetic Prandtl numbers $Pm=1,2,5$, an initial dipolar field, stress-free velocity boundaries, and insulating magnetic boundaries.

What would settle it

Repeat the simulations with much weaker viscosity (Ekman number from $10^{-5}$ down toward $10^{-7}$) at $Pm=1$; if zonal flows still appear at Lehnert number $10^{-3}$ rather than being suppressed, the claimed strong-field regime does not survive at stellar conditions.

Watch

Extended reading notes

Core claim

The paper's central claim is that the tidal response of a magnetised convective envelope splits into two regimes controlled by the strength of the large-scale poloidal field, measured by the Lehnert number $Le = B_0/(\sqrt{\mu\rho}R\Omega)$. For $Le \gtrsim 10^{-3}$ (at $Pm=1$), Maxwell stresses from the initial dipole suppress the zonal flow that otherwise builds up through nonlinear wave self-interactions; torsional Alfvén waves are excited, and the viscous dissipation settles near the linear hydrodynamic value. For $Le \lesssim 10^{-3}$, the zonal flow reaches the same amplitude as in the purely hydrodynamic case, and the dissipation matches the nonlinear hydrodynamic value. In this weak-field regime the zonal flow stretches the dipole into a toroidal field via the $\Omega$-effect, the quadrupolar tidal wave then restores a poloidal component, and for $Le \lesssim 6\times10^{-5}$ the Lorentz force becomes negligible so the magnetic energy scales as $Le^2$. The paper further claims that the transition shifts to smaller $Le$ as $Pm$ increases, and that with $Pm\sim 1$ neither regime is excluded for low-mass stars.

Load-bearing premise

The argument depends on the assumption that a simplified computer model of a convective shell, with an idealized initial magnetic field and much stronger friction than real stars have, behaves the same way as a real stellar or planetary convective envelope.

Editorial extensions

If this is right

  • Tidal dissipation in convective envelopes is magnetically controlled: at Lehnert numbers above about $10^{-3}$ and $Pm=1$, the nonlinear enhancement from zonal flows is absent, so dissipation returns to the linear hydrodynamic rate.
  • In the weak-field regime, the zonal flow winds the initial dipole into a toroidal field and the quadrupolar tidal flow then restores a poloidal component, so tides can restructure the large-scale magnetic field.
  • The transition between regimes is accompanied by torsional Alfvén waves that make the zonal flow and dissipation oscillate, and the critical Lehnert number shifts to lower values as the magnetic Prandtl number increases.
  • For turbulent magnetic Prandtl numbers near unity, both the zonal-flow-dominated and the zonal-flow-suppressed regimes can occur in low-mass stars for Lehnert numbers in the range $10^{-4}$ to $10^{-2}$.

Reading between the lines

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

  • If this regime switch holds at stellar parameters, stars with strong large-scale fields should circularize close binaries and hot-Jupiter orbits more slowly than pure hydrodynamics predicts, because the zonal-flow enhancement of tidal dissipation would be switched off.
  • A clean numerical test would initialize the same shell with a purely toroidal field at the same Lehnert number: the stress-balance mechanism in the paper predicts far weaker zonal-flow suppression than for a dipole, since the Lorentz force opposing the axisymmetric flow would be much smaller.
  • The weak-field branch's kinematic scaling, with magnetic energy proportional to $Le^2$ for $Le \lesssim 6\times10^{-5}$, suggests a diagnostic: time series of the toroidal-to-poloidal magnetic energy ratio should cleanly separate the two regimes in future simulations.
  • Because the critical Lehnert number decreases as $Pm$ increases, envelopes with turbulent magnetic Prandtl numbers near unity may sit close to the threshold, so modest changes in field strength could flip a star between linear and nonlinear tidal dissipation behavior.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. This SF2A conference proceeding presents the first nonlinear MHD simulations of tidally excited waves in an incompressible, rotating spherical shell with an initial dipolar magnetic field. The authors vary the Lehnert number Le and the magnetic Prandtl number Pm and report two regimes: for Le above roughly 10^-3, strong magnetic fields suppress tidally generated zonal flows, restoring tidal dissipation to linear hydrodynamic levels, while for weaker fields the zonal flow survives, drives nonlinear dissipation, and reshapes the magnetic field via the Omega-effect and the generation of toroidal fields. The paper concludes that both regimes may be relevant for low-mass stars when a turbulent Pm of order unity is used. The central two-regime claim is supported by the presented simulations, which are compared against hydrodynamic linear and nonlinear baselines, but several key numerical and physical parameters are not reported, limiting the reproducibility and robustness of the quantitative threshold.

Significance. If the central claim holds, this is a valuable step toward understanding how magnetic fields modify tidal dissipation in convective envelopes of stars and giant planets, with direct implications for orbital evolution models. The paper is explicitly exploratory and presents a clean qualitative picture: magnetic stresses can inhibit zonal flows and thereby change the dissipation rate from nonlinear back toward linear values. The use of direct numerical simulations, rather than a fitted model, is a strength, and the explicit comparison with hydrodynamic baselines (linear and nonlinear) helps anchor the interpretation. However, the significance is tempered by the absence of the tidal forcing amplitude, the lack of any convergence or resolution tests, and the speculative extrapolation from Ek=1e-5 to stellar Ek~1e-12. These issues do not necessarily invalidate the qualitative two-regime picture, but they prevent the reported critical Lehnert numbers from being used as quantitative predictions.

major comments (4)
  1. [Section 2, Eq. (2.1a); Section 3 (Lec values)] The effective tidal forcing ft in Eq. (2.1a) is never specified, in particular its amplitude and spatial structure. The paper's central quantitative result is the critical Lehnert number Lec (reported as ~4e-3, ~2e-3, ~1e-3 for Pm=1,2,5), but this threshold depends on a balance between Reynolds stresses (proportional to the square of the tidal forcing amplitude) that drive zonal flows and Maxwell stresses (proportional to Le^2) that suppress them. Without stating the forcing amplitude used in the simulations, the reported Lec is not reproducible, and the extrapolation to stellar parameters in the final paragraph of Section 3 is not justified. The authors must either report the forcing amplitude and show that the regime boundary is independent of it, or present the threshold as a function of both Le and forcing amplitude.
  2. [Section 3, Figs. 1 and 2] No error bars, uncertainty estimates, or convergence tests are provided for any of the reported quantities (Edr, Dnu, Lec). The two-regime claim rests on a small number of simulations (8 in Fig. 1, and only two Pm values shown in Fig. 2), and the absence of any resolution check means the reader cannot assess whether the sharp transition at Lec is physical or a numerical artifact. A single resolution test for representative weak- and strong-field simulations would substantially strengthen the claim.
  3. [Section 3, final paragraph and footnote] The prediction that both regimes are relevant for low-mass stars relies on extrapolating from Ek=1e-5 and Pm=1-5 to stellar Ek~1e-12 and Pm~1e-2 (microscopic) or ~1 (turbulent). The paper states that 'Lec to be much higher' for microscopic Pm but does not provide a scaling law or a theoretical argument for how Lec depends on Ek and Pm. This extrapolation is therefore speculative; either a scaling estimate should be given, or the claim should be softened to note that the relevance of the two regimes at stellar parameters is not yet established.
  4. [Section 3 and Conclusions (self-declared limitations)] The manuscript itself acknowledges that key supporting details are not included: the footnote on the torsional Alfven wave frequency states that 'more details will be given in a forthcoming article,' and the Conclusions state that the magneto-rotational instabilities observed at higher Pm are 'beyond the scope of this proceeding' and are the subject of an article in preparation. These statements indicate that the physical interpretation of the oscillations and instabilities is not fully substantiated within this paper, which limits the completeness of the two-regime explanation. The authors should at least outline the expected scaling of the torsional wave frequency with Le and Pm, or explicitly mark those aspects as preliminary.
minor comments (5)
  1. [Section 2] The 'effective tidal forcing ft' is not defined in this paper; it is only described as 'similar to Papers I & II.' Since this is a standalone proceedings paper, the amplitude of ft (or a reference to the exact definition with the value used) should be stated in the text or figure caption.
  2. [Fig. 1 caption] The bullets indicating the value of Edr at t≈10^4 are not explained in the caption; it would help to state explicitly that they denote the final-time values.
  3. [Fig. 2 caption and text] The 'grey zone' marking the transition between regimes is not quantitatively defined; please specify the range of Lep it covers, and define what 'powerful TOs' means in terms of an amplitude or frequency criterion.
  4. [Section 3, low-Le discussion] The term 'kinematic regime' (for Le ≲ 6·10^-5) is used but not defined; please clarify whether this means the Lorentz force is negligible in the momentum equation, as stated later in the same paragraph.
  5. [General presentation] The paper does not report the numerical resolution (e.g., spherical harmonic truncation or grid spacing) or the number of grid points used in the simulations. While proceedings often omit such details, including them would aid reproducibility and is particularly important given the exploratory nature of the study.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity; the MHD regime boundary is an emergent simulation outcome, with only minor contextual self-citations.

full rationale

The paper's central claim—that strong initial dipolar magnetic fields inhibit tidally generated zonal flows and restore linear-level viscous dissipation, while weak fields allow zonal flows to survive—is an output of direct numerical simulations, not a derivative of a fitted target. Equations (2.1) define the MHD model, and the Lehnert number is an input parameter, but the critical Lehnert number separating the two regimes is diagnosed from the simulations and is not imposed by the equations. The comparison dissipation rates are separate hydrodynamic simulation results, so the statement that strong-field cases match the linear rate and weak-field cases match the nonlinear rate is a measured outcome, not an identity. The only self-citations with any load are contextual: Papers I and II provide the hydrodynamic baseline and diagnostics, and Astoul et al. 2019 provides an estimated stellar Lehnert-number range for the final extrapolation. These do not substitute for the MHD calculation or assert its conclusion, and they are independently published results used as context. The unstated tidal-forcing amplitude is a reproducibility or robustness caveat, not a circular step, because the critical Lehnert number is not defined in terms of that amplitude. Overall, no significant circularity is present; the score reflects only minor contextual self-citations that are not load-bearing.

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

No fitted parameters are used to force a conclusion; the parameters listed are control parameters that determine the regime behavior. The paper relies on standard MHD equations and a public code, but the key domain assumptions are the simplified shell model and the extrapolation from Ek=1e-5, Pm O(1) to real stellar parameters. No new physical entities are invented.

free parameters (3)
  • Initial Lehnert number Le = 1e-5 to 1e-3 across runs
    Controls the initial magnetic field strength; chosen to sample both regimes, not fitted to match the conclusion, but the central claim depends on this grid.
  • Ekman number Ek = 1e-5
    Viscosity parameter chosen based on mixing-length theory (Ogilvie & Lin 2007); much larger than realistic stellar values, so the extrapolation to stars is a key assumption.
  • Magnetic Prandtl number Pm = 1, 2, 5
    Chosen to study the dependence of the regime transition on Ohmic diffusion; turbulent Pm is expected to be O(1), while microscopic Pm is much smaller, so the choice directly affects the stellar relevance claim.
assumptions (5)
  • domain assumption The incompressible, adiabatic, constant-density spherical shell with stress-free impenetrable velocity boundaries and current-free magnetic boundaries adequately represents a stellar or planetary convective envelope for tidal-wave studies.
    Model setup in Section 2; the paper acknowledges the wavelike velocity boundary is not the free-surface condition but argues it is adequate in a footnote.
  • domain assumption The effective tidal forcing ft, defined as in Astoul & Barker (2022, 2023), correctly excites the same tidal response as a real companion.
    Section 2 states ft is defined similarly to the authors' prior hydrodynamic papers; no independent derivation is given here.
  • domain assumption The initial dipolar magnetic field, defined exactly as in Lin & Ogilvie 2018, is a relevant proxy for the large-scale fields in low-mass stars.
    Section 2; the field decays after t=0, so the results depend on the chosen initial topology and amplitude.
  • domain assumption A constant effective turbulent viscosity and magnetic diffusivity with Ek=1e-5 and Pm around unity capture the unresolved convective transport in the simulations.
    Section 2 and Section 3; the authors note this is motivated by mixing-length theory but acknowledge that microscopic values differ by many orders of magnitude.
  • domain assumption MagIC, modified for tidal interactions, correctly solves the MHD system for the runs reported.
    Section 2; no resolution or convergence tests are reported in this proceedings paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Interactions between tidal flows and magnetic fields in stellar/planetary convective envelopes." pith.science (2026). https://pith.science/paper/5BYVJPP3

@misc{pith2026241116534,
  author       = {Pith},
  title        = {Pith review of: Interactions between tidal flows and magnetic fields in stellar/planetary convective envelopes},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5BYVJPP3}},
  note         = {Machine review of arXiv:2411.16534}
}
abstract

Stars and gaseous planets are magnetised objects but the influence of magnetic fields on their tidal responses and dissipation rates has not been well explored. We present the first exploratory nonlinear magnetohydrodynamic (MHD) simulations of tidally-excited waves in incompressible convective envelopes harbouring an initial dipolar magnetic field. Simulations with weak magnetic fields exhibit tidally-generated differential rotation in the form of zonal flows (like in the purely hydrodynamic case) that can modify tidal dissipation rates from prior linear predictions. Moreover, tidal waves and zonal flows affect the amplitude and structure of the magnetic field, notably through creation of toroidal fields via the $\Omega$-effect. In contrast, simulations with strong magnetic fields feature severely inhibited zonal flows, due to large-scale magnetic stresses, excitation of torsional waves, or magnetic instabilities. We predict that the different regimes observed for weak and strong magnetic fields may be both relevant for low-mass stars when using turbulent values of the magnetic Prandtl number.

Figures

Figures reproduced from arXiv: 2411.16534 by the authors.

Figure 1
Figure 1. Left: Energy in the differential rotation Edr versus the evolving poloidal Lehnert number Lep for 8 simulations having different initial Le (in different colours). The value of Edr at t ≈ 104 is indicated by a bullet and Edr for Le = 0 is shown by a horizontal blue line. Right: Magnetic energy M (solid lines) along with poloidal (dashed) and toroidal (dotted) components, all rescaled by Le2 , versus time Ωt, for sim… view at source ↗
Figure 2
Figure 2. Tidal viscous dissipation Dν versus evolving poloidal Lehnert number Lep. Hydrodynamical (Le = 0) linear and nonlinear tidal dissipation rates are shown in horizontal green and blue lines. Left: Pm = 1. Right: Pm = 2. 1. For high Lehnert numbers Lep > 10−3 , the tidally-generated zonal flow is destroyed early in the simula￾tions, as we can see in [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 9 canonical work pages

  1. [1]

    & Barker, A

    Astoul, A. & Barker, A. J. 2022, MNRAS, 516, 2913

    work page 2022

  2. [2]

    & Barker, A

    Astoul, A. & Barker, A. J. 2023, ApJ, 955, L23

    work page 2023

  3. [3]

    2019, A&A, 631, A111

    Astoul, A., Mathis, S., Baruteau, C., et al. 2019, A&A, 631, A111

    work page 2019

  4. [4]

    & Rieutord, M

    Baruteau, C. & Rieutord, M. 2013, Journal of Fluid Mechanics, 719, 47

    work page 2013

  5. [5]

    Brun, A. S. & Browning, M. K. 2017, Living Reviews in Solar Physics, 14, 4

    2017

  6. [6]

    J., Baruteau, C., & Ogilvie, G

    Favier, B., Barker, A. J., Baruteau, C., & Ogilvie, G. I. 2014, MNRAS, 439, 845

    work page 2014

  7. [7]

    Hori, K., Nilsson, A., & Tobias, S. M. 2023, Reviews of Modern Plasma Physics, 7, 5

    work page 2023

  8. [8]

    2015, A&A, 575, A106 K¨ apyl¨ a, P

    Jouve, L., Gastine, T., & Ligni` eres, F. 2015, A&A, 575, A106 K¨ apyl¨ a, P. J., Rheinhardt, M., Brandenburg, A., & K¨ apyl¨ a, M. J. 2020, A&A, 636, A93

    work page 2015

Show all 13 references
  1. [9]

    & Ogilvie, G

    Lin, Y. & Ogilvie, G. I. 2018, MNRAS, 474, 1644

    2018

  2. [10]

    Ogilvie, G. I. 2014, ARA&A, 52, 171

    2014

  3. [11]

    Ogilvie, G. I. & Lin, D. N. C. 2007, ApJ, 661, 1180

    2007

  4. [12]

    Spruit, H. C. 1999, A&A, 349, 189

    1999

  5. [13]

    2007, Phys

    Tilgner, A. 2007, Phys. Rev. Lett., 99, 194501 §Ek ∼ 10−12 for solar-like stars or ∼ 10−18 for giant gaseous planets

    2007

Pith tools

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