Hong-Yi Zhang, Heng Bian, Zhao-Yu Zuo

Magnetic white dwarfs can convert photons into axions in their strong magnetic fields, with the conversion probability modulating the light curve as the star rotates. However, this observable is degenerate with intrinsic stellar variability if the background is modeled too simplistically. We develop a controlled background-degeneracy framework that computes the axion-induced modulation from reconstructed stellar magnetic fields while fitting it simultaneously with a flexible Fourier model of the intrinsic light curve. Applying this framework to TESS observations of PG~1015+014 using two independent magnetic field reconstructions, we find that a sinusoidal stellar background can produce an apparent preference for nonzero axion-photon conversion. This preference is absorbed once the background light curve includes the second harmonic, indicating that higher-harmonic stellar variability is a leading degeneracy for precise photometric axion searches in magnetic white dwarfs. Interpreting the two-harmonic fit as a conservative baseline, we obtain competitive constraints for sub-$μ\mathrm{eV}$ axions. We further derive an analytic target-ranking estimate for other TESS magnetic white dwarfs, identifying systems where phase-resolved magnetic modeling would be most valuable for competitive axion probes.