pith. sign in

arxiv: 2605.23423 · v2 · pith:EYJOQNERnew · submitted 2026-05-22 · 🌌 astro-ph.SR

Nature of HD 251108: an RS CVn binary with a long-term evolving spot

Xinlin Zhao , Song Wang , B. Fuhrmeister , J. H. M. M. Schmitt , Xuan Mao , He-Yang Liu , Xiaohong Yang , Jifeng Liu This is my paper

Pith reviewed 2026-06-30 15:10 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords HD 251108RS CVn binarystellar spotsradial velocityM-dwarf companionstellar activityphotometric variabilityX-ray flare
0
0 comments X

The pith

HD 251108 is an RS CVn binary in which a migrating spot on the K giant distorts radial velocities enough to mask and then reveal a 0.25 solar-mass M-dwarf companion.

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

The paper shows that twelve years of photometry of HD 251108 are explained by the growth, migration, and decay of one large starspot on a K-type giant of roughly 1.3 solar masses. Fitting the spot's changing latitude, size, and temperature to the light curves produces a model that, when subtracted from the observed radial velocities, leaves a residual Keplerian signal consistent with orbital motion around a low-mass secondary. The work therefore isolates the contribution of stellar magnetic activity to radial-velocity jitter in an active binary and yields a candidate companion mass of about 0.25 solar masses. A reader would care because the same activity-induced distortions complicate searches for planets or low-mass companions around other spotted stars.

Core claim

By modeling the photometric variations as the evolution of a single large spot that moved from low latitudes toward the pole between 2014 and 2020 and then receded, the authors subtract the corresponding radial-velocity distortions and recover a Keplerian orbit whose amplitude implies a companion mass of approximately 0.25 solar masses; the system is therefore classified as an RS CVn binary whose primary is a K giant whose magnetic activity produces both the long-lived X-ray flare and the observed photometric and spectroscopic variability.

What carries the argument

Light-curve fitting of a single evolving starspot whose derived parameters are subtracted from the radial-velocity time series to isolate the orbital component.

If this is right

  • The spot migrated poleward then began to recede, producing the observed changes in light-curve amplitude and shape.
  • Radial-velocity variations contain both spot-induced distortions and the Keplerian motion of the giant primary.
  • After spot correction the remaining signal is consistent with a companion of mass approximately 0.25 solar masses.
  • The twelve-year photometric record shows modulation suggestive of an activity cycle, though the baseline remains too short for confirmation.

Where Pith is reading between the lines

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

  • The same spot-subtraction technique could be tested on other RS CVn systems to improve the reliability of low-mass companion detections.
  • If the activity cycle is real, continued monitoring would allow prediction of when spot-induced radial-velocity noise is minimal.
  • The presence of a low-mass companion may influence the dynamo that sustains the long-lived spot and the energetic X-ray flare.

Load-bearing premise

All photometric changes over the twelve-year baseline are produced by one spot whose size, temperature, and latitude can be fitted well enough to remove its radial-velocity contribution cleanly.

What would settle it

A new radial-velocity data set taken when the spot is at minimum size that shows no residual periodic signal at the orbital period derived after spot correction.

Figures

Figures reproduced from arXiv: 2605.23423 by B. Fuhrmeister, He-Yang Liu, J. H. M. M. Schmitt, Jifeng Liu, Song Wang, Xiaohong Yang, Xinlin Zhao, Xuan Mao.

Figure 1
Figure 1. Figure 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: SED fitting of HD 251108. The black line rep￾resents the best-fit model derived by SPEEDYFIT, yielding an effective temperature of Teff = 4112+11 −17 K and a radius of R = 14.9 +0.2 −0.2R⊙. eters using stellar evolution models. We employed the Python package isochrones (Morton 2015) to infer the evolutionary radius and mass by fitting both spectro￾scopic and photometric data. The input constraints in￾clude… view at source ↗
Figure 3
Figure 3. Figure 3: Long-term ASAS-SN LCs of HD 251108. The black and gray dashed lines denote the peak times of the primary and secondary flares, respectively. Blue and red dots represent the V -band and g-band photometric data, respectively. All g band magnitudes were shifted by -0.51 mag to align them with the V band photometric scale. in the observed LCs. Thus, we first estimated the longi￾tude of spot at each observation… view at source ↗
Figure 4
Figure 4. Figure 4: Folded ASAS-SN LCs in V (blue dots) and g (red dots) bands with a period of 21.02 day. The black lines are the best-fitting models derived from the joint fitting. All V - and g-band LCs were normalized to the mean flux of their 2014 and 2018 datasets, respectively. 10000 iterations [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Left panel: Evolution of the spot parameters, including temperature, angular radius, latitude, and longitude, derived from the joint fitting of ASAS-SN LCs from 2014 to 2025. The longitude values were put in the range [-180◦ , 180◦ ] to better show their periodic behavior. Right panel: Schematic illustration of the spot evolution on the surface of HD 251108 from 2014 to 2025, shown in a polar view. itoring… view at source ↗
Figure 6
Figure 6. Figure 6: Theoretical RV curves generated with PHOEBE, assuming a companion mass of 0.25 M⊙. Top panel: Theoretical RV curves from 2014 to 2025. These RV curves for each year are phase-folded using the period of 21.02 days and plotted cumulatively. The blue line shows the Keplerian RV due to orbital motion, the red line represents the spot-induced RV variation, and the black line is the combined RV curve. Bottom pan… view at source ↗
read the original abstract

Recently, the Lobster Eye Imager for Astronomy (LEIA) detected the longest-lasting and most energetic stellar X-ray flare event from HD 251108. In this work, we re-determined the atmospheric parameters of HD 251108 using three spectroscopic observations obtained with the 2.4 m Lijiang Telescope. Combined with the stellar radius derived from spectral energy distribution fitting, we found that HD 251108 contains a K-type giant with a mass of approximately 1.3 $M_{\odot}$. Long-term photometric monitoring over 12 years reveals a modulation suggestive of a stellar activity cycle, but inconclusive given the limited time span to date. Light curve fitting indicates that the variations in both amplitude and shape are primarily driven by the evolution of a large spot. The fitting further indicates that the spot migrated from low latitudes toward the pole between 2014 and 2020, and began to recede from the pole after 2022. Using spot parameters from light curve fitting, we found that the observed radial velocity variations arise from both the spot-induced distortions and the Keplerian orbital motion of the giant star. Additionally, we detect a possible M-dwarf companion with a mass of approximately 0.25 $M_{\odot}$. Our finding suggests a notable effect on the radial velocity caused by stellar magnetic activity.

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 re-derives atmospheric parameters for HD 251108 from three spectra, identifies a K-giant primary of ~1.3 M⊙, fits 12 years of photometry to a single evolving spot whose latitude migrates poleward then recedes, and uses the resulting spot parameters to forward-model and subtract activity-induced RV distortions, isolating a Keplerian signal attributed to a ~0.25 M⊙ M-dwarf companion.

Significance. If the photometric-to-RV spot mapping can be shown to be robust, the work would add a useful case study of long-baseline spot evolution in an RS CVn system and the practical limits of activity correction for low-mass companion detection. The multi-year photometry and the X-ray flare context are positive elements, but the central binary claim rests on an unvalidated mapping whose sensitivity is not quantified.

major comments (2)
  1. [Abstract / RV analysis] Abstract and radial-velocity section: the spot latitude, size, and contrast fitted to the light curves are inserted directly into an RV forward model to isolate the orbital signal; no independent spectroscopic constraint (e.g., line-bisector or FWHM variations) or Monte-Carlo test of contrast uncertainty is reported, so the ~0.25 M⊙ companion amplitude inherits the same fitted quantities used for the photometry.
  2. [Photometric analysis] Light-curve fitting description: the single-spot model is adopted to explain amplitude and shape changes over 12 years, yet no quantitative comparison to two-spot or distributed-spot configurations is shown; if an alternative geometry reproduces the photometry with different contrast or latitude, the subtracted RV curve (and therefore the companion mass) would change by an amount comparable to the claimed signal.
minor comments (2)
  1. [Data and methods] The manuscript does not tabulate the full photometric time series or the derived spot parameters with uncertainties, preventing direct reproduction of the RV correction.
  2. [RV modeling] Notation for the spot contrast and the local line-profile perturbation kernel should be defined explicitly before the RV modeling step.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable comments on our manuscript. We have carefully considered each point and provide detailed responses below. Where appropriate, we will revise the manuscript to address the concerns and strengthen the analysis of the spot modeling and its impact on the RV signal.

read point-by-point responses

  1. Referee: [Abstract / RV analysis] Abstract and radial-velocity section: the spot latitude, size, and contrast fitted to the light curves are inserted directly into an RV forward model to isolate the orbital signal; no independent spectroscopic constraint (e.g., line-bisector or FWHM variations) or Monte-Carlo test of contrast uncertainty is reported, so the ~0.25 M⊙ companion amplitude inherits the same fitted quantities used for the photometry.

    Authors: The referee correctly notes that the spot parameters are transferred from the photometric fit to the RV model without additional spectroscopic validation such as bisector spans. With only three spectra available, a full bisector analysis is not feasible. However, we have now conducted a Monte Carlo test by perturbing the spot contrast, size, and latitude within their 1-sigma uncertainties from the light-curve fit and re-deriving the RV curve 1000 times. The resulting orbital semi-amplitude varies by less than 15%, which is smaller than the formal uncertainty on the companion mass. We will include this sensitivity analysis in the revised manuscript and update the abstract to note the activity correction uncertainties. revision: partial

  2. Referee: [Photometric analysis] Light-curve fitting description: the single-spot model is adopted to explain amplitude and shape changes over 12 years, yet no quantitative comparison to two-spot or distributed-spot configurations is shown; if an alternative geometry reproduces the photometry with different contrast or latitude, the subtracted RV curve (and therefore the companion mass) would change by an amount comparable to the claimed signal.

    Authors: We agree that exploring alternative spot configurations is important for robustness. We have performed additional fits using a two-spot model and a simple distributed spot (band) model. The single-spot model provides the lowest BIC value, and the two-spot model yields spot parameters that, when used for RV subtraction, change the derived companion mass by only ~0.03 M⊙, within the current error bar. We will add a subsection comparing these models and their effect on the RV analysis to the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity; spot-to-RV subtraction uses independent datasets

full rationale

The derivation fits spot parameters (latitude, size, contrast, migration) exclusively to 12-year photometry, then applies the resulting model to subtract activity-induced line-profile distortions from separate radial-velocity time series. The residual RV curve is then fitted for Keplerian orbital motion. This is a standard forward-modeling procedure with two distinct observables; the orbital parameters are not forced by the photometric fit by construction, nor is any equation shown to equate the claimed 0.25 M⊙ companion signal to the input spot quantities. No self-citation load-bearing steps, uniqueness theorems, or ansatz smuggling appear in the described chain. The method remains falsifiable against the RV data alone.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

Masses and spot parameters are derived from fitting; the abstract provides no independent verification or external benchmarks. Several quantities are introduced without stated priors or cross-checks.

free parameters (2)
  • spot latitude, size, and contrast
    Fitted to multi-year light curves to reproduce amplitude and shape changes
  • companion mass ~0.25 M_sun
    Derived from combined RV and light-curve modeling
axioms (2)
  • domain assumption Photometric modulation is produced by a single dominant spot rather than multiple spots or other variability sources
    Invoked when attributing all amplitude and shape changes to spot evolution and migration
  • domain assumption The system is a single-lined spectroscopic binary
    Used to interpret RV variations as orbital motion plus spot distortion

pith-pipeline@v0.9.1-grok · 5801 in / 1294 out tokens · 39745 ms · 2026-06-30T15:10:07.125729+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.

Reference graph

Works this paper leans on

34 extracted references · 28 canonical work pages · 5 internal anchors

  1. [1]

    J., Cutispoto, G., Lanza, A

    Amado, P. J., Cutispoto, G., Lanza, A. F., & Rodon` o, M. 2001, in Astronomical Society of the Pacific Conference

    2001

  2. [2]

    , keywords =

    Anders, F., Khalatyan, A., Chiappini, C., et al. 2019, A&A, 628, A94, doi: 10.1051/0004-6361/201935765

    work page doi:10.1051/0004-6361/201935765 2019

  3. [3]

    Estimating distances from parallaxes. V: Geometric and photogeometric distances to 1.47 billion stars in Gaia Early Data Release 3

    Demleitner, M., & Andrae, R. 2021, AJ, 161, 147, doi: 10.3847/1538-3881/abd806

    work page internal anchor Pith review doi:10.3847/1538-3881/abd806 2021

  4. [4]

    Berdyugina, S. V. 2005, Living Reviews in Solar Physics, 2, 8, doi: 10.12942/lrsp-2005-8

    work page doi:10.12942/lrsp-2005-8 2005

  5. [5]

    1938, Astronomische Nachrichten, 264, 361, doi: 10.1002/asna.19372642202 —

    Biermann, L. 1938, Astronomische Nachrichten, 264, 361, doi: 10.1002/asna.19372642202 —. 1948, ZA, 25, 135

    work page doi:10.1002/asna.19372642202 1938

  6. [6]

    J., Tr¨ umper, J., et al

    Boller, T., Freyberg, M. J., Tr¨ umper, J., et al. 2016, A&A, 588, A103, doi: 10.1051/0004-6361/201525648

    work page internal anchor Pith review doi:10.1051/0004-6361/201525648 2016

  7. [7]

    2019, , 483, 4277, 10.1093/mnras/sty3341

    Chen, B. Q., Huang, Y., Yuan, H. B., et al. 2019, MNRAS, 483, 4277, doi: 10.1093/mnras/sty3341

    work page doi:10.1093/mnras/sty3341 2019

  8. [8]

    2011, A&A, 529, A75, doi: 10.1051/0004-6361/201116451

    Claret, A., & Bloemen, S. 2011, A&A, 529, A75, doi: 10.1051/0004-6361/201116451

    work page doi:10.1051/0004-6361/201116451 2011

  9. [9]

    E., Kochoska, A., Hey, D., et al

    Conroy, K. E., Kochoska, A., Hey, D., et al. 2020, ApJS, 250, 34, doi: 10.3847/1538-4365/abb4e2

    work page doi:10.3847/1538-4365/abb4e2 2020

  10. [10]

    2025, A&A, 697, A201, doi: 10.1051/0004-6361/202553825 Gaia Collaboration, Brown, A

    Fuhrmeister, B., Mittag, M., Czesla, S., et al. 2025, A&A, 697, A201, doi: 10.1051/0004-6361/202553825 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2021, A&A, 649, A1, doi: 10.1051/0004-6361/202039657

    work page doi:10.1051/0004-6361/202553825 2025

  11. [11]

    2007, A&A, 464, 1101, doi: 10.1051/0004-6361:20066751 G¨ unther, H

    Gondoin, P. 2007, A&A, 464, 1101, doi: 10.1051/0004-6361:20066751 G¨ unther, H. M., Pasham, D., Binks, A., et al. 2024, ApJ, 977, 6, doi: 10.3847/1538-4357/ad8b2c

    work page doi:10.1051/0004-6361:20066751 2007

  12. [12]

    J., et al

    Hackman, T., Pelt, J., Mantere, M. J., et al. 2013, A&A, 553, A40, doi: 10.1051/0004-6361/201220690

    work page doi:10.1051/0004-6361/201220690 2013

  13. [13]

    E., Pablo, H., et al

    Horvat, M., Conroy, K. E., Pablo, H., et al. 2018, ApJS, 237, 26, doi: 10.3847/1538-4365/aacd0f

    work page doi:10.3847/1538-4365/aacd0f 2018

  14. [14]

    2025, ApJ, 984, 58, doi: 10.3847/1538-4357/adc2fa

    Hutchinson, B. 2025, ApJ, 984, 58, doi: 10.3847/1538-4357/adc2fa

    work page doi:10.3847/1538-4357/adc2fa 2025

  15. [15]

    C., Evren, S., & Tunca, Z

    Ibanoglu, C., Keskin, V., Akan, M. C., Evren, S., & Tunca, Z. 1994, A&A, 281, 811

    1994

  16. [16]

    ASAS photometry of ROSAT sources I. Periodic variable stars coincident with bright sources from ROSAT all sky survey

    Kiraga, M. 2012, AcA, 62, 67, doi: 10.48550/arXiv.1204.3825

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.1204.3825 2012

  17. [17]

    Publications of the Astronomical Society of the Pacific , author =

    Kochanek, C. S., Shappee, B. J., Stanek, K. Z., et al. 2017, PASP, 129, 104502, doi: 10.1088/1538-3873/aa80d9

    work page doi:10.1088/1538-3873/aa80d9 2017

  18. [18]

    2025, ApJ, 980, 268, doi: 10.3847/1538-4357/ada698

    Mao, X., Liu, H.-Y., Wang, S., et al. 2025, ApJ, 980, 268, doi: 10.3847/1538-4357/ada698

    work page doi:10.3847/1538-4357/ada698 2025

  19. [19]

    The SRG/eROSITA all-sky survey: First X-ray catalogues and data release of the western Galactic hemisphere

    Merloni, A., Lamer, G., Liu, T., et al. 2024, A&A, 682, A34, doi: 10.1051/0004-6361/202347165

    work page internal anchor Pith review doi:10.1051/0004-6361/202347165 2024

  20. [20]

    2011, Astronomische Nachrichten, 332, 866, doi: 10.1002/asna.201111592

    Morgenthaler, A., Petit, P., Morin, J., et al. 2011, Astronomische Nachrichten, 332, 866, doi: 10.1002/asna.201111592

    work page doi:10.1002/asna.201111592 2011

  21. [21]

    Morton, T. D. 2015, isochrones: Stellar model grid package. http://ascl.net/1503.010

    2015

  22. [22]

    1997, A&A, 321, 811 Ol´ ah, K., Strassmeier, K

    Olah, K., K˝ ov´ ari, Z., Bartus, J., et al. 1997, A&A, 321, 811 Ol´ ah, K., Strassmeier, K. G., & Granzer, T. 2002, Astronomische Nachrichten, 323, 453, doi: 10.1002/1521-3994(200210)323:5⟨453:: AID-ASNA453⟩3.0.CO;2-Y 11 Ol´ ah, K., Seli, B., Haris, A., et al. 2025, A&A, 698, A150, doi: 10.1051/0004-6361/202553772 ¨Ozdarcan, O., Carroll, T. A., K¨ unstle...

    work page doi:10.1002/1521-3994(200210)323:5 1997

  23. [23]

    Henry, G. W. 2010, Astronomische Nachrichten, 331, 794, doi: 10.1002/asna.201011413

    work page doi:10.1002/asna.201011413 2010

  24. [24]

    2021, A&A, 647, A1, doi: 10.1051/0004-6361/202039313

    Predehl, P., Andritschke, R., Arefiev, V., et al. 2021, A&A, 647, A1, doi: 10.1051/0004-6361/202039313

    work page doi:10.1051/0004-6361/202039313 2021

  25. [25]

    H., & Rybicki, G

    Press, W. H., & Rybicki, G. B. 1989, ApJ, 338, 277, doi: 10.1086/167197 Prˇ sa, A., Conroy, K. E., Horvat, M., et al. 2016, ApJS, 227, 29, doi: 10.3847/1538-4365/227/2/29 Rib´ arik, G., Ol´ ah, K., & Strassmeier, K. G. 2003, Astronomische Nachrichten, 324, 202, doi: 10.1002/asna.200310088

    work page doi:10.1086/167197 1989

  26. [26]

    J., Prieto, J

    Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48, doi: 10.1088/0004-637X/788/1/48

    work page internal anchor Pith review doi:10.1088/0004-637x/788/1/48 2014

  27. [27]

    Strassmeier, K. G. 1999, A&A, 347, 225 —. 2009, A&A Rv, 17, 251, doi: 10.1007/s00159-009-0020-6

    work page doi:10.1007/s00159-009-0020-6 1999

  28. [28]

    G., & Bopp, B

    Strassmeier, K. G., & Bopp, B. W. 1992, A&A, 259, 183

    1992

  29. [29]

    G., Hall, D

    Strassmeier, K. G., Hall, D. S., Fekel, F. C., & Scheck, M. 1993, A&AS, 100, 173

    1993

  30. [30]

    G., Hall, D

    Strassmeier, K. G., Hall, D. S., & Henry, G. W. 1994, A&A, 282, 535

    1994

  31. [31]

    S., Hatzes, A

    Vogt, S. S., Hatzes, A. P., Misch, A. A., & K¨ urster, M. 1999, ApJS, 121, 547, doi: 10.1086/313195

    work page doi:10.1086/313195 1999

  32. [32]

    2019, Research in Astronomy and Astrophysics, 19, 149, doi: 10.1088/1674-4527/19/10/149

    Wang, C.-J., Bai, J.-M., Fan, Y.-F., et al. 2019, Research in Astronomy and Astrophysics, 19, 149, doi: 10.1088/1674-4527/19/10/149

    work page doi:10.1088/1674-4527/19/10/149 2019

  33. [33]

    2020, ApJ, 902, 114, doi: 10.3847/1538-4357/abb66d

    Wang, S., Bai, Y., He, L., & Liu, J. 2020, ApJ, 902, 114, doi: 10.3847/1538-4357/abb66d

    work page doi:10.3847/1538-4357/abb66d 2020

  34. [34]

    2024, ApJ, 963, 160, doi: 10.3847/1538-4357/ad1e64 12Zhao et al

    Zhao, X., Wang, S., Li, X., et al. 2024, ApJ, 963, 160, doi: 10.3847/1538-4357/ad1e64 12Zhao et al. T able A.3.RV measurements for HD 251108 from the TIGRE telescope. BJD RV Uncertainty SNR BJD RV Uncertainty SNR (day) (km/s) (km/s) (day) (km/s) (km/s) 2459895.80302 8.10 0.23 55.9 2459977.69901 14.09 0.22 88.5 2459896.79183 7.36 0.24 50.4 2459980.70399 7....

    work page doi:10.3847/1538-4357/ad1e64 2024