Pith. sign in

REVIEW 3 major objections 5 minor 31 references

3D Radiative MHD Simulations of Starspots II: Large-scale Structure

T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2412.16921 v1 pith:7U5IGNRN submitted 2024-12-22 astro-ph.SR

Tanayveer Singh Bhatia , Mayukh Panja , Robert Cameron , Sami Solanki This is my paper · ORCID
classification astro-ph.SR
keywords starspotsradiativemagnetohydrodynamicspenumbraEvershedflowMdwarfsstellarmagneticactivitylightcurvesMURa
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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).

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 (3)
  1. [Section 2 and Section 5] 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.
  2. [Appendix A] 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.
  3. [Section 4, Eq. (1)] 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.
minor comments (5)
  1. [Abstract] The abstract contains the typo 'the the previous starspot models'; it should read 'the previous starspot models'.
  2. [Section 4] The text says 'vertcal magnetic field' but should say 'vertical magnetic field'.
  3. [Section 4] The text refers to the 'M2V star' when the simulation is for an M0V star; the naming should be consistent with Table 1.
  4. [Section 2] The sentence 'The bin were split in opacity' should be 'The bins were split in opacity'.
  5. [Section 4 and References] The citation 'Smitha et. al., submitted ApJL' should be replaced with a proper reference consistent with the reference list entry 'Smitha et al. 2024'.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: simulation outputs are computed from stated initial/boundary conditions, not fitted to a target result.

full rationale

This paper is a numerical experiment, not a derivation whose conclusions reduce to its inputs. The authors prescribe an initial flux tube with a 3 kG core, a 2 kG ring to promote penumbra formation, and an alpha=1.5 top-boundary horizontal-field enhancement, and then solve the MHD equations. The reported umbral, penumbral, and quiet-star contrasts, thermodynamic stratifications, and Evershed-like flows are time-averaged diagnostics of those solutions, not parameters tuned to reproduce a particular observational data point. The cited self-work (Panja et al. 2020, 2021) supplies the flux-tube setup procedure and the ring idea, but the present paper does not claim to derive those setups as first-principles predictions; they are explicit modeling assumptions. The penumbra-promoting ring and top-boundary enhancement mean that the existence of a penumbra is influenced by construction, but the paper does not relabel that construction as a prediction: it states openly that the ring was added 'to promote the formation of a significant penumbra' and that the boundary enhancement follows Rempel (2012). No equation in the paper equates an output to an input by definition, and no fitted parameter is renamed as a prediction. Accordingly, there is no circular step of the kinds enumerated.

Assumptions & free parameters 6 free parameters · 4 assumptions · 0 invented entities

The simulations depend on several hand-chosen setup parameters (flux-tube field strengths, ring flux, alpha, opacity bin edges, intensity cut-offs, spot radius) and on domain assumptions about the monolithic flux-tube paradigm, the adequacy of 4-bin ODF radiative transfer, and convergence of the setup procedure. No new physical entities are postulated.

free parameters (6)
  • Main flux-tube magnetic field strengths = 3 kG photospheric, 30 kG bottom boundary
    Chosen to create a sunspot-like spot; affects field strength and pressure balance in the spot.
  • Penumbra-promoting ring flux = 2 kG photospheric, 10 kG bottom, thickness about half spot diameter
    Added to promote formation of a significant penumbra; directly influences the penumbral structure that is a central result.
  • Top-boundary horizontal-field enhancement alpha = 1.5
    Sets how much the horizontal field at the top boundary is enhanced relative to a potential field following Rempel (2012); this stabilizes the penumbra and could imprint the Evershed-like flow pattern.
  • Opacity bin edges = log10 tau = -0.5, -1.5, -3
    Bin-splitting for 4-bin ODF radiative transfer; standard but chosen and affects surface cooling and thus intensity contrasts.
  • Umbral/penumbral intensity cut-offs = G2V: 0.35/0.8, K2V: 0.5/0.9, M0V: 0.65/0.965
    Chosen by visual inspection to match filaments; directly determines the quoted umbral and penumbral contrasts and spot radii in Table 1 and Figures.
  • Initial spot radius = 0.3 Lx
    Sets spot size relative to the box; affects boundary effects and Wilson depression normalization.
assumptions (4)
  • domain assumption The monolithic flux-tube paradigm adequately represents real starspots
    Invoked in Section 5 ('all realistic radiative MHD simulations of sun/starspots work in the monolithic tube paradigm'); bundle-of-filaments and alternative geometries are not tested.
  • domain assumption MURaM MHD with 4-bin ODF radiative transfer is an adequate physical approximation for cool-star surface layers
    Radiative transfer is binned (Section 2) rather than full wavelength-resolved; this standard compromise affects the surface temperature and intensity contrasts.
  • domain assumption The spin-up, opacity iteration, and resolution doubling converge to a representative quasi-equilibrium state
    No convergence tests are shown; analysis covers one hour of stellar time (Section 2).
  • domain assumption The quiet-star reference region is unperturbed by the spot
    Appendix A uses a 0.2Lx by 0.2Ly corner square; no test that this region is insensitive to spot-induced flows.

how reviews work

0 comments
Cite this review

Pith. "Pith review of 3D Radiative MHD Simulations of Starspots II: Large-scale Structure." pith.science (2026). https://pith.science/paper/7U5IGNRN

@misc{pith2026241216921,
  author       = {Pith},
  title        = {Pith review of: 3D Radiative MHD Simulations of Starspots II: Large-scale Structure},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7U5IGNRN}},
  note         = {Machine review of arXiv:2412.16921}
}
read the original abstract

We compute realistic 3D radiative MHD near-surface models of starspots with substantial penumbrae on cool main-sequence stars using the MURaM simulation code. This work is an improvement on the the previous starspot models in a slab geometry. The umbra, penumbra and the quiet star for all starspots are distinct, not only in intensity and temperature, but also in thermodynamic and velocity structure. These models represent a significant step towards modelling contribution of starspots to stellar lightcurves.

Figures

Figures reproduced from arXiv: 2412.16921 by the authors.

Figure 1
Figure 1. Snapshots of surface (τ = 1) structure of bolometric intensity (first row), temperature (second row), magnetic field magnitude (third row), field inclination (fourth row), radial velocity (fifth row) and vertical velocity (sixth row) for the G2V (first column), K2V (second column) and M0V (third column) starspot. The rightmost (fourth) column shows the azimuthal average of the corresponding quantities in each row fo… view at source ↗
Figure 2
Figure 2. Time averages of azimuthally averaged relative change in temperature T /Tqs − 1 (first row), relative change in density ρ/ρqs − 1 (second row), vertical velocity (third row), radial velocity (fourth row), and magnetic field inclination (fifth row) for the G2V (first column, K2V (second column) and M0V (third column) starspot. The spot center corresponds to zero on the horizontal axis. The black contour lines near th… view at source ↗
Figure 3
Figure 3. Wilson depression, normalized by the average quiet star pressure scale height Hp at the τ = 1 level, plotted against normalized spot radius. The colors refer the same stellar types as in [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: The vertical velocity vz (top) and the radial velocity (bottom) averaged between ln p/p0 = 6 and 8. The error bars indicate the 1σ standard deviation in the averaging, from the data in [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: Top: Average entropy in the umbral (solid), penumbral (dashed) and quiet star (dotted) region plotted against pressure in units of the pressure scale height ln(p/p0) calculated from the pressure stratification in the quiet star region. Bottom: The temperature stratific…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 21 canonical work pages

  1. [1]

    L., Donahue, R

    Baliunas, S. L., Donahue, R. A., Soon, W. H., et al. 1995, ApJ, 438, 269

    1995

  2. [2]

    Baliunas, S. L. & Vaughan, A. H. 1985, ARA&A, 23, 379

    work page 1985

  3. [3]

    H., Reiners, A., & Schüssler, M

    Beeck, B., Cameron, R. H., Reiners, A., & Schüssler, M. 2013, A&A, 558, A48 Béky, B., Kipping, D. M., & Holman, M. J. 2014, MNRAS, 442, 3686

    work page 2013

  4. [4]

    Berdyugina, S. V. 2005, Living Reviews in Solar Physics, 2, 8 Article number, page 6 Tanayveer Singh Bhatia et al.: 3D Radiative MHD Simulations of Starspots II: Large-scale Structure

    work page 2005

  5. [5]

    S., Cameron, R

    Bhatia, T. S., Cameron, R. H., Solanki, S. K., et al. 2022, A&A, 663, A166

    2022

  6. [6]

    Borrero, J. M. & Ichimoto, K. 2011, Living Reviews in Solar Physics, 8, 4 Castellanos Durán, J. S., Lagg, A., & Solanki, S. K. 2021, A&A, 651, L1

    work page 2011

  7. [7]

    S., Leifer, S., et al

    Crass, J., Gaudi, B. S., Leifer, S., et al. 2021, arXiv e-prints, arXiv:2107.14291

    arXiv 2021

  8. [8]

    Dumusque, X., Boisse, I., & Santos, N. C. 2014, ApJ, 796, 132

    work page 2014

Show all 31 references
  1. [9]

    1909, MNRAS, 69, 454

    Evershed, J. 1909, MNRAS, 69, 454

    1909

  2. [10]

    A., Anglada-Escude, G., Arriagada, P., et al

    Fischer, D. A., Anglada-Escude, G., Arriagada, P., et al. 2016, PASP, 128, 066001

    2016

  3. [11]

    Gough, D. O. & Tayler, R. J. 1966, MNRAS, 133, 85

    1966

  4. [12]

    Hara, N. C. & Ford, E. B. 2023, Annual Review of Statistics and Its Application, 10, 623

    2023

  5. [13]

    2016, A&A, 586, A131

    Herrero, E., Ribas, I., Jordi, C., et al. 2016, A&A, 586, A131

    2016

  6. [14]

    Irwin, A. W. 2012, FreeEOS: Equation of State for stellar interiors cal- culations, Astrophysics Source Code Library, record ascl:1211.002 Jurčák, J., Schmassmann, M., Rempel, M., Bello González, N., &

    2012

  7. [15]

    2020, A&A, 638, A28

    Schlichenmaier, R. 2020, A&A, 638, A28

    2020

  8. [16]

    G., Lendl, M., Cubillos, P

    Juvan, I. G., Lendl, M., Cubillos, P. E., et al. 2018, A&A, 610, A15

    2018

  9. [17]

    Kurucz, R. L. 1993, in Astronomical Society of the Pacific Conference

    1993

  10. [18]

    2013, A&A, 557, A26 O’Neal, D., Neff, J

    Magic, Z., Collet, R., Asplund, M., et al. 2013, A&A, 557, A26 O’Neal, D., Neff, J. E., Saar, S. H., & Cuntz, M. 2004, AJ, 128, 1802

    2013

  11. [19]

    Panja, M., Cameron, R., & Solanki, S. K. 2020, ApJ, 893, 113

    2020

  12. [20]

    H., & Solanki, S

    Panja, M., Cameron, R. H., & Solanki, S. K. 2021, ApJ, 907, 102

    2021

  13. [21]

    V., Espinoza, N., Berdyugina, S

    Rackham, B. V., Espinoza, N., Berdyugina, S. V., et al. 2023, RAS Techniques and Instruments, 2, 148

    2023

  14. [22]

    2014, ApJ, 789, 132

    Rempel, M. 2014, ApJ, 789, 132

    2014

  15. [23]

    2015, ApJ, 814, 125

    Rempel, M. 2015, ApJ, 814, 125

    2015

  16. [24]

    & Schlichenmaier, R

    Rempel, M. & Schlichenmaier, R. 2011, Living Reviews in Solar Physics, 8, 3

    2011

  17. [25]

    Scharmer, G. B. 2009, Space Sci. Rev., 144, 229

    2009

  18. [26]

    2021, A&A, 656, A92

    Schmassmann, M., Rempel, M., Bello González, N., Schlichenmaier, R., & Jurčák, J. 2021, A&A, 656, A92

    2021

  19. [27]

    1958, Structure and evolution of the stars

    Schwarzschild, M. 1958, Structure and evolution of the stars. (Prince- ton University Press)

    1958

  20. [28]

    N., Shapiro, A

    Smitha, H. N., Shapiro, A. I., Witzke, V., et al. 2024, arXiv e-prints, arXiv:2411.14056

    2024 arXiv

  21. [29]

    Solanki, S. K. 2003, A&A Rev., 11, 153

    2003

  22. [30]

    Strassmeier, K. G. 2009, A&A Rev., 17, 251 van Noort, M., Lagg, A., Tiwari, S. K., & Solanki, S. K. 2013, A&A, 557, A24 Vögler, A., Shelyag, S., Schüssler, M., et al. 2005, A&A, 429, 335

    2009

  23. [31]

    I., Kostogryz, N

    Witzke, V., Shapiro, A. I., Kostogryz, N. M., et al. 2024, A&A, 681, A81 Article number, page 7 A&A proofs: manuscript no. main Appendix A: Intensity contours and averaging Fig. A.1. Snapshots of bolometric intensityIbol for G2V (top left), K2V (top right) and M0V (bottom left...

    2024

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.