{"id":"1cda2531-d1b0-4bad-9180-fadddf965095","arxiv_id":"2411.16534","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Nonlinear MHD simulations of tidally excited waves show that a strong initial magnetic dipole suppresses zonal flows and restores linear tidal dissipation, while a weak field allows flows to reshape the field and nonlinear dissipation.","lead":"Tidal forces on stars and planets can create powerful internal flows, and this paper adds magnetic fields to that picture for the first time in simulations. The researchers find two regimes: strong magnetic fields suppress the flows and return dissipation to old linear predictions, while weak fields let the flows grow and alter how tidal energy is lost.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Regime threshold likely depends on the unstated tidal forcing amplitude; without that parameter, the reported critical Lehnert number is not established as a universal condition.","rationale":"The reader's weakest assumption focuses on the extrapolation from idealized parameters (Ek=1e-5, Pm=1-5, incompressible shell) to real stellar conditions, which is indeed a valid concern. My stress-test identifies a more immediate and specific gap inside the simulations themselves: the tidal forcing amplitude is not reported, and the transition threshold is likely controlled by the ratio of Reynolds to Maxwell stresses, making Le_c a function of that amplitude. The paper's own Introduction states that zonal flows are particularly strong for high tidal amplitudes, so this is not an external critique but an internal control parameter that is missing. The central mechanism—strong fields inhibit zonal flows, weak fields allow them and reshape the field—is plausible and qualitatively supported by the figures, but the quantitative boundary (and hence the stellar prediction) cannot be assessed without the amplitude. This does not overturn the reader's conditional acceptance, as the work is exploratory and a detailed article is forthcoming; it reinforces that the conditional is warranted. The concrete test would directly resolve whether the discrepancy is severe or benign.","tokens_in":5079,"tokens_out":10105,"duration_ms":92300,"concrete_test":"Run a series of simulations at the same Le, Ek, and Pm but vary the tidal forcing amplitude by factors of 3 and 10 about the value used in the paper (requiring the authors to report that value). Determine whether the critical Le_p between the two regimes shifts roughly linearly with amplitude, as expected from the stress balance. If Le_c changes, the stated threshold is amplitude-dependent and the stellar extrapolation must specify the relevant tidal amplitude; if it does not, the ambiguity is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The two-regime claim hinges on a balance between Reynolds stresses (which generate zonal flows) and Maxwell stresses (which suppress them). The Reynolds stress is proportional to the square of the tidal forcing amplitude, while the Maxwell stress is proportional to Le^2, so the critical Lehnert number Le_c should vary with the tidal amplitude. Section 2 defines the effective tidal forcing ft 'in a similar way as in Papers I & II' but never states its amplitude; the paper's own Introduction notes that zonal flows are 'particularly strong for ... high tidal amplitudes' (Section 1). Thus the reported threshold (Le_c ≈ 1e-3 to 4e-3 for Pm = 1–5) is conditional on an unreported parameter. If the simulations used a large tidal amplitude to make the hydrodynamic zonal flow robust, then the inferred Le_c is an upper limit for realistic (smaller) amplitudes, and the paper's prediction that both regimes are relevant for low-mass stars (Section 3) could overstate the magnetic suppression of zonal flows. The missing amplitude also prevents an independent reproduction of the figures.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":5214,"tokens_out":3339,"duration_ms":32323,"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":[{"comment":"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.","section":"Section 2, Eq. (2.1a); Section 3 (Lec values)"},{"comment":"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.","section":"Section 3, Figs. 1 and 2"},{"comment":"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.","section":"Section 3, final paragraph and footnote"},{"comment":"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.","section":"Section 3 and Conclusions (self-declared limitations)"}],"minor_comments":[{"comment":"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.","section":"Section 2"},{"comment":"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.","section":"Fig. 1 caption"},{"comment":"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.","section":"Fig. 2 caption and text"},{"comment":"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.","section":"Section 3, low-Le discussion"},{"comment":"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.","section":"General presentation"}],"recommendation":"major_revision","confidential_remarks":"This is a short conference proceedings paper, so the scope of revision is limited. The main issue is the missing tidal forcing amplitude, which is central to the reproducibility and interpretation of the critical Lehnert number. The authors could address this by reporting the amplitude and showing that the regime boundary does not depend on it (or by parameterizing the threshold accordingly). The lack of convergence tests is concerning but perhaps acceptable for a proceedings if explicitly flagged; adding one resolution test would be a low-cost improvement. The paper's own caveats about forthcoming articles should be respected, but they should not be over-relied upon as a substitute for support within this manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Astoul and Barker report the first nonlinear MHD simulations of tidally excited waves in a rotating spherical shell with an initial dipolar field. The headline result—two regimes in which a magnetic field either kills the tidally driven zonal flow (Le above about 1e-3) or lets it survive and reshape the field (Le below)—is plausible and worth taking seriously. It extends their own hydrodynamic Papers I and II in a natural direction and goes beyond the linear MHD work of Lin & Ogilvie 2018.\n\nWhat they do well: they compare directly with hydrodynamic linear and nonlinear baselines; they track poloidal and toroidal magnetic energy separately, so the omega-effect and the quadrupolar feedback are visible in the diagnostics; and they vary Pm over 1, 2, 5, which gives a first look at how Ohmic diffusivity shifts the transition. For a proceedings paper, the level of physical interpretation is high. The conclusion that both regimes may matter for solar-mass stars is clearly flagged as a prediction, not a measurement.\n\nThe soft spots are real, but mostly inherited from the format. The tidal forcing amplitude ft is never stated; it is 'defined similarly to Papers I & II,' so the reported critical Le values are conditional on an unreported parameter. Since zonal flow generation scales with forcing amplitude while Maxwell stresses scale with Le^2, Le_c likely depends on amplitude, and the reader cannot reproduce the threshold from this text alone. There are no convergence or resolution tests and no error bars, and the extrapolation from Ek=1e-5, Pm=O(1) to stellar Ek~1e-12 with microscopic or turbulent Pm is exactly that—an extrapolation the authors acknowledge. None of this undermines the existence of two regimes; it just means the quantitative boundary is provisional.\n\nWho this is for: people working on tidal dissipation in low-mass stars and giant planets, star-planet orbital evolution, and MHD wave dynamics in convective envelopes. I would cite it, and I would send the full version to a serious referee. As a proceedings, it is a useful signal rather than a complete proof.","headline":"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.","tokens_in":5793,"tokens_out":2457,"would_cite":true,"duration_ms":23797,"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":"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.","keywords":["tidal interactions","magnetohydrodynamics","convective envelopes","zonal flows","Lehnert number","torsional Alfvén waves","tidal dissipation","low-mass stars"],"falsifier":"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.","tokens_in":4811,"feed_emoji":"🧲","tokens_out":14562,"duration_ms":125356,"temperature":0.7,"pith_summary":"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.","feed_headline":"Strong magnetic fields quench tidal zonal flows in simulations","feed_subtitle":"A dipolar field above a threshold restores linear tidal dissipation; weaker fields let nonlinear flows persist.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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}$."],"supporting_citations":[{"why":"Supplies the purely hydrodynamic nonlinear tidal simulations whose zonal-flow generation and dissipation rates this paper extends by adding magnetism.","marker":"Astoul & Barker 2022"},{"why":"Companion hydrodynamic study establishing how tidally generated differential rotation modifies dissipation, providing the nonlinear baseline.","marker":"Astoul & Barker 2023"},{"why":"Earlier global nonlinear tidal simulations that first demonstrated tidally driven zonal flows and their effect on dissipation.","marker":"Favier et al. 2014"},{"why":"Defines the initial dipolar field configuration and provides the linear MHD tidal framework used for comparison.","marker":"Lin & Ogilvie 2018"},{"why":"Supports the Maxwell-stress quenching of differential rotation and gives microscopic and turbulent Prandtl number estimates for stellar convection zones.","marker":"Brun & Browning 2017"},{"why":"Names and describes the $\\Omega$-effect by which differential rotation winds poloidal field into toroidal field.","marker":"Spruit 1999"},{"why":"Documents torsional Alfvén waves in rotating spherical shells, the phenomenon invoked to explain the zonal-flow oscillations.","marker":"Hori et al. 2023"},{"why":"Provides the timescale balance between zonal-flow winding and Alfvén propagation used to explain the regime transition.","marker":"Jouve et al. 2015"},{"why":"Supplies the estimate that turbulent magnetic Prandtl numbers are of order unity, which underlies the prediction that both regimes can occur in low-mass stars.","marker":"Käpylä et al. 2020"}],"fun_headline_variants":["Strong magnetic fields quench tidal zonal flows","Magnetic field strength flips tidal dissipation regime","Weak fields sustain tidal zonal flows, strong fields kill them","Magnetic fields decide whether tidal waves spin up flows","Strong magnetism restores linear tidal dissipation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Strong magnetic fields quench tidal zonal flows","Magnetic field strength flips tidal dissipation regime","Weak fields sustain tidal zonal flows, strong fields kill them","Magnetic fields decide whether tidal waves spin up flows","Strong magnetism restores linear tidal dissipation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000737,"raw_usage":{"total_tokens":3305,"prompt_tokens":967,"completion_tokens":2338,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":2265}},"tokens_in":583,"tokens_out":2338,"duration_ms":18222,"temperature":1.0,"reasoning_tokens":2265,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:59:48.888770+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"& Barker, A","cited_arxiv_id":null,"evidence_quote":"Supplies the purely hydrodynamic nonlinear tidal simulations whose zonal-flow generation and dissipation rates this paper extends by adding magnetism."},{"cited_title":"& Barker, A","cited_arxiv_id":null,"evidence_quote":"Companion hydrodynamic study establishing how tidally generated differential rotation modifies dissipation, providing the nonlinear baseline."},{"cited_title":"J., Baruteau, C., & Ogilvie, G","cited_arxiv_id":null,"evidence_quote":"Earlier global nonlinear tidal simulations that first demonstrated tidally driven zonal flows and their effect on dissipation."},{"cited_title":"& Ogilvie, G","cited_arxiv_id":null,"evidence_quote":"Defines the initial dipolar field configuration and provides the linear MHD tidal framework used for comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents torsional Alfvén waves in rotating spherical shells, the phenomenon invoked to explain the zonal-flow oscillations."},{"cited_title":"2015, A&A, 575, A106 K¨ apyl¨ a, P","cited_arxiv_id":null,"evidence_quote":"Provides the timescale balance between zonal-flow winding and Alfvén propagation used to explain the regime transition."}],"review_version":1}