Arxiv:2009.06972V2 [Astro-Ph.SR] 1 Oct 2020
Total Page:16
File Type:pdf, Size:1020Kb
MNRAS 000,1{19 (2020) Preprint 2 October 2020 Compiled using MNRAS LATEX style file v3.0 Asteroseismic Inference of Subgiant Evolutionary Parameters with Deep Learning Marc Hon1?, Earl P. Bellinger2;1;3, Saskia Hekker3;2, Dennis Stello1;2;4, and James S. Kuszlewicz3;2 1School of Physics, The University of New South Wales, Sydney NSW 2052, Australia 2Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, DK-8000 Aarhus C, Denmark 3Max-Planck-Institut fur¨ Sonnensystemforschung, Justus-von-Liebig-Weg 3, 37077 G¨ottingen, Germany 4Sydney Institute for Astronomy (SIfA), School of Physics, University of Sydney, NSW 2006, Australia Accepted XXX. Received YYY; in original form ZZZ ABSTRACT With the observations of an unprecedented number of oscillating subgiant stars ex- pected from NASA's TESS mission, the asteroseismic characterization of subgiant stars will be a vital task for stellar population studies and for testing our theories of stellar evolution. To determine the fundamental properties of a large sample of subgiant stars efficiently, we developed a deep learning method that estimates distributions of fundamental parameters like age and mass over a wide range of input physics by learn- ing from a grid of stellar models varied in eight physical parameters. We applied our method to four Kepler subgiant stars and compare our results with previously deter- mined estimates. Our results show good agreement with previous estimates for three of them (KIC 11026764, KIC 10920273, KIC 11395018). With the ability to explore a vast range of stellar parameters, we determine that the remaining star, KIC 10005473, is likely to have an age 1 Gyr younger than its previously determined estimate. Our method also estimates the efficiency of overshooting, undershooting, and microscopic diffusion processes, from which we determined that the parameters governing such pro- cesses are generally poorly-constrained in subgiant models. We further demonstrate our method's utility for ensemble asteroseismology by characterizing a sample of 30 Kepler subgiant stars, where we find a majority of our age, mass, and radius estimates agree within uncertainties from more computationally expensive grid-based modelling techniques. Key words: asteroseismology { stars: oscillations { stars: evolution { methods: data analysis 1 INTRODUCTION 1977). As the interiors of subgiants evolve over relatively arXiv:2009.06972v2 [astro-ph.SR] 1 Oct 2020 short timescales, the mixed mode behaviour of the star's Asteroseismology of solar-like oscillations is a powerful ap- oscillation spectrum also changes rapidly (e.g., Christensen- proach to measure ages of individual field stars. By prob- Dalsgaard et al. 1995). Hence, detailed measurements of sub- ing the stellar interior, asteroseismic measurements can re- giant mixed modes not only provide valuable diagnostics of veal structural changes that are indicators of stellar evo- the stellar interior (e.g., Deheuvels & Michel 2011; Benomar lution. This is especially the case for subgiant stars that et al. 2012; Benomar et al. 2014), but also yield precise stel- have begun to show mixed modes in their oscillation spec- lar age estimates (e.g., Deheuvels et al. 2014; Metcalfe et al. tra. These modes arise from the coupling of acoustic waves 2014; Li et al. 2017). that propagate in the stellar envelope with gravity (g-) waves that propagate near the core (Osaki 1975), and re- Owing to high-quality photometric observations from sult in perturbations to the near-uniform frequency spacing the Kepler space mission (Borucki et al. 2010), precise os- of acoustic (p-) modes (avoided crossings, Aizenman et al. cillation frequencies have been measured for subgiant stars (e.g., Appourchaux et al. 2012). Such measurements have en- abled the fundamental stellar parameters of subgiants to be ? E-mail: [email protected] determined using stellar modelling techniques (e.g., Metcalfe © 2020 The Authors 2 Hon et al. Figure 1. (a) General schematic of the deep neural network in this work. The network takes as input individual mode frequencies, xfreq, along with measurements from global seismic parameters and spectroscopic measurements, xobs, to predict the parameters describing a 10-dimensional Gaussian mixture distribution of stellar model parameters, y. These parameters are the mean (µ), deviations (σ), and mixture coefficient (π) of each Gaussian in the mixture. (b) Schematic of a mixture density network. The network maps input x, which is indicated by the neurons within the box, into conditional density p¹y j xº by predicting the shape parameters π¹xº; µ¹xº; σ¹xº for as many as k Gaussian functions, which are combined to form a mixture model as described in Equations4 and5. et al. 2010; Creevey et al. 2012; Doˇgan et al. 2013; Stokholm solutions. Grid interpolation methods mitigate the need for et al. 2019). Although only a small number of oscillating sub- a very fine grid of models; however they still struggle with giants were observed by Kepler, this number is expected to high computational complexity once additional dimensions be amplified by NASA's Transiting Exoplanet Survey Satel- are included in the search. A common alternative is to re- lite (TESS), where at least a few hundred oscillating sub- strict the search to only a few free parameters and use ap- giants are expected to be observed for a year (Campante proximations for other initial model parameters. These in- et al. 2016; Schofield et al. 2019). There will therefore be clude the adoption of a solar-calibrated value for the mixing further opportunities for studying subgiant stellar structure length parameter (αMLT), or the use of the Galactic enrich- and evolution along the subgiant branch. ment relation to estimate the initial helium abundance (Y0) Stellar models are necessary for inferring stellar ages but using the initial metal abundance (Z0). These assumptions the task of finding a model that best fits the observables may lead to underestimated uncertainties and/or systematic from a star is computationally demanding. Such a task is errors when inferring stellar properties from models. An ad- a non-linear, high-dimensional optimization problem, where ditional prohibiting factor in subgiant model searches is the the complex relations governing stellar structure and evolu- time-consuming calculation of non-radial modes for evolved tion (E) are sensitive to numerous input physical parameters stars, which makes it expensive for search methods that re- that are being optimized (P) such as the star's mass, initial quire either a large grid of models or the on-the-fly calcula- composition, and mixing parameters. Traditional optimiza- tion of stellar tracks (e.g., Metcalfe et al. 2009; Paxton et al. tion methods find a best-matching set of parameters (P∗) 2013). that best fits the observed properties of a star (O) by solv- ing the following: Bellinger et al.(2016, hereafter BA16) showed that E¹Pº − Oº2 these problems can be mitigated for main-sequence stars P∗ = arg min ; (1) by using machine learning to infer the parameters of stel- U2 P lar models from a given set of observables. Machine learning where U is the uncertainty from O. However, as the dimen- techniques, once trained, are able to statistically capture the sionality of P increases, the volume of the parameter space complex relations connecting observations to stellar models involved in the search increases exponentially. In an attempt at a fraction of the computational cost required for model to make stellar model searches tractable, traditional opti- grid searches. In other terms, machine learning algorithms mization methods typically deploy one or more of the fol- can learn to approximate the inverse relation E−1 between lowing strategies: lowering the model grid density, grid inter- model parameters and observed data. As a result, such al- polation (e.g., Rendle et al. 2019), or reducing the number gorithms output maximum likelihood estimates for P∗ by of initial model parameters that are explored in the search. computing E−1¹Oº. These algorithms have been shown by Lowering the grid density significantly reduces the number BA16 to be effective in the systematic age determination of models required to be generated, but comes at the cost of of all main-sequence stars within the high-quality Kepler parameter coverage that may result in finding sub-optimal LEGACY sample with an age precision closely comparable MNRAS 000,1{19 (2020) Subgiant properties with deep learning 3 to those inferred from traditional grid-based optimization the surface, which reduces a star's age at a given mean den- methods (Angelou et al. 2017; Bellinger et al. 2019a). sity (e.g., Miglio & Montalb´an 2005; Gai et al. 2009; Valle In this work, we seek to extend machine learning-based et al. 2015). Meanwhile, the undershooting parameter αunder stellar model inference towards subgiant stars using deep controls the inwards extent of the outer convective boundary learning. A major difference between our work and the BA16 of the stellar envelope and is often constrained to be equiv- study is the type of asteroseismic stellar age proxy used. The alent with αover. For exploratory purposes, we set αunder to observed oscillation frequency ratios r0;2, which are known be a free parameter. to be sensitive towards core hydrogen abundance (Roxburgh Despite much evidence in literature indicating the im- & Vorontsov 2003), are typically used as a stellar age proxy portance of these additional processes in stellar models (e.g., for main-sequence stars (e.g., Christensen-Dalsgaard 1984; Guzik & Cox 1993; Gruyters et al. 2013; Silva Aguirre et al. White et al. 2011; Bellinger & Christensen-Dalsgaard 2019). 2013), there remains significant uncertainty in both theory These ratios, however, are no longer effective age proxies for and observations regarding the nature and efficiency of such core hydrogen-depleted subgiant stars. Instead, observations processes. It is therefore common for modelling tasks to ei- of rapidly evolving mixed modes can be used to precisely ther disregard the parameters governing these additional constrain subgiant stellar ages.