pith. sign in

arxiv: 2606.29543 · v1 · pith:ZXJKS6JEnew · submitted 2026-06-28 · 🌌 astro-ph.SR

A Reproducible AAVSO Johnson-V Fourier Template for the Prototype Cepheid Delta Cephei

Zuhoor Elahi , Wafa Gull This is my paper

Pith reviewed 2026-06-30 01:56 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords Delta CepheiCepheid variableslight curveFourier analysisAAVSO observationsJohnson V photometrystellar pulsation
0
0 comments X

The pith

A third-order Fourier series provides a reproducible Johnson-V light curve template for Delta Cephei.

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

The paper constructs an empirical Fourier template from AAVSO Johnson-V observations of the prototype Cepheid Delta Cephei to serve as a benchmark for theoretical comparisons. It starts with 244 measurements spanning 355 days, rejects two outliers to leave 242 points, and phase-folds them using a period of 5.366531 days and an empirical maximum epoch. A third-order Fourier model is fitted and adopted after checking that orders four through six reduce the RMS residual by only 0.0012 magnitudes. The resulting template supplies explicit coefficients, amplitude, rise fraction, and asymmetry values. The work supplies an observational morphology target rather than a physical model of the pulsation.

Core claim

We present an empirical Fourier reconstruction of the observed Johnson-V light curve of Delta Cephei. Using an adopted period of 5.366531 d and 242 cleaned AAVSO points phased to an empirical bright-maximum epoch, a third-order Fourier model is adopted with A0 = 3.9031, A1 = 0.3434 mag, A2 = 0.1428 mag, and A3 = 0.0531 mag. This yields R21 = 0.4159, R31 = 0.1547, a full amplitude of 0.8544 mag, a rise fraction of 0.2885, and an asymmetry index of 0.4230. Higher orders add negligible improvement, so the third-order form is retained as the simplest adequate empirical template.

What carries the argument

The third-order Fourier series fit applied to the phase-folded AAVSO Johnson-V magnitudes of Delta Cephei.

If this is right

  • The N=3 template supplies a direct observational target for comparison with nonlinear pulsation calculations and synthetic photometry.
  • Bootstrap uncertainties on the Fourier coefficients permit quantitative tests of how well any given model reproduces the observed morphology.
  • The remaining non-random scatter in residuals indicates that observer-level effects must be accounted for in future data sets used with this template.

Where Pith is reading between the lines

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

  • The template can serve as a fixed reference shape against which light curves of other Cepheids can be compared to quantify morphological differences across the instability strip.
  • If radial-velocity data are later combined with this template, the resulting distance or radius estimates inherit the reported amplitude and asymmetry directly.

Load-bearing premise

The adopted period and bright-maximum epoch are accurate enough that phase folding does not distort the true light-curve shape after outlier removal.

What would settle it

A fresh set of Johnson-V observations of Delta Cephei, phased with the same period and epoch, that produces a significantly different amplitude, rise fraction, or asymmetry from the reported N=3 template.

Figures

Figures reproduced from arXiv: 2606.29543 by Wafa Gull, Zuhoor Elahi.

Figure 1
Figure 1. Figure 1: Phase-folded AAVSO Johnson-V light curve of Delta Cephei using Pobs = 5.366531 d. The data are repeated over two cycles for visual continuity. Blue points show the cleaned observations, and orange points show phase-binned medians. The magnitude axis is inverted so that brighter phases appear higher in the figure. The empirical template was fit with a Fourier series of the form V (ϕ) = A0 + X N k=1 Ak cos (… view at source ↗
Figure 2
Figure 2. Figure 2: Adopted empirical Fourier template for Delta Cephei. The blue points are the cleaned AAVSO Johnson-V observations and the solid curve is the adopted third-order Fourier template. The fit is intended as a low-order morphology template rather than a precision model for every individual observation [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Fourier harmonic amplitude spectrum for the adopted N = 3 empirical template. The nonzero A2 and A3 terms quantify the departure of the observed light curve from a pure sinusoid. fdecline = 1 − frise. (17) The asymmetry index is defined as Aasym = fdecline − frise fdecline + frise . (18) Since fdecline + frise = 1, this is equivalent to Aasym = 1 − 2frise. (19) For the adopted template, frise = 0.2885, (20… view at source ↗
Figure 4
Figure 4. Figure 4: Morphology summary for the adopted empirical Johnson-V Fourier template. The bright maximum corresponds to the minimum value of the magnitude curve, while the faint minimum corresponds to the maximum value. The adopted N = 3 template gives ∆V = 0.854439 mag, frise = 0.2885, and Aasym = 0.4230 [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: shows the residuals as a function of phase. The RMS residual is 0.096996 mag, while the median absolute deviation is 0.080869 mag. The median residual is 0.006955 mag, indicating no large global offset. The maximum absolute residual is 0.515884 mag. 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 Pulsation phase 0.5 0.4 0.3 0.2 0.1 0.0 0.1 0.2 0.3 Vobs VFourier m a g Fourier-template residuals, RMS = 0.0970 m… view at source ↗
Figure 6
Figure 6. Figure 6: shows the N = 3–6 templates over the phase-folded observations. The curves are visually very similar over most of the pulsation cycle. The higher-order fits slightly modify the extrema but do not change the main empirical morphology. We therefore adopt N = 3 for the main benchmark values. 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 Pulsation phase 3.4 3.6 3.8 4.0 4.2 4.4 Jo h nso n V m a g nitu d e Fourie… view at source ↗
Figure 7
Figure 7. Figure 7: Bootstrap distributions for the adopted N = 3 Fourier template. The distributions show the stability of the main empirical morphology diagnostics under resampling of the cleaned Johnson-V dataset. 13 [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Observer-level mean residuals for the adopted N = 3 Fourier template. Residuals are defined as Vobs − VFourier. The figure shows that the residual structure contains observer-dependent offsets at the several-hundredths to ∼ 0.1 mag level for the principal observer groups. 4. stellar oscillation models based on the GYRE code when coupled to an external light-curve or atmosphere calculation; 5. independent n… view at source ↗
read the original abstract

We present an empirical Fourier reconstruction of the observed Johnson-V light curve of the prototype Classical Cepheid Delta Cephei. The goal is not to infer a full physical stellar model but to establish a reproducible observed-light-curve benchmark for later comparison with nonlinear pulsation, synthetic photometry, Baade-Wesselink/SPIPS, GYRE-supported, and independent hydrodynamic calculations. Using an adopted period of Pobs = 5.366531 d, 244 AAVSO Johnson-V measurements were filtered to a cleaned sample of 242 points after rejecting two extreme outliers. The cleaned data span 355.09259 d and were phase folded using an empirical bright-maximum epoch of JD = 2460851.395800. We fit a low-order Fourier model to the phased light curve and adopt a third-order template as the preferred empirical morphology representation. The adopted N = 3 fit gives A0 = 3.9031, A1 = 0.3434 mag, A2 = 0.1428 mag, and A3 = 0.0531 mag, corresponding to R21 = 0.4159 and R31 = 0.1547. The template has a full Johnson-V amplitude of Delta V = 0.8544 mag, a rise fraction of frise = 0.2885, and an asymmetry index of Aasym = 0.4230. Bootstrap uncertainties are reported in the manuscript. Fourier orders N = 4-6 reduce the RMS residual by only about 0.0012 mag relative to the N = 3 model, so the third-order representation is retained as the simplest adequate empirical template. Observer-level residual diagnostics show that the remaining scatter is not purely random. This paper provides an observational morphology target rather than a physical explanation of the pulsation.

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 derives an empirical N=3 Fourier template for the Johnson-V light curve of Delta Cephei from 242 cleaned AAVSO observations spanning 355 d. Using an adopted period Pobs=5.366531 d and bright-maximum epoch JD=2460851.395800, the phased data are fit to yield A0=3.9031, A1=0.3434, A2=0.1428, A3=0.0531 mag (with R21=0.4159, R31=0.1547), Delta V=0.8544 mag, frise=0.2885, and Aasym=0.4230. The N=3 model is retained after comparing RMS residuals to higher orders, with bootstrap uncertainties provided; the template is positioned as a reproducible observational benchmark for theoretical comparisons rather than a physical model.

Significance. If the adopted phasing accurately represents the data without distortion, the work supplies a concrete, quantitative morphology target (with explicit coefficients and derived shape parameters) that can be directly compared against nonlinear pulsation models, synthetic photometry, Baade-Wesselink analyses, and hydrodynamic simulations. The explicit retention of N=3 on the basis of marginal RMS improvement and the reporting of bootstrap uncertainties strengthen its utility as a reproducible reference.

major comments (3)
  1. [Data preparation and phase-folding procedure] Data preparation section: the period Pobs=5.366531 d and epoch JD=2460851.395800 are adopted without re-derivation (e.g., via Lomb-Scargle or least-squares) from the AAVSO points or a sensitivity test for cumulative phase drift. Over the 355 d baseline (~66 cycles), even a 10^{-5} d error in P produces ~0.06-cycle drift that would systematically alter the folded morphology, the fitted amplitudes, R21, R31, frise, and Aasym.
  2. [Data preparation and phase-folding procedure] Data preparation section: the rejection of two extreme outliers (reducing 244 to 242 points) is stated without quantitative criteria such as a sigma threshold, magnitude limits, or explicit values of the rejected points. This directly affects whether the retained sample is representative for the central N=3 fit and derived parameters.
  3. [Fitting procedure and bootstrap uncertainties] Methods section on fitting and uncertainties: while bootstrap uncertainties are mentioned, the procedure (resampling method, number of iterations, handling of phase-folding) is not specified in sufficient detail to allow independent reproduction of the reported coefficient errors.
minor comments (2)
  1. [Abstract and Methods] The abstract and main text should explicitly state the quantitative outlier rejection threshold and the exact implementation of phase folding to support the reproducibility claim in the title.
  2. [Results] Notation for the asymmetry index Aasym and rise fraction frise should be defined with an equation or explicit formula in the text rather than only in the results.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the positive evaluation of the manuscript's significance and for the constructive comments aimed at improving reproducibility. We address each major comment point by point below.

read point-by-point responses

  1. Referee: Data preparation section: the period Pobs=5.366531 d and epoch JD=2460851.395800 are adopted without re-derivation (e.g., via Lomb-Scargle or least-squares) from the AAVSO points or a sensitivity test for cumulative phase drift. Over the 355 d baseline (~66 cycles), even a 10^{-5} d error in P produces ~0.06-cycle drift that would systematically alter the folded morphology, the fitted amplitudes, R21, R31, frise, and Aasym.

    Authors: The period and epoch were adopted as established literature values for Delta Cephei to maintain consistency with prior work. We acknowledge the validity of the phase-drift concern over the 355 d baseline. In the revised manuscript we will add an explicit sensitivity test by perturbing P by amounts up to ±10^{-5} d, re-phasing the data, and re-fitting to quantify any changes in the coefficients and shape parameters. revision: yes

  2. Referee: Data preparation section: the rejection of two extreme outliers (reducing 244 to 242 points) is stated without quantitative criteria such as a sigma threshold, magnitude limits, or explicit values of the rejected points. This directly affects whether the retained sample is representative for the central N=3 fit and derived parameters.

    Authors: We agree that the outlier rejection must be documented quantitatively for reproducibility. The revised data-preparation section will specify the exact criteria (including any sigma threshold or magnitude limits applied) and will list the observed magnitudes of the two rejected points. revision: yes

  3. Referee: Methods section on fitting and uncertainties: while bootstrap uncertainties are mentioned, the procedure (resampling method, number of iterations, handling of phase-folding) is not specified in sufficient detail to allow independent reproduction of the reported coefficient errors.

    Authors: We will expand the methods section to provide the missing details: the bootstrap uses sampling with replacement, 1000 iterations are performed, and each resampled dataset is re-phased using the same fixed period and epoch before refitting. revision: yes

Circularity Check

0 steps flagged

No circularity: empirical Fourier fit to AAVSO data presented explicitly as data reduction, not as a derived prediction.

full rationale

The manuscript states it fits a low-order Fourier series directly to the phase-folded AAVSO Johnson-V points after adopting an external period and epoch, then selects N=3 as the simplest adequate representation. The reported coefficients (A0, A1, A2, A3, R21, R31, etc.) are the explicit least-squares output; the paper makes no claim that these values are predicted from theory, prior models, or external benchmarks. No self-citations, uniqueness theorems, or ansatzes are invoked to justify the result. The procedure is a transparent data-reduction step whose output is definitionally the fit itself, with no reduction of a claimed derivation to its inputs.

Axiom & Free-Parameter Ledger

4 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumption that a low-order Fourier series adequately captures the periodic signal and on the accuracy of the adopted period, epoch, and data-cleaning choices. No new physical entities are introduced; the free parameters are the fitted coefficients and the discrete choices of order and folding parameters.

free parameters (4)
  • Fourier order N = 3
    Selected as the minimal order that adequately represents the data after checking higher orders improve RMS by only 0.0012 mag.
  • Period Pobs = 5.366531 d
    Adopted value used to phase-fold the observations.
  • Epoch JD = 2460851.395800
    Empirical bright-maximum time used for phase folding.
  • Fourier amplitudes A0-A3 = 3.9031, 0.3434, 0.1428, 0.0531
    Least-squares coefficients that define the adopted template.
axioms (2)
  • domain assumption A low-order Fourier series is sufficient to represent the periodic light variation of a Classical Cepheid.
    Invoked when the authors retain the N=3 model as the preferred representation.
  • domain assumption The AAVSO Johnson-V measurements, after removal of two extreme outliers, faithfully sample the true light curve.
    Basis for using the cleaned sample of 242 points.

pith-pipeline@v0.9.1-grok · 5873 in / 1568 out tokens · 63992 ms · 2026-06-30T01:56:07.271681+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

24 extracted references

  1. [1]

    Leavitt and Edward C

    Henrietta S. Leavitt and Edward C. Pickering. Periods of 25 variable stars in the small magellanic cloud.Harvard College Observatory Circular, 173:1–3, 1912

    1912

  2. [2]

    Madore and Wendy L

    Barry F. Madore and Wendy L. Freedman. The cepheid distance scale.Publications of the Astronomical Society of the Pacific, 103:933–957, 1991

    1991

  3. [3]

    Freedman, Barry F

    Wendy L. Freedman, Barry F. Madore, Brad K. Gibson, Laura Ferrarese, Daniel D. Kelson, Shoko Sakai, Jeremy R. Mould, Robert C. Kennicutt, Holland C. Ford, John A. Graham, John P. Huchra, Shaun M. G. Hughes, Garth D. Illingworth, Lucas M. Macri, and Peter B. Stetson. Final results from the hubble space telescope key project to measure the hubble constant.T...

    2001

  4. [4]

    Fritz Benedict, Barbara E

    G. Fritz Benedict, Barbara E. McArthur, Michael W. Feast, Thomas G. Barnes, Thomas E. Harrison, Richard J. Patterson, John W. Menzies, Jacob L. Bean, and Wendy L. Freedman. Hubble space telescope fine guidance sensor parallaxes of galactic cepheid variable stars: Period–luminosity relations.The Astronomical Journal, 133:1810–1827, 2007

    2007

  5. [5]

    Riess, Lucas M

    Adam G. Riess, Lucas M. Macri, Samantha L. Hoffmann, Dan Scolnic, Stefano Casertano, Alexei V. Filippenko, Brad E. Tucker, Mark J. Reid, David O. Jones, Jeffrey M. Silverman, Ryan Chornock, Peter Challis, Wenlong Yuan, Peter J. Brown, and Ryan J. Foley. A 2.4 percent determination of the local value of the hubble constant.The Astrophysical Journal, 826:56, 2016

    2016

  6. [6]

    Riess, Wenlong Yuan, Lucas M

    Adam G. Riess, Wenlong Yuan, Lucas M. Macri, Dan Scolnic, Dillon Brout, Stefano Casertano, David O. Jones, Yukei Murakami, Gagandeep S. Anand, Louise Breuval, Thomas G. Brink, Alexei V. Filippenko, Samantha Hoffmann, Saurabh W. Jha, W. D’Arcy Kenworthy, John Mackenty, Benjamin E. Stahl, and Weikang Zheng. A comprehensive measurement of the local value of ...

    2022

  7. [7]

    Fritz Benedict, Barbara E

    G. Fritz Benedict, Barbara E. McArthur, Laurence W. Fredrick, Thomas E. Harrison, Cather- ine L. Slesnick, Joseph Rhee, Richard J. Patterson, Michael F. Skrutskie, Otto G. Franz, Lawrence H. Wasserman, William H. Jefferys, Edmund Nelan, William van Altena, Peter J. Shelus, Paul D. Hemenway, Raynor L. Duncombe, D. Story, A. L. Whipple, and Arthur J. Bradle...

    2002

  8. [8]

    Ridgway, Guy Perrin, Michel Rieutord, Pierre Bordé, Theo ten Brummelaar, and Harold A

    Antoine Mérand, Pierre Kervella, Vincent Coudé du Foresto, Stephen T. Ridgway, Guy Perrin, Michel Rieutord, Pierre Bordé, Theo ten Brummelaar, and Harold A. McAlister. The projection factor of delta cephei: A calibration of the baade–wesselink method using the chara array. Astronomy and Astrophysics, 438:L9–L12, 2005

    2005

  9. [9]

    Delta cephei

    American Association of Variable Star Observers. Delta cephei. AAVSO Variable Star of the Season, 2026. URLhttps://www.aavso.org/vsots_delcep. Accessed 2026 June 16

    2026

  10. [10]

    Simon and A

    Norman R. Simon and A. S. Lee. The structural properties of cepheid light curves.The Astrophysical Journal, 248:291–297, 1981

    1981

  11. [11]

    Simon and Thomas J

    Norman R. Simon and Thomas J. Moffett. Classical cepheid light curves revisited.Publications of the Astronomical Society of the Pacific, 97:1078–1085, 1985

    1985

  12. [12]

    Antonello

    E. Antonello. Structural properties of the light curves of s-cepheids.Astronomy and Astrophysics, 169:149–154, 1986

    1986

  13. [13]

    G. K. Andreasen. Precise fourier decomposition parameters for classical cepheids.Astronomy and Astrophysics, 196:159–168, 1988

    1988

  14. [14]

    Kanbur and Chow-Choong Ngeow

    Shashi M. Kanbur and Chow-Choong Ngeow. On the use of principal component analysis in analysing cepheid light curves.Monthly Notices of the Royal Astronomical Society, 329:126–132, 2002

    2002

  15. [15]

    Kanbur, Sergei Nikolaev, Nial R

    Chow-Choong Ngeow, Shashi M. Kanbur, Sergei Nikolaev, Nial R. Tanvir, and Martin A. Hendry. Reconstructing a cepheid light curve with fourier techniques. i. the fourier expansion and interrelations.The Astrophysical Journal, 586:959–982, 2003

    2003

  16. [16]

    Chow-Choong Ngeow and Shashi M. Kanbur. Reconstructing cepheid light curve with fourier techniques. ii. catalogue of the fourier parameters for ogle lmc fundamental mode cepheids, 2004

    2004

  17. [17]

    Tanvir, Martin A

    Nial R. Tanvir, Martin A. Hendry, Alan Watkins, Shashi M. Kanbur, Leonid N. Berdnikov, and Chow-Choong Ngeow. Determination of cepheid parameters by light-curve template fitting. Monthly Notices of the Royal Astronomical Society, 363:749–759, 2005

    2005

  18. [18]

    Sukanta Deb and Harinder P. Singh. Light curve analysis of variable stars using fourier decomposition and principal component analysis.Astronomy and Astrophysics, 507:1729–1737, 2009

    2009

  19. [19]

    Kanbur, Harinder P

    Anupam Bhardwaj, Shashi M. Kanbur, Harinder P. Singh, Lucas M. Macri, and Chow-Choong Ngeow. On the variation of fourier parameters for galactic and large magellanic cloud cepheids at optical, near-infrared and mid-infrared wavelengths.Monthly Notices of the Royal Astronomical Society, 447:3342–3360, 2015. 18

    2015

  20. [20]

    Kanbur, Marcella Marconi, Marina Rejkuba, Harinder P

    Anupam Bhardwaj, Shashi M. Kanbur, Marcella Marconi, Marina Rejkuba, Harinder P. Singh, and Chow-Choong Ngeow. A comparative study of multiwavelength theoretical and observed light curves of cepheid variables.Monthly Notices of the Royal Astronomical Society, 466: 2805–2826, 2017

    2017

  21. [21]

    Szymański, Marcin Kubiak, Grzegorz Pietrzyński, Łukasz Wyrzykowski, Olaf Szewczyk, and Krzysztof Ulaczyk

    Igor Soszyński, Radosław Poleski, Andrzej Udalski, Michał K. Szymański, Marcin Kubiak, Grzegorz Pietrzyński, Łukasz Wyrzykowski, Olaf Szewczyk, and Krzysztof Ulaczyk. The optical gravitational lensing experiment. the ogle-iii catalog of variable stars. i. classical cepheids in the large magellanic cloud.Acta Astronomica, 58:163–185, 2008

    2008

  22. [22]

    Aavso international database, 2026

    American Association of Variable Star Observers. Aavso international database, 2026. URL https://www.aavso.org/aavso-international-database. AAVSO International Database observations of Delta Cephei, accessed 2026 June 16

    2026

  23. [23]

    Aavso data usage guidelines, 2025

    American Association of Variable Star Observers. Aavso data usage guidelines, 2025. URL https://www.aavso.org/data-usage-guidelines. Accessed 2026 June 16

    2025

  24. [24]

    David G. Turner. Visual observations of delta cephei: Time to update the finder chart.Journal of the American Association of Variable Star Observers, 2011. URLhttps://apps.aavso. org/jaavso/article/2711/. 19

    2011