Hydrodynamic models of pulsation period evolution in R Hydrae

Abstract

Pulsation period evolution during the helium-shell flash in the Mira variable R Hya is investigated using consistent stellar evolution and non-linear stellar pulsation computations. The initial and time-dependent inner boundary conditions for the equations of radiation hydrodynamics describing non-linear stellar oscillations were determined using a grid of TP-AGB model sequences with initial masses on the main sequence $1.5M_\odot\le M_\mathrm{ZAMS}\le 5.0M_\odot$ and the initial metallicity $Z=0.014$. The setup of initial conditions for hydrodynamic models corresponds to $\approx 100$ yr prior to the maximum of the helium-shell luminosity and ensures that the stellar envelope of the evolution model is under both hydrostatic and thermal equilibrium. Solution of the equations of hydrodynamics allowed us to determine the temporal variation of the pulsation period $\Pi(t)$ during $\approx 500$~yr. Within this time interval R Hya is a fundamental mode pulsator. The period temporal dependencies $\Pi(t)$ calculated for the AGB star models at the beginning of the third dredge-up phase and with masses $4.4M_\odot\le M\le 4.5M_\odot$ are in agreement with observational estimates of the period of R Hya obtained during last two centuries. The mean radius of R Hya pulsation models at the end of the XX century ($470 R_\odot < \bar{R} < 490 R_\odot$) agrees with observational estimates obtained using the interferometric angular diameter measurements.