Projection-factor systematics and the long-period calibration of the Galactic Cepheid period-luminosity relation from Baade-Wesselink-type (maximum-likelihood) distances

The maximum-likelihood (ML) calibration of the Galactic Cepheid period-luminosity (PL) relation by Lazovik & Rastorguev (2020), a Baade-Wesselink-type parallax-of-pulsation analysis of 44 Cepheids with multiphase effective-temperature measurements, adopted the empirical projection-factor law $p = 1.376 - 0.064\log P$ to convert observed radial velocities into pulsation velocities, and identified its dominant open systematic: the true period dependence of the projection factor, its possible cycle-to-cycle and phase-dependent variations, and their basis in atmospheric velocity gradients and line-formation physics are not quantitatively settled. The same study left the long-period ($P > 20$ days) end of the Galactic period-radius and PL relations constrained by only five calibrators, the recalibrated $\log T_{\rm eff}$-color and bolometric-correction relations unvalidated against phase dependence, the reddening system (including spatially varying extinction law $R_V$) in disagreement with an independent compilation, and the extension of the technique to Magellanic Cloud Cepheids undemonstrated. Later literature deepened rather than settled the p-factor question: infrared-surface-brightness calibrations find a steep period dependence ($p = 1.550 - 0.186\log P$), SPIPS (spectro-photo-interferometry of pulsating stars, a global-modeling implementation of the parallax of pulsation) modeling of 63 Cepheids with Gaia EDR3 parallaxes finds no period dependence at all ($p = 1.26 \pm 0.07$, dispersion about 0.15), and 2025 nonlinear pulsation-model fitting of Gaia DR3 Cepheids finds a constant $p = 1.22 \pm 0.05$. This problem, a physically grounded projection-factor prescription and a systematics-controlled, long-period-anchored calibration of the Cepheid distance scale from pulsation-based distances, remains open.

Background

Classical Cepheids are radially pulsating yellow supergiant stars whose mean absolute magnitude is a repeatable function of the pulsation period $P$, the period-luminosity (PL) or Leavitt relation. Because $P$ is easily measured while the luminosity is high, Cepheids calibrate the first rung of the extragalactic distance ladder and hence the local value of the Hubble constant $H_0$. One route to Cepheid luminosities that is independent of trigonometric parallaxes is the Baade-Wesselink (or parallax-of-pulsation) method: the star's changing angular size (inferred from photometry through color-temperature and bolometric-correction calibrations, or measured interferometrically) is compared with its changing linear radius, obtained by integrating the pulsation velocity over the pulsation cycle. The observable is the radial velocity $V_r$ from spectroscopic line shifts; converting it to the true pulsation velocity of the stellar surface requires the projection factor $p$, a number of order 1.2-1.5 that accounts for the projection of the surface velocity field onto the line of sight across the stellar disk and for the dynamics and line-formation physics of the extended pulsating atmosphere. It is conventionally decomposed into a geometric term $p_0$ (set by limb darkening), the atmospheric velocity gradient $f_{\rm grad}$ (spectral lines form at different atmospheric depths, which move at different velocities), and the offset $f_{o-g}$ between the optical photosphere and the corresponding mass elements. The interstellar reddening $E(B-V)$ (color excess due to dust along the line of sight) and the total-to-selective extinction ratio $R_V$ enter the same calibration because they set the intrinsic colors and magnitudes.

Lazovik & Rastorguev (2020) realized a maximum-likelihood (ML) version of the Baade-Wesselink method, building on the formulation introduced by Balona (1977) and its modernization by Rastorguev & Dambis (2011). The ML technique fits light, color, radial-velocity, and multiphase effective-temperature $T_{\rm eff}$ data simultaneously for the stellar radius, reddening, distance modulus, and absolute magnitude, using nonlinear calibrations of $\log T_{\rm eff}$ from the intrinsic color $(B-V)_0$ and of the bolometric correction that account for metallicity $[Fe/H]$ (heavy-element abundance) and surface gravity $\log g$. Applied to 44 Galactic Cepheids with multiphase spectroscopic temperatures, after restricting the pulsation phases to avoid data contaminated by atmospheric shock waves, it yielded the period-radius relation $\log R = 0.68\log P + 1.14$ and the visual-band PL relation $M_V = -2.67(\log P - 1) - 4.14$, adopting the then-standard empirical projection-factor law $p = 1.376 - 0.064\log P$ of Nardetto et al. (2007). A 0.09 mag zero-point offset relative to the Hubble Space Telescope (HST) parallax calibration was argued to correspond to a shift of $H_0$ from about 74 to about 71 km/s/Mpc, illustrating that this systematic is cosmologically consequential.

The projection factor was already contested before 2020. The infrared surface-brightness (IRSB) variant of the Baade-Wesselink method, which derives angular diameters from near-infrared colors via a surface-brightness-color relation (SBCR), yields a steep period dependence, $p = 1.550 - 0.186\log P$ (Storm et al. 2011), reinforced by the recalibration of Ngeow et al. (2012), $p = 1.462 - 0.172\log P$ with large scatter. Direct spectroscopy of the prototype delta Cep with HARPS-N (a high-resolution optical spectrograph on the Telescopio Nazionale Galileo) gave $p = 1.239 \pm 0.034 \pm 0.023$ for Gaussian line profiles (Nardetto et al. 2017), well below the Nardetto et al. (2007) law at that period, and SPIPS (spectro-photo-interferometry of pulsating stars, a global-modeling implementation of the parallax of pulsation) applied to Magellanic Cloud Cepheids gave the shallow relation $p = -0.08(\log P - 1.18) + 1.24$ with no detectable metallicity dependence (Gallenne et al. 2017). Theoretical predictions from hydrodynamical and static spherical model atmospheres are mutually inconsistent with these observational constraints (Neilson et al. 2012). After 2020, SPIPS modeling of 63 Galactic Cepheids with parallaxes from the Gaia mission's early third data release (EDR3; Trahin et al. 2021) and nonlinear convective pulsation-model fitting of 46 Cepheids with Gaia DR3 photometry and radial velocities (Molinaro et al. 2025) both find no significant period dependence (mean $p \approx 1.22$-1.26) but a large unexplained dispersion, while dedicated HARPS-N studies quantify strong method systematics: the choice of SBCR alone shifts $p$ by about 8%, and the method of deriving the radial-velocity curve by about 9% (Nardetto et al. 2023). For long-period Cepheids, circumstellar envelopes introduce further bias in the Baade-Wesselink quantities (Hocde et al. 2024, modeling Y Ophiuchi). Lazovik & Rastorguev (2020) themselves list what remains unresolved: the projection-factor law beyond the adopted empirical relation; the sparse long-period coverage (five calibrators with $P > 20$ days); the phase dependence and $\log g$/$[Fe/H]$ assumptions of the $\log T_{\rm eff}$ and bolometric-correction calibrations; the reddening-system disagreement with the compilation of Turner (2016), including metallicity-dependent and $R_V$-variable extinction; and the untested extension of the ML technique to Large and Small Magellanic Cloud Cepheids with multi-band photometry.

Problem Statement

Establish, and quantitatively justify, the projection-factor prescription used in Baade-Wesselink-type Cepheid distance determinations, and use it to complete a systematics-controlled calibration of the Galactic Cepheid period-radius and period-luminosity relations. Concretely, as left open by Lazovik & Rastorguev (2020): (1) determine the period dependence of the projection factor $p$ (or demonstrate its absence within quantified bounds), together with its possible cycle-to-cycle and phase-dependent variations and its physical origin in atmospheric velocity gradients and line-formation physics, so that the adoption of a contested empirical law such as $p = 1.376 - 0.064\log P$ no longer dominates the systematic error budget; (2) anchor the long-period end ($P > 20$ days) of the Galactic period-radius and PL relations with an adequate calibrator sample rather than extrapolation; (3) validate the $\log T_{\rm eff}((B-V)_0)$ and bolometric-correction calibrations, including their phase dependence and the assumption of constant $\log g$ and $[Fe/H]$ over the pulsation cycle; (4) resolve the reddening-system discrepancy between Cepheid-based and early-type reference-star-based determinations, including possible metallicity dependence and spatially varying extinction law $R_V$; and (5) demonstrate the pulsation-based calibration on Magellanic Cloud Cepheids with multi-band photometry as an extragalactic test. The central problem is (1) taken together with (2); (3)-(5) are related extensions that strengthen the calibration. Anchoring the long-period end requires pulsation-based Galactic calibrators at $P > 20$ days that span the long-period range and determine the long-period slope and zero point, rather than extrapolating the shorter-period relation. The objective is a pulsation-based Cepheid distance calibration whose zero point and systematics are controlled well enough to be compared with, and trusted alongside, purely geometric calibrations of the distance ladder. Determine, with uncertainties, whether the p-factor prescription changes the previously published Galactic PL zero point, and test the result against independent distances.

References

1. Lazovik, Y.; Rastorguev, A. (2020), 'Calibrating the Galactic Cepheid Period-Luminosity Relation from the Maximum-likelihood Technique', The Astronomical Journal 160, 136, DOI: 10.3847/1538-3881/aba627, https://doi.org/10.3847/1538-3881/aba627
2. Balona, L. A. (1977), 'Application of the Method of Maximum Likelihood to the Determination of Cepheid Radii', Monthly Notices of the Royal Astronomical Society 178, 231-243, DOI: 10.1093/mnras/178.2.231, https://doi.org/10.1093/mnras/178.2.231
3. Rastorguev, A.; Dambis, A. (2011), 'Classical Cepheids: a New version of the Baade-Becker-Wesselink method', Astrophysical Bulletin 66, 47-53, DOI: 10.1134/S1990341311010032, https://doi.org/10.1134/S1990341311010032
4. Nardetto, N.; Mourard, D.; Mathias, P.; Fokin, A.; Gillet, D. (2007), 'High-resolution spectroscopy for Cepheids distance determination', Astronomy & Astrophysics 471, 661-669, DOI: 10.1051/0004-6361:20066853, https://doi.org/10.1051/0004-6361:20066853
5. Turner, D. G. (2016), 'The Scale of Reddening for Classical Cepheid Variables', Revista Mexicana de Astronomía y Astrofísica 52, 223, arXiv:1603.02276, https://arxiv.org/abs/1603.02276
6. Storm, J.; Gieren, W.; Fouqué, P.; Barnes, T. G.; Pietrzyński, G.; Nardetto, N.; Weber, M.; Granzer, T.; Strassmeier, K. G. (2011), 'Calibrating the Cepheid period-luminosity relation from the infrared surface brightness technique', Astronomy & Astrophysics 534, A94, DOI: 10.1051/0004-6361/201117155, https://doi.org/10.1051/0004-6361/201117155
7. Ngeow, C.; Neilson, H.; Nardetto, N.; Marengo, M. (2012), 'Calibrating the projection factor for Galactic Cepheids', Astronomy & Astrophysics 543, A55, DOI: 10.1051/0004-6361/201218780, https://doi.org/10.1051/0004-6361/201218780
8. Neilson, H.; Nardetto, N.; Ngeow, C.; Fouqué, P.; Storm, J. (2012), 'Cepheid limb darkening, angular diameter corrections, and projection factor from static spherical model stellar atmospheres', Astronomy & Astrophysics 541, A134, DOI: 10.1051/0004-6361/201118550, https://doi.org/10.1051/0004-6361/201118550
9. Gallenne, A.; Kervella, P.; Mérand, A.; Pietrzyński, G.; Gieren, W.; Nardetto, N.; Trahin, B. (2017), 'Observational calibration of the projection factor of Cepheids', Astronomy & Astrophysics 608, A18, DOI: 10.1051/0004-6361/201731589, https://doi.org/10.1051/0004-6361/201731589 (arXiv title: '...IV. Period-projection factor relation of Galactic and Magellanic Cloud Cepheids')
10. Nardetto, N.; Poretti, E.; Rainer, M.; Fokin, A.; Mathias, P.; Anderson, R. I.; Gallenne, A.; Gieren, W.; Graczyk, D.; Kervella, P.; Mérand, A.; Mourard, D.; Neilson, H.; Pietrzyński, G.; Pilecki, B.; Storm, J. (2017), 'HARPS-N high spectral resolution observations of Cepheids I. The Baade-Wesselink projection factor of delta Cep revisited', Astronomy & Astrophysics 597, A73, DOI: 10.1051/0004-6361/201629400, https://doi.org/10.1051/0004-6361/201629400
11. Trahin, B.; Breuval, L.; Kervella, P.; Mérand, A.; Nardetto, N.; Gallenne, A.; Hocdé, V.; Gieren, W. (2021), 'Inspecting the Cepheid parallax of pulsation using Gaia EDR3 parallaxes', Astronomy & Astrophysics 656, A102, DOI: 10.1051/0004-6361/202141680, https://doi.org/10.1051/0004-6361/202141680 (arXiv title adds 'Projection factor and period-luminosity and period-radius relations')
12. Nardetto, N.; Gieren, W.; Storm, J.; Hocdé, V.; Pietrzyński, G.; Kervella, P.; Mérand, A.; Gallenne, A.; Graczyk, D.; Pilecki, B.; Poretti, E.; Rainer, M.; Zgirski, B.; Wielgórski, P.; Hajdu, G.; Górski, M.; Karczmarek, P.; Narloch, W.; Taormina, M. (2023), 'HARPS-N high spectral resolution observations of Cepheids II. The impact of the surface-brightness color relation on the Baade-Wesselink projection factor of η Aql', Astronomy & Astrophysics 671, A14, DOI: 10.1051/0004-6361/202245298, https://doi.org/10.1051/0004-6361/202245298
13. Hocdé, V.; Smolec, R.; Moskalik, P.; Singh Rathour, R.; Ziółkowska, O. (2024), 'Pulsation modeling of the Cepheid Y Ophiuchi with RSP/MESA', Astronomy & Astrophysics 683, A233, DOI: 10.1051/0004-6361/202348428, https://doi.org/10.1051/0004-6361/202348428 (arXiv title adds 'Impact of the circumstellar envelope and a high projection factor on Baade-Wesselink method')
14. Molinaro, R.; Marconi, M.; De Somma, G.; Ripepi, V.; Leccia, S.; Musella, I.; Sicignano, T.; Trentin, E.; Gatto, M. (2025), 'The model fitting of Gaia DR3 classical Cepheid light and radial velocity curves', Astronomy & Astrophysics 700, A212, DOI: 10.1051/0004-6361/202553855, https://doi.org/10.1051/0004-6361/202553855
15. Wielgórski, P.; Pietrzyński, G.; Gieren, W.; Zgirski, B.; Narloch, W.; Hajdu, G.; Storm, J.; Nardetto, N.; Kervella, P.; Pilecki, B.; Górski, M.; Smolec, R.; Salinas, R.; Graczyk, D.; Hocdé, V.; Karczmarek, P.; Taormina, M.; Pych, W.; Netzel, H.; Chini, R.; Hodapp, K.; Kałuszyński, M.; Pozo Nuñez, F.; Kotysz, K.; Moździerski, P.; Mikołajczyk, P.; Kołaczek-Szymański, P. (2026), 'Period-Luminosity Relations, projection factor and radii of Anomalous Cepheids', arXiv:2603.28472, https://arxiv.org/abs/2603.28472