Search

We study local magnetohydrodynamical instabilities of differential rotation in magnetized, stably stratified regions of stars and planets using a Cartesian Boussinesq model. We consider arbitrary latitudes and general shears (with gravity direction misaligned from this by an angle |$\phi$|⁠), to model radial (⁠|$\phi =0$|⁠), latitudinal (⁠|$\phi =\pm 90^\circ$|⁠), and mixed differential rotations, and study both non-diffusive [including magnetorotational instability (MRI) and Solberg–Høiland instability] and diffusive instabilities [including Goldreich–Schubert–Fricke (GSF) and MRI with diffusion]. These instabilities could drive turbulent transport and mixing in radiative regions, including the solar tachocline and the cores of red giant stars, but their dynamics are incompletely understood. We revisit linear axisymmetric instabilities with and without diffusion and analyse their properties in the presence of magnetic fields, including deriving stability criteria and computing growth rates, wave vectors, and energetics, both analytically and numerically. We present a more comprehensive analysis of axisymmetric local instabilities than prior work, exploring arbitrary differential rotations and diffusive processes. The presence of a magnetic field leads to stability criteria depending upon angular velocity rather than angular momentum gradients. We find MRI operates for much weaker differential rotations than the hydrodynamic GSF instability, and that it typically prefers much larger length-scales, while the GSF instability is impeded by realistic strength magnetic fields. We anticipate MRI to be more important for turbulent transport in the solar tachocline than the GSF instability when |$\phi \gt 0$| in the Northern (and vice versa in the Southern) hemisphere, though the latter could operate just below the convection zone when MRI is absent for |$\phi \lt 0$|⁠. [ABSTRACT FROM AUTHOR]

The inner and outer cores of neutron stars are believed to contain type-I and -II proton superconductors, respectively. The type-I superconductor exists in an intermediate state, comprising macroscopic flux-free and flux-containing regions, while the type-II superconductor is flux-free, except for microscopic, quantized flux tubes. Here, we show that the inner and outer cores are coupled magnetically, when the macroscopic flux tubes subdivide dendritically into quantized flux tubes, a phenomenon called flux branching. An important implication is that up to |${\sim} 10^{12} (r_1/10^6 \, {\rm cm}) \, {\rm erg}$| of energy are required to separate a quantized flux tube from its progenitor macroscopic flux tube, where |$r_1$| is the length of the macroscopic flux tube. Approximating the normal-superconducting boundary as sharp, we calculate the magnetic coupling energy between a quantized and macroscopic flux tube due to flux branching as a function of, |$f_1$|⁠ , the radius of the type-I inner core divided by the radius of the type-II outer core. Strong coupling delays magnetic field decay in the type-II superconductor. For an idealized inner core containing only a type-I proton superconductor and poloidal flux, and in the absence of ambipolar diffusion and diamagnetic screening, the low magnetic moments (⁠|${\lesssim} 10^{27} \, {\rm G \, cm^3}$|⁠) of recycled pulsars imply |$f_1 \lesssim 10^{-1.5}$|⁠. [ABSTRACT FROM AUTHOR]

Context. Gravito-inertial asteroseismology saw its birth from the 4-year-long light curves of rotating main-sequence stars assembled by the Kepler space telescope. High-precision measurements of internal rotation and mixing are available for about 600 stars of intermediate mass so far that are used to challenge the state-of-the-art stellar structure and evolution models. Aims. Our aim is to prepare for future large ensemble modelling of gravity-mode pulsators by relying on a new sample of such stars recently discovered from the third Data Release of the Gaia space mission and confirmed by space photometry from the TESS mission. This sample of potential asteroseismic targets is about 23 times larger than the Kepler sample. Methods. We use the effective temperature and luminosity inferred from Gaia to deduce evolutionary masses, convective core masses, radii, and ages for ∼14 000 gravity-mode pulsators classified as such from their nominal TESS light curves. We do so by constructing two dedicated grids of evolutionary models for rotating stars with input physics from the asteroseismic calibrations of Keplerγ Dor pulsators. These two grids consider the distribution of initial rotation velocities at the zero-age main sequence deduced from gravito-inertial asteroseismology, for two extreme values found for the metallicity of γ Dor stars deduced from spectroscopy ([M/H]=0.0 and −0.5). Results. We find the new gravity-mode pulsators to cover an extended observational instability region covering masses from about 1.3 M⊙ to about 9 M⊙. We provide their mass-luminosity and mass-radius relations, as well as convective core masses. Our results suggest that oscillations excited by the opacity mechanism occur uninterruptedly for the mass range above about 2 M⊙, where stars have a radiative envelope aside from thin convection zones in their excitation layers. Conclusions. Our evolutionary parameters for the sample of Gaia-discovered gravity-mode pulsators with confirmed modes by TESS offer a fruitful starting point for future TESS ensemble asteroseismology once a sufficient number of modes is identified in terms of the geometrical wave numbers and overtone for each of the pulsators. [ABSTRACT FROM AUTHOR]

The investigation of global, resonant oscillation modes in red giant stars offers valuable insights into their internal structures. In this study, we investigate in detail the information we can recover on the structural properties of core-helium burning (CHeB) stars by examining how the coupling between gravity- and pressure-mode cavities depends on several stellar properties, including mass, chemical composition, and evolutionary state. Using the structure of models computed with the stellar evolution code MESA, we calculate the coupling coefficient implementing analytical expressions, which are appropriate for the strong coupling regime and the structure of the evanescent region in CHeB stars. Our analysis reveals a notable anti-correlation between the coupling coefficient and both the mass and metallicity of stars in the regime M ≲ 1.8 M⊙, in agreement with Kepler data. We attribute this correlation primarily to variations in the density contrast between the stellar envelope and core. The strongest coupling is expected thus for red-horizontal branch stars, partially stripped stars, and stars in the higher-mass range exhibiting solar-like oscillations (M ≳ 1.8 M⊙). While our investigation emphasises some limitations of current analytical expressions, it also presents promising avenues. The frequency dependence of the coupling coefficient emerges as a potential tool for reconstructing the detailed stratification of the evanescent region. [ABSTRACT FROM AUTHOR]

Precise masses of red giant stars enable a robust inference of their ages, but there are cases where these age estimates are very precise but also very inaccurate. Examples are core-helium-burning (CHeB) stars that have lost more mass than predicted by standard single-star evolutionary models. Members of star clusters in the Kepler database represent a unique opportunity to identify such stars because they combine exquisite asteroseismic constraints with independent age information (members of a star cluster share a similar age and chemical composition). We focus on the single metal-rich (Z ≈ Z⊙) Li-rich low-mass CHeB star KIC4937011, which is a member of the open cluster NGC 6819 (turn-off mass of ≈1.6 M⊙, i.e. an age of ≈2.4 Gyr). This star has a lower mass by ≈1 M⊙ than expected for its age and metallicity, which might be explained by binary interactions or mass loss along the red giant branch (RGB). To infer formation scenarios for this object, we performed a Bayesian analysis by combining the binary stellar evolutionary framework BINARY_C V2.2.3 with the dynamic nested-sampling approach contained in the DYNESTY V2.1.1 package. We find that this star probably is the result of a common-envelope evolution (CEE) phase during the RGB stage of the primary star in which the low-mass (< 0.71 M⊙) main-sequence companion does not survive. The mass of the primary star at the zero-age main sequence is in the range [1.46, 1.71] M⊙, with a log-orbital period in the range [0.06, 2.4] log10(days). During the CEE phase, ≈1 M⊙ of material is ejected from the system, and the final star reaches the CHeB stage after helium flashes as if it were a single star with a mass of ≈0.7 M⊙, which is what we observe today. Although the proposed scenario is consistent with photometric and spectroscopic observations, a quantitative comparison with detailed stellar evolution calculations is needed to quantify the systematic skewness of the radius, luminosity, and effective temperature distributions towards higher values than observations. [ABSTRACT FROM AUTHOR]

Context. A central challenge in the field of stellar astrophysics lies in accurately determining the mass of stars, particularly when dealing with isolated ones. However, for pulsating white dwarf stars, the task becomes more tractable due to the availability of multiple approaches such as spectroscopy, asteroseismology, astrometry, and photometry, each providing valuable insights into the mass properties of white dwarf stars. Aims. Numerous asteroseismological studies of white dwarfs have been published, focusing on determining stellar mass using pulsational spectra and comparing it with spectroscopic mass, which uses surface temperature and gravity. The objective of this work is to compare these mass values in detail and, in turn, to compare them with the mass values derived using astrometric parallaxes or distances and photometry data from Gaia, employing astrometric and photometric methods. Methods. Our analysis involves a selection of pulsating white dwarfs with different surface chemical abundances that define the main classes of variable white dwarfs. We calculated their spectroscopic masses, compiled seismological masses, and determined astrometric masses. We also derived photometric masses, when possible. Subsequently, we compared all the sets of stellar masses obtained through these different methods. To ensure consistency and robustness in our comparisons, we used identical white dwarf models and evolutionary tracks across all four methods. Results. The analysis suggests a general consensus among the four methods regarding the masses of pulsating white dwarfs with hydrogen-rich atmospheres, known as DAV or ZZ Ceti stars, especially for objects with masses below approximately 0.75 M⊙, although notable disparities emerge for certain massive stars. For pulsating white dwarf stars with helium-rich atmospheres, called DBV or V777 Her stars, we find that astrometric masses generally exceed seismological, spectroscopic, and photometric masses. Finally, while there is agreement among the sets of stellar masses for pulsating white dwarfs with carbon-, oxygen-, and helium-rich atmospheres (designated as GW Vir stars), outliers exist, where mass determinations by various methods show significant discrepancies. Conclusions. Although a general agreement exists among different methodologies for estimating the mass of pulsating white dwarfs, significant discrepancies are prevalent in many instances. This shows the need to redo the determination of spectroscopic parameters and the parallax and/or improve asteroseismological models for many stars. [ABSTRACT FROM AUTHOR]

One way to understand planet formation is through studying the correlations between planet occurrence rates and stellar mass. However, measuring stellar mass in the red giant regime is very difficult. In particular, the spectroscopic masses of certain evolved stars, often referred to as 'retired A-stars', have been questioned in the literature. Efforts to resolve this mass controversy using spectroscopy, interferometry, and asteroseismology have so far been inconclusive. A recent ensemble study found a mass-dependent mass offset, but the result was based on only 16 stars. With NASA's Transiting Exoplanet SurveySatellite (TESS), we expand the investigation of the mass discrepancy to a total of 92 low-luminosity stars, synonymous with the retired A-stars. We measure their characteristic oscillation frequency, |$\mathrm{\nu }_{\mathrm{max}}$|⁠ , and the large frequency separation, |$\mathrm{\Delta \nu }$|⁠ , from their TESS photometric time series. Using these measurements and asteroseismic scaling relations, we derive asteroseismic masses and compare them with spectroscopic masses from five surveys, to comprehensively study the alleged mass-dependent mass offset. We find a mass offset between spectroscopy and seismology that increases with stellar mass. However, we note that adopting the seismic mass scale does not have a significant effect on the planet occurrence-mass-metallicity correlation for the so-called retired A-stars. We also report seismic measurements and masses for 157 higher luminosity giants (mostly helium-core-burning) from the spectroscopic surveys. [ABSTRACT FROM AUTHOR]

Aims. By adopting the recently empirically derived dependence of α-elements on [α/Fe] instead of the conventionally applied uniform one, we tested the agreement between stellar model predictions and observations for red giant branch (RGB) stars in the APO-K2 catalogue. We particularly focused on the biases in effective temperature scales and on the robustness of age estimations. Methods. We computed a grid of stellar models relying on the empirical scaling of α-elements, investigating the offset in effective temperature ΔT between these models and observations, using univariate analyses for both metallicity [Fe/H] and [α/Fe]. To account for potential confounding factors, we then employed a multivariate generalised additive model to study the dependence of ΔT on [Fe/H], [α/Fe], log g, and stellar mass. Results. The initial analysis revealed a negligible trend of ΔT with [Fe/H], in contrast with previous works in the literature, which adopt a uniform relation between the various α-elements and [α/Fe]. A slight ΔT difference of 25 K was detected between stars with high and low α-enhancement. Our multivariate analysis reveals a dependence of ΔT on both [Fe/H] and [α/Fe], and highlights a significant dependence on stellar mass. This suggests a discrepancy in how effective temperature scales with stellar mass in the models compared to observations. Despite differences in assumed chemical composition, our analysis, through a fortunate cancellation effect, yields ages that are largely consistent with recent studies of the same sample. Notably, our analysis identifies a 6% fraction of stars younger than 4 Ga within the high-α population. However, our analysis of the [C/N] ratio supports the possible origin of these stars as a result of mergers or mass transfer events. [ABSTRACT FROM AUTHOR]

Aims. We compared stellar radii derived from asteroseismic scaling relations with those estimated using two independent surface brightness-colour relations (SBCRs) combined with Gaia DR3 parallaxes. Methods. We cross-matched asteroseismic and astrometric data for over 6400 red giant branch (RGB) and red clump (RC) stars from the APO-K2 catalogue with the TESS Input Catalogue v8.2 to obtain precise V band magnitudes and E(B − V) colour excesses. We then adopted two different SBCRs from the literature to derive stellar radius estimates, denoted as Ra and Rb, respectively. We analysed the ratio of these SBCR-derived radii to the asteroseismic radius estimates, R, provided in the APO-K2 catalogue. Results. Both SBCRs exhibited good agreement with asteroseismic radius estimates. On average, Ra was overestimated by 1.2% with respect to R, while Rb was underestimated by 2.5%. For stars larger than 20 R⊙, SBCR radii are systematically lower than asteroseismic ones. The dispersion in the radius ratio was similar for the two methods (around 10%). The agreement with asteroseismic radii shows a strong dependence on the parallax. The dispersion is halved for stars with a parallax greater than 2.5 mas. In this subsample, Rb showed perfect agreement with R, while Ra remained slightly overestimated, by 3%. A trend with [Fe/H] was found at a level of 4% to 6% per dex. Additionally, a clear trend with asteroseismic mass is found. For stars less massive than about 0.95 M⊙, SBCR radii were significantly higher than asteroseismic ones, by about 6%. This overestimation correlated with the presence of extended helium cores in these stars' structures relative to their envelopes. Furthermore, radius ratios showed a dichotomous behaviour at higher masses, mainly due to the presence of several RC stars with SBCR radii significantly lower with respect to asteroseismology. This behaviour originates from a different response of asteroseismic scaling relations and SBCR to [α/Fe] abundance ratios for massive stars, both in RGB and RC phases, which is reported here for the first time. [ABSTRACT FROM AUTHOR]

The recent measurement of magnetic field strength inside the radiative interior of red giant stars has opened the way toward full 3D characterization of the geometry of stable large-scale magnetic fields. However, current measurements, which are limited to dipolar (ℓ = 1) mixed modes, do not properly constrain the topology of magnetic fields due to degeneracies on the observed magnetic field signature on such ℓ = 1 mode frequencies. Efforts focused toward unambiguous detections of magnetic field configurations are now key to better understand angular momentum transport in stars. We investigated the detectability of complex magnetic field topologies (such as the ones observed at the surface of stars with a radiative envelope with spectropolarimetry) inside the radiative interior of red giants. We focused on a field composed of a combination of a dipole and a quadrupole (quadrudipole) and on an offset field. We explored the potential of probing such magnetic field topologies from a combined measurement of magnetic signatures on ℓ = 1 and quadrupolar (ℓ = 2) mixed mode oscillation frequencies. We first derived the asymptotic theoretical formalism for computing the asymmetric signature in the frequency pattern for ℓ = 2 modes due to a quadrudipole magnetic field. To access asymmetry parameters for more complex magnetic field topologies, we numerically performed a grid search over the parameter space to map the degeneracy of the signatures of given topologies. We demonstrate the crucial role played by ℓ = 2 mixed modes in accessing internal magnetic fields with a quadrupolar component. The degeneracy of the quadrudipole compared to pure dipolar fields is lifted when considering magnetic asymmetries in both ℓ = 1 and ℓ = 2 mode frequencies. In addition to the analytical derivation for the quadrudipole, we present the prospect for complex magnetic field inversions using magnetic sensitivity kernels from standard perturbation analysis for forward modeling. Using this method, we explored the detectability of offset magnetic fields from ℓ = 1 and ℓ = 2 frequencies and demonstrate that offset fields may be mistaken for weak and centered magnetic fields, resulting in underestimating the magnetic field strength in stellar cores. We emphasize the need to characterize ℓ = 2 mixed-mode frequencies, (along with the currently characterized ℓ = 1 mixed modes), to unveil the higher-order components of the geometry of buried magnetic fields and to better constrain angular momentum transport inside stars. [ABSTRACT FROM AUTHOR]

Context. In the partially ionized material of stellar interiors, the strongest forces acting on electrons and ions are the Coulomb interactions between charges. The dynamics of the plasma as a whole depend on the magnitudes of the average electrostatic interactions and the average kinetic energies of the particles that constitute the stellar material. An important question is how these interactions of real gases are related to the observable stellar properties. Specifically, the relationships between rotation, magnetic activity, and the thermodynamic properties of stellar interiors are still not well understood. These connections are crucial for understanding and interpreting the abundant observational data provided by space-based missions, such as Kepler/K2 and TESS, and the future data from the PLATO mission. Aims. In this study, we investigate the electrostatic effects within the interiors of low-mass main sequence (MS) stars. Specifically, we introduce a global quantity, a global plasma parameter, which allows us to compare the importance of electrostatic interactions across a range of low-mass theoretical models (0.7 − 1.4 M⊙) with varying ages and metallicities. We then correlate the electrostatic properties of the theoretical models with the observable rotational trends on the MS. Methods. We use the open-source 1D stellar evolution code MESA to compute a grid of main-sequence stellar models. Our models span the log g − Teff space of a set of 66 Kepler main-sequence stars. Results. We identify a correlation between the prominence of electrostatic effects in stellar interiors and stellar rotation rates. The variations in the magnitude of electrostatic interactions with age and metallicity further suggest that understanding the underlying physics of the collective effects of plasma can clarify key observational trends related to the rotation of low-mass stars on the MS. These results may also advance our understanding of the physics behind the observed weakened magnetic braking in stars. [ABSTRACT FROM AUTHOR]

Context. The Kepler space mission provided high-quality light curves for more than 16 000 red giants. The global stellar oscillations extracted from these light curves carry information about the interior of the stars. Several hundred red giants were found to have low amplitudes in their dipole modes (i.e. they are suppressed dipole-mode stars). A number of hypotheses (involving e.g. a magnetic field, binarity, or resonant mode coupling) have been proposed to explain the suppression of the modes, yet none has been confirmed. Aims. We aim to gain insight into the mechanism at play in suppressed dipole-mode stars by investigating the mode properties (linewidths, heights, and amplitudes) of the radial oscillation modes of red giants with suppressed dipole modes. Methods. We selected from the literature suppressed dipole-mode stars and compared the radial-mode properties of these stars to the radial-mode properties of stars in two control samples of stars with typical (i.e. non-suppressed) dipole modes. Results. We find that the radial-mode properties of the suppressed dipole-mode stars are consistent with the ones in our control samples, and hence not affected by the suppression mechanism. Conclusions. From this we conclude that (1) the balance between the excitation and damping in radial modes is unaffected by the suppression, and by extrapolation the excitation of the non-radial modes is not affected either; and (2) the damping of the radial modes induced by the suppression mechanism is significantly less than the damping from turbulent convective motion, suggesting that the additional damping originates from the more central non-convective regions of the star, to which the radial modes are least sensitive. [ABSTRACT FROM AUTHOR]

Context. In order to improve the asteroseismic modelling efforts for red-giant stars, the numerical computation of theoretical individual oscillation modes for evolved red-giant models has to be made feasible. Aims. We aim to derive a method for circumventing the computational cost of computing oscillation spectra for models of red-giant stars with an average large frequency separation Δν < 15 μHz, thereby allowing asteroseismic investigations of giants utilising individual frequencies. Methods. The proposed truncated scanning method serves as a novel method detailing how the observable individual l = 0, 1, 2 frequencies of red giants may be computed on realistic timescales through so-called model truncation. By removing the innermost region of the stellar models, the g-mode influence on the oscillation spectra may be avoided. This then allows us to estimate the observable regions from the resulting pure p-mode oscillations. The appropriate observable frequency regions are subsequently scanned for the complete and non-truncated stellar model. The observable regions are determined by considering the limitations on observability from the internal mode coupling and damping, yielding consistent frequency spectra obtained at a much reduced computational cost. Results. The truncated scanning method proves that the individual frequencies of red-giant models can be obtained for a wide range of applications and research. The computational efficiency is improved by a factor of 10 or better. This means that the inclusion of l = 1, 2 individual frequencies is now a possibility in future asteroseismic modelling efforts of red-giant stars. Further potential avenues for improvements of this method are outlined for future pursuits. [ABSTRACT FROM AUTHOR]

Ultramassive white dwarfs with masses |$M\gtrsim 1.1\, {\rm M}_{\odot }$| probe extreme physics near the Chandrasekhar limit. Despite the rapid increase in observations, it is still unclear how many harbour carbon–oxygen (CO) versus oxygen–neon (ONe) cores. The origin of these white dwarfs and their strong magnetic fields – single stellar evolution or a stellar merger – is another open question. The steep mass–radius relation of the relativistic ultramassive white dwarfs shortens their crystallization time |$t_{\rm cryst}$|⁠ , such that the recently proposed crystallization dynamo mechanism may present an alternative to mergers in explaining the early appearance of magnetism in the observed population. However, the magnetic diffusion time from the convective dynamo to the white dwarf's surface delays the magnetic field's breakout time |$t_{\rm break}\gt t_{\rm cryst}$|⁠. We compute |$t_{\rm break}(M)$| for CO and ONe ultramassive white dwarfs and compare it to the local 40 pc volume-limited sample. We find that the breakout time from CO cores is too long to account for the observations. ONe crystallization dynamos remain a viable option, but their surrounding non-convective envelopes comprise only a few per cent of the total mass, such that |$t_{\rm break}$| is highly sensitive to the details of stellar evolution. [ABSTRACT FROM AUTHOR]

Strong metallicity-dependent winds dominate the evolution of core He-burning, classical Wolf–Rayet (cWR) stars, which eject both H and He-fusion products such as |$^{14}$| N, |$^{12}$| C, |$^{16}$| O, |$^{19}$| F, |$^{22}$| Ne, and |$^{23}$| Na during their evolution. The chemical enrichment from cWRs can be significant. cWR stars are also key sources for neutron production relevant for the weak s-process. We calculate stellar models of cWRs at solar metallicity for a range of initial Helium star masses (12–50  |$\rm M_{\odot }$|⁠), adopting recent hydrodynamical wind rates. Stellar wind yields are provided for the entire post-main sequence evolution until core O-exhaustion. While literature has previously considered cWRs as a viable source of the radioisotope |$^{26}$| Al, we confirm that negligible |$^{26}$| Al is ejected by cWRs since it has decayed to |$^{26}$| Mg or proton-captured to |$^{27}$| Al. However, in Paper I, we showed that very massive stars eject substantial quantities of |$^{26}$| Al, among other elements including N, Ne, and Na, already from the zero-age-main-sequence. Here, we examine the production of |$^{19}$| F and find that even with lower mass-loss rates than previous studies, our cWR models still eject substantial amounts of |$^{19}$| F. We provide central neutron densities (N |$_{n}$|⁠) of a 30  |$\rm M_{\odot }$| cWR compared with a 32  |$\rm M_{\odot }$| post-VMS WR and confirm that during core He-burning, cWRs produce a significant number of neutrons for the weak s-process via the |$^{22}$| Ne(⁠|$\alpha$|⁠ ,n) |$^{25}$| Mg reaction. Finally, we compare our cWR models with observed [Ne/He], [C/He], and [O/He] ratios of Galactic WC and WO stars. [ABSTRACT FROM AUTHOR]

One-dimensional (1D) stellar evolution models are widely used across various astrophysical fields, however they are still dominated by important uncertainties that deeply affect their predictive power. Among those, the merging of independent convective regions is a poorly understood phenomenon predicted by some 1D models but whose occurrence and impact in real stars remain very uncertain. Being an intrinsically multi-D phenomenon, it is challenging to predict the exact behaviour of shell mergers with 1D models. In this work, we conduct a detailed investigation of a multiple shell merging event in a 20 M |$_\odot$| star using 3D hydrodynamic simulations. Making use of the active tracers for composition and the nuclear network included in the 3D model, we study the merging not only from a dynamical standpoint but also considering its nucleosynthesis and energy generation. Our simulations confirm the occurrence of the merging also in 3D, but reveal significant differences from the 1D case. Specifically, we identify entrainment and the erosion of stable regions as the main mechanisms that drive the merging, we predict much faster convective velocities compared to the mixing-length theory velocities, and observe multiple burning phases within the same merged shell, with important effects for the chemical composition of the star, which presents a strongly asymmetric (dipolar) distribution. We expect that these differences will have important effects on the final structure of massive stars and thus their final collapse dynamics and possible supernova explosion, subsequently affecting the resulting nucleosynthesis and remnant. [ABSTRACT FROM AUTHOR]

The thermal evolution of isolated neutron stars is a key element in unravelling their internal structure and composition and establishing evolutionary connections among different observational subclasses. Previous studies have predominantly focused on one-dimensional or axisymmetric two-dimensional models. In this study, we present the thermal evolution component of the novel three-dimensional magnetothermal code MATINS (MAgneto-Thermal evolution of Isolated Neutron Star). MATINS employs a finite volume scheme and integrates a realistic background structure, along with state-of-the-art microphysical calculations for the conductivities, neutrino emissivities, heat capacity, and superfluid gap models. This paper outlines the methodology employed to solve the thermal evolution equations in MATINS , along with the microphysical implementation that is essential for the thermal component. We test the accuracy of the code and present simulations with non-evolving magnetic fields of different configurations (all with electrical currents confined to the crust and a magnetic field that does not thread the core), to produce temperature maps of the neutron star surface. Additionally, for a specific magnetic field configuration, we show one fully coupled evolution of magnetic field and temperature. Subsequently, we use a ray-tracing code to link the neutron star surface temperature maps obtained by MATINS with the phase-resolved spectra and pulsed profiles that would be detected by distant observers. This study, together with our previous article focused on the magnetic formalism, presents in detail the most advanced evolutionary code for isolated neutron stars, with the aim of comparison with their timing properties, thermal luminosities and the associated X-ray light curves. [ABSTRACT FROM AUTHOR]

As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate two-dimensional convection in the core at three different stages of evolution of a |$5\,\mathrm{ M}_{\odot }$| star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably stratified envelope, along with the overshooting length of convective plumes above the convective boundary. We find that the stratification in helium in evolved stars hinders radial motions and effectively shields the radiative envelope against plume penetration. This prevents convective overshooting from being an efficient mixing process in the radiative envelope. In addition, internal gravity waves are less excited in evolved models compared to the zero-age-main-sequence model, and are also more damped in the stratified region above the core. As a result, the wave power is several orders of magnitude lower in mid- and terminal-main-sequence models compared to zero-age-main-sequence stars. [ABSTRACT FROM AUTHOR]

We present time-dependent nova outburst models with optically thick winds for 1.2 and 1.35 |$\, M_{\odot }$| white dwarfs (WDs) with a mass-accretion rate of |$5 \times 10^{-9}\, M_{\odot }$|  yr−1 and for a 1.3 |$\, M_{\odot }$| WD with |$2 \times 10^{-9}\, M_{\odot }$|  yr−1. The X-ray flash occurs 11 d before the optical peak of the 1.2 |$\, M_{\odot }$| WD and 2.5 d before the peak of the 1.3 |$\, M_{\odot }$| WD. The wind mass-loss rate of the 1.2 |$\, M_{\odot }$| WD (1.3 |$\, M_{\odot }$| WD) reaches a peak of |$6.4 \times 10^{-5}\, M_{\odot }$| yr−1 (⁠|$7.4 \times 10^{-5}\, M_{\odot }$| yr−1) at the epoch of the maximum photospheric expansion with the lowest photospheric temperature of log  T ph (K) = 4.33 (4.35). The nuclear energy generated during the outburst is lost in the form of radiation (61% for the 1.2 |$\, M_{\odot }$| WD; 47% for the 1.3 |$\, M_{\odot }$| WD), gravitational energy of ejecta (39%; 52%), and kinetic energy of the wind (0.28%; 0.29%). We found an empirical relation for fast novae between the time to optical maximum from the outburst t peak and the expansion timescale τexp. With this relation, we are able to predict the time to optical maximum t peak from the ignition model (at t = 0) without following a time-consuming nova wind evolution. [ABSTRACT FROM AUTHOR]

We conducted a study on the thermal properties of stellar matter with the nuclear energy density functional BCPM. This functional is based on microscopic Brueckner–Hartree–Fock calculations and has demonstrated success in describing cold neutron stars. To enhance its applicability in astrophysics, we extended the BCPM equation of state to finite temperature for β-stable neutrino-free matter, taking into consideration the hot inner crust. Such an equation of state holds significant importance for hot compact objects, particularly those resulting from a binary neutron star merger event. Our exploration has shown that with increasing temperature, there is a fast decrease in the crust-core transition density, suggesting that for hot stars it is not realistic to assume a fixed value of this density. The microscopic calculations also reveal that the presence of nuclear clusters persists up to T = 7.21 MeV, identified as the limiting temperature of the crust. Above this threshold, the manifestation of clusters is not anticipated. Below this temperature, clusters within the inner crust are surrounded by uniform matter with varying densities, allowing for the distinction between the upper and lower transition density branches. Moreover, we computed mass–radius relations of neutron stars, assuming an isothermal profile for β-stable neutron star matter at various temperature values. Our findings highlight the significant influence of the hot inner crust on the mass–radius relationship, leading to the formation of larger and more inflated neutron stars. Consequently, under our prescription, the final outcome is a unified equation of state at finite temperature. [ABSTRACT FROM AUTHOR]

Aims. We aim to investigate the theoretical possibility of accurately determining the helium-to-metal enrichment ratio ΔY/ΔZ from precise observations of double-lined eclipsing binary systems. Methods. Using Monte Carlo simulations, we drew synthetic binary systems with masses between 0.85 and 1.00 M⊙ from a grid of stellar models. Both stars were sampled from a grid with ΔY/ΔZ = 2.0, with a primary star at 80% of its main-sequence evolution. Subsequently, a broader grid with ΔY/ΔZ from 1.0 to 3.0 was used in the fitting process. To account for observational uncertainties, two scenarios were explored: S1, with realistic uncertainties of 100 K in temperature and 0.1 dex in [Fe/H]; and S2, with halved uncertainties. We repeated the simulation at two baseline metallicities: [Fe/H] = 0.0 and −0.3. Results. The posterior distributions of ΔY/ΔZ revealed significant biases. The distributions were severely biased towards the edge of the allowable range in the S1 error scenario. The situation only marginally improved when considering the S2 scenario. The effect is due to the impact of changing ΔY/ΔZ in the stellar effective temperature and its interplay with [Fe/H] observational error, and it is therefore not restricted to the specific fitting method. Despite the presence of these systematic discrepancies, the age of the systems were recovered unbiased with 10% precision. Conclusions. Our findings indicate that the observational uncertainty in effective temperature and metallicity significantly hinders the accurate determination of the ΔY/ΔZ parameter from main-sequence binary systems. [ABSTRACT FROM AUTHOR]

We present new simulations investigating the impact of mass transfer on the asteroseismic signals of slowly pulsating B stars. We used MESA to simulate the evolution of a binary star system and GYRE to compute the asteroseismic properties of the accretor star. We show that, compared to a single star of the same final mass, a star that has undergone accretion (of non-enriched material) has a significantly different internal structure, which is evident in both the hydrogen abundance profile and the Brunt-Väisälä frequency profile. These differences result in significant changes in the observed period spacing patterns, implying that one may use this as a diagnostic to test whether a star's core has been rejuvenated as a result of accretion. We show that it is essential to consider the full multimodal posterior distributions when fitting stellar properties of mass-gainers to avoid drawing misleading conclusions. Even with these considerations, stellar ages will be significantly underestimated when assuming single star evolution for a mass-gainer. We find that future detectors with improved uncertainties would rule out single star models with the correct mass and central hydrogen fraction. Our proof of principle analysis demonstrates the need to further investigate the impact of binary interactions on stellar asteroseismic signals for a wide range of parameters, such as the initial mass, the amount of mass transferred, and the age of the accretor star at the onset of mass transfer. [ABSTRACT FROM AUTHOR]

Context. Numerical computations of stellar oscillations for models representative of B-type stars predict fewer modes to be excited than observations reveal from modern space-based photometric data. One shortcoming of state-of-the-art evolution models of B-type stars that may cause a lack of excited modes is the absence of microscopic diffusion in most such models. Aims. We investigate whether the inclusion of microscopic diffusion in stellar models of B-type stars, notably radiative levitation experienced by isotopes, leads to extra mode driving by the opacity mechanism compared to the case of models that do not include microscopic diffusion. Methods. We consider the case of slowly to moderately rotating stars and use non-rotating equilibrium models, while we account for (uniform) rotation in the computations of the pulsation frequencies. We calculate 1D stellar structure and evolution models with and without microscopic diffusion and examine the effect of radiative levitation on mode excitation, for both low-radial order pressure and gravity modes and for high-radial order gravity modes. As is common practice in asteroseismology, rotation is included in the pulsation computations according to the mode's frequency regime. For modes having frequencies below twice the rotation frequency, that is, modes in the sub-inertial regime, we adopt the traditional approximation of rotation. For modes in the super-inertial regime with frequency above twice the rotation frequency, rotation is treated perturbatively up to first order in the rotation. We consider macroscopic envelope mixing induced by internal gravity waves to compute the modes and study its effect on the surface abundances. Results. We find systematically more modes to be excited for the stellar models including microscopic diffusion compared to those without it, in agreement with observational findings of pulsating B-type dwarfs. Furthermore, the models with microscopic diffusion predict that excited modes occur earlier on in the evolution compared to modes without it. In order to maintain realistic surface abundances during the main sequence, we include macroscopic envelope mixing by internal gravity waves. Along with microscopic diffusion, such macroscopic envelope mixing ensures both more excited modes and surface abundances consistent with spectroscopic studies of B-type stars. Conclusions. While radiative levitation has so far largely been neglected in stellar evolution computations of B-type stars for computational convenience, it impacts mode excitation predictions for stellar models of such stars. We conclude that the process of radiative levitation is able to reduce the discrepancy between predicted and observed excited pulsation modes in B-type stars. [ABSTRACT FROM AUTHOR]

Aims. The most reliable sources for determining absolute stellar parameters are the double-lined eclipsing binary systems. Some of these systems also show apsidal motion, characterized by the variable log k2. This point grants us the possibility of investigating the stellar interior, specifically the degree of stellar mass concentration. The first studies carried out about four decades ago on this topic showed appreciable discrepancies not only with respect to the comparison between the observed absolute dimensions and their theoretical counterparts, but mainly concerning the degree of mass concentration through the analysis of their apsidal motions. Fortunately, this scenario has been gradually improving with the advances in the quality of observational techniques and advances in the input physics of the evolutionary stellar models (e.g. opacities, thermonuclear reactions, equations of state, numerical techniques). These new developments in the input physics has improved the comparison between observations and the values predicted by theory, including the apsidal motion rates. This progress has lead us to investigate second-order effects such as rotation and dynamic tides. In this paper we deal with the effects of rotation on the degree of mass concentration. Methods. The stellar models were calculated using the MESA package. The mass range studied here was 0.40–25.0 M⊙ for a generic chemical composition characterized by X = 0.70 and Z = 0.02. The present models were computed without taking into account core overshooting in order to highlight the effects of rotation. Each model was followed from the pre-main sequence until the central hydrogen content is of the order of or less than 0.03, covering the range of masses and evolutionary status of the double-lined eclipsing binaries showing apsidal motion. Regarding the calculation of the internal structure coefficients, we integrated the Radau equation using the fifth-order Runge-Kutta method and a tolerance level of 10−7. As an auxiliary tool, homology transformations have been used to explain, qualitatively, how the equation of state, thermonuclear reactions, and convection impact the degree of mass concentration of the models. Results. In our past studies on the effects of rotation on log k2 an average correction was used for all models together. For this a small range of masses (2.0, 7.0, 15.0 M⊙) has been used. In the present paper the corrections due to the effects of rotation in log k2 are presented for each mass individually and for three evolutionary phases. This point is particularly important, given that these corrections show a clear dependence on the mass and on the evolutionary status. Such corrections can be easily introduced into the theoretical calculation of apsidal motion rates through a linear equation. A typical correction due to the rotation effects is of the order of −0.03 in log k2. [ABSTRACT FROM AUTHOR]

Sharp structural variations induce specific signatures on stellar pulsations that can be studied to infer localised information on the stratification of the star. This information is key to improve our understanding of the physical processes that lead to the structural variations and how to model them. Here we revisit and extend the analysis of the signature of different types of buoyancy glitches in gravity-mode and mixed-mode pulsators presented in earlier works, including glitches with step-like, Gaussian-like, and Dirac-δ-like shapes. In particular, we provide analytical expressions for the perturbations to the periods and show that these can be reliably used in place of the expressions provided for the period spacings, with the advantage that the use of the new expressions does not require modes with consecutive radial orders to be observed. Based on a comparison with two limit cases and on simulated data, we further tested the accuracy of the expression for the Gaussian-like glitch signature whose derivation in an earlier work involved a significant approximation. We find that the least reliable glitch parameter inferred from fitting that expression is the amplitude, which can be up to a factor of two larger than the true amplitude, reaching this limit when the glitch is small. We further discuss the impact on the glitch signature of considering a glitch in the inner and outer half of the g-mode cavity, emphasising the break of symmetry that takes place in the case of mixed-mode pulsators. [ABSTRACT FROM AUTHOR]

Context. Modeling the convection process is a long-standing problem in stellar physics. To date, all ad hoc models have relied on a free parameter, α, (among others) that has no real physical justification and is therefore poorly constrained. However, a link exists between this free parameter and the entropy of the stellar adiabat. There are existing prescriptions, derived from 3D stellar atmospheric models, that treat entropy as a function of stellar atmospheric parameters (effective temperature, surface gravity, and chemical composition). This can offer sufficient constraints on α through the development of entropy-calibrated models. However, several questions have arisen as these models are increasingly used with respect to which prescription should be used and whether it ought to be used in its original form, along with the impacts of uncertainties on entropy-calibrated models. Aims. We aim to study the three existing prescriptions in detail and determine which of them demonstrate the most optimal performance and how it should be applied. Methods. We implemented the entropy-calibration method into the stellar evolution code (Cesam2k20) and performed comparisons with the Sun and the α Cen system. In addition, we used data from the CIFIST grid of 3D atmosphere models to evaluate the accuracy of the prescriptions. Results. Of the three entropy prescriptions currently available, we determined the one that has the best functional form for reproducing the entropies of the 3D models. However, the coefficients involved in this formulation must not be taken from the original paper because they were calibrated against a flawed set of entropies. We also demonstrate that the entropy obtained from this prescription should be corrected for the evolving chemical composition and for an entropy offset different between various EoS tables. This must be done following a precise procedure to ensure that the classical parameters obtained from the models are not strongly biased. Finally, we provide a data table with entropy of the adiabat of the CIFIST grid, along with the fits for these entropies. Conclusions. Thanks to a precise examination of entropy-calibrated modeling, we are able to offer our recommendations with respect to which adiabatic entropy prescription to use, how to correct it, and how to implement the method into a stellar evolution code. [ABSTRACT FROM AUTHOR]

The treatment of convection remains a major weakness in the modelling of stellar evolution with one-dimensional (1D) codes. The ever-increasing computing power makes now possible to simulate in three-dimensional (3D) part of a star for a fraction of its life, allowing us to study the full complexity of convective zones with hydrodynamics codes. Here, we performed state-of-the-art hydrodynamics simulations of turbulence in a neon-burning convective zone, during the late stage of the life of a massive star. We produced a set of simulations varying the resolution of the computing domain (from 1283 to 10243 cells) and the efficiency of the nuclear reactions (by boosting the energy generation rate from nominal to a factor of 1000). We analysed our results by the mean of Fourier transform of the velocity field, and mean-field decomposition of the various transport equations. Our results are in line with previous studies, showing that the behaviour of the bulk of the convective zone is already well captured at a relatively low resolution (2563), while the details of the convective boundaries require higher resolutions. The different boosting factors used show how various quantities (velocity, buoyancy, abundances, and abundance variances) depend on the energy generation rate. We found that for low boosting factors, convective zones are well mixed, validating the approach usually used in 1D stellar evolution codes. However, when nuclear burning and turbulent transport occur on the same time-scale, a more sophisticated treatment would be needed. This is typically the case when shell mergers occur. [ABSTRACT FROM AUTHOR]

Aims. The most reliable sources for determining absolute stellar parameters are the double-lined eclipsing binary systems. Some of these systems also show apsidal motion, characterized by the variable log k2. This point grants us the possibility of investigating the stellar interior, specifically the degree of stellar mass concentration. The first studies carried out about four decades ago on this topic showed appreciable discrepancies not only with respect to the comparison between the observed absolute dimensions and their theoretical counterparts, but mainly concerning the degree of mass concentration through the analysis of their apsidal motions. Fortunately, this scenario has been gradually improving with the advances in the quality of observational techniques and advances in the input physics of the evolutionary stellar models (e.g. opacities, thermonuclear reactions, equations of state, numerical techniques). These new developments in the input physics has improved the comparison between observations and the values predicted by theory, including the apsidal motion rates. This progress has lead us to investigate second-order effects such as rotation and dynamic tides. In this paper we deal with the effects of rotation on the degree of mass concentration. Methods. The stellar models were calculated using the MESA package. The mass range studied here was 0.40–25.0 M⊙ for a generic chemical composition characterized by X = 0.70 and Z = 0.02. The present models were computed without taking into account core overshooting in order to highlight the effects of rotation. Each model was followed from the pre-main sequence until the central hydrogen content is of the order of or less than 0.03, covering the range of masses and evolutionary status of the double-lined eclipsing binaries showing apsidal motion. Regarding the calculation of the internal structure coefficients, we integrated the Radau equation using the fifth-order Runge-Kutta method and a tolerance level of 10−7. As an auxiliary tool, homology transformations have been used to explain, qualitatively, how the equation of state, thermonuclear reactions, and convection impact the degree of mass concentration of the models. Results. In our past studies on the effects of rotation on log k2 an average correction was used for all models together. For this a small range of masses (2.0, 7.0, 15.0 M⊙) has been used. In the present paper the corrections due to the effects of rotation in log k2 are presented for each mass individually and for three evolutionary phases. This point is particularly important, given that these corrections show a clear dependence on the mass and on the evolutionary status. Such corrections can be easily introduced into the theoretical calculation of apsidal motion rates through a linear equation. A typical correction due to the rotation effects is of the order of −0.03 in log k2. [ABSTRACT FROM AUTHOR]

Mature neutron stars are expected to exhibit gravity g modes due to stratification caused by a varying matter composition. These modes are affected by nuclear reactions, leading to complex (damped) mode frequencies and the suppression of high-order g modes. This is in contrast with the common non-dissipative analysis which leads to an infinite g-mode spectrum. Focusing on the transition between the fast- and slow-reaction regimes, we examine the precise impact of nuclear reactions on the g-mode spectrum. The general framework for the analysis is presented along with sample numerical results for a matter model based on the BSk21 equation of state with a suitable parametrization for the reaction rates. [ABSTRACT FROM AUTHOR]

Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light-curve observations. We present the oscillation properties of high-resolution hydrodynamic simulations of a |$25\,\,{\rm{M}_\odot }$| star performed with the PPMstar code. The model star has a convective core mass of |$\approx 12\,\,{\rm{M}_\odot }$| and approximately half of the envelope simulated. From this simulation, we extract light curves from several directions, average them over each hemisphere, and process them as if they were real photometric observations. We show how core convection excites waves with a similar frequency as the convective time-scale in addition to significant power across a forest of low and high angular degree l modes. We find that the coherence of these modes is relatively low as a result of their stochastic excitation by core convection, with lifetimes of the order of 10s of days. Thanks to the still significant power at higher l and this relatively low coherence, we find that integrating over a hemisphere produces a power spectrum that still contains measurable power up to the Brunt–Väisälä frequency. These power spectra extracted from the stable envelope are qualitatively similar to observations, with the same order of magnitude yet lower characteristic frequency. This work further shows the potential of long-duration, high-resolution hydrodynamic simulations for connecting asteroseismic observations to the structure and dynamics of core convection and the convective boundary. [ABSTRACT FROM AUTHOR]

Context. The theory of rotational and chemical evolution is incomplete, thereby limiting the accuracy of model-dependent stellar mass and age determinations. The γ Doradus (γ Dor) pulsators are excellent points of calibration for the current state-of-the-art stellar evolution models, as their gravity modes probe the physical conditions in the deep stellar interior. Yet, individual asteroseismic modelling of these stars is not always possible because of insufficient observed oscillation modes. Aims. This paper presents a novel method to derive distributions of the stellar mass, age, core-boundary mixing efficiency, and initial rotation rates for γ Dor stars. Methods. We computed a grid of rotating stellar evolution models covering the entire γ Dor instability strip. We then used the observed distributions of the luminosity, effective temperature, buoyancy travel time, and near-core rotation frequency of a sample of 539 stars to assign a statistical weight to each of our models. This weight is a measure of how likely the combination of a specific model is. We then computed weighted histograms to derive the most likely distributions of the fundamental stellar properties. Results. We find that the rotation frequency at zero-age main sequence follows a normal distribution, peaking at around 25% of the critical Keplerian rotation frequency. The probability-density function for extent of the core-boundary mixing zone, given by a factor of fCBM times the local pressure scale height (assuming an exponentially decaying parameterisation), decreases linearly with increasing fCBM. Conclusions. Converting the distribution of fractions of critical rotation at the zero-age main sequence to units of d−1, we find most F-type stars start the main sequence with a rotation frequency between 0.5 d−1 and 2 d−1. Regarding the core-boundary mixing efficiency, we find that it is generally weak in this mass regime. [ABSTRACT FROM AUTHOR]

Context. Modeling the surface brightness distribution of stars is of prime importance to interpret the large amount of available interferometric, spectropolarimetric, or photometric observations. Beyond stellar physics, this is also a prerequisite to characterize exoplanets or our Galaxy. Nevertheless, this remains quite challenging for cool stars as it requires one to model the magnetohydrodynamic turbulence that develops in their convective envelope. Aims. In Paper I (Raynaud, R., Rieutord, M., Petitdemange, L., Gastine, T., & Putigny, B. 2018, A&A, 609, A124), the effect of the Coriolis acceleration on the surface heat flux has been studied by means of hydrodynamic simulations. In this paper, we aim to investigate the additional effect of dynamo magnetic fields that can be generated in the thick convective envelopes of cool stars. We focus on an envelope thickness that is representative of either a ∼0.35 M⊙ M dwarf, a young red giant star or a pre-main sequence star. Methods. We performed a parametric study using numerical magnetohydrodynamic simulations of anelastic convection in thick rotating spherical shells. The stratification in density ranges from a few tens to a few hundreds. The setup assumes a constant entropy jump between the inner and outer layers to force convection, with stress-free boundary conditions for the velocity field. The magnetic Prandtl number was systematically varied in order to vary the magnetic field intensity. For each model, we computed the azimuthally and temporally averaged surface distribution of the heat flux, and examined the leading-order effect of the magnetic field on the obtained latitudinal luminosity profile. Results. We identify three different regimes. Close to the onset of convection, while the first unstable modes tend to convey heat more efficiently near the equator, magnetic fields are shown to generally enhance the mean heat flux close to the polar regions (and the tangent cylinder). By progressively increasing the Rayleigh number, the development of a prograde equatorial jet was previously shown to make the equator darker when no magnetic field is taken into account. For moderate Rayleigh numbers, magnetic fields can instead inverse the mean pole-equator brightness contrast (which means going from a darker to a brighter equator when a dynamo sets in) and finally induce a similar regime to that found close to the onset of convection. For more turbulent models with larger Rayleigh numbers, magnetic fields alternatively tend to smooth out the brightness contrast. This general behavior is shown to be related to the quenching of the surface differential rotation by magnetic fields and remains valid regardless of the magnetic morphology. Conclusions. Mean global trends regarding the impact of rotation and magnetic fields on the surface brightness distribution of cool stars are theoretically depicted and need to be tested by future observations. This work opens the door to more detailed theoretical studies including the effect of nonaxisymmetric and time-variable surface features associated with magnetic activity. [ABSTRACT FROM AUTHOR]

Aims. We investigated the theoretical biases affecting the asteroseismic grid-based estimates of stellar mass, radius, and age in the presence of a mismatch between the heavy element mixture of observed stars and stellar models. Methods. We performed a controlled simulation adopting a stellar effective temperature, [Fe/H], an average large frequency spacing, and a frequency of maximum oscillation power as observational constraints. Synthetic stars were sampled from grids of stellar models computed with different [α/Fe] values from 0.0 to 0.4. The mass, radius, and age of these objects were then estimated by adopting a grid of models with a fixed [α/Fe] value of 0.0. The experiment was repeated assuming different sets of observational uncertainties. In the reference scenario, we adopted an uncertainty of 1.5% in seismic parameters, 50 K in effective temperature, and 0.05 dex in [Fe/H]. A higher uncertainty in the atmospheric constraints was also adopted in order to explore the impact on the precision of the observations of the estimated stellar parameters. Results. Our Monte Carlo experiment showed that estimated parameters are biased up to 3% in mass, 1.5% in radius, and 4% in age when the reference uncertainty scenario was adopted. These values correspond to 45%, 48%, and 16% of the estimated uncertainty in the stellar parameters. These non-negligible biases in mass and radius disappear when adopting larger observational uncertainties because of the possibility of the fitting algorithm exploring a wider range of possible solutions. However, in this scenario, the age is significantly biased by −8%. Finally, we verified that the stellar mass, radius, and age can be estimated with a high accuracy by adopting a grid with the incorrect value of [α/Fe] if the metallicity [Fe/H] of the target is adjusted to match the Z in the fitting grid. In this scenario, the maximum bias in the age was reduced to 1.5%. [ABSTRACT FROM AUTHOR]

Gravitational-wave (GW) observations are revealing the population of compact objects from a new angle. Yet their stellar progenitors remain uncertain because few observational clues on their progenitors exist. Theoretical models typically assume that the progenitor evolution can be approximated with single-star models. We explore how binary evolution affects the pre-supernova (SN) structure of stars, and the resulting distribution of compact object remnants. We focus on the differences in the core properties of single stars and of donor stars that transfer their outer layers in binary systems and become binary-stripped. We show that the final structures of binary-stripped stars that lose their outer layers before the end of core helium burning are systematically different compared to single stars. As a result, we find that binary-stripped stars tend to explode more easily than single stars and preferentially produce neutron stars and fewer black holes, with consequences for GW progenitors. [ABSTRACT FROM AUTHOR]

In recent years, it has been discovered that massive stars commonly exhibit a non-coherent form of variability in their light curves referred to as stochastic low frequency (SLF) variability. Various physical mechanisms can produce SLF variability in such stars, including stochastic gravity waves excited at the interface of convective and radiative regions, dynamic turbulence generated in the near-surface layers, and clumpy winds. Gravity waves in particular are a promising candidate for explaining SLF variability as they can be ubiquitously generated in main sequence stars owing to the presence of a convective core, and because they provide the large-scale predominantly tangential velocity field required to explain macroturbulence in spectral line fitting. Here, I provide an overview of the methods and results of studying SLF variability in massive stars from time series photometry and spectroscopy. [ABSTRACT FROM AUTHOR]

Our knowledge of massive star evolution is limited by uncertainties linked with multi-dimensional processes taking place in stellar interiors. Important examples are convective boundary mixing (CBM) and entrainment, which are implemented in 1D stellar evolution models assuming simplified prescriptions. 3D hydrodynamics models can improve these prescriptions by studying realistic multi-D processes for a short timerange (minutes or hours). In these proceedings, we present results coming from a new set of high-resolution hydrodynamics simulations of a neon-burning shell in a massive star, and discuss how the entrainment law can be calibrated from 3D models and then used to improve 1D stellar evolution prescriptions. [ABSTRACT FROM AUTHOR]

Although red clump stars function as reliable standard candles, their surface characteristics (i.e. Teff, log g , and [Fe/H]) overlap with those of red giant branch stars, which are not standard candles. Recent results have revealed that spectral features containing carbon (e.g. CN molecular bands) carry information correlating with the 'gold-standard' asteroseismic classifiers that distinguish red clump from red giant branch stars. However, the underlying astrophysical processes driving the correlation between these spectroscopic and asteroseismic quantities in red giants remain inadequately explored. This study aims to enhance our understanding of this 'spectro-seismic' effect by refining the list of key spectral features predicting red giant evolutionary state. In addition, we conduct further investigation into those key spectral features to probe the astrophysical processes driving this connection. We employ the data-driven The Cannon algorithm to analyse high-resolution (R ∼ 80 000) Veloce Rosso spectra from the Anglo-Australian Telescope for 301 red giant stars (where asteroseismic classifications from the TESS mission are known for 123 of the stars). The results highlight molecular spectroscopic features, particularly those containing carbon (e.g. CN), as the primary indicators of the evolutionary states of red giant stars. Furthermore, by investigating CN isotopic pairs (that is, 12C14N and 13C14N), we find suggestions of statistically significant differences in the reduced equivalent widths of such lines, suggesting that physical processes that change the surface abundances and isotopic ratios in red giant stars, such as deep mixing, are the driving forces of the 'spectro-seismic' connection of red giants. [ABSTRACT FROM AUTHOR]

Convection in stars and planets must be maintained against viscous and Ohmic dissipation. Here, we present the first systematic investigation of viscous dissipation in simulations of rotating, density-stratified plane layers of convection. Our simulations consider an anelastic ideal gas, and employ the open-source code Dedalus. We demonstrate that when the convection is sufficiently vigorous, the integrated dissipative heating tends towards a value that is independent of viscosity or thermal diffusivity, but depends on the imposed luminosity and the stratification. We show that knowledge of the dissipation provides a bound on the magnitude of the kinetic energy flux in the convection zone. In our non-rotating cases with simple flow fields, much of the dissipation occurs near the highest possible temperatures, and the kinetic energy flux approaches this bound. In the rotating cases, although the total integrated dissipation is similar, it is much more uniformly distributed (and locally balanced by work against the stratification), with a consequently smaller kinetic energy flux. The heat transport in our rotating simulations is in good agreement with results previously obtained for 3D Boussinesq convection, and approaches the predictions of diffusion-free theory. [ABSTRACT FROM AUTHOR]

We study the surface magnetic field fluctuations due to radial oscillations as a viable cause for the microstructures of the radio pulsar pulse patterns. The electrical conductivity of matter in the outer layer of the crust of a neutron star (NS) plays a crucial role in the resulting surface magnetic field if we assume that the magnetic field is confined to this layer. This outer layer has a rapidly varying matter density that changes the microphysics of the material, affecting the electrical conductivity at every stage of the density change. In this study, the varying electrical conductivity in this rapidly varying density regime of the outer layer of the NS crust – from ∼1011 to about 104 g cm−3 – has been used to calculate the surface magnetic field using the induction equation. A finite effect of the strong gravitational field at the NS surface has also been taken into account. The equations have been solved in MATLAB using the method of lines. Any minor radial fluctuation due to stellar oscillation, in particular the radial oscillations, leads to a fluctuation of the electrical conductivity in the outer layer of the crust. This leads to fluctuations in the surface magnetic field with a frequency equal to the frequency of the stellar oscillation. We find that not only is the variation of the surface magnetic field substantial, but also it does not remain constant throughout the lifetime of the NS. [ABSTRACT FROM AUTHOR]

A convective dynamo operating during the crystallization of white dwarfs is one of the promising channels to produce their observed strong magnetic fields. Although the magnitude of the fields generated by crystallization dynamos is uncertain, their timing may serve as an orthogonal test of this channel's contribution. The carbon–oxygen cores of |$M\approx 0.5-1.0\, {\rm M}_{\odot }$| white dwarfs begin to crystallize at an age t cryst ∝  M −5/3, but the magnetic field is initially trapped in the convection zone – deep inside the CO core. Only once a mass of m cryst has crystallized, the convection zone approaches the white dwarf's helium layer, such that the magnetic diffusion time through the envelope shortens sufficiently for the field to break out to the surface, where it can be observed. This breakout time is longer than t cryst by a few Gyr, scaling as t break ∝  t cryst f −1/2, where f ≡ 1 − m cryst/ M depends on the white dwarf's initial C/O profile before crystallization. The first appearance of strong magnetic fields B ≳ 1 MG in volume-limited samples approximately coincides with our numerically computed t break(M) – potentially signalling crystallization dynamos as a dominant magnetization channel. However, some observed magnetic white dwarfs are slightly younger, challenging this scenario. The dependence of the breakout process on the white dwarf's C/O profile implies that magnetism may probe the CO phase diagram, as well as uncertainties during the core helium burning phase in the white dwarf's progenitor, such as the 12C(α, γ)16O nuclear reaction. [ABSTRACT FROM AUTHOR]

The inner structure of core helium burning (CHeB) stars remains uncertain due to the yet unknown nature of mixing at the boundary of their cores. Large convective cores beyond a bare Schwarzschild model are favoured both from theoretical arguments and from asteroseismological constraints. However, the exact nature of this extra mixing, and in particular the possible presence of semiconvective layers, is still debated. In this work, we approach this problem through a new avenue by performing the first full-sphere 3D hydrodynamics simulations of the interiors of CHeB stars. We use the ppmstar explicit gas dynamics code to simulate the inner 0.45 |$\, {\rm M}_{\odot }$| of a 3 M⊙ CHeB star. Simulations are performed using different Cartesian grid resolutions (7683, 11523, and 17283) and heating rates. We use two different initial states, one based on mesa s's predictive mixing scheme (which significantly extends the core beyond the Schwarzschild boundary) and one based on the convective premixing approach (which exhibits a semiconvective interface). The general behaviour of the flow in the convective core and in the stable envelope (where internal gravity waves are observed) is consistent with our recent simulations of core convection in massive main-sequence stars, and so are the various luminosity scaling relations. The semiconvective layers are dominated by strong internal gravity waves that do not produce measurable species mixing, but overshooting motions from the convective core gradually homogenize the semiconvective interface. This process can possibly completely erase the semiconvective layers, which would imply that CHeB stars do not harbour a semiconvection zone. [ABSTRACT FROM AUTHOR]

Strong magnetic fields are expected to significantly modify the pulsation frequencies of waves propagating in the cores of red giants or in the radiative envelopes of intermediate- and high-mass main-sequence stars. We calculate the g-mode frequencies of stars with magnetic dipole fields which are aligned with their rotational axes, treating both the Lorentz and Coriolis forces non-perturbatively. We provide a compact asymptotic formula for the g-mode period spacing and universally find that strong magnetism decreases this period spacing substantially more than is predicted by perturbation theory. These results are validated with explicit numerical mode calculations for realistic stellar models. The approach we present is highly versatile: once the eigenvalues λ of a certain differential operator are pre-computed as a function of the magnetogravity and rotational frequencies (in units of the mode frequency), the non-perturbative impact of the Coriolis and Lorentz forces is understood under a broad domain of validity and is readily incorporated into asteroseismic modelling. [ABSTRACT FROM AUTHOR]

This is the fourth paper of our new release of the BaSTI (a Bag of Stellar Tracks and Isochrones) stellar model and isochrone library. Following the updated solar-scaled, α-enhanced, and white dwarf model libraries, we present here α-depleted ([α/Fe] = −0.2) evolutionary tracks and isochrones, suitable to study the α-depleted stars discovered in Local Group dwarf galaxies and in the Milky Way. These calculations include all improvements and updates of the solar-scaled and α-enhanced models, and span a mass range between 0.1 and 15 M⊙, and 21 metallicities between [Fe/H] = −3.20 and +0.45 with a helium-to-metal enrichment ratio Δ Y /Δ Z  = 1.31, homogeneous with the solar-scaled and α-enhanced models. The isochrones – available in several photometric filters – cover an age range between ∼20 Myr and 14.5 Gyr, including the pre-main-sequence phase. We have compared our isochrones with independent calculations of α-depleted stellar models, available for the same α-element depletion adopted in the present investigation. We have also discussed the effect of an α-depleted heavy element distribution on the bolometric corrections in different wavelength regimes. Our α-depleted evolutionary tracks and isochrones are publicly available at our BaSTI website. [ABSTRACT FROM AUTHOR]

The theory of nested figures of equilibrium, expanded in Papers I and II, is investigated in the limit where the number of layers of the rotating body is infinite, enabling to reach full heterogeneity. In the asymptotic process, the discrete set of equations becomes a differential equation for the rotation rate. In the special case of rigid rotation (from centre to surface), we are led to an integro-differential equation (IDE) linking the ellipticity of isopycnic surfaces to the equatorial mass-density profile. In contrast with most studies, these equations are not restricted to small flattenings, but are valid for fast rotators as well. We use numerical solutions obtained from the self-consistent-field method to validate this approach. At small ellipticities (slow rotation), we fully recover Clairaut's equation. Comparisons with Chandrasekhar's perturbative approach and with Roberts' work based on virial equations are successful. We derive a criterion to characterize the transition from slow to fast rotators. The treatment of heterogeneous structures containing mass-density jumps is proposed through a modified IDE. [ABSTRACT FROM AUTHOR]

We present the results of our analysis of 122 light-curves from 50 Wolf-Rayet (WR) stars using a red+white noise analysis, where we compare the fitted red noise features with stellar parameters to assess the presence of correlations with stellar parameters. A significant correlation between the amplitude of variability α0 and v∞ was found for the whole sample, along with several other correlations satisfying the Spearman-Rank p <0.001 criterion for both He-burning and WNh stars. Our results are compatible with several plausible processes that can have an influence on the level of variability in the winds of these stars, including a subsurface convection zone and core-generated internal gravity waves. [ABSTRACT FROM AUTHOR]

The current angular momentum (AM) transport models fail to reproduce asterosieimic observations. One of the best candidates to explain this discrepancy is the magnetic field in radiative zones with its various possible topologies, for instance axisymmetric toroidal magnetic field. If such azimuthal field is strong enough, the Tayler's instability could occur which induces a magnetic torque that allows a very efficient transport of AM and could trigger dynamo action in radiative layers. If such important field does not emerge at the surface, spectropolarimetry is blind. In this case, the only way to detect and characterise the field is by using magneto-asteroseismology. It consists in searching for the characteristic signatures of magnetic field in the observed frequency spectra of stellar oscillations. [ABSTRACT FROM AUTHOR]

The galactic binary star LB-1 contains a recently stripped B-type star. Comparing its properties to detailed binary star models shows that tidal braking and magnetic torques lead to low surface rotational velocities in the stripped donors after Roche-lobe overffiow. Models without magnetic torques cannot reproduce the observed low surface rotation. [ABSTRACT FROM AUTHOR]

We present a first 3D magnetohydrodynamic (MHD) simulation of oxygen, neon, and carbon shell burning in a rapidly rotating |$16\hbox{-}\mathrm{M}_\odot$| core-collapse supernova progenitor. We also run a purely hydrodynamic simulation for comparison. After |$\mathord \approx 180\mathrm{s}$| (⁠|$\mathord \approx$| 15 and 7 convective turnovers, respectively), the magnetic fields in the oxygen and neon shells achieve saturation at 1011 and 5 × 1010 G. The strong Maxwell stresses become comparable to the radial Reynolds stresses and eventually suppress convection. The suppression of mixing by convection and shear instabilities results in the depletion of fuel at the base of the burning regions, so that the burning shell eventually move outward to cooler regions, thus reducing the energy generation rate. The strong magnetic fields efficiently transport angular momentum outwards, quickly spinning down the rapidly rotating convective oxygen and neon shells and forcing them into rigid rotation. The hydrodynamic model shows complicated redistribution of angular momentum and develops regions of retrograde rotation at the base of the convective shells. We discuss implications of our results for stellar evolution and for the subsequent core-collapse supernova. The rapid redistribution of angular momentum in the MHD model casts some doubt on the possibility of retaining significant core angular momentum for explosions driven by millisecond magnetars. However, findings from multidimensional models remain tentative until stellar evolution calculations can provide more consistent rotation profiles and estimates of magnetic field strengths to initialize multidimensional simulations without substantial numerical transients. We also stress the need for longer simulations, resolution studies, and an investigation of non-ideal effects. [ABSTRACT FROM AUTHOR]

The observation of an unprecedented number of solar-like oscillating subgiant (SG) stars by the Kepler and TESS missions is crucial for the asteroseismic characterization of these stars, stellar population studies, and the study of stellar evolution theories. Owing to these missions, the fundamental parameters of the solar-like oscillating stars are precisely calculated from the evolution codes using the observed oscillation frequencies. Herein, we considered four solar-like oscillating SG stars. We obtained the fundamental parameters of the SG stars by constructing interior models using asteroseismic and non-asteroseismic observed parameters. The interior models of the four SG stars are constructed using the Modules for Experiments in Stellar Astrophysics code to effectively determine the fundamental properties. Using this method, the four solar-like oscillating SG stars are found to have masses and radii within the ranges of 1.16–1.75 M⊙ and 2.26–3.17 R⊙, respectively. The estimation accuracy of the typical asteroseismic radius, mass, and age is increased by fitting the observed and model reference frequencies. The typical uncertainties of the mass and radius are 3–4 |${{\ \rm per\ cent}}$| and 1–2 |${{\ \rm per\ cent}}$|⁠ , respectively. Furthermore, the observed l  = 1 frequencies, which showed a mixed mode for the first time, were also fitted to the models. Information regarding the gravity and density of the helium core was obtained by examining the mixed modes. Moreover, new asteroseismic methods for determining the age of SG stars are developed for the first time in this study. [ABSTRACT FROM AUTHOR]

G 29 − 38 (TIC 422526868) is one of the brightest (V  = 13.1) and closest (d  = 17.51 pc) pulsating white dwarfs with a hydrogen-rich atmosphere (DAV/ZZ Ceti class). It was observed by the TESS spacecraft in sectors 42 and 56. The atmosphere of G 29 − 38 is polluted by heavy elements that are expected to sink out of visible layers on short time-scales. The photometric TESS data set spans ∼51 d in total, and from this, we identified 56 significant pulsation frequencies, that include rotational frequency multiplets. In addition, we identified 30 combination frequencies in each sector. The oscillation frequencies that we found are associated with g -mode pulsations, with periods spanning from ∼ 260 to ∼ 1400 s. We identified rotational frequency triplets with a mean separation δνℓ = 1 of 4.67 μHz and a quintuplet with a mean separation δνℓ = 2 of 6.67 μHz, from which we estimated a rotation period of about 1.35 ± 0.1 d. We determined a constant period spacing of 41.20 s for ℓ = 1 modes and 22.58 s for ℓ = 2 modes. We performed period-to-period fit analyses and found an asteroseismological model with M ⋆/M⊙ = 0.632 ± 0.03, |$T_{\rm eff}=11\, 635\pm 178$| K, and log  g  = 8.048 ± 0.005 (with a hydrogen envelope mass of M H ∼ 5.6 × 10−5 M ⋆), in good agreement with the values derived from spectroscopy. We obtained an asteroseismic distance of 17.54 pc, which is in excellent agreement with that provided by Gaia (17.51 pc). [ABSTRACT FROM AUTHOR]