No Arabic abstract
Aims: We investigate from a theoretical perspective if space asteroseismology can be used to distinguish between different thermal structures and shapes of the near-core mixing profiles for different types of coherent oscillation modes in massive stars with convective cores, and if this capacity depends on the evolutionary stage of the models along the main sequence. Methods: We compute 1D stellar structure and evolution models for four different prescriptions of the mixing and temperature gradient in the near-core region. Their effect on the frequencies of dipole prograde gravity modes in both Slowly Pulsating B and $beta$ Cep stars is investigated, as well as for pressure modes in $beta$ Cep stars. Results: A comparison between the mode frequencies of the different models at various stages during the main sequence evolution reveals that they are more sensitive to a change in temperature gradient than to the exact shape of the mixing profile in the near-core region. Depending on the duration of the observed light curve, one can distinguish between either just the temperature gradient, or also between the shapes of the mixing coefficient. The relative frequency differences are in general larger for more evolved models, and are largest for the higher-frequency pressure modes in $beta$ Cep stars. Conclusions:In order to unravel the core boundary mixing and thermal structure of the near-core region, one must have asteroseismic masses and radii with $sim 1%$ relative precision for hundreds of stars.
Convective boundary mixing (CBM) is ubiquitous in stellar evolution. It is a necessary ingredient in the models in order to match observational constraints from clusters, binaries and single stars alike. We compute `effective overshoot measures that reflect the extent of mixing and which can differ significantly from the input overshoot values set in the stellar evolution codes. We use constraints from pressure modes to infer the CBM properties of Kepler and CoRoT main-sequence and subgiant oscillators, as well as in two radial velocity targets (Procyon A and $alpha$ Cen A). Collectively these targets allow us to identify how measurement precision, stellar spectral type, and overshoot implementation impact the asteroseismic solution. With these new measures we find that the `effective overshoot for most stars is in line with physical expectations and calibrations from binaries and clusters. However, two F-stars in the CoRoT field (HD 49933 and HD 181906) still necessitate high overshoot in the models. Due to short mode lifetimes, mode identification can be difficult in these stars. We demonstrate that an incongruence between the radial and non-radial modes drives the asteroseismic solution to extreme structures with highly-efficient CBM as an inevitable outcome. Understanding the cause of seemingly anomalous physics for such stars is vital for inferring accurate stellar parameters from TESS data with comparable timeseries length.
Yearslong time series of high-precision brightness measurements have been assembled for thousands of stars with telescopes operating in space. Such data have allowed astronomers to measure the physics of stellar interiors via nonradial oscillations, opening a new avenue to study the stars in the Universe. Asteroseismology, the interpretation of the characteristics of oscillation modes in terms of the physical properties of the stellar interior, brought entirely new insights in how stars rotate and how they build up their chemistry throughout their evolution. Data-driven space asteroseismology delivered a drastic increase in the reliability of computer models mimicking the evolution of stars born with a variety of masses and metallicities. Such models are critical ingredients for modern physics as a whole, because they are used throughout various contemporary and multidisciplinary research fields in space science, including the search for life outside the solar system, archaeological studies of the Milky Way, and the study of single and binary supernova progenitors, among which are future gravitational wave sources. The specific role and potential of asteroseismology for those modern research fields are illustrated. The review concludes with current limitations of asteroseismology and highlights how they can be overcome with ongoing and future large infrastructures for survey astronomy combined with new theoretical research in the era of high-performance computing. This review presents results obtained through major community efforts over the past decade. These breakthroughs were achieved in a collaborative and inclusive spirit that is characteristic of the asteroseismology community. The reviews aim is to make this research field accessible to graduate students and readers coming from other fields of physics, with incentives to join future applications in this domain of astrophysics.
The size of convective cores remains uncertain, despite its substantial influence on stellar evolution, and thus on stellar ages. The seismic modeling of young subgiants can be used to obtain indirect constraints on the core structure during main sequence, thanks to the high probing potential of mixed modes. We selected the young subgiant KIC10273246, observed by Kepler, based on its mixed-mode properties. We thoroughly modeled this star, with the aim of placing constraints on the size of its main sequence convective core. We first extracted the parameters of the oscillation modes of the star using the full Kepler data set. To overcome the challenges posed by the seismic modeling of subgiants, we proposed a method which is specifically tailored for subgiants with mixed modes and consists in a nested optimization. We then applied this method to perform a detailed seismic modeling of KIC10273246. We obtained models that show good statistical agreements with the observations, both seismic and non-seismic. We showed that including core overshooting in the models significantly improves the quality of the seismic fit, optimal models being found for $alpha_{mathrm{ov}} = 0.15$. Higher amounts of core overshooting strongly worsen the agreement with the observations and are thus firmly ruled out. We also found that having access to two g-dominated mixed modes in young subgiants allows us to place stronger constraints on the gradient of molecular weight in the core and on the central density. This study confirms the high potential of young subgiants with mixed modes to investigate the size of main-sequence convective cores. It paves the way for a more general study including the subgiants observed with Kepler, TESS, and eventually PLATO.
The Gaia M-dwarf gap is a significant under-density of stars observed near $M_G = 10.2$ in a color-magnitude diagram for stars within 200 pc of the Sun. It has been proposed that the gap is the manifestation of structural instabilities within stellar interiors due to non-equilibrium $^{3}$He fusion prior to some stars becoming fully convective. To test this hypothesis, we use Dartmouth stellar evolution models, MARCS model atmospheres, and simple stellar population synthesis to create synthetic $M_G$-($G_{rm BP} - G_{rm RP})$ color-magnitude diagrams. We confirm that the proposed $^{3}$He instability is responsible for the appearance of the M-dwarf gap. Our synthetic gap shows qualitatively similar features to the observed gap including: its vertical extent in $M_G$, its slope in the color-magnitude diagram, and its relative prominence at bluer colors as compared to redder colors. Furthermore, corresponding over-densities of stars above the gap are reproduced by the models. While qualitatively similar, the synthetic gap is approximately 0.2 magnitudes bluer and, accounting for this color offset, 0.16 magnitudes brighter than the observed gap. Our results reveal that the Gaia M dwarf gap is sensitive to conditions within cores of M dwarf stars, making the gap a powerful tool for testing the physics of M dwarf stars and potentially using M dwarfs to understand the local star formation history.
Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with the PROMPI code. We performed a mean field analysis of the simulations within the Reynolds-averaged Navier-Stokes framework. Both the vertical RMS velocity within the convective region and the bulk Richardson number of the boundaries are found to scale with the driving luminosity as expected from theory. The positions of the convective boundaries were estimated through the composition profiles across them, and the strength of convective boundary mixing was determined by analysing the boundaries within the framework of the entrainment law. We find that the entrainment is approximately inversely proportional to the bulk Richardson number. Although the entrainment law does not encompass all the processes occurring at boundaries, our results support the use of the entrainment law to describe convective boundary mixing in 1D models, at least for the advanced phases. The next steps and challenges ahead are also discussed.