Mira variables in the Milky Way's nuclear stellar disc: discovery and classification

TL;DR

This study identifies and classifies 1782 Mira variable candidates in the Milky Way's nuclear stellar disc using multi-epoch infrared data, revealing insights into its star formation history and stellar properties.

Contribution

It presents a new large sample of Mira variables in the Galactic center, employing advanced light curve modeling and multi-band photometry to analyze their properties and implications.

Findings

Majority are consistent with O-rich chemistry.

Most short-period stars are at the Galactic Centre distance.

Many long-period variables are dusty with high mass-loss rates.

Abstract

The properties of the Milky Way's nuclear stellar disc give crucial information on the epoch of bar formation. Mira variables are promising bright candidates to study the nuclear stellar disc, and through their period-age relation dissect its star formation history. We report on a sample of $1782$ Mira variable candidates across the central $3\times3\,\mathrm{deg}^2$ of the Galaxy using the multi-epoch infrared VISTA Variables in Via Lactea (VVV) survey. We describe the algorithms employed to select candidate variable stars and then model their light curves using periodogram and Gaussian process methods. By combining with WISE, 2MASS and other archival photometry, we model the multi-band light curves to refine the periods and inspect the amplitude variation between different photometric bands. The infrared brightness of the Mira variables means many are too bright and missed by VVV. However, our sample follows a well-defined selection function as expected from artificial star tests. The multi-band photometry is modelled using stellar models with circumstellar dust that characterise the mass loss rates. We demonstrate how $\gtrsim90$ per cent of our sample is consistent with O-rich chemistry. Comparison to period-luminosity relations demonstrates that the bulk of the short period stars are situated at the Galactic Centre distance. Many of the longer period variables are very dusty, falling significantly under the O-rich Magellanic Cloud and solar neighbourhood period-luminosity relations and exhibit high mass-loss rates of $\sim2.5\times10^{-5}M_\odot\,\mathrm{yr}^{-1}$. The period distribution appears consistent with the nuclear stellar disc forming $\gtrsim8\,\mathrm{Gyr}$ ago although it is not possible to disentangle the relative contributions of the nuclear stellar disc and the contaminating bulge.