{"id":"8edd4926-cf6f-40e1-a9b1-63a713095552","arxiv_id":"2412.16921","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"3D MURaM simulations of starspots on G2V, K2V, and M0V stars reveal distinct umbral, penumbral, and quiet-region structures with spectral-type-dependent contrasts and Evershed-like flows.","lead":"This paper presents the first 3D radiative magnetohydrodynamic simulations of starspots with penumbrae on cool main-sequence stars, covering G2V, K2V, and M0V types. The simulations show that spot umbra, penumbra, and quiet star regions differ in intensity, temperature, flow, and magnetic structure, providing a more realistic basis for modeling stellar light curves and exoplanet observations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The penumbra and Evershed flow may be imposed by the alpha=1.5 top boundary condition and the 2 kG penumbra-promoting ring, not self-organized; no sensitivity test is reported.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the penumbra-promoting ring and the alpha=1.5 horizontal-field enhancement are imposed in Section 2 and could be responsible for the penumbral filaments and radial outflows. I agree with that assessment. The strongest claim is not merely geometric; it is a claim about realistic starspot structure, and realism requires that the penumbral structure and flows persist when the penumbra-promoting inputs are removed or varied. The paper contains no such sensitivity study, and the authors themselves describe the models as exploratory in Section 5. I do not count the absence of code or data as the primary issue, since MURaM sunspot simulations are an established method; the load-bearing issue is specifically the unquantified imprint of the boundary and initialization choices. The proposed alpha/ring variation is computationally feasible and directly tests whether the central claim survives. If that test showed strong dependence, the verdict would need to move toward REJECT for the stated strongest claim; without the test, the reader's CONDITIONAL verdict is exactly right. I therefore set verdict_should_be to UNCHANGED, meaning my read does not move the reader's conditional verdict.","tokens_in":9323,"tokens_out":4636,"duration_ms":44511,"concrete_test":"For the G2V case, rerun with alpha=1.0 (potential-field top boundary) and with the penumbra-promoting ring removed, holding all other parameters fixed. Measure the azimuthally averaged penumbral intensity contrast, penumbral area, field inclination, and radial velocity at tau=1 over the same 1-hour interval. If removing the ring and alpha reduces penumbral area or Evershed flow to noise, or changes the intensity contrast by more than about 20%, the distinct large-scale structure is largely imposed by the setup; if the profiles are statistically unchanged, the self-organized interpretation is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the simulations produce realistic, distinct umbral, penumbral, and quiet-star regions with Evershed-like flows—requires that these structures emerge from convective MHD self-organization. The most load-bearing assumption is that the imposed initial/boundary conditions in Section 2 do not directly manufacture the penumbra. The setup places a 2 kG photospheric ring around the 3 kG flux tube, with ring thickness about half the tube diameter, explicitly 'to promote the formation of a significant penumbra,' and enhances the horizontal field at the top boundary via alpha=1.5 following Rempel (2012). Penumbral filaments, inclined fields, and radial outflows are then diagnosed exactly in the region where this ring and boundary enhancement act. If these imposed conditions are controlling, the abstract's claim that umbra, penumbra, and quiet star are distinct in thermodynamic and velocity structure is a property of the chosen boundary-value problem rather than a robust stellar result. The paper reports no runs varying alpha, ring flux, resolution, or box size, and Section 5 labels the models 'somewhat exploratory.' Solar MURaM sunspot simulations use a similar boundary treatment, but for the Sun independent observations anchor the calibration; for K2V/M0V the same tuning is unvalidated. The penumbral contrasts and Evershed amplitudes in Table 1 and Figure 1 should therefore be treated as provisional until dependence on these imposed parameters is quantified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 3D radiative MHD simulations of starspots on G2V, K2V, and M0V stars using the MURaM code. The spots are initialized as monolithic flux tubes with an added 2 kG penumbra-promoting ring and an alpha=1.5 top-boundary horizontal-field enhancement. The authors report that all three spots develop a distinct umbra, penumbra, and quiet-star region, with umbral and penumbral intensity contrasts decreasing toward cooler stars and a persistent radially outward Evershed-like flow in the penumbra. Thermodynamic and velocity stratifications are compared across spectral types. The paper is an extension of the slab-geometry models of Panja et al. (2020).","tokens_in":9596,"tokens_out":5943,"duration_ms":51871,"significance":"If the results are robust, the paper provides a valuable first look at 3D starspot structure across spectral types, with implications for starspot contrasts in light curves and RV signals. The use of the established MURaM code, the detailed description of the numerical setup, and the explicit comparison with the slab-geometry predecessor are strengths. However, the quantitative conclusions rest on imposed boundary conditions and subjective intensity cut-offs, so the current results are best viewed as exploratory.","major_comments":[{"comment":"The penumbral structure and Evershed-like flow may be imposed by the initial and boundary conditions, specifically the 2 kG penumbra-promoting ring and the alpha=1.5 top-boundary horizontal-field enhancement. No sensitivity runs varying alpha, ring flux, or ring thickness are reported, and Section 5 itself labels the models 'somewhat exploratory.' Because the abstract's central claim is that the three regions are distinct in thermodynamic and velocity structure, the authors should either add such tests or clearly qualify the results as dependent on the imposed parameters.","section":"Section 2 and Section 5"},{"comment":"The umbral and penumbral boundaries are defined by intensity cut-offs (Iu/Iqs and Ip/Iqs) chosen from histograms and visual inspection. This makes the contrast ratios in Table 1 and the normalized radii used throughout Figs. 1-4 dependent on subjective choices. A quantitative, reproducible boundary definition (e.g., based on magnetic field inclination or on a fixed intensity threshold) would strengthen the paper.","section":"Appendix A"},{"comment":"Equation (1) states that the requirement for convection is nabla - nabla_ad - B_z^2/(B_z^2+4*pi*Gamma*p) < 0. Setting B_z=0 then implies nabla - nabla_ad < 0, which is the Schwarzschild stability condition, not the convection condition. The inequality sign appears reversed; the correct Gough-Tayler criterion should require nabla - nabla_ad greater than the magnetic term. The discussion of umbral stability in Section 4 relies on this equation, so it needs to be corrected.","section":"Section 4, Eq. (1)"}],"minor_comments":[{"comment":"The abstract contains the typo 'the the previous starspot models'; it should read 'the previous starspot models'.","section":"Abstract"},{"comment":"The text says 'vertcal magnetic field' but should say 'vertical magnetic field'.","section":"Section 4"},{"comment":"The text refers to the 'M2V star' when the simulation is for an M0V star; the naming should be consistent with Table 1.","section":"Section 4"},{"comment":"The sentence 'The bin were split in opacity' should be 'The bins were split in opacity'.","section":"Section 2"},{"comment":"The citation 'Smitha et. al., submitted ApJL' should be replaced with a proper reference consistent with the reference list entry 'Smitha et al. 2024'.","section":"Section 4 and References"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of A&A. The main revision should focus on establishing how much of the penumbral structure and Evershed-like flow is a consequence of the imposed ring flux and alpha=1.5 boundary condition; a sensitivity test would substantially strengthen the conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this is the first paper to show 3D radiative MHD starspots in circular geometry for cool main-sequence stars, with a genuine penumbra and Evershed-like outflow. It extends the earlier slab models of Panja et al. to G2V, K2V, and M0V spectral types, and the qualitative trends (lower contrast, weaker flows for cooler stars) match the slab results and the solar MURaM sunspot simulations. That is a real step forward, and it gives exoplanet people concrete structures to test against transit and RV observations.\n\nWhat the paper does well: the setup is described in enough detail to reproduce, the MURaM code is established, and the analysis is thorough—thermodynamic stratification, Wilson depression, convective ring patterns, velocity structure. The authors are also honest: they call the models 'somewhat exploratory' and note the monolithic flux-tube paradigm. The appendix admits the umbral/penumbral boundaries are chosen from histograms and visual inspection. That sort of transparency deserves credit.\n\nThe soft spots are real but not fatal. The 2 kG ring and the alpha=1.5 top-boundary field enhancement are explicitly there 'to promote the formation of a significant penumbra.' The penumbral filaments and radial outflows are then diagnosed in exactly that region, and there are no sensitivity runs varying alpha, ring flux, resolution, or box size. So the quantitative contrasts and flow amplitudes in Table 1 should be treated as provisional until the dependence on those imposed parameters is quantified. The stress-test concern lands on reading the paper; it is a legitimate caveat, not a takedown. For solar spots, observations anchor the calibration; for K2V/M0V stars, there is no equivalent anchor, so the self-organization claim is weaker. The intensity-threshold segmentation is also subjective, though it likely does not change the qualitative picture.\n\nWho this is for: stellar activity modelers, exoplanet RV/transit people, and anyone using starspot models. It is a solid simulation paper that deserves referee time, not desk rejection. My recommendation: send it to A&A or ApJ, and ask the authors to add at least one perturbation test (e.g., lower alpha or no ring) and to quantify how much the penumbral structure depends on the boundary conditions. Even without that, the paper is worth publishing—just with the caveats clearly stated.","headline":"First 3D circular-geometry MURaM starspot models with real penumbrae for G2V/K2V/M0V stars; worth a serious referee, but the penumbra is partly nursed by hand-tuned boundary conditions and no sensitivity runs exist.","tokens_in":10178,"tokens_out":1411,"would_cite":true,"duration_ms":14226,"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":"Three-dimensional radiative magnetohydrodynamic simulations produce starspots with distinct umbra, penumbra, and quiet-star regions on G2V, K2V, and M0V stars, with contrast and flow speeds decreasing toward cooler stars.","keywords":["starspots","radiative magnetohydrodynamics","penumbra","Evershed flow","M dwarfs","stellar magnetic activity","light curves","MURaM"],"falsifier":"Run the same three simulations with the top-boundary horizontal field set to a potential-field extrapolation (alpha=0) and without the penumbral flux ring; if the penumbral filaments and the persistent outward radial flow disappear or change character, the structure is boundary-forced rather than self-organized. A complementary check is to compare the simulated umbral and penumbral intensity contrasts for the K2V and M0V stars against contrasts measured from high-precision photometry or spectroscopy of spotted late-type stars.","tokens_in":9092,"feed_emoji":"🌟","tokens_out":7428,"duration_ms":63200,"temperature":0.7,"pith_summary":"Using three-dimensional radiative magnetohydrodynamic simulations, this paper builds self-consistent models of starspots on a Sun-like G2V star, a K2V star, and an M0V star, each showing a dark umbra, a filamentary penumbra, and a distinct quiet-star region. It argues that these three components differ not only in brightness and temperature but also in their thermodynamic and velocity structure, and that the differences become smaller as the star gets cooler. The significance, if the models are right, is that the widespread shortcut of treating a starspot as a uniform dark circle with a single effective temperature is systematically inadequate, and that realistic spot structure matters for interpreting exoplanet transit spectra, light curves, and radial-velocity measurements.","feed_headline":"Cooler stars have fainter, slower starspots in 3D models","feed_subtitle":"Simulations of spots on Sun-like, K, and M stars show distinct umbra, penumbra, and quiet-star structure.","key_machinery":"The machinery is a set of three-dimensional radiative magnetohydrodynamic simulations run with the MURaM code, in which a monolithic vertical flux tube with a 3 kG photospheric field is inserted into a relaxed stellar convection zone, ringed by an extra 2 kG flux annulus intended to promote penumbra formation, and anchored at the top boundary by an enhanced horizontal field set through a parameter alpha=1.5. The penumbral filaments, inclined fields, and outward flows emerge from convection in strongly inclined magnetic fields, while the near-surface suppression of convection in the umbra is interpreted through the Gough-Tayler criterion, which modifies the Schwarzschild criterion when a strong vertical magnetic field is present. This setup yields stable, long-lived spots whose structure is then azimuthally and temporally averaged to extract the large-scale radial profiles.","core_discovery":"The paper's central claim is that, for the first time, three-dimensional radiative magnetohydrodynamic models of spots on cool main-sequence stars reproduce the essential sunspot phenomenology outside the Sun: a dark umbra, a penumbra made of bright and dark filaments with inclined magnetic fields, and a persistent radially outward flow resembling the Evershed effect. Across the G2V, K2V, and M0V cases, the umbral and penumbral intensity contrasts relative to the quiet star decrease with decreasing effective temperature, while the radial profiles of magnetic field strength and inclination collapse onto nearly the same shape when scaled by spot radius. The temperature and density stratification inside the umbral trunk departs strongly from the quiet-star stratification, signaling magnetic suppression of convection, and the mean flows include ring-like convective patterns whose strength also decreases toward cooler stars. The authors conclude that spot structure, not just spot brightness, has to be part of realistic models of stellar variability and radial-velocity signals.","pith_inferences":["A natural testable extension is to scale spot size to the stellar radius rather than the pressure scale height; because pressure scale height grows more slowly than radius, larger, more realistic spots on M dwarfs might reverse or modify the contrast trends found here.","The Gough-Tayler interpretation implies that the depth at which umbral convection shuts off depends on field strength and local gas pressure; one could search for a systematic relation between umbral brightness and spot size or field strength, as observed on the Sun.","The ringed convective pattern around the spot, with the extra inner flow pair most evident for the M0V star, suggests that the spot's thermal 'moat' changes with spectral type; time-resolved high-precision photometry of spotted M dwarfs might reveal brightness rings or moat-related variability that the current models predict on small scales."],"forward_implications":["The intensity and temperature contrast between spot and quiet star falls with decreasing effective temperature, so a single contrast recipe cannot be applied to all spectral types.","The umbral, penumbral, and quiet-star stratifications are thermodynamically distinct, meaning 1D radiative-equilibrium spot models misrepresent the structure underlying spot spectra.","Spot flows are fast enough to contribute to radial-velocity signals depending on the spot's position on the stellar disk, with hotter stars showing larger flow speeds.","The models provide realistic spot contrasts and limb-darkening behavior that can replace uniform-circle or blackbody approximations in light-curve and transit-spectrum analyses."],"supporting_citations":[{"why":"Supplies the slab-geometry starspot models and flux-tube initialization procedure that this work extends to three-dimensional circular geometry.","marker":"Panja et al. (2020)"},{"why":"Introduces the penumbra-promoting ring of magnetic flux used in the initial condition.","marker":"Panja et al. (2021)"},{"why":"Provides the alpha-enhanced horizontal field at the top boundary that stabilizes the simulated penumbra.","marker":"Rempel (2012)"},{"why":"Gives the sunspot simulation whose convective flow structure and downflows are compared with the present results.","marker":"Rempel (2015)"},{"why":"Documents the version of the MURaM code used for the radiative MHD computations.","marker":"Rempel (2014)"},{"why":"Describes the original MURaM code on which the simulations are based.","marker":"Vögler et al. (2005)"},{"why":"Supplies the criterion for convection in the presence of a vertical magnetic field used to interpret the umbral stratification.","marker":"Gough & Tayler (1966)"},{"why":"Demonstrates distinct convective regimes in umbra, penumbra, and quiet sun using the same criterion, which this paper applies to stellar spots.","marker":"Schmassmann et al. (2021)"}],"fun_headline_variants":["3D starspot models reveal sunspot-like structure on cool stars","Simulated starspots on cool stars show umbra, penumbra, and Evershed flow","3D MHD simulations capture spot structure on G, K, M stars","First 3D starspot models show penumbrae on cool main-sequence stars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the simulated penumbral filaments, inclined fields, and radial outflow arise self-consistently from convection in inclined magnetic fields, rather than being forced by the imposed initial flux-tube ring or by the alpha=1.5 horizontal-field enhancement at the top boundary.","fun_headline_variants_meta":{"raw":{"variants":["3D starspot models reveal sunspot-like structure on cool stars","Simulated starspots on cool stars show umbra, penumbra, and Evershed flow","3D MHD simulations capture spot structure on G, K, M stars","First 3D starspot models show penumbrae on cool main-sequence stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001038,"raw_usage":{"total_tokens":4308,"prompt_tokens":824,"completion_tokens":3484,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":440,"completion_tokens_details":{"reasoning_tokens":3397}},"tokens_in":440,"tokens_out":3484,"duration_ms":21162,"temperature":1.0,"reasoning_tokens":3397,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:58:01.432844+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same three simulations with the top-boundary horizontal field set to a potential-field extrapolation (alpha=0) and without the penumbral flux ring; if the penumbral filaments and the persistent outward radial flow disappear or change character, the structure is boundary-forced rather than self-organized. A complementary check is to compare the simulated umbral and penumbral intensity contrasts for the K2V and M0V stars against contrasts measured from high-precision photometry or spectroscopy of spotted late-type stars.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the slab-geometry starspot models and flux-tube initialization procedure that this work extends to three-dimensional circular geometry."},{"cited_title":"H., & Solanki, S","cited_arxiv_id":null,"evidence_quote":"Introduces the penumbra-promoting ring of magnetic flux used in the initial condition."},{"cited_title":"2015, ApJ, 814, 125","cited_arxiv_id":null,"evidence_quote":"Gives the sunspot simulation whose convective flow structure and downflows are compared with the present results."},{"cited_title":"2021, A&A, 656, A92","cited_arxiv_id":null,"evidence_quote":"Demonstrates distinct convective regimes in umbra, penumbra, and quiet sun using the same criterion, which this paper applies to stellar spots."}],"review_version":1}