Stellar Evolution Models Research Articles

The Global Asteroseismology Project Proof of Concept: Asteroseismology of Massive Stars with Continuous Ground-based Observations

Massive (≳8M ⊙) stars are the progenitors of many astrophysical systems, yet key aspects of their structure and evolution are poorly understood. Asteroseismology has the potential to solve these open puzzles; however, sampling both the short period pulsations and long period beat patterns of massive stars poses many observational challenges. Ground-based single-site observations require years or decades to discern the main oscillation modes. Multisite campaigns were able to shorten this time span, but have not been able to scale up to population studies on sample of objects. Space-based observations can achieve both continuous sampling and observe large numbers of objects; however, most lack the multiband data that is often necessary for mode identification and removing model degeneracies. Here, we develop and test a new ground-based observational strategy for discerning and identifying the main oscillation modes of a massive star in a few months, in a way that can be scaled to large samples. We do so using the Las Cumbres Observatory—a unique facility consisting of robotic, homogeneous telescopes operating as a global network, overcoming most of the challenges of previous multisite efforts, but presenting new challenges which we tailor our strategy to address. This work serves as the proof of concept for the Global Asteroseismology Project, which aims to move massive star asteroseismology from single-objects to bulk studies, unleashing its full potential in constraining stellar structure and evolution models. This work also demonstrates the ability of the Las Cumbres Observatory to perform multisite continuous observations for various science goals.

Modelling stellar variability in archival HARPS data: I - Rotation and activity properties with multidimensional Gaussian processes

ABSTRACT Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-metre per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m s–1 in Sun-like stars. To overcome this, a multidimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S + LEAF 2, enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to (a) measure rotation periods for 49 slow rotators in our sample, (b) explore the impact of these results on the spin-down of middle-aged late F, G, and K stars, and (c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations programme, and help fine-tune current stellar structure and evolution models.

Recent results and future perspectives with solid targets at LUNA

The stellar evolution and chemical make-up of the Universe are determined by nuclear reactions occurring in a wide variety of stellar sites. Precise determinations of the cross sections of these reactions are crucial for the calculation of reaction rates and for the development of stellar evolution models. The Laboratory for Underground Nuclear Astrophysics (LUNA) collaboration has been at the forefront of the direct measurement of nuclear reactions at the low energies of astrophysical interest for the last 35 years. The many significant results achieved at LUNA have been made possible due to the low background conditions uniquely available thanks to its location deep underground at the Laboratori Nazionali del Gran Sasso. Another key aspect of these successes is due to the experience of the LUNA collaboration in the production and characterization of a variety of solid targets used in reaction measurements. In this review, the main production techniques of solid targets are described, as well as the common methods adopted for target degradation monitoring. We also present the results of recent measurements using these targets and the future plans of the LUNA collaboration for measurements using solid targets at the LUNA400 kV and the new Ion Beam Facility (IBF) 3.5 MV are also presented.

Angular momentum transport near convective-core boundaries of Gamma Doradus stars

Recent asteroseismic studies have revealed that the convective core of γ Doradus stars rotates faster than their radiative interior. We study the development of differential rotation near the convective core to test angular momentum transport processes that are typically adopted in stellar evolution models. Models that only include the advection of angular momentum by meridional circulation and shear instabilities cannot reproduce current rotational constraints, irrespective of the initial conditions. The latest formulation of internal magnetic fields based on the Tayler instability is indeed able to reproduce the internal rotation rate of post-main sequence stars; however, it appears too efficient during the main sequence and has thus been disfavoured. A less efficient version of the same transport process can simultaneously reproduce the rotation rate of the convective core, the rotation rate in radiative regions as probed by gravity-modes, and the surface rotational velocities of γ Doradus stars. Our work suggests that there are additional physical processes apart from internal magnetic fields at work in the stellar interiors of post-main sequence stars.

HD 214220 (BD+56 2813): An eclipsing binary with its primary component at the end of the main sequence

Context. Analysis of the light curves of an eclipsing binary allows one to derive the absolute dimensions of the system. This in turn yields information on the radii of the components, which allows the stars to be accurately placed on the Hertzsprung-Russell diagram and their evolutionary phase to be interpreted via comparisons to tracks of stellar evolution models. Aims. I aim to derive the stellar and system parameters of HD 214220. Methods. I measured the epochs of three primary and three secondary minima of the eclipsing binary HD 214220 from 2019 to 2022 from photometric fluxes obtained by the TESS satellite. I modeled the light curve and the velocity amplitudes, which were obtained by the Gaia satellite, with the software PHOEBE. Results. HD 214220 is an eclipsing binary system with an orbital period of P = 43.14 d, eclipse depths of 17% and 13%, and masses of 2.49 M⊙ and 2.42 M⊙. The sum of the radii is R1 + R2 ≈ 8.5 R⊙, and the temperatures of the components are similar, with a ratio of T2/T1 ≈ 1.03. Conclusions. By consulting stellar evolution models, I find that the primary component has ended core hydrogen burning and is potentially in the contraction phase, prior to shell burning.

S-Process Enriched Evolved Binaries in the Galaxy and the Magellanic Clouds

Abstract Post-asymptotic giant branch stars (post-AGB) in binary systems, with typical orbital periods between $\sim\!100$ to $\sim$ 1 000 days, result from a poorly understood interaction that terminates their precursory AGB phase. The majority of these binaries display a photospheric anomaly called ‘chemical depletion’, thought to arise from an interaction between the circumbinary disc and the post-AGB star, leading to the reaccretion of pure gas onto the star, devoid of refractory elements due to dust formation. In this paper, we focus on a subset of chemically peculiar binary post-AGBs in the Galaxy and the Magellanic Clouds (MCs). Our detailed stellar parameter and chemical abundance analysis utilising high-resolution optical spectra from VLT+UVES revealed that our targets span a $T_{\rm eff}$ of 4 900–7 250 K and [Fe/H] of −0.5 - −1.57 dex. Interestingly, these targets exhibit a carbon ([C/Fe] ranging from 0.5 - 1.0 dex, dependant on metallicity) and s-process enrichment ( $\textrm{[s/Fe]}\,\geq\!1$ dex) contrary to the commonly observed chemical depletion pattern. Using spectral energy distribution (SED) fitting and period–luminosity–colour (PLC) relation methods, we determine the luminosity of the targets (2 700–8 300 $\rm L_{\odot}$ ), which enables confirmation of their evolutionary phase and estimation of initial masses (as a function of metallicity) (1–2.5 $\textrm{M}_{\odot}$ ). In conjunction with predictions from dedicated ATON stellar evolutionary models, our results indicate a predominant intrinsic enrichment of carbon and s-process elements in our binary post-AGB targets. We qualitatively rule out extrinsic enrichment and inherited s-process enrichment from the host galaxy as plausible explanations for the observed overabundances. Our chemically peculiar subset of intrinsic carbon and s-process enriched binary post-AGBs also hints at potential variation in the efficiency of chemical depletion between stars with C-rich and O-rich circumbinary disc chemistries. However, critical observational studies of circumbinary disc chemistry, along with specific condensation temperature estimates in C-rich environments, are necessary to address gaps in our current understanding of disc-binary interactions inducing chemical depletion in binary post-AGB systems.

Black Hole Ultracompact X-Ray Binaries as Galactic Low-frequency Gravitational Wave Sources: The He Star Channel

Black hole (BH) ultracompact X-ray binaries (UCXBs) are potential Galactic low-frequency gravitational wave (GW) sources. As an alternative channel, BH UCXBs can evolve from BH+He star binaries. In this work, we perform a detailed stellar evolution model for the formation and evolution of BH UCXBs evolving from the He star channel to diagnose their detectability as low-frequency GW sources. Our calculations found that some nascent BH+He star binaries after the common-envelope (CE) phase could evolve into UCXB-LISA sources with a maximum GW frequency of ∼5 mHz, which can be detected in a distance of 10 kpc (or 100 kpc). Once BH+He star systems become UCXBs through mass transfer, they would emit X-ray luminosities of ∼1038 erg s−1, making them ideal multimessenger objects. If the initial He-star masses are ≥0.7 M ⊙, those systems are likely to experience two Roche lobe overflows, and the X-ray luminosity can reach a maximum of 3.5 × 1039 erg s−1 in the second mass-transfer stage. The initial He-star masses and initial orbital periods of progenitors of Galactic BH UCXB-LISA sources are in the range of 0.32–2.9 M ⊙ and 0.02–0.19 days, respectively. Nearly all BH+He star binaries in the above parameter space can evolve into GW sources whose chirp masses can be accurately measured. Employing a population synthesis simulation, we predict the birthrate and detection number of Galactic BH UCXB-LISA sources evolving from the He star channel are R = 2.2 × 10−6 yr−1 and 33 for an optimistic CE parameter, respectively.

Is there a black hole in the center of the Sun?

There is probably not a black hole in the center of the Sun. Despite this detail, our goal in this work to convince the reader that this question is interesting and that work studying stars with central black holes is well motivated. If primordial black holes exist then they may exist in sufficiently large numbers to explain the dark matter in the universe. While primordial black holes may form at almost any mass, the asteroid-mass window between 10−16−10−10M⊙\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$10^{-16} - 10^{-10}\\ \ extrm{M}_{\\odot }$\\end{document} remains a viable dark matter candidate and these black holes could be captured by stars upon formation. Such a star, partially powered by accretion luminosity from a microscopic black hole in its core, has been called a ‘Hawking star.’ Stellar evolution of Hawking stars is highly nontrivial and requires detailed stellar evolution models, which were developed in our recent work. We present here full evolutionary models of solar mass Hawking stars using two accretion schemes: one with a constant radiative efficiency, and one that is new in this work that uses an adaptive radiative efficiency to model the effects of photon trapping.

Dynamical Masses for the Hyades Binary System vB 120

We report spectroscopic observations of vB 120 (HD 30712), a 5.7 yr astrometric–spectroscopic binary system in the Hyades cluster. We combine our radial velocities with others from the literature, and with existing speckle interferometry measurements, to derive an improved 3D orbit for the system. We infer component masses of M 1 = 1.065 ± 0.018 M ☉ and M 2 = 1.008 ± 0.016 M ☉, and an orbital parallax of 21.86 ± 0.15 mas, which we show to be more accurate than the parallax from Gaia DR3. This is the ninth binary or multiple system in the Hyades with dynamical mass determinations, and one of the examples with the highest precision. An analysis of its spectral energy distribution yields the absolute radii of the stars, R 1 = 0.968 ± 0.012 R ☉ and R 2 = 0.878 ± 0.013 R ☉, and effective temperatures of 5656 ± 56 K and 5489 ± 60 K for the primary and secondary, respectively. A comparison of these properties with the predictions of current stellar evolution models for the known age and metallicity of the cluster shows only minor differences.

ALMA detection of CO rotational line emission in red supergiant stars of the massive young star cluster RSGC1

Context. The fate of stars largely depends on the amount of mass lost during the end stages of evolution. For single stars with an initial mass between ∼8–30 M⊙, most mass is lost during the red supergiant (RSG) phase, when stellar winds deplete the H-rich envelope. However, the RSG mass-loss rate (Ṁ) is poorly understood theoretically, and so stellar evolution models rely on empirically derived mass-loss rate prescriptions. However, it has been shown that these empirical relations differ largely, with differences up to 2 orders of magnitude. Aims. We aim to derive a new mass-loss rate prescription for RSGs that is not afflicted with some uncertainties inherent in preceding studies. Methods. We have observed CO rotational line emission towards a sample of RSGs in the open cluster RSGC1 that all are of a similar initial mass. The ALMA CO(2–1) line detections allowed us to retrieve the gas mass-loss rates (ṀCO). In contrast to mass-loss rates derived from the analysis of dust spectral features (ṀSED), the data allowed us a direct determination of the wind velocity and no uncertain dust-to-gas correction factor was needed. Results. Five RSGs in RSGC1 have been detected in CO(2–1). The retrieved ṀCO values are systematically lower than ṀSED. Although only five RSGs in RSGC1 have been detected, the data allow us to propose a new mass-loss rate relation for M-type red supergiants with effective temperatures between ∼3200 and 3800 K that is dependent on the luminosity and initial mass, and that is valid during the phase where nuclear burning determines the evolution along the RSG branch. The new mass-loss rate relation is based on the new ṀCO values for the RSGs in RSGC1 and on prior ṀSED values for RSGs in four clusters, including RSGC1. The new Ṁ-prescription yields a good prediction for the mass-loss rate of some well-known Galactic RSGs that are observed in multiple CO rotational lines, including α Ori, μ Cep and VX Sgr. Moreover, there are indications that a stronger, potentially eruptive, mass-loss process is occurring during some fraction of the RSG lifetime, suggesting that RSGs might experience a phase change in mass loss leading to the wind mass-loss rate dominating the RSG evolution at that stage. Conclusions. Implementing a lower mass-loss rate in evolution codes for massive stars has important consequences as to the nature of their end-state. A reduction of the RSG mass-loss rate implies that quiescent RSG mass loss is not enough to strip a single star’s hydrogen-rich envelope. Upon core collapse such single stars would explode as RSGs. Mass-loss rates of order ∼6 times higher would be needed to strip the H-rich envelope and produce a Wolf-Rayet star while evolving back to the blue side of the Hertzsprung–Russell diagram. Future observations of a larger sample of RSGs in open clusters should allow a more stringent determination of the ṀCO–luminosity relation and a sharper diagnostic as to when the phase change in mass loss is occurring.

Atomic diffusion and mixing in old stars - VIII. Chemical abundance variations in the globular cluster M4 (NGC 6121).

Variations in chemical abundances with evolutionary phase have been identified among stars in globular and open clusters with a wide range of metallicities. In the metal-poor clusters, these variations compare well with predictions from stellar structure and evolution models considering the internal diffusive motions of atoms and ions, collectively known as atomic diffusion, when moderated by an additional mixing process with a fine-tuned efficiency. We present here an investigation of these effects in the Galactic globular cluster NGC 6121 (M4) ([Fe/H]= -1.13) through a detailed chemical abundance analysis of 86 stars using high-resolution ESO Very Large Telescope (VLT) Fibre Large Array Multi Element Spectrograph (FLAMES) spectroscopy. The stars range from the main-sequence turnoff point (TOP) to the red giant branch (RGB) just above the bump. We identify C-N-O and Mg-Al-Si abundance anticorrelations, and confirm the presence of a bimodal population differing by 1 dex in nitrogen abundance. The composition of the second-generation stars imply pollution from both massive (20-40 [Formula: see text]) and asymptotic giant branch stars. We find evolutionary variations in chemical abundances between the TOP and RGB, which are robust to uncertainties in stellar parameters and modelling assumptions. The variations are weak, but match predictions well when employing efficient additional mixing. Without correcting for Galactic production of lithium, we derive an initial lithium abundance 2.63±0.10, which is marginally lower than the predicted primordial big-bang nucleosynthesis value.

RR Lyrae from binary evolution: abundant, young, and metal-rich

ABSTRACT RR Lyrae are a well-known class of pulsating horizontal branch stars widely used as tracers of old, metal-poor stellar populations. However, mounting observational evidence shows that a significant fraction of these stars may be young and metal-rich. Here, through detailed binary stellar evolution modelling, we show that all such metal-rich RR Lyrae can be naturally produced through binary interactions. Binary companions of these RR Lyrae stars formed through binary interactions partly strip their progenitor’s envelopes during a preceding red giant phase. As a result, stripped horizontal branch stars become bluer than their isolated stellar evolution counterparts and thus end up in the instability strip. In contrast, in the single evolution scenario, the stars can attain such colours only at large age and low metallicity. While binary-made RR Lyrae can possess any ages and metallicities, their Galactic population is relatively young (1 – $9\, {\rm Gyr}$) and dominated by the thin disc and the bulge. We show that Galactic RR Lyrae from binary evolution are produced at rates compatible with the observed metal-rich population and have consistent G-band magnitudes, Galactic kinematics, and pulsation properties. Furthermore, these systems dominate the RR Lyrae population in the solar neighbourhood. We predict that all metal-rich RR Lyrae have an A, F, G, or K-type companion with a long orbital period ($P \gtrsim 1000\, {\rm d}$). Observationally characterizing the orbital periods and masses of such stellar companions will provide valuable new constraints on mass and angular momentum-loss efficiency for Sun-like accretors and the nature of RR Lyrae populations.

Surface chemical composition of single WNh stars

Context. Wolf–Rayet (WR) stars of the WNh category contain a significant fraction of hydrogen at their surface. They can be hydrogen-burning, very massive stars or stars in a post-main sequence phase of evolution. Also, WNh stars are sometimes not included in population synthesis models. Aims. We aim to better characterise the properties of single WNh stars in the Galaxy and the Magellanic Clouds. In particular, we want to constrain their surface chemistry beyond the hydrogen content by determining the helium, carbon, and nitrogen surface abundances. Methods. We perform a spectroscopic analysis of 22 single WNh stars. We fit their ultraviolet and/or optical spectra using synthetic spectra computed with the code CMFGEN. We determine the main stellar parameters (temperature, luminosity, mass-loss rates) and the surface H, He, C, and N mass fractions. We investigate the ability of current evolutionary models to reproduce all parameters at the same time. Results. We find that all WNh stars show the signatures of CNO-cycle material at their surface: they are carbon-depleted and nitrogen-rich. A clear trend of higher nitrogen content at higher metallicity is observed, as expected. The amount of hydrogen (X) varies significantly from one star to another, independently of luminosity. Values of X larger than 0.4 are not exceptional. The majority of Galactic WNh stars can be explained by evolutionary models, provided sufficient fine-tuning of the input parameters of evolutionary calculations. At lower metallicity, most stars escape predictions from evolutionary models. This has been noted in the literature but constraints on the surface nitrogen content exacerbate this severe issue. Conclusions. Our study highlights the need to refine the treatment of WR stars in both stellar evolution and population synthesis models.

Orbits and dynamical masses for the active Hyades multiple system HD 284163

ABSTRACT We report near-infrared long-baseline interferometric observations of the Hyades multiple system HD 284163, made with the Center for High Angular Resolution Astronomy array, as well as almost 43 yr of high-resolution spectroscopic monitoring at the Center for Astrophysics. Both types of observations resolve the 2.39 d inner binary, and also an outer companion in a 43.1 yr orbit. Our observations, combined with others from the literature, allow us to solve for the 3D inner and outer orbits, which are found to be at nearly right angles to each other. We determine the dynamical masses of the three stars (good to better than 1.4 per cent for the inner pair), as well as the orbital parallax. The secondary component (0.5245 ± 0.0047 M⊙) is now the lowest mass star with a dynamical mass measurement in the cluster. A comparison of these measurements with current stellar evolution models for the age and metallicity of the Hyades shows good agreement. All three stars display significant levels of chromospheric activity, consistent with the classification of HD 284163 as an RS CVn object. We present evidence that a more distant fourth star is physically associated, making this a hierarchical quadruple system.

Planetary Engulfment Prognosis within the ρ CrB System

Exoplanets have been detected around stars at various stages of their lives, ranging from young stars emerging from formation to the latter stages of evolution, including white dwarfs and neutron stars. Post-main-sequence stellar evolution can result in dramatic, and occasionally traumatic, alterations to the planetary system architecture, such as tidal disruption of planets and engulfment by the host star. The ρ CrB system is a particularly interesting case of advanced main-sequence evolution, due to the relative late age and brightness of the host star, its similarity to solar properties, and the harboring of four known planets. Here, we use stellar evolution models to estimate the expected trajectory of the stellar properties of ρ CrB, especially over the coming 1.0–1.5 billion yr as it evolves off the main sequence. We show that the inner three planets (e, b, and c) are engulfed during the red giant phase and asymptotic giant branch, likely destroying those planets via either evaporation or tidal disruption at the fluid-body Roche limit. The outer planet, planet d, is briefly engulfed by the star several times toward the end of the asymptotic giant branch, but the stellar mass loss and subsequent changing planetary orbit may allow the survival of the planet into the white dwarf phase of the stellar evolution. We discuss the implications of this outcome for similar systems and describe the consequences for planets that may lie within the habitable zone of the system.

CD –27°11535: Evidence for a Triple System in the β Pictoris Moving Group

We present new spatially resolved astrometry and photometry of the CD –27°11535 system, a member of the β Pictoris moving group consisting of two resolved K-type stars on a ∼20 yr orbit. We fit an orbit to relative astrometry measured from NIRC2, GPI, and archival NaCo images, in addition to literature measurements. However, the total mass inferred from this orbit is significantly discrepant from that inferred from stellar evolutionary models using the luminosity of the two stars. We explore two hypotheses that could explain this discrepant mass sum: a discrepant parallax measurement from Gaia due to variability, and the presence of an additional unresolved companion to one of the two components. We find that the ∼20 yr orbit could not bias the parallax measurement, but that variability of the components could produce a large-amplitude astrometric motion, an effect that cannot be quantified exactly without the individual Gaia measurements. The discrepancy could also be explained by an additional star in the system. We jointly fit the astrometric and photometric measurements of the system to test different binary and triple architectures for the system. Depending on the set of evolutionary models used, we find an improved goodness of fit for a triple system architecture that includes a low-mass (M = 0.177 ± 0.055 M ⊙) companion to the primary star. Further studies of this system will be required in order to resolve this discrepancy, either by refining the parallax measurement with a more complex treatment of variability-induced astrometric motion or by detecting a third companion.

A Lithium Depletion Age for the Carina Association

The dispersed remnants of stellar nurseries, stellar associations, provide unparalleled samples of coeval stars critical for studies of stellar and planetary formation and evolution. The Carina Stellar Association is one of the closest stellar associations to Earth, and yet measurements of its age have varied from 13 to 45 Myr. We aim to update the age of Carina using the lithium depletion boundary (LDB) method. We obtain new measurements of the Li 6708 Å absorption feature in likely members using optical spectra from the Goodman High Throughput Spectrograph on SOAR and NRES on LCO. We detect the depletion boundary at M K ≃ 6.8 (M5). This age is consistent within uncertainties across six different models, including those that account for magnetic fields and spots. We also estimate the age through analysis of the group’s overall variability, and by comparing the association members’ color–magnitude diagram to stellar evolutionary models using a Gaussian Mixture Model, recovering ages consistent with the LDB. Combining these age measures we obtain an age for the Carina association of Myr. The resulting age agrees with the older end of previous age measurements and is consistent with the lithium depletion age for the neighboring Tucana-Horologium moving group.

Inferring Coupling Strengths of Mixed-mode Oscillations in Red Giant Stars Using Deep Learning

Asteroseismology is a powerful tool that may be applied to shed light on stellar interiors and stellar evolution. Mixed modes, behaving like acoustic waves in the envelope and buoyancy modes in the core, are remarkable because they allow for probing the radiative cores and evanescent zones of red giant stars. Here, we have developed a neural network that can accurately infer the coupling strength, a parameter related to the size of the evanescent zone, of solar-like stars in ∼5 ms. In comparison with existing methods, we found that only ∼43% of inferences were in agreement with a difference less than 0.03 in a sample of ∼1700 Kepler red giants. To understand the origin of these differences, we analyzed a few of these stars using independent techniques such as the Monte Carlo Markov Chain method and echelle diagrams. Through our analysis, we discovered that these alternate techniques are supportive of the neural-net inferences. We also demonstrate that the network can be used to yield estimates of coupling strength and period spacing in stars with structural discontinuities. Our findings suggest that the rate of decline in the coupling strength in the red giant branch is greater than previously believed. These results are in closer agreement with calculations of stellar-evolution models than prior estimates, further underscoring the remarkable success of stellar evolution theory and computation. Additionally, we show that the uncertainty in measuring period spacing increases rapidly with diminishing coupling strength.

Probing the physics in the core boundary layers of the double-lined B-type binary KIC 4930889 from its gravito-inertial modes

Context. Stellar evolution models of B-type stars are still uncertain in terms of internal mixing properties, notably in the area between the convective core and the radiative envelope. This impacts age determination of such stars in addition to the computation of chemical yields produced at the end of their life. Aims. We investigated the thermal and chemical structure and rotation rate in the near-core boundary layer of the double-lined B-type binary KIC 4930889 from its four-year Kepler light curve, ground-based spectroscopy, and Gaia astrometry. Methods. We computed grids of 1D stellar structure and evolution models for different mixing profiles and prescriptions of the temperature gradient in the near-core region. We examined the preferred prescription and the near-core rotation rate using 22 prograde dipole modes detected by Kepler photometry of KIC 4930889. We employed a Mahalanobis distance merit function and considered various nested stellar model grids, rewarding goodness of fit but penalising model complexity. Results. We were able to constrain the near-core rotation rate of the pulsator to Ωrot = 0.73−0.06+0.02 d−1. Furthermore, we found a preference for either an exponentially decaying mixing profile in the near-core region or absence of additional near-core mixing, but found no preference among the various options for the temperature gradient in this region. The frequency (co)variances of our theoretical predictions are much larger than the errors on the observed frequencies. This forms the main limitation on further constraining the individual parameters of our models. A combination of spectroscopic, astrometric, binary, and asteroseismic information was used to achieve these constraints. Additionally, non-adiabatic pulsation computations of our best models indicate a need for opacity enhancements to accurately reproduce the observed mode excitation. Conclusions. The eccentric close binary system KIC 4930889 proves to be a promising target to investigate additional physics in close binaries by developing new modelling methods with the capacity to include the effect of tidal interactions for full exploitation of all detected oscillation modes.

The SN 2023ixf Progenitor in M101. I. Infrared Variability

Observational evidence points to a red supergiant (RSG) progenitor for SN 2023ixf. The progenitor candidate has been detected in archival images at wavelengths (≥0.6 μm) where RSGs typically emit profusely. This object is distinctly variable in the infrared (IR). We characterize the variability using pre-explosion mid-IR (3.6 and 4.5 μm) Spitzer and ground-based near-IR (JHK s ) archival data jointly covering 19 yr. The IR light curves exhibit significant variability with rms amplitudes in the range 0.2–0.4 mag, increasing with decreasing wavelength. From a robust period analysis of the more densely sampled Spitzer data, we measure a period of 1091 ± 71 days. We demonstrate using Gaussian process modeling that this periodicity is also present in the near-IR light curves, thus indicating a common physical origin, which is likely pulsational instability. We use a period–luminosity relation for RSGs to derive a value of M K = −11.58 ± 0.31 mag. Assuming a late M spectral type, this corresponds to at T eff = 3200 K and to at T eff = 3500 K. This gives an independent estimate of the progenitor’s luminosity, unaffected by uncertainties in extinction and distance. Assuming the progenitor candidate underwent enhanced dust-driven mass loss during the time of these archival observations, and using an empirical period–luminosity–based mass-loss prescription, we obtain a mass-loss rate of around (2–4) × 10−4 M ⊙ yr−1. Comparing the above luminosity with stellar evolution models, we infer an initial mass for the progenitor candidate of 20 ± 4 M ⊙, making this one of the most massive progenitors for a Type II SN detected to date.