{"id":"3623cc1d-324c-4b53-8d0b-3d29061c7eba","arxiv_id":"2009.10473","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using Paper I's rotating convection scalings, the paper derives gravito-inertial wave flux ratios for two excitation channels and identifies a critical convective Rossby number for enhanced sub-inertial waves.","lead":"This paper estimates how much gravito-inertial wave energy rotating convection can generate in stars and planets. It finds that sub-inertial waves can be strongly enhanced below a critical rotation rate set by the stiffness of the convective-radiative boundary.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The monomodal delta-function source makes the predicted flux peaks at Roc≈1/√5 and the order-of-magnitude sub-inertial enhancements fragile; the paper acknowledges the idealization but does not quantify sensitivity to a broad spectrum or to damping.","rationale":"I read the paper as a self-consistent analytic extension of Paper I's heat-flux-maximized convection model to GIW excitation, not as a claim of immediate quantitative agreement with any particular star. Its nonrotating limits reproduce the Press and Lecoanet-Quataert scalings (Eq. 41, Appendix A), and the formal algebra of Sections 4-5 is laid out in enough detail to be checked. The stress-test question is whether the headline mechanism, sub-inertial wave flux enhanced by orders of magnitude at low-to-moderate convective Rossby number, is robust. The weakest link is the source model and the undamped cavity response: Section 4 states that the source is a delta function in Fourier space, Section 5 uses a single-mode velocity field, and the discrete resonances in Figure 6 are exactly the kind of prediction that can be an artifact of that idealization. Because the authors explicitly defer damping to future work, the flux they compute is an energy-injection rate rather than an observable wave flux; that is an acceptable first step, but it means the quantitative peaks should not be treated as measured amplitudes. A single source-broadening test would show whether the Roc≈1/√5 feature survives. I do not see an internal inconsistency that would require rejection; the reader's CONDITIONAL verdict is the right disposition, and my concern is the same one the reader identified, so no change in verdict is needed.","tokens_in":27390,"tokens_out":22489,"duration_ms":196786,"concrete_test":"Recompute the Reynolds-stress flux of Eq. (65) for the same parameter choices as Figure 6, replacing the monochromatic velocity field (Eqs. 69-70) with a spectrum of modes whose frequency and horizontal wavenumber are broadened around (ω_c, k_⊥) by, say, a Lorentzian with half-width Δω/ω_c=0.3 and a corresponding width in k, while keeping the Paper I rms velocity v(Roc) fixed. If the peak flux at Roc≈1/√5 is smeared by more than an order of magnitude or moves to a different Roc, the delta-function source is load-bearing. A useful companion check is to repeat the same calculation with a small damping rate γ/ω≈0.1 added in the amplitude equation (50), which would also show whether the undamped cavity resonances in Figure 6 survive.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims are the interfacial peak below Roc,crit (Eq. 43) and the Reynolds-stress enhancements of Figure 6. Both inherit the same single-scale model: Section 4 explicitly states that the interfacial source permits only a single horizontal spatial scale and a single time scale, and Section 5 uses the monochromatic velocity field of Eqs. (69)-(70), whose products force only the discrete Fourier component ω=2ω_c, k_⊥=k_n/2. The wave response itself is computed as undamped normal modes of a finite box with Dirichlet boundaries at z=0 and z=L (Eqs. 60-63); no damping appears in the amplitude solution (Eq. 50), so the response at resonance is formally divergent, and the plotted peaks are selected cavity resonances. Real rotating convection has a broad k-ω spectrum, finite correlation times, and damping; each of these will broaden and reduce resonant enhancements. Because the paper does not test the sensitivity of Figure 6 to replacing the delta-function source with a measured or synthetic spectrum, the 'orders of magnitude' sub-inertial enhancement is not yet established as a robust prediction. The limitations are listed in Section 6, but no quantitative bound is given.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript (Paper II of a series) uses the heat-flux-maximized monomodal convection model of Paper I to estimate gravito-inertial wave (GIW) fluxes driven by rotating convection in stellar and planetary interiors. Section 4 generalizes the Press (1981) interfacial pressure-excitation calculation to GIWs, yielding Eq. (43) for the flux ratio Fz/F0 as a function of the convective Rossby number, interface stiffness S, normalized frequency σ, and geometric angles; a peak near the lower cutoff for Roc below Roc,crit = sqrt(6)/(5πS) is reported. Section 5 instead solves the forced Poincaré equation with Reynolds-stress sources, using the amplitude formula Eq. (50) and exact box eigenfunctions for a discontinuous Brunt-Väisälä profile. This leads to flux expressions Eqs. (65)-(72) and Figure 6, which show discrete sub-inertial resonances with enhancements of orders of magnitude near Roc ~ 1/sqrt(5), as well as super-inertial peaks. The nonrotating limits recover the Press and Lecoanet-Quataert scalings. The paper closes with caveats about the single-scale convection model, neglected damping, and the omission of global-scale and plume-like flows.","tokens_in":27659,"tokens_out":15224,"duration_ms":142129,"significance":"If the results are robust, this paper provides the first quantitative, rotation-dependent prescriptions for GIW excitation amplitudes across a wide range of stellar convective Rossby numbers, with direct implications for asteroseismic mode visibility and angular-momentum transport in rotating early-type stars and pre-main-sequence low-mass stars. The strengths are substantial: the derivations are fully analytic, no parameter in the wave-flux ratios is fitted to observations or simulations, the nonrotating limits correctly recover the Press (1981) and Lecoanet-Quataert (2013) scalings, and the f-plane treatment retains nontraditional Coriolis components. The final predictions are falsifiable against global simulations, laboratory experiments, and asteroseismic data. The main reservation is that the quantitative peaks and 'orders-of-magnitude' enhancements are not yet shown to survive the two most obvious idealizations of the model: a monochromatic delta-function convective source and zero wave damping.","major_comments":[{"comment":"The central quantitative claim of order-of-magnitude sub-inertial flux enhancements rests on an undamped, monochromatic forcing. The source in Eqs. (69)-(70) contains a single time frequency omega_c and a single horizontal wavenumber k_perp, so the time and horizontal integrals impose omega = 2 omega_c and k_perp = k_n/2, as stated after Eq. (71); Eq. (50) contains no damping term. In an undamped cavity, resonant forcing produces formally divergent or secular responses, and the finite peaks in Figure 6 are selected cavity resonances. The manuscript acknowledges the neglect of damping in Section 5 and lists related limitations in Section 6, but it gives no quantitative sensitivity estimate. I request at minimum a test with a finite-width source spectrum and a damping rate in the amplitude equation, and a demonstration that the Roc ~ 1/sqrt(5) peaks survive such broadening; as it stands, the 'orders-of-magnitude enhancement' is not established as a robust prediction.","section":"§5, Eqs. (50), (65)-(72), Fig. 6"},{"comment":"The interfacial flux ratio inherits the Press (1981) delta-function source assumption, which Section 4 explicitly states 'permits only a single horizontal spatial scale ... and a single time scale'. The critical-Rossby peak at Roc,crit = sqrt(6)/(5*pi*S) and the decrease below it are therefore predictions of a single-scale model. The authors should quantify the effect of a finite spectral width of the convection, for instance by convolving the right-hand side of Eq. (43) with a plausible k-omega spectrum, or at least give a quantitative criterion for when the monomodal approximation is valid. Without such a test, the location and existence of the interfacial peak are not robust to the assumed source spectrum.","section":"§4, Eq. (43), Figs. 2-3"},{"comment":"As printed, the denominator of the amplitude formula is not the same as the normalization coefficient defined just below it. Equation (50) contains (f^2 - N^2)/(f^2 - omega^2) + ((omega^2 + f^2) f_ss^2)/(f^2 - omega^2)^2 in the denominator, whereas Eq. (51) defines c_n^2 with (omega^2 - N^2)/(f^2 - omega^2) + (omega^2 f_ss^2)/(f^2 - omega^2)^2. Since the absolute flux level in Figure 6 is controlled by this denominator, the authors should clarify which inner product results from projecting Eq. (49) onto the conjugate mode and correct any typographical inconsistencies between the two expressions.","section":"§5, Eqs. (50) and (51)"}],"minor_comments":[{"comment":"The transition from Eq. (39) to Eq. (40) uses k_perp ~ omega_c/v and changes the scaling from v^4 to v^3; please state explicitly that the velocity in the final expression is the rms convective velocity v(Roc) from Paper I and check the notation consistency.","section":"§4, Eq. (40)"},{"comment":"The symbol S is used both for the interface stiffness and for the convective source term in Section 5; this is confusing because both appear in close proximity (e.g., Eq. (44) and Eq. (52)). Please rename one of them.","section":"Table 2"},{"comment":"The Heaviside function H that confines the convection to the convective region is introduced in the text but never defined explicitly; please define H(z) and its support clearly.","section":"§5, Eq. (56)"},{"comment":"Typographical errors such as 'semi-analytically been modelled' in the Introduction (second paragraph) and 'tantilizing' in Section 6 (third paragraph) should be corrected.","section":"Introduction and §6"},{"comment":"The caption states 'S = NR/N0' as the stiffness ratio, but the figure panel descriptions refer to S = NR/N0 with N0 only defined implicitly through Paper I; please define N0 at first use and make the normalization of the y-axis in each panel explicit.","section":"Fig. 6 caption"},{"comment":"The discrete nature of the plotted flux values is mentioned in the text but is not apparent in the figure; it would improve readability to mark the discrete eigenvalues k_n used for each curve or to state explicitly that points are connected only for visualization.","section":"Fig. 6"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThis is a serious analytic extension of Augustson & Mathis's rotating convection model to gravito-inertial wave excitation. It gives the field its first quantitative, parameter-dependent flux scalings for both interfacial and Reynolds-stress driving, including a critical convective Rossby number below which the interfacial sub-inertial flux drops and a predicted enhancement near Roc ~ 1/sqrt5. The nonrotating limits recover the Press and Lecoanet-Quataert scalings, and nothing is fitted to target wave fluxes. The self-citation is heavy but justified: this paper builds directly on Paper I and Mathis et al. (2014), and the wave eigenmode calculation is independent.\n\nThe soft spots are where the stress-test note lands. The convective source is a single monomodal sinusoid, so the Fourier-space source is a delta function. That forces the wave response into discrete cavity resonances with no damping, so the peaks in Figure 6 are selected resonant eigenmodes. Real rotating convection has a broad spectrum, finite correlation times, and damping, all of which will broaden and reduce those peaks. The authors list these limitations in Section 6, but they do not quantify the sensitivity. That is not a fatal objection—the calculation is internally consistent and the limiting cases are credible—but the quantitative amplitudes, especially the \"orders of magnitude\" sub-inertial enhancements, should be treated as qualitative predictions.\n\nThe other modeling choices—Malkus heat-flux maximization, discontinuous Brunt-Väisälä profile—are recognized approximations. The Malkus closure is not something I would defend as a physical principle, but it is not a hidden circularity; the results are analytic consequences of stated assumptions.\n\nWho gets value: anyone modeling wave-induced angular momentum transport or interpreting GIW asteroseismology in rotating early-type stars. It deserves a serious referee; I would send it out rather than desk reject. The referee should press for sensitivity tests to a broad convective spectrum and to damping, and for a clear statement that Figure 6's peaks are not secure amplitudes. Even with those caveats, this is a useful contribution.","headline":"First quantitative rotation scalings for gravito-inertial wave excitation, but the delta-function source and undamped resonances mean the peak fluxes should be read as qualitative predictions rather than measured amplitudes.","tokens_in":28213,"tokens_out":2397,"would_cite":true,"duration_ms":27358,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Rotation sets the strength of gravito-inertial waves leaving convective zones, with a critical spin rate dividing boosted from suppressed regimes.","keywords":["gravito-inertial waves","rotating convection","convective Rossby number","wave excitation","Reynolds stresses","interfacial driving","asteroseismology","internal wave transport"],"falsifier":"Direct numerical simulations of rotating convection in a two-layer stratified domain could measure the radiated gravito-inertial wave flux as a function of $Ro_c$ and frequency and check whether the flux peaks near $Ro_c \\approx 1/\\sqrt{5}$ and whether the sub-inertial interfacial flux turns over at $Ro_{c,\\mathrm{crit}} = \\sqrt{6}/(5\\pi S)$; observationally, one could search for a spin-rate-dependent enhancement of gravito-inertial mode amplitudes among stars of similar mass and age.","tokens_in":27169,"feed_emoji":"🌊","tokens_out":10565,"duration_ms":90003,"temperature":0.7,"pith_summary":"Rotation sets the strength of the gravito-inertial waves that leave stellar and planetary convection zones and transport angular momentum, chemicals, and energy through stable regions. Using the single-eddy, heat-flux-maximized convection model of Paper I, the paper computes two excitation channels: pressure-driven waves at the convective-radiative interface and Reynolds-stress-driven waves distributed through the convection zone. The central quantitative claim is that the flux relative to nonrotating gravity waves depends on the convective Rossby number $Ro_c$, the interfacial stiffness $S$, and the wave frequency, with a critical value $Ro_{c,\\mathrm{crit}} = \\sqrt{6}/(5\\pi S)$; below it, interfacially excited sub-inertial waves lose flux and frequency range, while near $Ro_c \\approx 1/\\sqrt{5}$ Reynolds-stress-driven sub-inertial modes can be boosted by orders of magnitude while staying below the convective flux. If these regimes hold, they give asteroseismology a rotation-sensitive probe of convection and imply that wave-driven transport peaks at specific evolutionary stages.","feed_headline":"Rotation can amplify internal-wave flux by orders of magnitude","feed_subtitle":"A rotating-convection model places the boost at a narrow spin range near convective Rossby number 1/√5.","key_machinery":"The load-bearing object is the heat-flux-maximized monomodal convection model from Paper I: convection is represented as a single vertical sinusoid with one horizontal wavenumber $k_\\perp$ and one convective frequency $\\omega_c$, with velocity amplitude $v(Ro_c)$ fixed by maximizing the heat flux at fixed maximum flux. This single eddy makes the wave source a delta function in Fourier space, so the wave response can be written in closed form. Two wave calculational frameworks are attached to it: interfacial excitation, which matches the turbulent pressure to the gravito-inertial wave pressure perturbation across the boundary, and volumetric Reynolds-stress excitation, which solves the linearized Boussinesq wave equation (the Poincar\\'e equation, retaining the full Coriolis acceleration) in a two-layer domain with a discontinuous Brunt-V\\\"ais\\\"al\\\"a jump, using exact eigenmodes and dispersion relations for sub-inertial and super-inertial waves. The ratio $F_z/F_0$ cancels many unknown amplitudes and isolates the rotation and stiffness dependence.","core_discovery":"The paper's central discovery is that rotation leaves a strong, parameter-dependent fingerprint on the flux of gravito-inertial waves, and that the fingerprint is different for the two main excitation mechanisms. For interfacial driving, the wave flux relative to nonrotating gravity waves is approximately $(v/v_0)^3$ times a rotational factor (Equation 43); sub-inertial waves are enhanced as $Ro_c$ decreases only down to $Ro_{c,\\mathrm{crit}} = \\sqrt{6}/(5\\pi S)$, after which the flux and the propagative frequency band both shrink, while super-inertial waves are appreciably enhanced only for $Ro_c$ near unity. For Reynolds-stress driving, solving the forced Poincar\\'e equation with the monomodal convection source yields sub-inertial modes whose flux can exceed the nonrotating gravity-wave flux by orders of magnitude near $Ro_c \\approx 1/\\sqrt{5}$, with the exact pattern depending on interface stiffness through resonance; the super-inertial flux peaks near the same Rossby number and asymptotes to the nonrotating case. In all computed regimes the gravito-inertial wave flux remains below the total convective flux. The same parameter map places the predicted peaks in regimes that low-mass stars cross during pre-main-sequence evolution and near the base of their convective envelopes, so rotation-dependent wave transport and seismic signatures are expected to vary along stellar evolution.","pith_inferences":["Our inference: a realistic broad-band convection spectrum would smear the delta-function source, so the sharp resonant peaks in the Reynolds-stress flux likely become broader and weaker, and the precise value $1/\\sqrt{5}$ should be interpreted as a model-dependent marker rather than a universal scale.","Our inference: because the paper computes energy injection but not damping, the transported flux could be substantially lower than the injected flux for resonantly excited modes, so asteroseismic inferences need a companion damping model to compare with observed amplitudes.","Our inference: replacing the discontinuous Brunt-V\\\"ais\\\"al\\\"a jump with a smooth stratification may trap sub-inertial waves in regions of weak stratification, moving the energy deposition away from the interface and changing mixing profiles.","Our inference: plumes and global-scale shear excluded by the local monomodal model may dominate interfacial excitation; a plume-based or global rotating convection model could shift the predicted peak to lower or higher Rossby numbers."],"forward_implications":["Below $Ro_{c,\\mathrm{crit}}=\\sqrt{6}/(5\\pi S)$, interfacially driven sub-inertial gravito-inertial waves lose both energy flux and propagative frequency range, so wave-driven angular-momentum transport should weaken in the most rapidly rotating stars and planets.","Near $Ro_c\\approx 1/\\sqrt{5}$, Reynolds-stress-driven sub-inertial modes can carry orders of magnitude more energy than nonrotating gravity waves while remaining below the convective flux, making their period spacings attractive asteroseismic targets.","Super-inertial waves are enhanced only for $Ro_c$ near unity, so their seismic visibility should track a narrow range of rotation rates.","The integrated interfacial flux decreases with increasing rotation overall, with a modest peak near $\\Omega/\\Omega_b \\approx 0.091$ for a stiffness $S=10^3$, implying a rotation window for interfacial wave transport.","Because low-mass stars with typical initial rotation pass through the predicted enhancement regimes during pre-main-sequence evolution, wave-induced transport and mixing should be stronger at those ages."],"supporting_citations":[{"why":"Establishes the interfacial pressure-matching method and the nonrotating gravity-wave flux $F_0$ against which all rotating fluxes are normalized.","marker":"Press (1981)"},{"why":"Paper I, whose heat-flux-maximized monomodal convection model provides $v(Ro_c)$, $k_\\perp$, and $\\omega_c$ for both excitation channels.","marker":"Augustson & Mathis (2019)"},{"why":"Provides the rotating-convection scaling laws from which the Paper I model is derived.","marker":"Stevenson (1979)"},{"why":"Supplies the gravito-inertial energy-flux and pressure-perturbation relations used to build the interfacial flux ratio.","marker":"André et al. (2017)"},{"why":"Establishes the sub-inertial/super-inertial classification and the forced Poincaré equation used in the Reynolds-stress calculation.","marker":"Mathis et al. (2014)"},{"why":"Provides the nonrotating Reynolds-stress excitation formalism and the $F_z \\propto F_c/S$ limit that the rotating calculation generalizes.","marker":"Lecoanet & Quataert (2013)"},{"why":"Supplies the modal expansion of the Poincaré equation with nontraditional Coriolis terms used for the wave eigenfunctions.","marker":"Gerkema & Shrira (2005)"},{"why":"Supplies the stochastic-excitation amplitude formalism that underlies the Reynolds-stress flux computation.","marker":"Goldreich & Kumar (1990)"}],"fun_headline_variants":["Rotation boosts internal wave flux at specific spins","Wave flux peaks near Rossby 1/√5 for stress-driven modes","Interfacial wave flux peaks at a stiffness-dependent spin","Sub-inertial waves gain from rotation, but not all spins","Gravito-inertial wave flux shows sharp rotational dependence"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the convection that drives the waves can be represented by one sinusoidal eddy with a single horizontal scale and a single frequency; if the broad spectrum, large-scale shear, and plume-like events of real rotating convection contribute substantially, the location and even the existence of the predicted flux peaks are not guaranteed.","fun_headline_variants_meta":{"raw":{"variants":["Rotation boosts internal wave flux at specific spins","Wave flux peaks near Rossby 1/√5 for stress-driven modes","Interfacial wave flux peaks at a stiffness-dependent spin","Sub-inertial waves gain from rotation, but not all spins","Gravito-inertial wave flux shows sharp rotational dependence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000458,"raw_usage":{"total_tokens":2380,"prompt_tokens":1109,"completion_tokens":1271,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":725,"completion_tokens_details":{"reasoning_tokens":1186}},"tokens_in":725,"tokens_out":1271,"duration_ms":12203,"temperature":1.0,"reasoning_tokens":1186,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-27T23:21:45.311440+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Direct numerical simulations of rotating convection in a two-layer stratified domain could measure the radiated gravito-inertial wave flux as a function of $Ro_c$ and frequency and check whether the flux peaks near $Ro_c \\approx 1/\\sqrt{5}$ and whether the sub-inertial interfacial flux turns over at $Ro_{c,\\mathrm{crit}} = \\sqrt{6}/(5\\pi S)$; observationally, one could search for a spin-rate-dependent enhancement of gravito-inertial mode amplitudes among stars of similar mass and age.","supporting_citations":[{"cited_title":", archivePrefix = \"arXiv\", eprint =","cited_arxiv_id":null,"evidence_quote":"Establishes the sub-inertial/super-inertial classification and the forced Poincaré equation used in the Reynolds-stress calculation."},{"cited_title":"Journal of Fluid Mechanics , year = 2005, month = apr, volume = 529, pages =","cited_arxiv_id":null,"evidence_quote":"Supplies the modal expansion of the Poincaré equation with nontraditional Coriolis terms used for the wave eigenfunctions."}],"review_version":1}