skip to main content

Title: Modelling stars with Gaussian Process Regression: augmenting stellar model grid

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 on to 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 characterize 1000 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.

; ; ; ; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 5597-5610
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. Magnetic fields influence the formation and evolution of stars and impact the observed stellar properties. magnetic A-type stars (Ap stars) are a prime example of this. Access to precise and accurate determinations of their stellar fundamental properties, such as masses and ages, is crucial to understand the origin and evolution of fossil magnetic fields. Aims. We propose using the radii and luminosities determined from interferometric measurements, in addition to seismic constraints when available, to infer fundamental properties of 14 Ap stars préviously characterised. Methods. We used a grid-based modelling approach, employing stellar models computed with the CESTAM stellar evolution code, and the parameter search performed with the AIMS optimisation method. The stellar model grid was built using a wide range of initial helium abundances and metallicities in order to avoid any bias originating from the initial chemical composition. The large frequency separations (Δ ν ) of HR 1217 (HD 24712) and α Cir (HD 128898), two rapidly oscillating Ap stars of the sample, were used as seismic constraints. Results. We inferred the fundamental properties of the 14 stars in the sample. The overall results are consistent within 1 σ with previous studies, however, the stellar masses inferred in thismore »study are higher. This trend likely originates from the broader range of chemical compositions considered in this work. We show that the use of Δ ν in the modelling significantly improves our inferences, allowing us to set reasonable constraints on the initial metallicity which is, otherwise, unconstrained. This gives an indication of the efficiency of atomic diffusion in the atmospheres of roAp stars and opens the possibility of characterising the transport of chemical elements in their interiors.« less

    Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the mesa software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3–60 M⊙), surface equatorial magnetic field strength (0–50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly rotating, nitrogen-enriched (‘Group 2’) stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this gridmore »of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of (i) revisiting the derived stellar parameters of known magnetic stars, and (ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes.

    « less

    We calculate the α-enhancement ratio [α/Fe] for the Mapping Nearby Galaxies at APO (MaNGA) Stellar Library (MaStar) while also fitting for the fundamental atmospheric parameters effective temperature, surface gravity, and metallicity – Teff, log g, [Fe/H]. This approach builds upon a previous catalogue of stellar parameters, whereby only the fundamental atmospheric parameters are fit with solar-scaled models. Here, we use the same Markov Chain Monte Carlo method with the additional free parameter [α/Fe]. Using the full spectral fitting code pPXF, we are able to fit multiple lines sensitive to [α/Fe] for a more robust measurement. Quality flags based on the convergence of the sampler, errors in [α/Fe] and a cut in the χ2 of the model fit are used to clean the final catalogue, returning 17 214 spectra and values in the range of −0.25 < [α/Fe] < 0.48. Comparing our calculated [α/Fe] with literature values reveals a degeneracy in cool stars with log g ≥ ∼4; this comparison is then used to create an alternative and calibrated parameter set. We also plot the final catalogue in an [Fe/H] versus [α/Fe] diagram and recover the expected result of increasing [α/Fe] with decreasing [Fe/H] for Milky Way disc-halo stars. We applymore »our method to a subsample of spectra of uniform resolution and higher signal to noise that finds that our results are independent of this higher signal to noise. In the context of stellar population models, we are able to cover a parameter space for the creation of intermediate to old age models at solar-scaled [α/Fe], high [Fe/H] and enhanced [α/Fe], low [Fe/H].

    « less
  4. Abstract

    Bayesian optimization (BO) has been leveraged for guiding autonomous and high-throughput experiments in materials science. However, few have evaluated the efficiency of BO across a broad range of experimental materials domains. In this work, we quantify the performance of BO with a collection of surrogate model and acquisition function pairs across five diverse experimental materials systems. By defining acceleration and enhancement metrics for materials optimization objectives, we find that surrogate models such as Gaussian Process (GP) with anisotropic kernels and Random Forest (RF) have comparable performance in BO, and both outperform the commonly used GP with isotropic kernels. GP with anisotropic kernels has demonstrated the most robustness, yet RF is a close alternative and warrants more consideration because it is free from distribution assumptions, has smaller time complexity, and requires less effort in initial hyperparameter selection. We also raise awareness about the benefits of using GP with anisotropic kernels in future materials optimization campaigns.

  5. Abstract Accurate stellar ages are essential for our understanding of the star formation history of the Milky Way and Galactic chemical evolution, as well as to constrain exoplanet formation models. Gyrochronology, a relationship between stellar rotation and age, appears to offer a reliable age indicator for main-sequence (MS) stars over the mass range of approximately 0.6–1.3 M ⊙ . Those stars lose their angular momentum due to magnetic braking and as a result their rotation speeds decrease with age. Although current gyrochronology relations have been fairly well tested for young MS stars with masses greater than 1 M ⊙ , primarily in young open clusters, insufficient tests exist for older and lower mass MS stars. Binary stars offer the potential to expand and fill in the range of ages and metallicity over which gyrochronology can be empirically tested. In this paper, we demonstrate a Monte Carlo approach to evaluate gyrochronology models using binary stars. As examples, we used five previously published wide binary pairs. We also demonstrate a Monte Carlo approach to assess the precision and accuracy of ages derived from each gyrochronology model. For the traditional Skumanich models, the age uncertainties are σ age /age = 15%–20% for starsmore »with B − V = 0.65 and σ age /age = 5%–10% for stars with B − V = 1.5 and rotation period P ≤ 20 days.« less