Modelling stars with Gaussian Process Regression: Augmenting Stellar Model Grid
Authors:
Tanda Li,
Guy R. Davies,
Alexander J. Lyttle,
Warrick H. Ball,
Lindsey M. Carboneau,
Rafael A. Garcia
Abstract:
Grid-based modelling is widely used for estimating stellar parameters. However, stellar model grid is sparse because of the computational cost. This paper demonstrates an application of a machine-learning algorithm using the Gaussian Process (GP) Regression that turns a sparse model grid onto a continuous function. We train GP models to map five fundamental inputs (mass, equivalent evolutionary ph…
▽ More
Grid-based modelling is widely used for estimating stellar parameters. However, stellar model grid is sparse because of the computational cost. This paper demonstrates an application of a machine-learning algorithm using the Gaussian Process (GP) Regression that turns a sparse model grid onto a continuous function. We train GP models to map five fundamental inputs (mass, equivalent evolutionary phase, initial metallicity, initial helium fraction, and the mixing-length parameter) to observable outputs (effective temperature, surface gravity, radius, surface metallicity, and stellar age). We test the GP predictions for the five outputs using off-grid stellar models and find no obvious systematic offsets, indicating good accuracy in predictions.As a further validation, we apply these GP models to characterise 1,000 fake stars. Inferred masses and ages determined with GP models well recover true values within one standard deviation. An important consequence of using GP-based interpolation is that stellar ages are more precise than those estimated with the original sparse grid because of the full sampling of fundamental inputs.
△ Less
Submitted 16 February, 2022;
originally announced February 2022.
Hierarchically modelling Kepler dwarfs and subgiants to improve inference of stellar properties with asteroseismology
Authors:
Alexander J. Lyttle,
Guy R. Davies,
Tanda Li,
Lindsey M. Carboneau,
Ho-Hin Leung,
Harry Westwood,
William J. Chaplin,
Oliver J. Hall,
Daniel Huber,
Martin B. Nielsen,
Sarbani Basu,
Rafael A. García
Abstract:
With recent advances in modelling stars using high-precision asteroseismology, the systematic effects associated with our assumptions of stellar helium abundance ($Y$) and the mixing-length theory parameter ($α_\mathrm{MLT}$) are becoming more important. We apply a new method to improve the inference of stellar parameters for a sample of Kepler dwarfs and subgiants across a narrow mass range (…
▽ More
With recent advances in modelling stars using high-precision asteroseismology, the systematic effects associated with our assumptions of stellar helium abundance ($Y$) and the mixing-length theory parameter ($α_\mathrm{MLT}$) are becoming more important. We apply a new method to improve the inference of stellar parameters for a sample of Kepler dwarfs and subgiants across a narrow mass range ($0.8 < M < 1.2\,\mathrm{M_\odot}$). In this method, we include a statistical treatment of $Y$ and the $α_\mathrm{MLT}$. We develop a hierarchical Bayesian model to encode information about the distribution of $Y$ and $α_\mathrm{MLT}$ in the population, fitting a linear helium enrichment law including an intrinsic spread around this relation and normal distribution in $α_\mathrm{MLT}$. We test various levels of pooling parameters, with and without solar data as a calibrator. When including the Sun as a star, we find the gradient for the enrichment law, $ΔY / ΔZ = 1.05^{+0.28}_{-0.25}$ and the mean $α_\mathrm{MLT}$ in the population, $μ_α= 1.90^{+0.10}_{-0.09}$. While accounting for the uncertainty in $Y$ and $α_\mathrm{MLT}$, we are still able to report statistical uncertainties of 2.5 per cent in mass, 1.2 per cent in radius, and 12 per cent in age. Our method can also be applied to larger samples which will lead to improved constraints on both the population level inference and the star-by-star fundamental parameters.
△ Less
Submitted 11 June, 2021; v1 submitted 10 May, 2021;
originally announced May 2021.
