pith. sign in

arxiv: 2511.08481 · v2 · submitted 2025-11-11 · 🌌 astro-ph.SR

Stellar cycle variability in Mount Wilson stars and dynamo models: Rotation rate and dynamo number dependency

Suyog Garg , Bidya Binay Karak , Rohan B. Mandrai This is my paper

Pith reviewed 2026-05-17 23:25 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords stellar magnetic cyclescycle variabilityrotation ratedynamo numberRossby numberdynamo modelsMount Wilson stars
0
0 comments X

The pith

The variability of stellar magnetic cycles decreases with increasing rotation rate and dynamo number.

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

The paper analyzes magnetic cycle data from 81 stars observed at Mount Wilson to determine how cycle variability depends on rotation. It reports a modest positive correlation between variability and rotation period, and negative correlations with the inverse squared Rossby number and the squared logarithm of the cycle-to-rotation period ratio. Dynamo models show the same pattern of decreasing variability with faster rotation or larger dynamo number. The work also identifies better scaling relations for the dynamo number from the models.

Core claim

Our results demonstrate that the stellar magnetic cycle variability decreases with the increase of the rotation rate or the dynamo number. Analysis of Mount Wilson data reveals correlations with rotation period, Rossby number inverse square, and cycle duration ratio. Dynamo models confirm these trends and suggest revised exponents for dynamo number measures.

What carries the argument

Stellar magnetic cycle variability measured from time series, related to rotation period, Rossby number, and mean cycle period over rotation period.

If this is right

  • Variability decreases as rotation rate increases.
  • Variability decreases as dynamo number increases.
  • Dynamo models exhibit similar variability trends to observations.
  • Ro to the power -0.6 and cycle-to-rotation ratio to the power 0.6 serve as good dynamo number measures.

Where Pith is reading between the lines

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

  • These trends could help predict cycle behavior in stars of different ages.
  • More complete long-term observations might confirm or refine the correlations found.
  • The results suggest a general property of stellar dynamos where faster rotation leads to more stable cycles.
  • Applications may extend to assessing magnetic activity on stars with planets.

Load-bearing premise

The limited time series data accurately reflect the intrinsic cycle variability without significant biases from short observation periods.

What would settle it

A clear example of a rapidly rotating star with high cycle variability in extended monitoring would challenge the main finding.

Figures

Figures reproduced from arXiv: 2511.08481 by Bidya Binay Karak, Rohan B. Mandrai, Suyog Garg.

Figure 1
Figure 1. Figure 1: Time series plots of MWO raw S-index data and the respective quasi-Planck fittings for the seven additional GK25 stars considered in this study. Two stars, HD 81809 and HD 166620, that have relatively well defined cycles, and which were studied by Garg et al. (2019), are also shown for cycle comparison. Properties of the stars are tabulated in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Scatter plot between the parametrized variability (Equation (2)) and the inverse squared Rossby number for the GK19 set of stars. The Spearman rank correlation r is also provided. in the activity can be quantified in terms of the quantity var(R), as defined in Section 2. We now want to investigate how this stellar variability changes with the rotation rate of the stars. It is expected that with increasing … view at source ↗
Figure 4
Figure 4. Figure 4: Scatter plot between the parametrized variability and the rotational period for the B95+ set of stars. The variability class data is obtained from [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: Stellar cycle variability as obtained from three different dynamo models as a function of the stellar rotation period. The dynamo models are discussed in detail in Section 3, while the variability is defined in Equation (3). stars of similar type with the same mass of 1M⊙, and thus our stellar model data represents only a small sample of the observed stars, which have varying structures. Nevertheless, our … view at source ↗
read the original abstract

Similar to the solar cycle, the magnetic cycles of other solar-type stars are also variable. How the variability of the stellar cycle changes with the rotation rate or the dynamo number is a valuable information for understanding the stellar dynamo process. We examine the variability in the stellar magnetic cycles by studying 81 stars from the data of the Mount Wilson Observatory, which started observations in 1966. For 28 stars, we have time series data available till 2003, while for others, the data are limited till 1995. We specifically explore how the variability changes with respect to three rotation-related parameters. We find a modest positive correlation between the variability and the stellar rotation period. In addition, we find suggestive negative correlations between the variability and the inverse squared Rossby number ($Ro^{-2}$), and the ratio of the mean cycle duration and rotation period ($\log \, (\langle P_{\rm cyc} \rangle / P_{\rm rot})^2$). Variability computed from the magnetic field of stellar dynamo models also show similar trends. Finally, inspired by previous studies, we examine dynamo number scaling in our model data and find that $Ro^{-0.6}$ (instead of $Ro^{-2}$ as suggested in the linear $\alpha \Omega$ dynamo theory) and $(\langle P_{\rm cyc} \rangle /P_{\rm rot})^{0.6}$ (instead of $\log \, (\langle P_{\rm cyc} \rangle / P_{\rm rot})^2$ as predicted in previous observations) are a good measure of the dynamo number. In conclusion, our results demonstrate that the stellar magnetic cycle variability decreases with the increase of the rotation rate or the dynamo number.

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

3 major / 2 minor

Summary. The manuscript analyzes S-index time series from the Mount Wilson Observatory for 81 solar-type stars to study how magnetic cycle variability depends on rotation rate and dynamo number. It reports a modest positive correlation between variability and rotation period, suggestive negative correlations with Ro^{-2} and log((<P_cyc>/P_rot)^2), and similar trends in dynamo models. The authors propose adjusted scalings Ro^{-0.6} and (<P_cyc>/P_rot)^{0.6} as improved dynamo number measures and conclude that cycle variability decreases with increasing rotation rate or dynamo number.

Significance. If the trends prove robust, the work supplies observational and modeling constraints on how stellar cycle variability scales with rotation and dynamo parameters, aiding refinement of stellar dynamo theory. The dual use of Mount Wilson data and numerical models is a positive feature, though the modest correlation strengths and data limitations reduce the overall impact.

major comments (3)
  1. [Observational data and variability measurement] For 53 of the 81 stars the time series terminate in 1995, yielding at most ~29 years of coverage. With typical cycle lengths of 5–20 yr this supplies only 1–3 cycles per star, rendering robust quantification of cycle-to-cycle variability (amplitude scatter, period jitter, envelope modulation) unreliable; any excess scatter or bias in this subsample directly weakens the reported correlations with rotation period and dynamo proxies.
  2. [Dynamo number scaling analysis] The adjusted exponents -0.6 for Ro and 0.6 for (P_cyc/P_rot) are obtained by fitting the authors’ own dynamo-model outputs rather than independent theoretical predictions or external benchmarks. This renders the claim that these scalings constitute a “good measure” of the dynamo number circular, as the functional form is tuned to reproduce the very trends being demonstrated.
  3. [Statistical results and correlation analysis] The correlations are characterized only as “modest” and “suggestive,” yet the manuscript provides no detailed error propagation, bootstrap or Monte-Carlo significance tests, or comparison against null models that account for heterogeneous baseline lengths and possible secular trends or calibration drifts.
minor comments (2)
  1. [Abstract and §1] The abstract and introduction should explicitly state the fraction of stars with data extending only to 1995 versus 2003 and discuss the implications for cycle coverage.
  2. [Methods and notation] Notation such as <P_cyc> and the precise definition of cycle variability (e.g., standard deviation of amplitudes or periods) should be introduced earlier and used consistently.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thoughtful and constructive comments, which help us clarify limitations and strengthen the statistical presentation of our results. We respond to each major comment below, indicating planned revisions where appropriate.

read point-by-point responses

  1. Referee: For 53 of the 81 stars the time series terminate in 1995, yielding at most ~29 years of coverage. With typical cycle lengths of 5–20 yr this supplies only 1–3 cycles per star, rendering robust quantification of cycle-to-cycle variability (amplitude scatter, period jitter, envelope modulation) unreliable; any excess scatter or bias in this subsample directly weakens the reported correlations with rotation period and dynamo proxies.

    Authors: We agree that the limited number of observed cycles for the 53 stars with data ending in 1995 represents a genuine limitation for quantifying cycle-to-cycle variability. The manuscript already notes the differing baseline lengths. To address the concern, we will add an explicit subsample analysis restricted to the 28 stars with coverage through 2003 and verify that the reported correlations with rotation period and dynamo proxies remain qualitatively unchanged (though with larger uncertainties). We will also expand the discussion section to quantify the typical number of cycles per star and discuss how this affects the robustness of the variability measures. revision: yes

  2. Referee: The adjusted exponents -0.6 for Ro and 0.6 for (P_cyc/P_rot) are obtained by fitting the authors’ own dynamo-model outputs rather than independent theoretical predictions or external benchmarks. This renders the claim that these scalings constitute a “good measure” of the dynamo number circular, as the functional form is tuned to reproduce the very trends being demonstrated.

    Authors: The referee is correct that the specific exponents were determined by optimizing the correlation strength within our own set of dynamo models. We will revise the relevant section to make clear that these values are empirical adjustments motivated by the mismatch between linear αΩ theory and both our models and the Mount Wilson trends, rather than new first-principles predictions. The manuscript will now frame them as practical proxies that improve the observed correlation with cycle variability, while acknowledging the data-driven nature of the choice. revision: yes

  3. Referee: The correlations are characterized only as “modest” and “suggestive,” yet the manuscript provides no detailed error propagation, bootstrap or Monte-Carlo significance tests, or comparison against null models that account for heterogeneous baseline lengths and possible secular trends or calibration drifts.

    Authors: We accept that the current statistical analysis is insufficiently rigorous. In the revised manuscript we will add bootstrap resampling (with 10,000 iterations) to obtain confidence intervals and p-values for all reported correlations. We will also construct null models that preserve the heterogeneous observation lengths and test against randomized cycle amplitudes to assess the influence of baseline differences and potential secular drifts. These additions will be described in a new subsection on statistical methods. revision: yes

Circularity Check

1 steps flagged

Dynamo number scalings fitted to own model data reduce central claim to data-driven adjustment

specific steps
  1. fitted input called prediction [Abstract (final paragraph before conclusion)]
    "Finally, inspired by previous studies, we examine dynamo number scaling in our model data and find that Ro^{-0.6} (instead of Ro^{-2} as suggested in the linear αΩ dynamo theory) and (⟨P_cyc⟩ /P_rot)^{0.6} (instead of log (⟨P_cyc⟩ / P_rot)^2 as predicted in previous observations) are a good measure of the dynamo number."

    The exponents -0.6 and 0.6 are obtained by inspecting trends within the authors' own stellar dynamo model runs. This fitted scaling is then adopted as the operative definition of dynamo number to interpret both the model and observational variability results, so the reported decrease of variability with dynamo number is partly a restatement of the fit performed on the same data.

full rationale

The paper reports observational correlations of cycle variability with rotation period and two initial dynamo proxies (Ro^{-2} and log(P_cyc/P_rot)^2). It then uses its own dynamo models to calibrate new exponents (Ro^{-0.6} and (P_cyc/P_rot)^{0.6}) and presents these as improved measures of the dynamo number. Because the calibration is performed on the same model outputs that are later shown to exhibit the variability trend, the support for the headline claim that variability decreases with dynamo number is partly enforced by the fitting step rather than arising independently. The observational data and basic rotation-period correlation remain external, so the circularity is partial rather than total.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard dynamo theory assumptions plus data-derived scalings. No new entities are postulated. The adjusted exponents appear fitted rather than derived from first principles.

free parameters (2)
  • exponent -0.6 for Ro in dynamo number
    Chosen to fit the model data better than the linear theory value of -2.
  • exponent 0.6 for (P_cyc/P_rot) in dynamo number
    Fitted to model data instead of the previously suggested log-squared form.
axioms (2)
  • domain assumption Variability metric computed from Mount Wilson time series accurately captures intrinsic dynamo fluctuations.
    Invoked when correlating variability with rotation parameters without detailed discussion of measurement uncertainties or cycle coverage.
  • domain assumption Dynamo models reproduce the relevant stellar regimes.
    Used to confirm observational trends without specifying model parameters or validation against independent data.

pith-pipeline@v0.9.0 · 5624 in / 1463 out tokens · 25610 ms · 2026-05-17T23:25:33.375429+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    L., Donahue, R

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

    work page 1995

  2. [2]

    Barnes, S. A. 2007, ApJ, 669, 1167

    work page 2007

  3. [3]

    C., Wright, J

    Baum, A. C., Wright, J. T., Luhn, J. K., & Isaacson, H. 2022, AJ, 163, 183

    work page 2022

  4. [4]

    B., Usoskin, I., & Weisshaar, E

    Biswas, A., Karak, B. B., Usoskin, I., & Weisshaar, E. 2023, SSRv, 219, 19 B¨ohm-Vitense, E. 2007, ApJ, 657, 486 Boro Saikia, S., Marvin, C. J., Jeffers, S. V ., et al. 2018, A&A, 616, A108

    work page 2023

  5. [5]

    2020, Living Reviews in Solar Physics, 17, 4

    Charbonneau, P. 2020, Living Reviews in Solar Physics, 17, 4

    work page 2020

  6. [6]

    2005, ApJ, 619, 613

    Charbonneau, P., St-Jean, C., & Zacharias, P. 2005, ApJ, 619, 613

    work page 2005

  7. [7]

    R., & Latour, J

    Durney, B. R., & Latour, J. 1977, Geophysical and Astrophysical Fluid Dynamics, 9, 241

    work page 1977

  8. [8]

    2017, PhD thesis, Montana State University, Bozeman,

    Egeland, R. 2017, PhD thesis, Montana State University, Bozeman,

    work page 2017

  9. [9]

    B., Egeland, R., Soon, W., & Baliunas, S

    Garg, S., Karak, B. B., Egeland, R., Soon, W., & Baliunas, S. 2019, ApJ, 886, 132

    work page 2019

  10. [10]

    Gilman, P. A. 1980, in IAU Colloq. 51: Stellar Turbulence, ed. D. F. Gray & J. L. Linsky, V ol. 114, 19

    work page 1980

  11. [11]

    H., Wilson, R

    Hathaway, D. H., Wilson, R. M., & Reichmann, E. J. 1994, Solar Physics, 151, 177

    work page 1994

  12. [12]

    B., & Kitchatinov, L

    Hazra, G., Jiang, J., Karak, B. B., & Kitchatinov, L. 2019, ApJ, 884, 35 10 GARGET AL. Table 3.Cycle fit parameters for the additional GK25 stars. Star cyct start Pcyc a b c t 0 χ2 red HD 1835 1* 1967.805 - 0.000001 261.22799 1.11764 1843.470 11.6 2 1972.638 4.67−2.293e+05−11.17197−8.342e+10 2027.498 14.9 3 1977.305 4.5 0.00039 11.94459 0.84901 1967.274 1...

    work page arXiv 2019

  13. [13]

    Karak, B. B. 2023, Living Reviews in Solar Physics, 20, 3

    work page 2023

  14. [14]

    B., & Cameron, R

    Karak, B. B., & Cameron, R. 2016, ApJ, 832, 94

    work page 2016

  15. [15]

    B., & Choudhuri, A

    Karak, B. B., & Choudhuri, A. R. 2011, MNRAS, 410, 1503

    work page 2011

  16. [16]

    B., Kitchatinov, L

    Karak, B. B., Kitchatinov, L. L., & Choudhuri, A. R. 2014, ApJ, 791, 59

    work page 2014

  17. [17]

    B., Mandal, S., & Banerjee, D

    Karak, B. B., Mandal, S., & Banerjee, D. 2018, ApJ, 866, 17

    work page 2018

  18. [18]

    B., & Miesch, M

    Karak, B. B., & Miesch, M. 2017, ApJ, 847, 69

    work page 2017

  19. [19]

    L., & Olemskoy, S

    Kitchatinov, L. L., & Olemskoy, S. V . 2011, MNRAS, 411, 1059 —. 2012, MNRAS, 423, 3344

    work page 2011

  20. [20]

    L., & R¨udiger, G

    Kitchatinov, L. L., & R¨udiger, G. 1999, A&A, 344, 911

    work page 1999

  21. [21]

    Krause, F., & R¨adler, K. H. 1980, Mean-field magnetohydrodynamics and dynamo theory (Oxford: Pergamon Press)

    work page 1980

  22. [22]

    Kumar, P., Biswas, A., & Karak, B. B. 2022, MNRAS, 513, L112

    work page 2022

  23. [23]

    B., & Vashishth, V

    Kumar, P., Karak, B. B., & Vashishth, V . 2021, ApJ, 913, 65

    work page 2021

  24. [24]

    L., & Brueckner, G

    Lean, J. L., & Brueckner, G. E. 1989, ApJ, 337, 568

    work page 1989

  25. [25]

    K., Wright, J

    Luhn, J. K., Wright, J. T., Henry, G. W., Saar, S. H., & Baum, A. C. 2022, ApJL, 936, L23

    work page 2022

  26. [26]

    Nandy, D., & Martens, P. C. H. 2007, Advances in Space Research, 40, 891

    work page 2007

  27. [27]

    Vaughan, A. H. 1984, ApJ, 279, 763 Ol´ah, K., K˝ov´ari, Z., Petrovay, K., et al. 2016, A&A, 590, A133

    work page 1984

  28. [28]

    2004, PASP, 116, 699

    Paxton, B. 2004, PASP, 116, 699

    work page 2004

  29. [29]

    2020, Living Reviews in Solar Physics, 17, 2

    Petrovay, K. 2020, Living Reviews in Solar Physics, 17, 2

    work page 2020

  30. [30]

    Rengarajan, T. N. 1984, ApJL, 283, L63

    work page 1984

  31. [31]

    H., & Brandenburg, A

    Saar, S. H., & Brandenburg, A. 2002, Astronomische Nachrichten, 323, 357 Schr¨oder, K. P., Mittag, M., Hempelmann, A., Gonz´alez-P´erez, J. N., & Schmitt, J. H. M. M. 2013, A&A, 554, A50

    work page 2002

  32. [32]

    P., Wright, J

    Shah, S. P., Wright, J. T., Isaacson, H., Howard, A. W., & Curtis, J. L. 2018, ApJL, 863, L26

    work page 2018

  33. [33]

    1972, ApJ, 171, 565

    Skumanich, A. 1972, ApJ, 171, 565

    work page 1972

  34. [34]

    H., Baliunas, S

    Soon, W. H., Baliunas, S. L., & Zhang, Q. 1993, ApJL, 414, L33 —. 1994, SoPh, 154, 385

    work page 1993

  35. [35]

    2023, Journal of Space Weather and Space Climate, 13, 32

    Thibeault, C., Miara, L., & Charbonneau, P. 2023, Journal of Space Weather and Space Climate, 13, 32

    work page 2023

  36. [36]

    B., & Kitchatinov, L

    Vashishth, V ., Karak, B. B., & Kitchatinov, L. 2023, MNRAS, 522, 2601

    work page 2023

  37. [37]

    A., Gregory, S

    Vidotto, A. A., Gregory, S. G., Jardine, M., et al. 2014, MNRAS, 441, 2361

    work page 2014

  38. [38]

    1935, Astronomische Mitteilungen der Eidgen¨ossischen Sternwarte Zurich, 14, 105

    Waldmeier, M. 1935, Astronomische Mitteilungen der Eidgen¨ossischen Sternwarte Zurich, 14, 105

    work page 1935

  39. [39]

    J., Drake, J

    Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48

    work page 2011