{"id":"e704753c-ddba-4e0f-87c9-5789db963a47","arxiv_id":"2607.19730","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Global simulations show a prescribed latitudinal entropy gradient, in thermal-wind balance, can sustain the Sun's observed differential rotation without strong convective Reynolds stresses.","lead":"Does the Sun's differential rotation require turbulent convection, or can a background heat imbalance alone keep the equator spinning faster than the poles? This paper's simulations show that a latitude-dependent thermal forcing — the same thermal-wind balance that drives atmospheric winds — can produce and sustain solar-like rotation even without vigorous convection.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The prescribed, solar-derived entropy gradient Θ_B is the unmodeled load-bearing premise; simulations show dynamical consistency but not that such a gradient is sustained in the Sun.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: the prescribed latitudinal entropy gradient has no modeled physical origin in the Sun, so solar applicability is conditional on an ingredient the paper explicitly defers. I agree with that assessment. The manuscript is honest and well-scoped as a dynamical-viability study, but the central claim for the Sun cannot be stronger than the support for the sustained Θ_B. The reader's CONDITIONAL verdict is appropriate. The proposed persistence test would quantify how much of the maintained differential rotation is directly due to the nudging, and thus how much weight falls on the unmodeled source. No change to the reader's verdict is needed.","tokens_in":12245,"tokens_out":5741,"duration_ms":73096,"concrete_test":"Turn off the Newtonian cooling after a statistically steady solar-like state is reached (set τ→∞ in Eq. 3, or equivalently hold Θ_B only in the initial condition) and evolve the same run for several hundred years, measuring the decay rate of ⟨Ω⟩'s equator-pole contrast and of the latitudinal Θ' gradient. If the solar-like state decays on a timescale shorter than the solar cycle, the maintained state is entirely a product of the externally prescribed Θ_B, and the solar scenario needs an independent, sustained entropy-gradient source; if it persists, the baroclinic state is self-sustaining and the concern is substantially reduced.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is solar applicability: that a background latitudinal entropy gradient in thermal wind balance can generate and maintain the observed differential rotation without strong convective Reynolds stresses. The load-bearing premise is that such a gradient is present and sustained in the real convection zone. The paper does not model or justify this premise. In Appendix A Eq. (A1), Θ_B is constructed by integrating thermal wind balance from the helioseismic Ω profile, so the target itself encodes the observed rotation. In Eq. (3), Newtonian cooling with τ ≈ 2.4–4.9 yr continuously relaxes Θ' toward this target. The simulations therefore demonstrate that a baroclinic state designed to match the Sun is dynamically compatible with nearly solar Ω, not that the Sun possesses the entropy gradient. The authors acknowledge this: §1 states the gradient is 'prescribed as a controlled forcing... rather than to model its physical origin,' and §4 leaves the origin for future work. Moreover, Section 3 reports that only weakly convective cases reproduce the observed radial increase at low latitudes and that Reynolds stresses contribute to maintaining that increase, qualifying the 'without relying on Reynolds stresses' headline even within the model. None of this is internally inconsistent; but if convection in the Sun mixes away any such gradient, the proposed mechanism does not operate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates whether solar differential rotation can be sustained by a prescribed large-scale latitudinal entropy gradient rather than by Reynolds stresses from vigorous convection. Global anelastic EULAG-MHD simulations of a spherical shell (0.72–0.96 R_sun) are run for 20 combinations of polytropic index (stable, adiabatic, weakly unstable) and two Newtonian-cooling timescales, with a target potential-temperature profile Θ_B derived from helioseismic differential rotation via thermal-wind balance. Time-averaged results show a fast equator and slower poles across most runs, with radially aligned mid-latitude contours; torque analysis attributes the maintenance primarily to meridional circulation, with Reynolds stresses playing a smaller or auxiliary role. The authors conclude that baroclinic forcing can maintain solar-like differential rotation without strong convective Reynolds stresses, offering a possible resolution to the convective conundrum, while acknowledging that the physical origin of the imposed entropy gradient is not modeled.","tokens_in":12510,"tokens_out":8147,"duration_ms":93415,"significance":"The paper has genuine strengths: the equations, boundary conditions, grid, and parameter choices are fully specified; the parameter sweep is broad; and the simulations are evolved long enough for robust temporal averages. The torque decomposition and thermal-wind diagnostics are appropriate and clearly presented. If the claims are scoped carefully, the study is a useful dynamical proof-of-concept: it shows that a baroclinic state designed to be in thermal-wind balance is nonlinearly compatible with a fast-equator rotation profile across a range of stratifications. However, the central solar-application claim is weakened by the construction of Θ_B from the very differential-rotation profile the paper explains, and by the fact that the entropy gradient is externally maintained rather than shown to exist or persist in the Sun. The paper is therefore more convincing as a controlled dynamics study than as an explanation of the origin of solar differential rotation.","major_comments":[{"comment":"The target profile Θ_B is constructed by integrating the thermal-wind balance using the observed helioseismic Ω profile (Eqs. A2–A3), and Eq. (3) relaxes the simulation toward this target. The solar-like Ω in the runs is therefore at least partly imposed: the forcing encodes the very observable being explained. The paper states explicitly that the physical origin of the entropy gradient is not modeled (§1, §4), but the abstract and conclusion present the result as an 'alternative scenario' for the solar rotation profile. As it stands, this is a demonstration of dynamical consistency, not a demonstration that the Sun possesses a sustained entropy gradient of the required amplitude and structure. Please narrow the claims accordingly, or provide a quantitative estimate of the required Θ_B and a discussion of whether such a gradient is consistent with observational constraints and with the r","section":"Appendix A, Eq. (A1); §4"},{"comment":"Eq. (3) relaxes Θ′ toward Θ_B, but the boundary conditions immediately below require Θ′=0 at r=0.72 and 0.96 R_sun. The target Θ_B shown in Fig. 2 has latitudinal structure at these radii (its latitudinal average is zero, but pointwise values are not), so the homogeneous Dirichlet condition is incompatible with the nudging target at the boundaries. This will create boundary layers in Θ′ and corresponding baroclinic torques. Please clarify how the boundary condition is imposed in the code, whether the source term is suppressed at the boundary cells, and whether the reported time averages—including Fig. 3 and the torque integrals of Fig. 4—are sensitive to this treatment.","section":"§2, boundary conditions after Eq. (3)"},{"comment":"The paper's headline claim that solar-like differential rotation is maintained 'without relying on Reynolds stresses' is internally qualified. §3 states that 'weakly convective models ... are the only ones capable of reproducing the observed radial increase of Ω at low latitudes' and that 'the Reynolds stresses are responsible for torques capable of maintaining the radial increase of the rotation rate.' The radial increase at low latitudes is part of the observed solar profile, so the abstract and conclusions overstate the result. Please recast the claim to distinguish the baroclinically driven fast-equator/pole contrast from the low-latitude radial shear, which in the successful simulations is maintained with assistance from weak convective Reynolds stresses.","section":"§3, first paragraph; §4"}],"minor_comments":[{"comment":"The symbols C and B are used before they are defined; define them explicitly or insert a sentence introducing them.","section":"Eq. (10)"},{"comment":"The phrase 'wavelength agnostic diffusivity' is vague. Since the Fourier-space form −τ^{-1}Θ̂′ is given, please explain the intended physical interpretation more concretely.","section":"§2, paragraph on nudging"},{"comment":"Spelling is inconsistent: 'superabadicity' vs. 'superadiabaticity' and 'heliosseismic' vs. 'helioseismic'. Please check the whole text.","section":"Throughout"},{"comment":"The concluding paragraph ends with a stray 'a' after '...viable alternative to convection-dominated paradigms and contributes toward resolving long-standing discrepancies.'","section":"§4"},{"comment":"Please state the units and color-scale range for the potential-temperature panel; currently the reader cannot estimate the amplitude of Θ_B.","section":"Fig. 2 caption"},{"comment":"Entries in the second column like '1.5−1×10−4' are easy to misread as 1.5×10−4. Use explicit arithmetic (e.g., '1.4999') or add a space.","section":"Table 1"},{"comment":"In the reference for Stefan et al. 2026, 'Guerreo' appears to be a typo for 'Guerrero'. Please check.","section":"References"},{"comment":"No explicit viscosity or diffusivity appears in Eqs. (1)–(3). If the runs rely on EULAG's numerical dissipation as an implicit subgrid model, please state this explicitly, since it is relevant for interpreting the Reynolds-stress torque decomposition.","section":"§2, numerical model"}],"recommendation":"major_revision","confidential_remarks":"The dynamical-consistency result appears credible and the parameter study is well executed. The main issues are the circular construction of the forcing and the overstatement in the abstract/conclusions relative to what is actually demonstrated. The boundary-condition incompatibility at the shells also needs a clear answer, as it could affect the torque balance. I recommend major revision rather than rejection; the central idea is worth publishing after the claims are properly scoped and the technical point is addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Honest, well-scoped numerical experiment. The new thing is not the thermal-wind idea — Balbus, Gunderson, Hester and others have been there — but the nonlinear spherical-shell demonstration. Across 20 runs spanning stable, adiabatic, and weakly unstable stratifications, the authors show that a prescribed latitudinal entropy gradient in thermal wind balance drives a steady, fast-equator/slower-pole rotation profile with mid-latitude radial contours, and their torque decomposition attributes the maintenance to meridional circulation rather than Reynolds stresses. That is real evidence for dynamical viability, and the diagnostics are careful: 456 years of evolution, 76-year averages, explicit equations and boundary conditions.\n\nThe soft spot is the load-bearing premise. The entropy gradient is not modeled; it is imposed as Θ_B, and that target is constructed from the very helioseismic Ω the paper aims to explain (Appendix A). The Newtonian cooling nudges toward it, so a solar-like outcome is partly inherited. The authors are open about this — they call it a controlled forcing and defer the origin — but it means the paper demonstrates consistency, not that the Sun possesses such a gradient. If convection mixes the gradient away, the mechanism fails. Their own results also show fragility: only weakly convective cases reproduce the observed radial increase at low latitudes, and strongly superadiabatic cases collapse to cylindrical contours. And the headline claim about not relying on Reynolds stresses is qualified even within the model, since §3 says Reynolds stresses help maintain that radial increase in the weakly convective cases.\n\nMinor but real: no error bars on the time averages, no quantitative match metric against helioseismic Ω, and no released code or data. All addressable.\n\nMy overall take agrees with the reader's conditional verdict: this is a useful proof-of-concept that the baroclinic scenario is dynamically plausible, not a demonstration that it operates in the Sun. The paper is honest, internally consistent, and clearly scoped. It deserves a serious referee; the authors have done the work to make a referee's job straightforward. I'd want the entropy-gradient origin acknowledged more forcefully as a limitation, but it already is. Worth engaging with.","headline":"A careful spherical-shell demonstration that a prescribed baroclinic forcing can sustain solar-like differential rotation without strong convection, but the solar scenario rests on an unmodeled entropy gradient.","tokens_in":13146,"tokens_out":2441,"would_cite":true,"duration_ms":24755,"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":"The Sun's differential rotation—fast equator, slow poles—can be sustained by a small latitudinal entropy gradient in thermal wind balance, without needing vigorous convection or Reynolds stresses. The authors demonstrate this through 20 glo","keywords":["solar differential rotation","baroclinic forcing","thermal wind balance","convective conundrum","latitudinal entropy gradient","global hydrodynamic simulations","gyroscopic pumping","solar convection zone"],"falsifier":"A direct helioseismic measurement of the latitudinal entropy gradient in the convection zone, showing it to be absent or far smaller than the thermal-wind value needed to balance the observed rotation, would falsify the mechanism. Alternatively, detection of Reynolds stress torques dominating over meridional-circulation torques in the Sun (as deduced from inversions of flow correlations) would contradict the paper's central claim.","tokens_in":11994,"feed_emoji":"☀️","tokens_out":1566,"duration_ms":18457,"temperature":0.7,"pith_summary":"This paper argues that the Sun's observed differential rotation (equator rotating faster than poles, with radially aligned mid-latitude contours) can be generated and maintained without requiring strong turbulent convection and its associated Reynolds stresses. Instead, the authors propose that a background latitudinal entropy gradient, held in thermal wind balance, is enough to drive the flow. They test this idea with a suite of global hydrodynamic simulations spanning stable, adiabatic, and weakly unstable stratifications, finding that all produce a solar-like rotation profile. The result offers a way around the 'convective conundrum', where mixing-length theory predicts much stronger convection than helioseismology observes.","feed_headline":"Solar rotation may be set by heat, not convection","feed_subtitle":"Simulations show a latitudinal entropy gradient alone can sustain the fast equator and slow poles, easing the convective conundrum.","key_machinery":"The central object is the prescribed latitudinal potential-temperature gradient Θ_B, built by integrating the thermal wind balance relation (Eq. A1) against the helioseismic rotation profile. The stratification is enforced through a Newtonian cooling term in the potential-temperature equation (Eq. 3) with a long timescale (2.4–4.9 yr). The dynamical balance is expressed in the thermal wind equation: the centrifugal term C = r sinθ ∂Ω²/∂z is balanced against the baroclinic term B = (g/(Θ₀ r)) ∂Θ′/∂θ. The small residual C−B drives gyroscopic pumping, producing meridional circulation cells that redistribute angular momentum and sustain the differential rotation without the need for turbulent st","core_discovery":"The central claim is that solar-like differential rotation can be generated and maintained in the convection zone without relying on Reynolds stresses from vigorous turbulent convection, provided a large-scale latitudinal entropy gradient exists and is sustained in thermal wind balance. In a rotating spherical shell with anelastic equations and a nudged potential-temperature profile, all simulations—whether stably stratified, adiabatic, or weakly superadiabatic—produce a rotation profile with a faster equator and slower poles. The torque analysis shows that the meridional circulation, rather than Reynolds stresses, is the main carrier of angular momentum in these runs; the Reynolds stresses","pith_inferences":["A testable corollary of the paper's logic is that the observed solar rotation profile implies the existence of a sustained latitudinal entropy gradient of roughly the magnitude of Θ_B plotted in Fig. 2; if helioseismic or other measurements ever rule out such a gradient over multi-year timescales, the mechanism would be in conflict.","The paper's preference for weak or stable stratifications suggests that the Sun's interior may contain subadiabatic layers (as in the Deardorff-layer picture); an inference is that the location and sharpness of such layers would control the shape of the rotation profile and could be constrained by the observed radial behavior near the tachocline.","Extending this hydrodynamical framework to include magnetic fields, as the authors suggest, would test whether the baroclinically maintained rotation can survive Lorentz forces; one might predict that the magnetically active cycle modifies the entropy gradient and hence the rotation profile, providing a link to the observed torsional oscillations.","The paper does not prescribe the origin of the entropy gradient; an inference is that any convective or tachocline-driven process that creates a large-scale latitudinal entropy variation (e.g., rotational modulation of heat transport or entropy rain) would suffice, meaning the mechanism is robust to the exact source as long as the gradient persists."],"forward_implications":["If correct, the solar convection zone can be much less turbulent than mixing-length theory predicts, easing the convective conundrum.","The observed differential rotation profile becomes a signature of large-scale thermodynamic balance rather than a direct measure of convective Reynolds stresses.","The same baroclinic mechanism could operate in other solar-type stars, giving a common explanation for their rotation profiles.","Models of the solar dynamo, which rely on differential rotation to stretch magnetic fields, may need to account for a rotation profile that is thermodynamically constrained rather than convectively driven.","The work suggests that the radial tilt of the rotation contours at mid-latitudes emerges naturally from gyroscopic pumping of the baroclinic imbalance, offering a direct explanation for this long-puzzling feature."],"fun_headline_variants":["Heat, not turbulence, sets Sun's rotation","Solar differential rotation from entropy gradients alone","Baroclinic heat flow spins the Sun, not convection","Latitudinal heat gradient can power Sun's differential spin","Thermal wind balance may drive Sun's rotation without convection"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The Sun's convection zone must actually maintain a persistent large-scale latitudinal entropy gradient, in thermal wind balance with the observed rotation, over timescales of years; if convection mixes that gradient away faster than it can be replenished, the proposed mechanism will not operate in the real Sun.","fun_headline_variants_meta":{"raw":{"variants":["Heat, not turbulence, sets Sun's rotation","Solar differential rotation from entropy gradients alone","Baroclinic heat flow spins the Sun, not convection","Latitudinal heat gradient can power Sun's differential spin","Thermal wind balance may drive Sun's rotation without convection"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000507,"raw_usage":{"total_tokens":2247,"prompt_tokens":620,"completion_tokens":1627,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":364,"completion_tokens_details":{"reasoning_tokens":1552}},"tokens_in":364,"tokens_out":1627,"duration_ms":12428,"temperature":1.0,"reasoning_tokens":1552,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T11:52:37.079950+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct helioseismic measurement of the latitudinal entropy gradient in the convection zone, showing it to be absent or far smaller than the thermal-wind value needed to balance the observed rotation, would falsify the mechanism. Alternatively, detection of Reynolds stress torques dominating over meridional-circulation torques in the Sun (as deduced from inversions of flow correlations) would contradict the paper's central claim.","supporting_citations":[],"review_version":1}