Variable Star Light Curves & Phase Folding

Physical & Observational Principles

Because ground-based observations are interrupted by daylight, weather, and seasonal solar conjunctions, raw brightness measurements m(t_i) of a periodic variable star are sparsely distributed across hundreds or thousands of cycles. Once period-search algorithms—such as the Lomb–Scargle periodogram for unevenly spaced sinusoidal data, Phase Dispersion Minimization (PDM), or Box Least Squares (BLS) for flat-bottomed eclipses—identify the fundamental period P, astronomers fold all timestamps onto a unit phase interval φ = ((t − T₀) / P) − floor((t − T₀) / P). By convention in the General Catalogue of Variable Stars (GCVS), the zero-phase epoch T₀ is defined at primary maximum light for pulsating stars (Cepheids, RR Lyrae, Miras) and at primary mid-eclipse minimum light for eclipsing binaries (Algol, Beta Lyrae, W Ursae Majoris systems).

The resulting phase-folded light curve acts as a physical fingerprint distinguishing extrinsic geometric variables from intrinsic radial pulsators. In detached Algol-type (EA) eclipsing binaries like Beta Persei (Algol, P = 2.867328 days), brightness stays nearly constant across most of the orbit before plunging into a sharp primary minimum when the cooler K-type subgiant occults the hot B8 V primary, followed at φ = 0.50 (for circular orbits) by a shallow secondary minimum. In contrast, classical δ Cephei pulsators (driven by the κ-mechanism opacity valve in the doubly ionized helium He II zone) exhibit a rapid, steep rise to maximum light followed by a slower, asymmetric decline spanning ΔV ≈ 0.6–1.2 magnitudes.

In asymptotic giant branch Mira variables (like Omicron Ceti, P ≈ 331.96 days), pulsation coupled with titanium oxide (TiO) molecular band formation—which dims optical V light heavily at minimum temperature—produces visual amplitudes of ΔV = 6 to 9 magnitudes even though bolometric infrared luminosity varies by only ~1 magnitude. Tracking Observed minus Calculated maxima or minima (O − C = T_obs − (T₀ + E × P)) across cycle count E directly measures secular stellar evolution and mass transfer rates (dP/dt).

Worked Calculation: Phase Folding an Observation of Delta Cephei and Applying Leavitt's Law

1. Delta Cephei has GCVS elements T₀ = BJD 2,436,075.445 (maximum light) and pulsation period P = 5.366270 days.
2. Suppose we observe Delta Cephei at t = BJD 2,461,000.000: elapsed time is Δt = 2,461,000.000 − 2,436,075.445 = 24,924.555 days.
3. Divide by P: Δt / P = 24,924.555 / 5.366270 = 4,644.6703 cycles ⇒ cycle integer E = 4,644 and fractional pulsation phase φ = 0.670 (on the descending branch, ~67% of the way from maximum toward minimum).
4. Using the Milky Way Cepheid Period–Luminosity relation M_V = −2.43(log₁₀(5.3663) − 1) − 4.05 = −2.43(−0.2703) − 4.05 = −3.39 mag, which with mean <V> = +3.95 and A_V ≈ 0.23 gives distance modulus μ_0 = 7.11 mag (d ≈ 265–295 pc).

Common Systematic Pitfalls & Misconceptions

• Confusing a 1-day sidereal sampling alias (f_alias = f_true ± 1.00274 day⁻¹) in ground-based periodograms with the star's true physical pulsation frequency.
• Folding an eclipsing binary light curve at half its true orbital period (P/2) when primary and secondary eclipses have nearly equal depths, which superimposes odd and even eclipses on top of each other.

Benchmark Celestial Objects & Interactive Calculators

Authoritative Standards & Sources

• General Catalogue of Variable Stars (GCVS) & AAVSO VSX (v5.1 (2026 update)): Samus, N. N., et al. 2017, Astronomy Reports, 61, 80 (General Catalogue of Variable Stars: Version GCVS 5.1) — Sternberg Astronomical Institute / AAVSO
• ESA Gaia Archive — Data Release 3 (DR3 (2022-06-13)): Gaia Collaboration, Vallenari, A., et al. 2023, A&A, 674, A1 (Gaia Data Release 3) — European Space Agency / Gaia DPAC