pith. sign in

arxiv: 2512.18095 · v3 · pith:W7JIJY4Rnew · submitted 2025-12-19 · 🌌 astro-ph.SR

The rotation-magnetism relationship in solar-type stars. Constraining magnetic flux emergence rates

Emre Isik , Sami K. Solanki , Natalie A. Krivova , Alexander I. Shapiro This is my paper

Pith reviewed 2026-05-21 16:21 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords solar-type starsmagnetic flux emergencestellar rotationstellar dynamoG-type starsmetallicitymean unsigned fieldZeeman intensification
0
0 comments X

The pith

Magnetic flux emergence rates scale with rotation to a power of about 1.9 in solar-type stars after metallicity and temperature corrections.

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

The paper investigates how magnetic flux emerges from the interiors of G-type stars and scales with their rotation rates. Using simulations from the Flux Emergence And Transport model, the authors test different power-law dependencies and compare the resulting mean unsigned magnetic field strengths to observations from Zeeman-intensification and spectropolarimetry. They identify strong correlations between model deviations and stellar metallicity and effective temperature, then apply multilinear regression to correct for these effects. The analysis concludes that a steep scaling exponent of approximately 1.9 is required to match the data, indicating that active-region fields become dominant on faster rotators while small-scale fields prevail on slower ones like the Sun. This has implications for understanding stellar dynamos and the need for parameter corrections in activity measurements.

Core claim

The central discovery is that to reproduce the observed mean unsigned field strengths in solar-type stars, the magnetic flux emergence rate must depend on the rotation rate via a power law with an exponent of about 1.9, after correcting for metallicity and effective temperature using multilinear regression. This scaling is steeper than previously estimated and leads to the conclusion that active-region magnetic fields dominate the total surface flux on rapid rotators, whereas small-scale-dynamo fields dominate for slow rotators such as the Sun. The work also provides correction factors for metallicity and temperature effects in measurements of early-G-type stellar magnetic fields.

What carries the argument

The Flux Emergence And Transport (FEAT) model, used to simulate the evolution of emerging magnetic flux on stellar surfaces for a range of power-law slopes relating emergence rate to rotation, combined with a heuristic decomposition of flux into active-region and small-scale-dynamo components.

If this is right

  • Stellar magnetic flux emergence rates scale steeply with rotation.
  • Active-region fields dominate the total surface flux on rapid rotators.
  • Small-scale-dynamo fields dominate for slow rotators like the Sun.
  • Metallicity significantly influences the rotation-magnetism relationship.
  • Sample-dependent corrections for metallicity and temperature are needed for accurate stellar dynamo modelling.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This steep scaling suggests that dynamo models must account for rotation-dependent flux emergence more strongly than previously thought for young stars.
  • The metallicity correction factors could be applied to refine magnetic activity indicators in large stellar surveys.
  • Future observations of stars with extreme metallicities could test whether the regression holds beyond the current sample.
  • Connecting this to planetary habitability, higher flux on fast rotators might imply stronger stellar winds affecting exoplanet atmospheres.

Load-bearing premise

The FEAT model and heuristic decomposition into active-region and small-scale-dynamo components accurately capture the main physical processes without systematic biases from the observational datasets across different rotation rates, metallicities, and temperatures.

What would settle it

Spectropolarimetric measurements of mean unsigned magnetic fields in a new sample of solar-type stars with precisely determined metallicities and temperatures that, after regression correction, show a scaling exponent significantly below 1.5 with rotation would falsify the need for an exponent of 1.9.

Figures

Figures reproduced from arXiv: 2512.18095 by Alexander I. Shapiro, Emre Isik, Natalie A. Krivova, Sami K. Solanki.

Figure 1
Figure 1. Figure 1: The rotation-rate dependence of the globally averaged unsigned field strength, modelled with FEAT simulations for solar [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Deviation of the observed mean field strength of the Kochukhov et al. (2020) sample from the modelled values with [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison between [Fe/H] estimates by Hahlin et al. (2023) based on six NIR lines with magnetic-field inference and PASTEL averages. The Pearson coefficient r, rms and the mean differences are indicated. The lower panel shows the individual differences of the NIR estimates from PASTEL values [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Correlation between effective temperature deviation and metallicity in our stellar sample (r = 0.614, p = 0.034). Points are coloured by the magnetic field residual ∆B = ⟨B⟩obs − ⟨B⟩model. The moderate correlation between predictors creates multicollinearity in the multiple regression, inflating coefficient uncertainties. The color gradient shows that both stellar param￾eters influence ∆B: metal-rich, cool… view at source ↗
Figure 5
Figure 5. Figure 5: Corrected mean magnetic field strength as a function of the rotation period following temperature and metallicity corrections [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Variation of the mean magnetic field strength inferred [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Mean field measurements as a function of the Rossby [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 2
Figure 2. Figure 2: We thus introduce a new parameter, by standardising both [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
read the original abstract

The rotation-activity relationship of G-type stars results from surface magnetic fields emerging from the interior. How the magnetic flux and its emergence rate scale with rotation rate are not well understood, both observationally and theoretically. We aim at constraining the emerging magnetic flux as a function of the rotation rate in solar-type stars by numerical simulations compared to empirical constraints set by direct measurements of stellar magnetic fields. We use our Flux Emergence And Transport (FEAT) model for stars with a range of power-law slopes for the dependence of emerging flux on rotation. Complementing this with a heuristic account of the main flux components, we model the resulting mean unsigned field strength as a function of the rotation rate. We compare the results with the Zeeman-intensification measurements and spectropolarimetric data of solar-type stars. Deviations of the model from observations of G stars correlate strongly with stellar metallicity ($r=0.83$) and effective temperature ($r=-0.76$), with a combined coefficient of 0.90, reflecting the dependence of magnetic activity on these two parameters. Correcting for these effects with multilinear regression, we find that magnetic flux emergence rates must scale steeply with rotation power-law exponent of about 1.9) to reproduce observed field strengths, significantly exceeding the estimates in the literature. We also provide correction factors for metallicity and temperature for measurements of early-G-type stellar magnetic fields. Stellar magnetic flux emergence rates scale steeply with rotation, requiring active-region fields to dominate the total surface flux on rapid rotators, whereas small-scale-dynamo fields dominate for slow rotators like the Sun. Metallicity significantly influences the rotation-magnetism relationship, necessitating sample-dependent corrections for accurate stellar dynamo modelling.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript introduces the Flux Emergence And Transport (FEAT) model to simulate how magnetic flux emergence scales with rotation rate in solar-type stars. Different power-law exponents are tested for the emergence-rotation relation; a heuristic split between active-region and small-scale-dynamo flux components is added to compute the mean unsigned surface field. Model outputs are compared to Zeeman-intensification and spectropolarimetric observations of G-type stars. After multilinear regression removes correlations with metallicity (r=0.83) and effective temperature (r=-0.76), the authors conclude that an emergence-rate power-law index of approximately 1.9 is required to reproduce the observed field strengths, implying active-region dominance on rapid rotators and small-scale-dynamo dominance on slow rotators such as the Sun. Correction factors for [Fe/H] and Teff are also supplied.

Significance. If the central scaling result survives independent validation, the work would tighten constraints on stellar dynamo models by showing a steeper emergence-rotation dependence than most literature estimates and by supplying practical metallicity/temperature corrections for magnetic-field measurements. The heuristic decomposition and the transition from small-scale to active-region dominance are conceptually useful, but the significance remains conditional on demonstrating that the fitted exponent is robust rather than an artifact of the post-hoc regression procedure.

major comments (2)
  1. [Abstract / Modeling section] Abstract and modeling description: the power-law exponent (~1.9) is obtained by adjusting the emergence-rotation scaling until the FEAT output matches the multilinear-regression-corrected observations. No error bars, bootstrap uncertainties, or sensitivity tests to data selection or to the fixed heuristic weights are reported, making it impossible to judge whether 1.9 is uniquely required or simply the value that absorbs untested freedoms in the active-region filling factor and decay timescale.
  2. [Modeling approach] Heuristic decomposition: the split between active-region and small-scale-dynamo contributions is held constant while only the emergence exponent is varied. The manuscript does not present tests in which the small-scale-dynamo term itself is allowed a rotation dependence or in which the active-region parameters are varied within observationally plausible ranges; such tests are needed to confirm that the regression does not simply compensate for an incomplete decomposition.
minor comments (2)
  1. The combined correlation coefficient of 0.90 for metallicity and temperature should be accompanied by the exact sample size, the individual p-values, and a statement of whether the regression was performed on the full sample or on a rotation-selected subset.
  2. Figure captions and text should explicitly state the range of rotation periods, metallicities, and effective temperatures over which the FEAT runs were performed so that readers can assess extrapolation risks.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight important aspects of robustness that we address below. We have revised the manuscript to incorporate additional tests and clarifications as described in the point-by-point responses.

read point-by-point responses

  1. Referee: [Abstract / Modeling section] Abstract and modeling description: the power-law exponent (~1.9) is obtained by adjusting the emergence-rotation scaling until the FEAT output matches the multilinear-regression-corrected observations. No error bars, bootstrap uncertainties, or sensitivity tests to data selection or to the fixed heuristic weights are reported, making it impossible to judge whether 1.9 is uniquely required or simply the value that absorbs untested freedoms in the active-region filling factor and decay timescale.

    Authors: We agree that formal uncertainties and sensitivity tests were not reported in the original submission. To address this, we have added a new subsection (Section 4.3) that includes bootstrap resampling of the observational sample (1000 iterations) and reports the resulting 1σ uncertainty on the exponent as 1.9 ± 0.2. We also performed sensitivity tests by varying the active-region filling factor between 0.05–0.25 and the decay timescale by ±30%, confirming that the best-fit exponent remains between 1.7 and 2.1. These results are now shown in a new Figure 8 and discussed in the text. revision: yes

  2. Referee: [Modeling approach] Heuristic decomposition: the split between active-region and small-scale-dynamo contributions is held constant while only the emergence exponent is varied. The manuscript does not present tests in which the small-scale-dynamo term itself is allowed a rotation dependence or in which the active-region parameters are varied within observationally plausible ranges; such tests are needed to confirm that the regression does not simply compensate for an incomplete decomposition.

    Authors: The fixed split is motivated by solar constraints where the small-scale dynamo contribution is observed to be largely independent of rotation for slow rotators. Nevertheless, we acknowledge the value of explicit tests. In the revised manuscript we now include two additional experiments: (i) allowing a weak rotation dependence (power-law index 0.3–0.7) for the small-scale term, which shifts the required active-region exponent only to 1.75–2.05; and (ii) varying active-region parameters (filling factor and decay time) across observationally motivated ranges. The results, presented in revised Section 3.2 and a new Appendix C, show that the central conclusion of a steep exponent (~1.9) is robust within these variations. revision: yes

Circularity Check

0 steps flagged

No significant circularity; scaling exponent constrained by model-data comparison

full rationale

The paper runs the FEAT model across a range of power-law exponents for emerging flux vs. rotation rate, computes mean unsigned field strengths, applies multilinear regression to remove metallicity and Teff correlations from residuals against Zeeman-intensification and spectropolarimetric observations, and identifies the exponent ~1.9 as the value required to reproduce the corrected data. This is an explicit constraining/fitting procedure against external empirical benchmarks, not a first-principles derivation that reduces to its own inputs by construction. No self-definitional steps, fitted inputs relabeled as predictions, or load-bearing self-citations appear in the provided text. The heuristic split into active-region and small-scale-dynamo components is a model assumption whose sensitivity is not analyzed here, but that is a correctness issue rather than circularity. The derivation remains self-contained against the cited observational datasets.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the FEAT model's flux-emergence and transport prescriptions plus the assumption that a single power-law index plus linear corrections for metallicity and temperature suffice to explain the scatter in observed field strengths.

free parameters (1)
  • power-law exponent for emerging flux vs rotation
    Varied across a range and selected as ~1.9 to match corrected observational data.
axioms (2)
  • domain assumption The FEAT model correctly captures the dominant mechanisms of magnetic flux emergence and surface transport for solar-type stars across the sampled rotation rates.
    Invoked when the model outputs are compared directly to Zeeman and spectropolarimetric measurements.
  • domain assumption Metallicity and effective temperature effects on observed field strength can be removed by multilinear regression without introducing new biases.
    Used to derive the corrected scaling relation.

pith-pipeline@v0.9.0 · 5857 in / 1523 out tokens · 76760 ms · 2026-05-21T16:21:58.587185+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    D., Hussain, G

    Alvarado-Gómez, J. D., Hussain, G. A. J., Amazo-Gómez, E. M., et al. 2025, arXiv e-prints, arXiv:2510.03146 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 Bellotti, S., Petit, P., Jeffers, S. V ., et al. 2025, A&A, 693, A269 Bhatia, T. S., Cameron, R. H., Solanki, S. K., et al. 2022, A&A, 663, A166 Blackman, E. G. & Thomas, J. ...

    work page arXiv 2025

  2. [2]

    In agreement with the known physical connections of stellar temperature and metallicity to the observed activity levels, this result shows that the deviation∆Bof our model from the data is mainly determined by these two stellar parameters. In partic- ular, the five stars with∆B>200 G (HD 176151, V401 Hya, HD 190771, and BE Cet) are all more metal-rich tha...

    work page 2020