Abstract
We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (
MESA
). The new
auto
_
diff
module implements automatic differentiation ...in
MESA
, an enabling capability that alleviates the need for hard-coded analytic expressions or finite-difference approximations. We significantly enhance the treatment of the growth and decay of convection in
MESA
with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron-degenerate ignition events. We strengthen
MESA
’s implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in
MESA
, we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator-split nuclear burning mode. We close by discussing major updates to
MESA
’s software infrastructure that enhance source code development and community engagement.
We present a method for measuring internal stellar structure based on asteroseismology that we call "inversions for agreement." The method accounts for imprecise estimates of stellar mass and radius ...as well as the relatively limited oscillation mode sets that are available for distant stars. By construction, the results of the method are independent of stellar models. We apply this method to measure the isothermal sound speeds in the cores of the solar-type stars 16 Cyg A and B using asteroseismic data obtained from Kepler observations. We compare the asteroseismic structure that we deduce against best-fitting evolutionary models and find that the sound speeds in the cores of these stars exceed those of the models.
ABSTRACT
We present two approaches to determine the dynamical stability of a hierarchical triple-star system. The first is an improvement on the Mardling–Aarseth stability formula from 2001, where we ...introduce a dependence on inner orbital eccentricity and improve the dependence on mutual orbital inclination. The second involves a machine learning approach, where we use a multilayer perceptron (MLP) to classify triple-star systems as ‘stable’ and ‘unstable’. To achieve this, we generate a large training data set of 106 hierarchical triples using the N-body code MSTAR. Both our approaches perform better than previous stability criteria, with the MLP model performing the best. The improved stability formula and the machine learning model have overall classification accuracies of $93{{\ \rm per\ cent}}$ and $95{{\ \rm per\ cent}}$ respectively. Our MLP model, which accurately predicts the stability of any hierarchical triple-star system within the parameter ranges studied with almost no computation required, is publicly available on Github in the form of an easy-to-use python script.
ABSTRACT
The period of pulsation and the structure of the light curve for Cepheid and RR Lyrae variables depend on the fundamental parameters of the star: mass, radius, luminosity, and effective ...temperature. Here, we train artificial neural networks on theoretical pulsation models to predict the fundamental parameters of these stars based on their period and light-curve structure. We find significant improvements to estimates of these parameters made using light-curve structure and period over estimates made using only the period. Given that the models are able to reproduce most observables, we find that the fundamental parameters of these stars can be estimated up to 60 per cent more accurately when light-curve structure is taken into consideration. We quantify which aspects of light-curve structure are most important in determining fundamental parameters, and find, for example, that the second Fourier amplitude component of RR Lyrae light curves is even more important than period in determining the effective temperature of the star. We apply this analysis to observations of hundreds Cepheids in the Large Magellanic Cloud and thousands of RR Lyrae in the Magellanic Clouds and Galactic bulge to produce catalogues of estimated masses, radii, luminosities, and other parameters of these stars. As an example application, we estimate Wesenheit indices and use those to derive distance moduli to the Magellanic Clouds of μLMC,CEP = 18.688 ± 0.093, μLMC,RRL = 18.52 ± 0.14, and μSMC,RRL = 18.88 ± 0.17 mag.
The goal of stellar evolution theory is to predict the structure of stars throughout their lifetimes. Usually, these predictions can be assessed only indirectly, for example by comparing predicted ...and observed effective temperatures and luminosities. Thanks now to asteroseismology, which can reveal the internal structure of stars, it becomes possible to compare the predictions from stellar evolution theory to actual stellar structures. In this work, we present an inverse analysis of the oscillation data from the solar-type star KIC 6225718, which was observed by the Kepler space observatory during its nominal mission. As its mass is about 20% greater than solar, this star is predicted to transport energy by convection in its nuclear-burning core. We find significant differences between the predicted and actual structure of the star in the radiative interior near to the convective core. In particular, the predicted sound speed is higher than observed in the deep interior of the star, and too low at a fractional radius of 0.25 and beyond. The cause of these discrepancies is unknown, and is not remedied by known physics in the form of convective overshooting or elemental diffusion.
Abstract
The theoretical oscillation frequencies of even the best asteroseismic models of solar-like oscillators show significant differences from observed oscillation frequencies. Structure ...inversions seek to use these frequency differences to infer the underlying differences in stellar structure. While used extensively to study the Sun, structure inversion results for other stars have so far been limited. Applying sound speed inversions to more stars allows us to probe stellar theory over a larger range of conditions, as well as look for overall patterns that may hint at deficits in our current understanding. To that end, we present structure inversion results for 12 main-sequence solar-type stars with masses between 1 and 1.15
M
⊙
. Our inversions are able to infer differences in the isothermal sound speed in the innermost 30% by radius of our target stars. In half of our target stars, the structure of our best-fit model fully agrees with the observations. In the remainder, the inversions reveal significant differences between the sound speed profile of the star and that of the model. We find five stars where the sound speed in the core of our stellar models is too low and one star showing the opposite behavior. For the two stars in which our inversions reveal the most significant differences, we examine whether changing the microphysics of our models improves them and find that changes to nuclear reaction rates or core opacities can reduce, but do not fully resolve, the differences.
Astronomy is in an era where all-sky surveys are mapping the Galaxy. The plethora of photometric, spectroscopic, asteroseismic, and astrometric data allows us to characterize the comprising stars in ...detail. Here we quantify to what extent precise stellar observations reveal information about the properties of a star, including properties that are unobserved, or even unobservable. We analyze the diagnostic potential of classical and asteroseismic observations for inferring stellar parameters such as age, mass, and radius from evolutionary tracks of solar-like oscillators on the lower main sequence. We perform rank correlation tests in order to determine the capacity of each observable quantity to probe structural components of stars and infer their evolutionary histories. We also analyze the principal components of classic and asteroseismic observables to highlight the degree of redundancy present in the measured quantities and demonstrate the extent to which information of the model parameters can be extracted. We perform multiple regression using combinations of observable quantities in a grid of evolutionary simulations and appraise the predictive utility of each combination in determining the properties of stars. We identify the combinations that are useful and provide limits to where each type of observable quantity can reveal information about a star. We investigate the accuracy with which targets in the upcoming TESS and PLATO missions can be characterized. We demonstrate that the combination of observations from GAIA and PLATO will allow us to tightly constrain stellar masses, ages, and radii with machine learning for the purposes of Galactic and planetary studies.
ABSTRACT Owing to the remarkable photometric precision of space observatories like Kepler, stellar and planetary systems beyond our own are now being characterized en masse for the first time. These ...characterizations are pivotal for endeavors such as searching for Earth-like planets and solar twins, understanding the mechanisms that govern stellar evolution, and tracing the dynamics of our Galaxy. The volume of data that is becoming available, however, brings with it the need to process this information accurately and rapidly. While existing methods can constrain fundamental stellar parameters such as ages, masses, and radii from these observations, they require substantial computational effort to do so. We develop a method based on machine learning for rapidly estimating fundamental parameters of main-sequence solar-like stars from classical and asteroseismic observations. We first demonstrate this method on a hare-and-hound exercise and then apply it to the Sun, 16 Cyg A and B, and 34 planet-hosting candidates that have been observed by the Kepler spacecraft. We find that our estimates and their associated uncertainties are comparable to the results of other methods, but with the additional benefit of being able to explore many more stellar parameters while using much less computation time. We furthermore use this method to present evidence for an empirical diffusion-mass relation. Our method is open source and freely available for the community to use.6
Solar Evolution Models with a Central Black Hole Bellinger, Earl P.; Caplan, Matt E.; Ryu, Taeho ...
Astrophysical journal/The Astrophysical journal,
12/2023, Letnik:
959, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Abstract
Hawking proposed that the Sun may harbor a primordial black hole (BH) whose accretion supplies some of the solar luminosity. Such an object would have formed within the first 1 s after the ...Big Bang with the mass of a moon or an asteroid. These light BHs are a candidate solution to the dark matter problem, and could grow to become stellar-mass BHs if captured by stars. Here we compute the evolution of stars having such a BH at their center. We find that such objects can be surprisingly long-lived, with the lightest BHs having no influence over stellar evolution, while more massive ones consume the star over time to produce a range of observable consequences. Models of the Sun born about a BH whose mass has since grown to approximately 10
−6
M
⊙
are compatible with current observations. In this scenario, the Sun would first dim to half its current luminosity over a span of 100 Myr as the accretion starts to generate enough energy to quench nuclear reactions. The Sun would then expand into a fully convective star, where it would shine luminously for potentially several gigayears with an enriched surface helium abundance, first as a sub-subgiant star, and later as a red straggler, before becoming a subsolar-mass BH. We also present results for a range of stellar masses and metallicities. The unique internal structures of stars harboring BHs may make it possible for asteroseismology to discover them, should they exist. We conclude with a list of open problems and predictions.
ABSTRACT
With the observations of an unprecedented number of oscillating subgiant stars expected 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 learning 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 determined 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 processes 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.