{"id":"80b7776c-3d3b-4fea-a3d5-ba0ab8705ae8","arxiv_id":"2411.13925","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A semi-analytical prescription for angular momentum transport by nonlinearly breaking internal gravity waves in stellar radiation zones, adapted from the Earth's atmospheric saturation model.","lead":"Internal gravity waves inside stars can crash and redistribute the star's rotation, just as similar waves do in Earth's atmosphere. This paper adapts an atmospheric saturation model to stellar interiors and derives a new formula for the angular momentum deposited by such breaking waves.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (44) applies min to signed fluxes, so it fails to cap angular-momentum transport for prograde modes; a sign-preserving magnitude cap is needed.","rationale":"The strongest claim is that Eq. (44) provides a ready-to-implement parametrization that limits the angular-momentum flux to the minimum of radiative damping and nonlinear saturation. For that claim to hold, the mathematical operation must actually cap the flux. In Eq. (44) the three signed fluxes are all negative for prograde modes (ε = -1, m < 0, F_E < 0), and the pointwise min over negative numbers returns the most negative value, i.e., the largest magnitude. Consequently the cap is never applied in the regime where it is needed. This is an internal inconsistency, independent of the physical approximations: even granting the saturation model of §3.1–3.2, Eq. (44) does not deliver what the text says. It is more load-bearing than the reader's radiative-diffusion concern about Eq. (25), which is a physical approximation that could be tested or debated with simulations; the min operation is a clear mathematical defect. A sign-preserving cap fixes it, and the rest of the derivation (saturated flux formulas, damping factor, modal sum) can stand. The paper is therefore still conditionally acceptable, but the conditions must include correcting Eq. (44) and explicitly stating the sign convention. The reader's identified weakness is real but secondary to this central formula error.","tokens_in":17797,"tokens_out":18791,"duration_ms":189204,"concrete_test":"Choose a solar-type test case (ε = -1), one prograde mode (m < 0), and a radius r where the quasi-adiabatic flux magnitude exceeds the saturation flux magnitude. Evaluate the three arguments inside the min in Eq. (44) with representative stellar values (ρ, N, K, l, m, ω, Ω). Verify that min returns the damped flux (largest magnitude) rather than the saturation value, and that a sign-preserving cap sign(F)·min(|F|, |F_CWB|, |F_SWB|) returns the saturation value. Then recompute the sum over modes: if the net flux changes between the two operations, Eq. (44) does not correctly implement the stated saturation condition.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central prescription Eq. (44) is a pointwise minimum of three signed flux expressions. For each mode, all three terms have the same sign: the factors -m/ω (first term) and -m ε ρ (the two saturation terms) make the flux positive for retrograde waves and negative for prograde waves under the convention of §3.3. The minimum of signed numbers selects the most negative value, i.e., the largest magnitude. Thus whenever the damped quasi-adiabatic flux is more negative than the convective or shear saturation fluxes (exactly the regime where the wave is supposed to break), min returns the damped flux, not the saturation limit. The stated cap in the text ('the flux cannot exceed its value in the saturated regime') is therefore not enforced for prograde modes. The correct operation is a sign-preserving cap, e.g., sign(F_damp) · min(|F_damp|, |F_CWB|, |F_SWB|), or equivalently max in the negative-sign branch. This is not a minor typo: it changes the sign-selected angular-momentum transport and would directly affect the output of any stellar evolution code implementing Eq. (44). The physical saturation model of §3.1–3.2 may be accepted; the mathematical form of the central equation does not implement it.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript adapts a saturation scheme originally developed for internal gravity waves in the Earth's atmosphere to deep spherical stellar radiation zones. Starting from the JWKB description of low-frequency IGWs, it derives saturated vertical and horizontal velocities for convective breaking (Eqs. 27-28) and for shear-driven breaking without and with thermal diffusion (Eqs. 36 and 40), from which it computes energy and angular-momentum fluxes (Eqs. 31-32, 37-38, 41-42). The central prescription, Eq. (44), is a sum over modes of the minimum of the quasi-adiabatic radiatively damped flux and the two nonlinearly saturated fluxes. A rotating, magnetised extension is given in Section 4 with Hough functions replacing spherical harmonics and culminating in Eq. (52). The paper claims to provide the first prescription for angular-momentum transport by nonlinearly breaking IGWs in stellar radiation zones.","tokens_in":18066,"tokens_out":11543,"duration_ms":118385,"significance":"If Eq. (44) can be corrected and the saturation criterion properly justified, the paper would supply a ready-to-implement expression for a missing ingredient in one-dimensional stellar evolution codes. The derivation is internally consistent and traceable: the saturated amplitudes and fluxes follow from stated instability criteria, no parameters are fitted to stellar data (the only free parameter is the critical Richardson number Ri_c), and the geophysical saturation scheme has been independently validated against stratospheric balloon measurements. The algebraic steps are explicit, and the physical discussion of the dependence on l, m, N, rho, K, and the Doppler-shifted frequency is useful. The main weaknesses are that the central formula as written does not implement the saturation statement it claims, and the convective saturation criterion is imported without a quantitative treatment of thermal diffusion in stellar conditions.","major_comments":[{"comment":"The pointwise minimum in Eq. (44) is taken over signed fluxes, but for a fixed mode all three bracketed terms have the same sign: each is proportional to -m*epsilon times a positive factor, where the energy-flux sign epsilon is carried by the first term as well as by the two saturation terms, assuming the usual branch with hatomega > 0. Consequently, for the sign branch in which the flux is negative (e.g. prograde waves, m < 0, in the solar-type case epsilon = -1), min selects the most negative term, i.e. the largest magnitude. In precisely the regime where nonlinear breaking should cap the transport, the formula returns the quasi-adiabatic damped flux rather than the smaller saturated flux. The statement preceding Eq. (44) that the flux cannot exceed its value in the saturated regime is therefore not enforced for prograde modes. This is not a typographical issue: it changes the sign-selected angular-momentum transport and would affect any stellar-evolution implementation. The correct operation is a sign-preserving magnitude cap, for example F = sign(F_damped) * min(|F_damped|, |F_CWB|, |F_SWB|), or equivalently a branchwise max for the negative-sign branch; the same repair is needed in Eq. (52).","section":"Section 3.3, Eq. (44); also Section 4, Eq. (52)"},{"comment":"The saturation criterion for convective breaking uses the adiabatic heat-transport equation Dt T' + Gamma ur = 0 and identifies |dr T'| ~ kr |T'| = Gamma. In stellar radiation zones the Prandtl number is very small but the thermal diffusivity K is large, and heat diffusion is retained, not neglected, in the quasi-adiabatic damping treatment of Section 2.2. It is therefore not self-evident that the temperature perturbation reaches the gradient threshold before being diffusively damped or phase-shifted. The derivation of the saturated velocity |ur|_sat = hatomega / kr, and hence of Eqs. (27)-(32) and the corresponding terms in Eq. (44), needs an explicit timescale comparison (e.g. hatomega versus K kr^2) or a derivation starting from the linearised heat equation including K. If the diffusion-modified saturation differs, the saturated amplitudes and fluxes must be revised. The conclusion's later admission that the local criterion is an approximation does not by itself settle this quantitative point.","section":"Section 3.1.1, Eqs. (24)-(27)"}],"minor_comments":[{"comment":"There are several typographical errors: 'convectivelly-excited' in Section 1, 'adress' in Section 3.1, 'ballons' in Section 3.1.2, 'loose' for 'lose' in Section 3.1.2, and 'beyound' in Section 5.","section":"Throughout"},{"comment":"The typesetting of the damping integrand appears corrupted: 'K N N 2 T' should presumably be 'K N N_T^2' (or the intended combination of N and N_T). Please correct the displayed equations.","section":"Eqs. (14) and (21)"},{"comment":"The sign convention for the Doppler-shifted frequency hatomega is not stated. The manuscript says prograde waves have m < 0 and retrograde waves have m > 0, but it does not state whether hatomega is positive for the propagating waves considered. Since every term in Eq. (44) changes sign with hatomega, the convention should be made explicit.","section":"Section 2.1 and Eq. (44)"},{"comment":"The conclusion that the local saturation criterion is 'an approximation of the reality' once turbulence develops is useful, but it should be stated at the point where the criterion is introduced, near Eq. (25), with a brief indication of the expected validity domain in stellar parameters.","section":"Section 3.1.1 and Section 5"},{"comment":"The angular average in Eq. (51) is written with the notation <Theta^2>_theta, which averages only over colatitude, while the flux is described as a horizontal average over the sphere. Please clarify how the azimuthal dependence of the Hough functions is removed or normalised.","section":"Section 4, Eq. (51)"}],"recommendation":"major_revision","confidential_remarks":"To the editor: I agree with the reader's conditional assessment. The saturation-flux algebra is internally consistent, but Eq. (44) as written fails to implement its own cap because of the signed-min issue, and the convective saturation criterion needs a quantitative justification under stellar thermal-diffusion conditions. Both issues are load-bearing for the claimed prescription and are correctable in revision. If addressed, the paper would be a valuable contribution to the IGW angular-momentum transport literature and is suitable for A&A."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the one-line: Mathis has done a real service by adapting the atmospheric saturation scheme of Lott and co-workers to spherical stellar radiation zones, and the resulting closed-form fluxes are exactly what 1D stellar evolution codes have been missing. The derivations are careful, the JWKB treatment is standard, and the geophysical validation of the parent model is a point in its favor. Worth reading closely.\n\nWhat's new: the saturated vertical velocity from convective overturning (Eq. 27) and from vertical shear, with and without heat diffusion (Eqs. 36, 40), and the corresponding energy and angular momentum fluxes (Eqs. 31-32, 37-38, 41-42). The extension to rotating/magnetized stars via TARM in Sec. 4 is a bonus. I checked the algebra on the main steps and it is internally consistent. The paper also reads honestly, with limitations stated.\n\nThe soft spot that matters: Eq. (44) takes the pointwise minimum of three signed fluxes. Since all three terms carry the same sign for a given mode, min for prograde (negative) fluxes returns the most negative, i.e., the largest magnitude. So in exactly the regime where the wave breaks, the damping flux wins and the saturation cap is not applied. The text says the flux cannot exceed the saturated value, but the formula does not implement that for prograde waves. The fix is straightforward: apply min to absolute values and restore the sign, e.g. sign(F_damp) min(|F_damp|, |F_CWB|, |F_SWB|). This is not cosmetic; it changes the sign-selected transport in any code that implements Eq. (44).\n\nA second concern, less central: the convective saturation criterion in Eq. (25) is imported from the atmosphere without discussing radiative diffusion in stellar interiors (Prandtl number ~1e-6). If diffusion erodes the temperature perturbation before overturning, the saturated amplitudes and fluxes in Sec. 3.1 could shift. The shear case with diffusion is treated, but the convective case is not.\n\nAlso, the paper would be stronger with a comparison to existing numerical simulations (e.g., Barker & Ogilvie, Rogers) or asteroseismic constraints. The 'none has been proposed' claim is slightly overstrong, but not by much.\n\nWho should read it: anyone building or using angular-momentum transport prescriptions in stellar evolution, and geophysicists interested in spherical extensions of wave-breaking parameterizations. It deserves a serious referee, with the expectation of a revision to fix the signed-min issue and to address the diffusion caveat.","headline":"A genuinely useful, mostly sound prescription for angular-momentum transport by breaking internal gravity waves, but the central cap equation uses min on signed fluxes and so fails for prograde modes.","tokens_in":18579,"tokens_out":5000,"would_cite":true,"duration_ms":48650,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives the first prescription for angular momentum transport by breaking internal gravity waves in stellar radiation zones, expressed as a minimum of three fluxes.","keywords":["internal gravity waves","stars: rotation","stars: evolution","methods: analytical","angular momentum transport","wave breaking","saturation scheme","asteroseismology"],"falsifier":"Run a direct numerical simulation of internal gravity waves in a spherical shell with stellar-like stratification and Prandtl number $10^{-6}$, measure the rms vertical velocity and angular momentum flux at the onset of breaking, and compare with $|\\hat{u}_r|_{\\mathrm{sat}}=\\hat{\\omega}/k_r$ and with Eq. (44): if breaking starts at a noticeably different amplitude or deposits momentum at a different radius, the saturation cap is wrong.","tokens_in":17609,"feed_emoji":"🌊","tokens_out":14496,"duration_ms":98833,"temperature":0.7,"pith_summary":"This paper supplies the piece of physics that stellar evolution codes have been missing: what happens when internal gravity waves in stably stratified stellar radiation zones break because they overturn the stratification or because their own vertical shear goes unstable. The central claim is that each wave mode's angular momentum flux is the minimum of three values: the usual radiatively damped flux, the flux at convective-breaking saturation, and the flux at shear-breaking saturation. The prescription adapts a saturation model developed for the Earth's atmosphere, and the paper argues it brings stellar wave–mean-flow interaction to the same level of realism as atmospheric global circulation models. If correct, it gives stellar evolution codes a ready-to-implement term that could help close the gap between predicted and asteroseismically observed angular momentum extraction, which is about two orders of magnitude when waves and magnetic fields are neglected. The paper also extends the formula to rotating, magnetised stars, where the Coriolis and Lorentz forces reduce the convective-breaking flux.","feed_headline":"New formula caps spin lost to breaking waves in stars","feed_subtitle":"Asteroseismology sees stars losing spin far faster than models predict; this adds the missing wave-breaking term.","key_machinery":"The carrying machinery is a saturation scheme adapted from atmospheric physics. For convective breaking, the wave amplitude is capped where the radial gradient of the temperature perturbation, $|\\partial_r T'|\\simeq k_r|T'|$, equals the background gradient scale $\\Gamma=T N_T^2/(g\\delta)$; this yields the saturated vertical velocity $\\hat{\\omega}/k_r$ and, through the wave polarisation relation, the saturated energy flux (Eq. 31) and angular momentum flux (Eq. 32). For shear breaking, the cap is set by the Richardson criterion $N^2/|d\\langle \\mathbf{u}_h\\cdot\\mathbf{u}_h^*\\rangle/dr|^2\\le Ri_c\\,\\beta_{l,m}$ with $\\beta_{l,m}=m^2/[l(l+1)]$, giving saturated fluxes that, when heat diffusion is included, grow with the heat diffusivity $K$. The final formula puts the three transport channels into one minimum, so the limiting mechanism sets the deposit. In the rotating, magnetised generalisation, spherical harmonics are replaced by Hough functions and the dispersion relation becomes $k_r=(N/\\omega_M)k_h$ with $\\omega_M^2=\\omega_s^2-m^2\\omega_A^2$, which shifts the same machinery to magneto-gravito-inertial waves.","core_discovery":"The paper's central result is Eq. (44): the total angular momentum flux transported by internal gravity waves at radius $r$ is the sum over modes of the minimum of three fluxes, namely the standard quasi-adiabatic flux $-(m/\\hat{\\omega})\\,F_{E;l,m}(r_0)\\,e^{-\\tau_{l,m}(r)}$, the convective-breaking flux $F^{\\mathrm{CWB}}_{J;l,m}$, and the shear-breaking flux $F^{\\mathrm{SWB}}_{J;l,m}$. The convective-breaking flux follows from the saturation rule that the wave's temperature gradient reaches the background gradient, $|\\partial_r T'|\\simeq k_r|T'|=\\Gamma$, which fixes the saturated radial velocity to $|\\hat{u}_{r;l,m}|_{\\mathrm{sat}}=\\hat{\\omega}/k_r=\\hat{\\omega}^2/(N k_h)$. The shear-breaking flux follows from the Richardson criterion for the wave's own vertical shear, with the spherical-geometry factor $\\beta_{l,m}=m^2/[l(l+1)]$, and with heat diffusion included through a Péclet number. The result is that breaking deposits angular momentum preferentially where the Doppler-shifted frequency $\\hat{\\omega}$ tends to zero, at high latitudinal degree $l$, at large $N/\\hat{\\omega}$, and where density is low. The paper states this is the first prescription for angular momentum transport triggered by nonlinear breaking of internal gravity waves in stellar radiation zones.","pith_inferences":["The same min-of-three structure could be turned into a diagnostic: in a stellar evolution model, the radius where the saturated flux first caps the radiatively damped flux should coincide with the shallowest gradient in an asteroseismically inferred rotation profile, so disagreement would localise where the saturation model needs revision.","Because the convective-breaking flux scales with $\\rho$ while the heat-diffusive shear-breaking flux has a strong dependence on $K$, the parameterisation could plausibly be extended to chemical transport by replacing the angular momentum flux with a tracer flux along the same saturation curves.","The dependence of the rotating, magnetised version on $\\omega_M$ predicts an angular momentum 'dead zone' at the magnetic critical layer; comparing rotation profiles of stars with known magnetic fields against that prediction would provide a test independent of the hydrodynamic case."],"forward_implications":["The prescription implies that internal gravity waves in stars can break before reaching critical layers, so the saturation cap, not radiative damping alone, sets where angular momentum is deposited.","The breaking flux grows with latitudinal degree $l$, with $N/\\hat{\\omega}$, and as density drops, and it is deposited most efficiently where the Doppler-shifted frequency $\\hat{\\omega}$ vanishes; retrograde waves extract angular momentum from the radiative cores of solar-type stars while prograde waves deposit it, with the opposite signs in early-type stars.","The formula is analytically ready to be implemented in one-dimensional stellar structure and evolution codes, raising wave-driven angular momentum transport in stars to the same level of description as atmospheric global circulation models.","In rotating, magnetised stars the Coriolis and Lorentz forces reduce the convective-breaking flux, and the transport vanishes at the magnetic critical layer where $\\omega_M=0$.","Near the low-density surfaces of early-type stars, the enhanced breaking-driven momentum deposit may participate in driving matter ejection, as in active Be stars."],"supporting_citations":[{"why":"Supplies the saturation scheme for convectively excited breaking internal waves whose saturated velocity and flux expressions are adapted here to spherical stellar geometry.","marker":"Lott et al. (2012)"},{"why":"Provides the wave-mean-flow parameterisation for convectively excited breaking waves that the stellar version generalises to deep spherical shells.","marker":"Lott & Guez (2013)"},{"why":"Reports in-situ stratospheric balloon measurements of wave-induced Reynolds stresses, used as observational validation of the atmospheric scheme being adapted.","marker":"Lott et al. (2023)"},{"why":"Introduces the convective saturation criterion that the wave's temperature gradient reaches the background temperature gradient.","marker":"Lindzen (1981)"},{"why":"Establishes convective instability as a nonlinear breaking mechanism for internal gravity waves whose threshold the saturation condition encodes.","marker":"Sutherland (2001)"},{"why":"Identifies internal gravity wave breaking in stars and the criterion that the wave velocity overcomes the phase velocity, which the saturated amplitudes satisfy.","marker":"Press (1981)"},{"why":"Provides the standard quasi-adiabatic radiative-damping flux and the angular momentum transport equation that the final minimum formula builds on.","marker":"Zahn et al. (1997)"},{"why":"Derives the Richardson criterion for vertical shear instability of internal gravity waves in spherical geometry, including the spherical-geometry factor depending on the harmonic orders.","marker":"Alvan et al. (2013)"},{"why":"Gives the rotating, magnetised formalism with Hough functions, the frequency omega_M, and the thermal damping rate used for the Coriolis and Lorentz generalisation.","marker":"Mathis & de Brye (2012)"}],"fun_headline_variants":["Breaking waves in stars: new spin transport formula","Star spin loss from wave breaking now predictable","Formula for stellar wave breaking fills spin gap","Geophysical wave model adapted to stellar spin","First wave-breaking spin transport prescription for stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the saturation condition imported from atmospheric physics: a wave breaks when the radial gradient of its temperature perturbation, $k_r|T'|$, reaches the background temperature gradient $\\Gamma$, fixing the saturated velocity to $\\hat{\\omega}/k_r$; this assumes the temperature perturbation survives to that point even though stellar radiation zones have extremely strong heat diffusion (Prandtl number $\\sim10^{-6}$) that could soften it before breaking.","fun_headline_variants_meta":{"raw":{"variants":["Breaking waves in stars: new spin transport formula","Star spin loss from wave breaking now predictable","Formula for stellar wave breaking fills spin gap","Geophysical wave model adapted to stellar spin","First wave-breaking spin transport prescription for stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000345,"raw_usage":{"total_tokens":2011,"prompt_tokens":1183,"completion_tokens":828,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":799,"completion_tokens_details":{"reasoning_tokens":761}},"tokens_in":799,"tokens_out":828,"duration_ms":72008,"temperature":1.0,"reasoning_tokens":761,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:44:22.854998+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a direct numerical simulation of internal gravity waves in a spherical shell with stellar-like stratification and Prandtl number $10^{-6}$, measure the rms vertical velocity and angular momentum flux at the onset of breaking, and compare with $|\\hat{u}_r|_{\\mathrm{sat}}=\\hat{\\omega}/k_r$ and with Eq. (44): if breaking starts at a noticeably different amplitude or deposits momentum at a different radius, the saturation cap is wrong.","supporting_citations":[{"cited_title":"2012, Geophys","cited_arxiv_id":null,"evidence_quote":"Supplies the saturation scheme for convectively excited breaking internal waves whose saturated velocity and flux expressions are adapted here to spherical stellar geometry."},{"cited_title":"& Guez, L","cited_arxiv_id":null,"evidence_quote":"Provides the wave-mean-flow parameterisation for convectively excited breaking waves that the stellar version generalises to deep spherical shells."},{"cited_title":"2023, Journal of Geophysical Research (Atmospheres), 128, e2022JD037585","cited_arxiv_id":null,"evidence_quote":"Reports in-situ stratospheric balloon measurements of wave-induced Reynolds stresses, used as observational validation of the atmospheric scheme being adapted."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the convective saturation criterion that the wave's temperature gradient reaches the background temperature gradient."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes convective instability as a nonlinear breaking mechanism for internal gravity waves whose threshold the saturation condition encodes."},{"cited_title":"P., Talon, S., & Matias, J","cited_arxiv_id":null,"evidence_quote":"Provides the standard quasi-adiabatic radiative-damping flux and the angular momentum transport equation that the final minimum formula builds on."}],"review_version":1}