Dramatically improved data from observatories like the CoRoT and Kepler spacecraft have recently facilitated nonlinear time series analysis and phenomenological modeling of variable stars, including the search for strange (aka fractal) or chaotic dynamics. We recently argued [Lindner et al., Phys. Rev. Lett. 114 (2015) 054101] that the Kepler data includes "golden" stars, whose luminosities vary quasiperiodically with two frequencies nearly in the golden ratio, and whose secondary frequencies exhibit power-law scaling with exponent near -1.5, suggesting strange nonchaotic dynamics and singular spectra. Here we use a series of phenomenological models to make plausible the connection between golden stars and fractal spectra. We thereby suggest that at least some features of variable star dynamics reflect universal nonlinear phenomena common to even simple systems.

The unprecedented light curves of the Kepler space telescope document how the brightness of some stars pulsates at primary and secondary frequencies whose ratios are near the golden mean, the most irrational number. A nonlinear dynamical system driven by an irrational ratio of frequencies generically exhibits a strange but nonchaotic attractor. For Kepler's "golden" stars, we present evidence of the first observation of strange nonchaotic dynamics in nature outside the laboratory. This discovery could aid the classification and detailed modeling of variable stars.

Learned et. al. proposed that a sufficiently advanced extra-terrestrial civilization may tickle Cepheid and RR Lyrae variable stars with a neutrino beam at the right time, thus causing them to trigger early and jogging the otherwise very regular phase of their expansion and contraction. This would turn these stars into beacons to transmit information throughout the galaxy and beyond. The idea is to search for signs of phase modulation (in the regime of short pulse duration) and patterns, which could be indicative of intentional, omnidirectional signaling. We have performed such a search among variable stars using photometric data from the Kepler space telescope. In the RRc Lyrae star KIC 5520878, we have found two such regimes of long and short pulse durations. The sequence of period lengths, expressed as time series data, is strongly auto correlated, with correlation coefficients of prime numbers being significantly higher ($p=99.8$\%). Our analysis of this candidate star shows that the prime number oddity originates from two simultaneous pulsation periods and is likely of natural origin. Simple physical models elucidate the frequency content and asymmetries of the KIC 5520878 light curve. Despite this SETI null result, we encourage testing other archival and future time-series photometry for signs of modulated stars. This can be done as a by-product to the standard analysis, and even partly automated.