Abstract:Stellar evolution theory predicts that asymptotic giant branch (AGB) stars undergo a series of short thermal pulses that significantly change their luminosity and mass on timescales of hundreds to thousands of years. These pulses are confirmed observationally by the existence of the short-lived radioisotope technetium in the spectra of some of these stars, but other observational consequences of thermal pulses are subtle and may only be detected over many years of observations. Secular changes in these stars resulting from thermal pulses can be detected as measurable changes in period if the star is undergoing Mira pulsations. It is known that a small fraction of Mira variables exhibit large secular period changes, and the detection of these changes among a larger sample of stars could therefore be useful in evolutionary studies of these stars. The American Association of Variable Star Observers (AAVSO) International Database currently contains visual data for over 1500 Mira variables. Light curves for these stars span nearly a century in some cases, making it possible to study the secular evolution of the pulsation behavior on these timescales. In this paper we present the results of our study of period change in 547 Mira variables using data from the AAVSO. We use wavelet analysis to measure the period changes in individual Mira stars over the span of available data. By making linear fits to the period versus time measurements, we determine the average rates of period change, d ln P/dt, for each of these stars. We find nonzero d ln P/dt at the 2 σ significance level in 57 of the 547 stars, at the 3 σ level in 21 stars, and at the level of 6 σ or greater in eight stars. The latter eight stars have been previously noted in the literature, and our derived rates of period change largely agree with published values. The largest and most statistically significant d ln P/dt are consistent with the rates of period change expected during thermal pulses on the AGB. A number of other stars exhibit nonmonotonic period change on decades-long timescales, the cause of which is not yet known. In the majority of stars, the period variations are smaller than our detection threshold, meaning the available data are not sufficient to unambiguously measure slow evolutionary changes in the pulsation period. It is unlikely that more stars with large period changes will be found among heretofore well-observed Mira stars in the short term, but continued monitoring of these and other Mira stars may reveal new and serendipitous candidates in the future.

The Mira Distance Ladder

The Mira Distance to M101 and a 4% Measurement of H0

The period--luminosity relation for Mira variables in the Milky Way using Gaia DR3: a further distance anchor for $H_0$

The Mira Distance to M101 and a 4% Measurement of H 0

Simultaneous inference of periods and period-luminosity relations for Mira variable stars

Near-infrared narrow-band photometry of M-giant and Mira stars: models meet observations

Models of Mira variables of the Large Magellanic Cloud

A BayeSN Distance Ladder: $H_0$ from a consistent modelling of Type Ia supernovae from the optical to the near infrared

The impact of third dredge-up on the mass loss of Mira variables

Interplay between pulsation, mass loss, and third dredge-up: More about Miras with and without technetium

Period -- mass-loss rate relation of Miras with and without technetium

The M33 Synoptic Stellar Survey. III. Miras and LPVs in griJHKS

Mira Variables Based on Emission Lines in LAMOST DR9

Unveiling VVV/WISE Mira variables on the far side of the Galactic disk: Distances, kinematics and a new extinction law

A MultiWavelength Study of the Symbiotic Mira HM Sge with SOFIA and HST

A Very Long Baseline Array Trigonometric Parallax for RR Aql and the Mira Period–Luminosity Relation

Secular Evolution in Mira Variable Pulsations

Connecting Infrared Surface Brightness Fluctuation Distances to Type Ia Supernova Hosts: Testing the Top Rung of the Distance Ladder