{"id":"2f1b5e68-6d02-493a-960b-1bd42110bc0d","arxiv_id":"2502.07006","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"3D radiative hydrodynamic simulations of a 1.47-solar-mass F star show rotation-induced radius decrease, differential rotation, meridional flows, and latitude-dependent roll-like convection.","lead":"Computer simulations of a rotating F-type star show that faster rotation shrinks the simulated photosphere, slows surface flows, and organizes convection into elongated rolls that vary with latitude. The results offer a peek at how rotation and turbulent convection interact in stars that are heavier than the Sun.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The f-plane model cannot support the claimed global radius decrease and von Zeipel gravity-darkening attribution, because the omitted centrifugal force dominates rotational deformation.","rationale":"The paper's simulations are carefully constructed and likely capture genuine local effects of the Coriolis force on turbulent convection in a shallow F-star convection zone: the development of differential rotation, meridional flows, and roll-like structures follows from the f-plane equations and is internally plausible. I therefore do not reject the entire manuscript. However, the abstract and conclusions elevate these local results into global stellar properties: a decrease of the stellar radius by about 29/58 km and a 'gravity-darkening effect.' These claims are not derivable from independent Cartesian boxes with periodic horizontal boundaries and no centrifugal force. The concrete estimate of Omega^2 R^3 / (G M) ~ 0.028 shows that the centrifugal force, if included, would alter the equatorial radius by ~10^4 km, dwarfing the quoted tens of kilometers and reversing the sign of the radius change. Thus the f-plane approximation is not a minor simplification for the radius and gravity-darkening claims; it omits the dominant physical mechanism. The internal inconsistency between Section 4.3's 40 km photosphere shift and the abstract's 29 km at the equator further weakens the quantitative headline. The von Zeipel attribution is also questionable because the simulated pole is cooler than the equator, whereas classical gravity darkening makes the pole hotter. These issues are correctable by reframing the results as local rotation-convection coupling and by removing or rigorously qualifying the global radius and gravity-darkening statements. That is why the verdict should remain conditional on major revision rather than moving to accept or reject outright.","tokens_in":19722,"tokens_out":11711,"duration_ms":108756,"concrete_test":"Compute the rotating-structure prediction for this star using a 1D Roche or ESTER model with CESAM inputs and P=1 day, and compare the equator-to-pole effective-temperature difference and equatorial radius change with the paper's values. If the predicted pole-equator temperature difference is positive (about 50 K for a radiative envelope, from Delta T/T = (1/4) Omega^2 R^3 / (G M)) and the equatorial radius change is about +10^4 km rather than -29 km, the f-plane interpretation of 'radius decrease' and 'gravity darkening' is refuted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and Section 7 claim a stellar radius decrease of about 29 km at the equator and 58 km at higher latitudes, plus a 'gravity-darkening effect,' but these global statements cannot be derived from the f-plane Cartesian boxes described in Section 2, which explicitly neglect surface curvature and centrifugal force. For this 1.47 M_sun, R=1.456 R_sun star at P=1 day, Omega^2 R^3 / (G M) is approximately 0.028, so the centrifugal contribution to the equatorial radius is of order (1/2) Omega^2 R^3 / (G M) times R, roughly 1.4e9 cm or about 14,000 km, five orders of magnitude larger than the quoted tens of kilometers and of opposite sign. A fixed box with periodic horizontal boundaries cannot capture this global deformation. Section 4.3 also reports a 'downward shift of the photosphere by about 40 km,' while the abstract and conclusions quote 29 km at the equator and 58 km at higher latitudes; these numbers are not reconciled. Furthermore, the simulated photospheric temperature pattern (equator hotter by 90 K at 30 degrees and 200 K at 60 degrees) is opposite in sign to the classical von Zeipel prediction, which gives a hotter pole for radiative envelopes, so labeling this 'gravitational darkening' is unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a series of 3D radiative hydrodynamic simulations of the outer layers of a 1.47 solar-mass F-type star using the StellarBox code. The simulations are Cartesian f-plane domains placed at latitudes of 0, 30, and 60 degrees, with imposed rotation periods of 1 and 14 days. The authors report that rotation modifies the convection structure and large-scale dynamics, including differential rotation, meridional flows, roll-like convective patterns, shifts of ionization zones, a photosphere/radius decrease of tens of kilometers, and a latitude-dependent surface temperature pattern that they attribute to gravity darkening. The comparison baseline is the authors' earlier non-rotating simulation of the same star.","tokens_in":19994,"tokens_out":4018,"duration_ms":40288,"significance":"If the global interpretations were supported, the simulations would provide novel quantitative predictions for how rotation alters the structure and dynamics of a shallow-envelope F-type star, with implications for stellar evolution, angular momentum transport, and activity. The numerical experiments are self-contained and use a realistic equation of state, chemical composition, radiative transfer, and compressible hydrodynamics, which are strengths relative to anelastic global models. The reported large-scale flow structures are falsifiable predictions for future observations or global simulations. However, the paper's central global claims—the stellar radius decrease and the gravity-darkening attribution—are not supported by the f-plane local-box setup, which explicitly neglects surface curvature and the centrifugal force.","major_comments":[{"comment":"The claim that the stellar radius decreases by about 29 km at the equator and about 58 km at higher latitudes for P_rot = 1 day is not supported by the model. Section 2 states that the simulations use the f-plane approximation and that 'the effects of surface curvature and centrifugal force are neglected.' A local Cartesian box with periodic horizontal boundaries cannot capture a global hydrostatic radius change. For this star at P = 1 day, Omega^2 R^3/(GM) ~ 0.028, so the centrifugal contribution to the equatorial radius is of order several thousand km, several orders of magnitude larger and opposite in sign to the quoted tens of kilometers. The authors should either remove the global radius language, recast it explicitly as a local photosphere-height shift within the box, or supplement the local simulations with a global stellar-structure calculation that includes the centrifugal force.","section":"Abstract; Section 2; Section 7 bullet list"},{"comment":"The attribution of the latitude-dependent photosphere temperature decrease to 'gravitational darkening' (von Zeipel 1924) is unsupported. The simulations are independent boxes at different latitudes, with no centrifugal force and no global surface geometry, so they cannot self-consistently produce the effective-gravity variation that defines gravity darkening. Moreover, the simulated pattern has a hotter equator (90 K hotter than 30 degrees and about 200 K hotter than 60 degrees for P_rot = 1 day), whereas the classical von Zeipel law for a radiative envelope would predict a cooler equator under centrifugal deformation. The manuscript should reframe this result as a rotation-induced, latitude-dependent photosphere temperature variation in local models, without invoking the gravity-darkening mechanism.","section":"Section 7; Abstract"},{"comment":"The reported magnitude of the radial/radius change is internally inconsistent. Section 4.3 states that rotation causes 'a downward shift of the photosphere by about 40 km for the 1-day period of rotation,' while the abstract and the Section 7 bullet list quote 29 km at the equator and 58 km at higher latitudes. These numbers do not agree, and the paper does not describe how the radius change was defined or measured (e.g., optical-depth surface versus a fixed pressure level). The authors should specify the diagnostic used and present a single, consistent set of values, or explain the discrepancy explicitly.","section":"Section 4.3; Abstract; Section 7"},{"comment":"The large-scale-flow profiles (differential rotation, meridional flows) are averaged over one hour in time, yet the paper presents them as robust properties of the convection zone and uses them to support the formation of persistent structures such as rolls and a tachocline-like shear layer. For a convection zone 28.5 Mm deep with strong downdrafts and large-scale rolls extending 60-80 Mm, one hour may be too short for converged statistics, especially for the weaker meridional flows (around 0.2-1.35 km/s). The authors should provide convergence tests or longer time averages to demonstrate that the reported profiles are not dominated by transients or by the finite horizontal box size.","section":"Figures 4-6; Section 4"}],"minor_comments":[{"comment":"In the Figure 4 caption, '60 o (equator, black curves)' should presumably read '0 o (equator, black curves)'; the same typo appears in the Figure 5 caption.","section":"Section 2, Figure 4 caption"},{"comment":"The subgrid-scale model coefficients C_C = C_S = 0.01 are stated to have been 'initially determined for modeling solar convection.' The paper should justify or discuss the transferability of these coefficients to a 1.47 M_sun F star, since this is an implicit modeling assumption.","section":"Section 2"},{"comment":"The inset in Figure 10b is referenced in the text ('see the inset plot in Fig. 10b'), but the inset is not visible in the reproduced figure; the authors should ensure the final figure includes it or adjust the reference.","section":"Section 5"},{"comment":"The roll-like structures are inferred from averages over the meridional plane. Because the horizontal box size is only about 5.8 degrees, the authors should explicitly discuss how the finite azimuthal extent may affect the measured roll lengths (the paper notes this for the longest rolls but could extend the caveat to the reported roll properties at all latitudes).","section":"Section 6"}],"recommendation":"major_revision","confidential_remarks":"The core local simulations appear competently executed, but the global claims in the abstract and Section 7 go beyond what an f-plane Cartesian model can support. The most direct fix is to reframe the radius and gravity-darkening statements as local photosphere-height and temperature variations and to add a clear statement of the model's limitations. The inconsistent radius numbers in Section 4.3 versus the abstract/conclusions also need to be reconciled. If the authors are unwilling to restrict the global interpretations, the paper may not be suitable for publication in its current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper is a solid numerical experiment in local convection dynamics, but the abstract and conclusions overclaim what f-plane boxes can say about a global stellar radius and gravity darkening. That part needs to be cut or rewritten before this is citable.\n\nWhat's genuinely new: this is the first 3D realistic RHD simulation of a rotating 1.47 M_sun F star, with StellarBox's radiative transfer, realistic EOS, and three latitudes (0, 30, 60) at P_rot = 1 and 14 days. The finding that rotation organizes convection into retrograde-propagating roll-like structures, and that their scale shrinks with latitude, is a clear result. The self-generated differential rotation and meridional flows, with return flows at the base, are also plausible within the modeled domain. Those parts are worth reporting.\n\nThe soft spots are real and they hit the paper's central claims. The model uses periodic Cartesian boxes on the f-plane, and Section 2 explicitly says the centrifugal force and surface curvature are neglected. You can't then claim a global radius decrease of 29-58 km and call the latitude-dependent temperature pattern gravitational darkening. The centrifugal term is order Omega^2 R^3/(GM) ~ 0.03; the resulting equatorial bulge is five orders of magnitude larger than 29 km and has the opposite sign. Even setting that aside, their simulated photosphere is hotter at the equator than at 60 degrees, the opposite of von Zeipel's prediction for an envelope like this. So the 'gravity-darkening effect' label is unsupported. Also, Section 4.3 quotes a 40 km photosphere shift while the abstract says 29 km at the equator and 58 km at higher latitudes; the three numbers are never reconciled. There are also no error bars on the radius/temperature changes, and the time averaging is only one hour.\n\nThe simulations themselves are not wasted. The local dynamics - the rolls, the shear layer, the meridional flow structure - are physically sensible and clearly described. I'd ask the authors to rework the paper around what the model can actually deliver: a local, f-plane study of rotation-convection coupling. Drop the global radius and gravity-darkening language, report the photosphere shifts as local diagnostics, and add some quantification of the averaging and boundary influences. Then it's a useful contribution.\n\nBottom line: send it to peer review, but with a referee report that demands those changes. The overclaims are fixable; the numerical work is too expensive to throw away.","headline":"A capable f-plane simulation study whose headline claims about global radius and gravitational darkening outrun the model; the local dynamics are the real story.","tokens_in":20506,"tokens_out":4269,"would_cite":false,"duration_ms":35006,"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":"Rotation measurably changes a 1.47-solar-mass F-type star: fast spin shrinks its radius by tens of kilometers, shifts its ionization zones, and drives differential rotation, meridional flows, and gravity darkening, according to 3D…","keywords":["stellar rotation","stellar convective zones","gravity darkening","hydrodynamical simulations","radiative transfer simulations","F-type stars","differential rotation","meridional circulation"],"falsifier":"Run a global spherical simulation of the same 1.47 solar-mass star at a 1-day rotation period with centrifugal force and curvature included: if the photosphere is not lower by about 29 km at the equator and about 58 km at high latitudes relative to the non-rotating model, or if high latitudes are not cooler than the equator, the radius-decrease and gravity-darkening interpretations fail. A direct observational check would be interferometric or asteroseismic radius and latitudinal brightness mapping of a similar fast-rotating F star.","tokens_in":19541,"feed_emoji":"🌟","tokens_out":10584,"duration_ms":91022,"temperature":0.7,"pith_summary":"The paper aims to establish that rotation, through the Coriolis force, changes both the dynamics and the thermodynamic structure of a shallow outer convection zone in a 1.47-solar-mass F-type main-sequence star. In a series of 3D radiative hydrodynamic simulations run in local Cartesian boxes at the equator, 30 degrees, and 60 degrees latitude, with rotation periods of 1 and 14 days, fast 1-day rotation produces a stellar radius decrease of about 29 km at the equator and 58 km at higher latitudes, shifts the hydrogen and helium ionization zones toward the photosphere, and cools the photosphere at high latitudes relative to the equator in the gravity-darkening sense. The models also self-consistently develop differential rotation, meridional flows, a tachocline, and latitude-dependent roll-like convection, all absent or weaker without rotation. A sympathetic reader cares because these signatures are potentially observable and because they show that even a very shallow convection zone (about 2.8 percent of the stellar radius) can sustain large-scale rotational coupling.","feed_headline":"Rotation shrinks an F-type star's radius by up to 58 km","feed_subtitle":"3D simulations show fast spin also cools high latitudes and reorganizes convection into roll-like flows.","key_machinery":"The central machinery is a set of 3D radiative hydrodynamic simulations in the f-plane approximation: a constant rotation vector appropriate to each latitude is imposed on a Cartesian box, with periodic horizontal boundaries, no surface curvature, and no centrifugal force. Each domain spans about 5 percent of the stellar radius in depth (50.5 Mm), from the upper radiative zone through the whole 28.5 Mm convection zone into a low atmosphere, with 102.4 Mm horizontal extent and roughly 100 km horizontal resolution. The code uses a realistic equation of state and chemical composition, time-dependent radiative transfer in four spectral bins with long-characteristics ray tracing, and a compressible Smagorinsky subgrid-scale model; initial conditions come from a stellar-evolution model of the same star. The argument depends on comparing horizontally and temporally averaged azimuthal and meridional velocities, temperature, density, energy fluxes, and vorticity against the same models without rotation, so that the mean flows, roll structures, ionization-zone shifts, and radius changes can be attributed to the imposed rotation.","core_discovery":"On its own terms, the paper's central discovery is that the Coriolis force measurably couples to the full depth of a shallow stellar convection zone. In the fast-rotation case ($P_{\\rm rot}=1$ day), the mean azimuthal flow is slower than the imposed rotation near the surface, by up to 3.1 km/s at the equator at about 3.3 Mm depth, and faster than the imposed rotation in the lower convection zone and overshoot layer; meridional flows reach about 1.35 km/s northward at 60 degrees latitude, with return flows near the base of the convection zone. These shearing flows organize convection into large roll-like structures at the equator that become smaller and shorter-lived toward the poles. Thermodynamically, rotation lowers the mean temperature throughout the convection zone, with photospheric deviations of 7.5 to 9.8 percent relative to the non-rotating model and a poleward temperature decrease of 90 to 202 K for the 1-day period and 220 to 240 K for the 14-day period relative to the equator; the radius decrease of roughly 29 km at the equator and 58 km at higher latitudes is inferred from the downward shift of the photosphere, and the latitude-dependent temperature pattern is attributed to gravity darkening. The slower 14-day rotation reproduces the same qualitative behaviors with weaker amplitudes.","pith_inferences":["A natural extension the paper does not pursue is to repeat the runs in a global spherical geometry that includes centrifugal force; the radius-change and gravity-darkening interpretations would be confirmed only if the global run shows the same sign and magnitude of effects.","The 14-day runs show a larger photosphere temperature contrast relative to the equator (220 to 240 K) than the 1-day runs (90 to 202 K), a non-monotonic trend the paper does not discuss; probing intermediate rotation periods would show whether gravity darkening really weakens at faster rotation.","Synthetic observables, such as intensity maps, limb-darkened profiles, and interferometric visibilities, could be generated from these simulations; comparing them with observations of rapidly rotating F stars would test the gravity-darkening and radius-shift signatures directly.","The simulated tachocline and overshoot in a shallow convection zone imply the ingredients for a thin-shell stellar dynamo; adding magnetic fields to the models would test whether the roll-like flows and meridional circulation can sustain one."],"forward_implications":["For a star rotating near a 1-day period, stellar radius determinations that ignore rotation would be off by tens of kilometers, a shift that matters for precise asteroseismic and eclipsing-binary radii.","The predicted pole-to-equator photosphere temperature contrast of roughly 90 to 200 K for 1-day rotation implies that brightness, color, and spectral-line measurements of inclined fast rotators cannot assume a latitude-independent surface temperature.","Differential rotation and meridional flows can develop in a convection zone only 28.5 Mm thick, so shallow-convection F stars should be treated as capable of sustaining large-scale angular-momentum transport and shear layers, not just small-scale turbulence.","Roll-like convective patterns imply anisotropic heat and momentum transport, so mean-field models of F-star convection zones should include a latitude-dependent roll contribution rather than isotropic turbulent diffusion.","The simulated shift of hydrogen and helium ionization zones under rotation changes the adiabatic gradient profile, which would alter acoustic-mode frequencies and therefore asteroseismic structure inversions."],"supporting_citations":[{"why":"Supplies the non-rotating baseline model of the same 1.47 solar-mass star and the initial structure; all rotation effects in this paper are measured as deviations from it.","marker":"Kitiashvili et al. 2016"},{"why":"Describes the code that performs the first-principles radiative hydrodynamic simulations used here.","marker":"Wray et al. 2018"},{"why":"Establishes the method on solar convection and identifies the leptocline and ionization-zone fine structure used to interpret meridional-flow and ionization-zone correlations.","marker":"Kitiashvili et al. 2023"},{"why":"Provides the classical gravity-darkening relation by which the simulated latitude-dependent photosphere temperature decrease is interpreted.","marker":"von Zeipel 1924"},{"why":"Defines the convective and kinetic energy flux diagnostics used to quantify how rotation suppresses energy transport.","marker":"Nordlund & Stein 2001"},{"why":"Supplies the stellar-evolution initial conditions for the 1.47 solar-mass model and its convection-zone depth.","marker":"Morel & Lebreton 2008"},{"why":"Provides the subgrid-scale turbulence closure used to represent unresolved small-scale dissipation and transport.","marker":"Smagorinsky 1963"},{"why":"Provides the long-characteristics radiative transfer method used for the angle-dependent four-bin radiation transport.","marker":"Feautrier 1964"}],"fun_headline_variants":["Fast spin shrinks F-star radius by 58 km","Rotation slows surface flow, speeds deep flow in F-star","Coriolis forces shape roll-like convection in F-star","Fast rotation restructures convection, shrinks radius","Rotation drives 1.35 km/s meridional flows in F-star"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that flat, horizontally periodic local boxes that neglect surface curvature and the centrifugal force can stand in for the whole star, so that a local photosphere shift is read as a stellar radius decrease and latitude-dependent photosphere cooling as gravity darkening; if curvature or centrifugal effects are essential to those effects, the central claims overstate the model.","fun_headline_variants_meta":{"raw":{"variants":["Fast spin shrinks F-star radius by 58 km","Rotation slows surface flow, speeds deep flow in F-star","Coriolis forces shape roll-like convection in F-star","Fast rotation restructures convection, shrinks radius","Rotation drives 1.35 km/s meridional flows in F-star"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000888,"raw_usage":{"total_tokens":3919,"prompt_tokens":1119,"completion_tokens":2800,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":735,"completion_tokens_details":{"reasoning_tokens":2718}},"tokens_in":735,"tokens_out":2800,"duration_ms":19756,"temperature":1.0,"reasoning_tokens":2718,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T14:04:45.717242+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a global spherical simulation of the same 1.47 solar-mass star at a 1-day rotation period with centrifugal force and curvature included: if the photosphere is not lower by about 29 km at the equator and about 58 km at high latitudes relative to the non-rotating model, or if high latitudes are not cooler than the equator, the radius-decrease and gravity-darkening interpretations fail. A direct observational check would be interferometric or asteroseismic radius and latitudinal brightness mapping of a similar fast-rotating F star.","supporting_citations":[{"cited_title":"A., Bensassy, K., Kitiashvili, I","cited_arxiv_id":null,"evidence_quote":"Describes the code that performs the first-principles radiative hydrodynamic simulations used here."},{"cited_title":"1964, Comptes Rendus Academie des Sciences (serie non specifiee), 258, 3189","cited_arxiv_id":null,"evidence_quote":"Provides the long-characteristics radiative transfer method used for the angle-dependent four-bin radiation transport."}],"review_version":1}