Stellar Evolution Research Articles

Predicting gravitational wave signals from BPASS White Dwarf Binary and Black Hole Binary populations of a Milky Way-like galaxy model for LISA

Abstract Galactic white dwarf binaries (WDBs) and black hole binaries (BHBs) will be gravitational wave (GW) sources for LISA. Their detection will provide insights into binary evolution and the evolution of our Galaxy through cosmic history. Here, we make predictions of the expected WDB and BHB population within our Galaxy. We combine predictions of the compact remnant binary populations expected by stellar evolution from the detailed Binary Population and Spectral Synthesis (BPASS) code, with a Milky Way analogue galaxy model from the Feedback In Realistic Environment (FIRE) simulations. We use PhenomA and legwork to simulate LISA observations. Both packages make similar predictions that on average four Galactic BHBs and 673 Galactic WDBs are above the signal-to-noise ratio (SNR) threshold of 7 after a four-year mission. We compare these predictions to earlier results using the Binary Star Evolution (BSE) code with the same FIRE model galaxy. We find that BPASS predicts a few more LISA observable Galactic BHBs and a twentieth of the Galactic WDBs. The differences are due to the different physical assumptions that have gone into the binary evolution calculations. These results indicate that the expected population of compact binaries that LISA will detect depends very sensitively on the binary population synthesis models used and thus observations of the LISA population will provide tight constraints on our modelling of binary stars. Finally, from our synthetic populations we have created mock LISA signals that can be used to test and refine data processing methods of the eventual LISA observations.

Radiative neutron capture reaction rates for stellar nucleosynthesis

AbstractThere is a high demand for nuclear data in multidisciplinary subject like nuclear astrophysics. The two areas of nuclear physics which are most clearly related to one another are stellar evolution and nucleosynthesis. The necessity for nuclear data for astrophysical applications puts experimental methods as well as reliability and predicative ability of current nuclear models to the test. Despite recent, considerable advances, there are still significant issues and mysteries. Only a few characteristics of nuclear astrophysics are covered in the current work which include $$^{52}$$ 52 Fe(n,$$\gamma $$ γ )$$^{53}$$ 53 Fe, $$^{53}$$ 53 Fe(n,$$\gamma $$ γ )$$^{54}$$ 54 Fe, $$^{54}$$ 54 Fe(n,$$\gamma $$ γ )$$^{55}$$ 55 Fe, $$^{55}$$ 55 Fe(n,$$\gamma $$ γ )$$^{56}$$ 56 Fe, $$^{56}$$ 56 Fe(n,$$\gamma $$ γ )$$^{57}$$ 57 Fe, $$^{57}$$ 57 Fe(n,$$\gamma $$ γ )$$^{58}$$ 58 Fe and $$^{58}$$ 58 Fe(n,$$\gamma $$ γ )$$^{59}$$ 59 Fe reactions which are important in stellar nucleosynthesis. The reaction rates are calculated using nuclear statistical model. These rates are subsequently fitted to polynomials of temperature T$$_9$$ 9 in order to facilitate calculations for stellar nucleosynthesis.

Exploring the dynamical evolution of Cepheid multiplicity in star clusters and its implications for B-star multiplicity at birth

Classical Cepheid variable stars provide a unique probe of binary evolution in intermediate-mass stars over the course of several tens to hundreds of Myr. In addition, understanding binary evolution with the inclusion of cluster dynamics is desirable for obtaining a more complete picture of these stars, especially as they play a vital role in distance determinations. We studied the binary and multiple properties of Cepheids, assuming that all mid-B stars form in binaries inside star clusters. We also estimated the birth multiplicity of mid-B stars by comparing the observed multiplicity statistics of Cepheids with models based on particular assumptions. The clusters were modelled with the code, including synthetic stellar and binary evolutionary tracks. The Cepheids were identified from their position on the Hertzsprung-Russell diagram. The dynamical cluster environment results in a higher binary fraction among the Cepheids that remain in star clusters ($ 60$<!PCT!>) than among the Cepheids which have escaped to the field ($ 35$<!PCT!>). The fraction of Cepheids in triples ($ 30$<!PCT!> and $ 10$<!PCT!> in clusters and field, respectively) follows the same trend. In clusters, the binary, triple, and multiple fraction decreases with increasing cluster mass. More massive clusters have binaries of shorter orbital periods than lower mass clusters and field Cepheids. Mergers are very common with $ 30$<!PCT!> of mid-B stars not evolving to Cepheids because of the interaction with their companion. Approximately $40$ <!PCT!> of Cepheids have merged with their companion, and the merger event impacts stellar evolution, so that $ 25$<!PCT!> of all Cepheids occur at an age by more than $40$<!PCT!> different than what would be expected from their mass and the current cluster age; the expected age of Cepheids can differ from the age of their host cluster. Our models predict that one in five Cepheids is the result of a merger between stars with mass below the lower mass limit for Cepheids; in clusters, these objects occur substantially later than expected from their mass. Approximately $10$<!PCT!> of binary Cepheids have a different companion from the zero-age main sequence (ZAMS) one, and $ 3$ to $5$<!PCT!> of all Cepheids have a compact companion ($ 0.15$ <!PCT!> of all Cepheids are accompanied by a black hole). The binary fraction derived from our simulations (42<!PCT!>) underestimates the observed binary Cepheid fraction by approximately a factor of 2. This suggests that the true multiplicity fraction of B-stars at birth could be substantially larger than unity and, thus, that mid-B stars may typically form in triple and higher order systems.

Formation and Evolution of Special Stellar Types: A Comprehensive Study of Special Stars.

This article presents a comprehensive study of the formation, evolution, and distinctive characteristics of special stellar types, including protostars, super-giants, pulsars, red giants, etc. By exploring the life cycles of these stars, we gain a deeper understanding of the nuclear reactions, gravitational influences, and electromagnetic spectra that govern their development and eventual demise. Through a synthesis of observational data, theoretical models, and recent advancements in astrophysics, this research aims to elucidate the intricate processes driving the life cycles of these stellar phenomena. The findings contribute to the broader knowledge of stellar dynamics, offering insights into the mechanisms of stellar evolution and the fundamental forces shaping our universe.

The Evolution of Massive Binary Stars

Massive stars play a major role in the evolution of their host galaxies and serve as important probes of the distant Universe. It has been established that the majority of massive stars reside in close binaries and interact with their companion stars during their lifetimes. Such interactions drastically alter their life cycles and complicate our understanding of their evolution, but are also responsible for the production of interesting and exotic interaction products. ▪Extensive observation campaigns with well-understood detection sensitivities have enabled the conversion of observed properties into intrinsic characteristics, facilitating a direct comparison to theory.▪Studies of large samples of massive stars in our Galaxy and the Magellanic Clouds have unveiled new types of interaction products, providing critical constraints on the mass transfer phase and the formation of compact objects.▪The direct detection of gravitational waves has revolutionized the study of stellar mass compact objects, providing a new window to study massive star evolution. Their formation processes are, however, still unclear. The known sample of compact object mergers will increase by orders of magnitude in the coming decade, which is vastly outgrowing the number of stellar-mass compact objects detected through electromagnetic radiation.

White dwarf eccentricity fluctuation and dissipation by AGB convection

ABSTRACT Millisecond pulsars with white dwarf companions have typical eccentricities $e\sim 10^{-6}{\!-\!}10^{-3}$. The eccentricities of helium white dwarfs are explained well by applying the fluctuation–dissipation theorem to convective eddies in their red giant progenitors. We extend this theory to more massive carbon–oxygen (CO) white dwarfs with asymptotic giant branch (AGB) progenitors. Due to the radiation pressure in AGB stars, the dominant factor in determining the remnant white dwarf’s eccentricity is the critical residual hydrogen envelope mass $m_{\rm env}$ required to inflate the star to giant proportions. Using a suite of mesa stellar evolution simulations with $\Delta m_{\rm c}=10^{-3}\, {\rm M}_{\odot }$ core-mass intervals, we resolved the AGB thermal pulses and found that the critical $m_{\rm env}\propto m_{\rm c}^{-6}$. The resulting eccentricity $e\sim 3\times 10^{-3}$ is almost independent of the remnant CO white dwarf’s mass $m_{\rm c}$. Nearly all of the measured eccentricities lie below this robust theoretical limit, indicating that the eccentricity is damped during the common-envelope inspiral that follows the unstable Roche lobe overflow of the AGB star. Specifically, we focused on white dwarfs with median masses $m_{\rm c}\gt 0.6\, {\rm M}_{\odot }$. These massive white dwarfs begin their inspiral with practically identical orbital periods and eccentricities, eliminating any dependence on the initial conditions. For this sub-sample, we find an empirical relation $e\propto P^{3/2}$ between the final period and eccentricity that is much tighter than previous studies – motivating theoretical work on the eccentricity evolution during the common envelope phase. The eccentricities of lower mass CO white dwarfs may be explained by alternative formation channels.

Suggested magnetic braking prescription derived from field complexity fails to reproduce the cataclysmic variable orbital period gap

Magnetic wind braking drives the spin-down of low-mass stars and the evolution of most interacting binary stars. A magnetic braking prescription that was claimed to reproduce both the period distribution of cataclysmic variables (CVs) and the evolution of the rotation rates of low-mass stars is based on a relation between the angular momentum loss rate and magnetic field complexity. The magnetic braking model based on field complexity has been claimed to predict a detached phase that could explain the observed period gap in the period distribution of but has never been tested in detailed models of evolution. Here we fill this gap. We incorporated the suggested magnetic braking law in MESA and simulated the evolution of for different initial stellar masses and initial orbital periods. We find that the prescription for magnetic braking based on field complexity fails to reproduce observations of The predicted secondary star radii are smaller than measured, and an extended detached phase that is required to explain the observed period gap (a dearth of non-magnetic with periods sim 2$ and $ sim 3$ hours ) is not predicted. Proposed magnetic braking prescriptions based on a relation between the angular momentum loss rate and field complexity are too weak to reproduce the bloating of donor stars in derived from observations and, in contrast to previous claims, do not provide an explanation for the observed period gap. The suggested steep decrease in the angular momentum loss rate does not lead to detachment. Stronger magnetic braking prescriptions and a discontinuity at the fully convective boundary are needed to explain the evolution of close binary stars that contain compact objects. The tension between braking laws derived from the spin-down of single stars and those required to explain and other close binaries containing compact objects remains.

Towards a universal analytical model for Population III star formation: interplay between feedback and fragmentation

ABSTRACT JWST has brought us new insights into Cosmic Dawn with tentative detection of the unique signatures of metal-free Population III (Pop III) stars, such as strong He II emission, extremely blue ultraviolet spectrum, and enhanced nitrogen abundance. Self-consistent theoretical predictions of the formation rates, sites, and masses of Pop III stars are crucial for interpreting the observations, but are challenging due to complex physical processes operating over the large range of length-scales involved. One solution is to combine analytical models for the small-scale star formation process with cosmological simulations that capture the large-scale physics such as structure formation, radiation backgrounds, and baryon-dark matter streaming motion that regulate the conditions of Pop III star formation. We build an analytical model to predict the final masses of Pop III stars/clusters from the properties of star-forming clouds, based on the key results of small-scale star formation simulations and stellar evolution models. Our model for the first time considers the interplay between feedback and fragmentation and covers different modes of Pop III star formation ranging from ordinary small ($\sim\!{10{-}2000}\ \rm M_\odot$) clusters in molecular-cooling clouds to massive ($\gtrsim\!{10^{4}}\ \rm M_\odot$) clusters containing supermassive ($\sim\!{10^{4}{-}3}\times 10^{5}\ \rm M_\odot$) stars under violent collapse of atomic-cooling clouds with large gas accretion rates of $\gtrsim\!{0.1}\ \rm M_\odot \ yr^{-1}$. As an example, the model is applied to the Pop III star-forming clouds in the progenitors of typical haloes hosting high-z luminous quasars ($M_{\rm h}\sim 10^{12}\ \rm M_\odot$ at $z\sim 6$), which shows that formation of Pop III massive clusters is common ($\sim\!{20{-}70}{{\ \rm per\ cent}}$) in such biased ($\sim\!{4}\sigma$) regions, and the resulting heavy black hole seeds from supermassive stars can account for a significant fraction of observed luminous ($\gtrsim\!{10^{46}}\ \rm erg\ s^{-1}$) quasars at $z\sim 6$.

Anomalously low-mass core-He-burning star in NGC 6819 as a post-common-envelope phase product

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 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 Li-rich low-mass CHeB star KIC4937011, which is a member of the open cluster NGC 6819 (turn-off mass of $ i.e. an age of $ 2.4$ Gyr). This star has a lower mass by $ 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$ main-sequence companion does not survive. The mass of the primary star at the zero-age main sequence is in the range with a log-orbital period in the range days )$. During the CEE phase, $ 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 $ 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.

Twenty-three new Heartbeat Star systems discovered based on TESS data

ABSTRACT Heartbeat stars (HBSs) are ideal astrophysical laboratories to study the formation and evolution of binary stars in eccentric orbits and the internal structural changes of their components under strong tidal action. We discover 23 new HBSs based on TESS photometric data. The orbital parameters, including orbital period, eccentricity, orbital inclination, argument of periastron, and epoch of periastron passage of these HBSs, are derived by using a corrected version of Kumar et al.'s model based on the Markov Chain Monte Carlo (MCMC) method. The preliminary results show that these HBSs have orbital periods in the range from 2.7 to 20 d and eccentricities in the range from 0.08 to 0.70. The eccentricity-period relation of these objects shows a positive correlation between eccentricity and period and also shows the existence of orbital circularization. The Hertzsprung–Russell diagram shows that the HBSs are not all located in a particular area. The distribution of the derived parameters suggests a selection bias within the TESS survey towards HBSs with shorter periods. These objects are a very useful source to study the structure and evolution of eccentricity orbit binaries and to extend the TESS HBS catalog.

Younger age for the oldest magnetic white dwarfs

Abstract Sufficiently old white dwarfs cool down through a convective envelope that directly couples their degenerate cores to the surface. Magnetic fields may inhibit this convection by stiffening the criterion for convective instability. We consistently implemented the modified criterion in the stellar evolution code mesa, and computed the cooling of white dwarfs as a function of their mass and magnetic field B. In contrast to previous estimates, we find that magnetic fields can significantly change the cooling time t even if they are relatively weak $B^2\ll 8\pi P$, where P is the pressure at the edge of the degenerate core. Fields B ≳ 1 MG open a radiative window that decouples the core from the convective envelope, effectively lowering the luminosity to that of a fully radiative white dwarf. We identified a population of observed white dwarfs that are younger by Δt ∼ Gyr than currently thought due to this magnetic inhibition of convective energy transfer – comparable to the cooling delay due to carbon–oxygen phase separation. In volume-limited samples, the frequency and strength of magnetic fields increase with age. Accounting for magnetic inhibition is therefore essential for accurate cooling models for cosmic chronology and for determining the origin of the magnetic fields.

The rates and host galaxies of pair-instability supernovae through cosmic time: predictions from BPASS and IllustrisTNG

ABSTRACT Pair-instability supernovae (PISNe) have long been predicted to be the final fates of near-zero-metallicity very massive stars ($Z \lt Z_\odot /3$, M$_\mathrm{ZAMS} \gtrsim 140\, \text{M}_\odot$). However, no definite PISN has been observed to date, leaving theoretical modelling validation open. To investigate the observability of these explosive transients, we combine detailed stellar evolution models for PISNe formation, computed from the binary population and spectral synthesis code suite, bpass, with the star formation history of all individual computational elements in the Illustris-TNG simulation. This allows us to compute comic PISN rates and predict their host galaxy properties. Of particular importance is that IllustrisTNG galaxies do not have uniform metallicities throughout, with metal-enriched galaxies often harbouring metal-poor pockets of gas where PISN progenitors may form. Accounting for the chemical inhomogeneities within these galaxies, we find that the peak redshift of PISNe formation is $z=3.5$ instead of the value of $z=6$ when ignoring chemical inhomogeneities within galaxies. Furthermore, the rate increases by an order of magnitude from 1.9 to 29 PISN Gpc$^{-3}$ yr$^{-1}$ at $z=0$, if the chemical inhomogeneities are considered. Using state-of-the-art theoretical PISN light curves, we find an observed rate of 13.8 (1.2) visible PISNe per year for the Euclid-Deep survey, or 83 (7.3) over the 6-yr lifetime of the mission when considering chemically inhomogeneous (homogenous) systems. Interestingly, only 12 per cent of helium PISN progenitors are sufficiently massive to power a superluminous supernova event, which can potentially explain why PISN identification in time-domain surveys remains elusive and progress requires dedicated strategies.

Impacts of stellar wind and supernovae on star cluster formation: Origins of extremely high N/O ratios and multiple stellar populations

Abstract We study metal enrichment originating from stellar wind and supernovae in low-metallicity clouds by performing three-dimensional radiation hydrodynamics simulations. We find that metals ejected from stellar wind are accumulated, leading to subsequent star formation in the nitrogen-enriched gas. During this early phase, the ${\rm N/O}$ ratios are similar to observed nitrogen-enriched galaxies (${\rm [N/O]}\gtrsim 0.5$). Then, once supernovae occur, the ${\rm N/O}$ ratios decrease significantly. If the duration of star formation is comparable to the time-scale of supernovae, the mass fraction of nitrogen-enriched stars reaches half the mass of star clusters. We suggest that the mass of the star cluster needs to exceed $\sim \!10^6$ M$_{\odot }$ to have multiple populations due to stellar wind, considering the condition for massive star cluster formation and the timescales of stellar evolution.

Exploring the stellar rotation of early-type stars in the LAMOST medium-resolution survey

Context. Stellar rotation significantly shapes the evolution of massive stars, yet the interplay of mass and metallicity remains elusive, limiting our capacity to construct accurate stellar evolution models and to better estimate the impact of rotation on the chemical evolution of galaxies. Aims. Our goal is to investigate how mass and metallicity influence the rotational evolution of A-type stars on the main sequence (MS). We seek to identify deviations in rotational behaviors that could serve as new constraints for existing stellar models. Methods. Using the LAMOST Median-Resolution Survey Data Release 9, we derived stellar parameters for a population of 104 752 A-type stars. Our study focused on the evolution of surface rotational velocities and their dependence on mass and metallicity in 84 683 “normal” stars. Results. Normalizing surface rotational velocities to zero age main sequence (ZAMS) values revealed a prevailing evolutionary profile from 1.7 to 4.0 M⊙. This profile features an initial rapid acceleration until t/tMS = 0.25 ± 0.1 and potentially a second acceleration peak near t/tMS = 0.55 ± 0.1 for stars heavier than 2.5 M⊙, followed by a steady decline and a “hook” feature at the end. Surpassing theoretical expectations, the initial acceleration likely stems from a concentrated distribution of angular momentum at the ZAMS, resulting in a prolonged increase in speed. A transition phase for stars with 2.0 < M/M⊙ < 2.3 emerged, a region where evolutionary tracks remain uncertain. Stellar expansion primarily drives the spin down in the latter half of the MS, accompanied by significant influence from inverse meridional circulation. The inverse circulation becomes more efficient at lower metallicities, explaining the correlation of the slope of this deceleration phase with metallicity from –0.3 dex up to 0.1 dex. The metal-poor subsample (−0.3 dex < [M/H]< − 0.1 dex) starts with lower velocities at the ZAMS, suggesting that there is a metallicity-dependent mechanism that removes angular momentum during star formation. The proportion of fast rotators decreases with an increase in metallicity, up to log(Z/Z⊙)∼ − 0.2, a trend consistent with observations of OB-type stars found in the Small and Large Magellanic Clouds.

The impact of asteroseismically calibrated internal mixing on nucleosynthetic wind yields of massive stars

Context. Asteroseismology gives us the opportunity to look inside stars and determine their internal properties, such as the radius and mass of the convective core. Based on these observations, estimations can be made for the amount of the convective boundary mixing and envelope mixing of such stars and for the shape of the mixing profile in the envelope. However, these results are not typically included in stellar evolution models. Aims. We aim to investigate the impact of varying convective boundary mixing and envelope mixing in a range based on asteroseismic modelling in stellar models up to the core collapse, both for the stellar structure and for the nucleosynthetic yields. In this first study, we focus on the pre-explosive evolution and we evolved the models to the final phases of carbon burning. This set of models is the first to implement envelope mixing based on internal gravity waves for the entire evolution of the star. Methods. We used the MESA stellar evolution code to simulate stellar models with an initial mass of 20 M⊙ from zero-age main sequence up to a central core temperature of 109 K, which corresponds to the final phases of carbon burning. We varied the convective boundary mixing, implemented as ‘step-overshoot’, with the overshoot parameter (αov) in the range 0.05−0.4. We varied the amount of envelope mixing (log(Denv/cm2 s−1)) in the range 0−6 with a mixing profile based on internal gravity waves. To study the nucleosynthesis taking place in these stars in great detail, we used a large nuclear network of 212 isotopes from 1H to 66Zn. Results. Enhanced mixing according to the asteroseismology of main-sequence stars, both at the convective core boundary and in the envelope, has significant effects on the nucleosynthetic wind yields. This is especially the case for 36Cl and 41Ca, whose wind yields increase by ten orders of magnitude compared to those of the models without enhance envelope mixing. Our evolutionary models beyond the main sequence diverge in yields from models based on rotational mixing, having longer helium-burning lifetimes and lighter helium-depleted cores. Conclusions. We find that the asteroseismic ranges of internal mixing calibrated from core hydrogen-burning stars lead to similar wind yields as those resulting from the theory of rotational mixing. Adopting the seismic mixing levels beyond the main sequence, we find earlier transitions to radiative carbon burning compared to models based on rotational mixing because they have lower envelope mixing in that phase. This influences the compactness and the occurrence of shell mergers, which may affect the supernova properties and explosive nucleosynthesis.

The evolution of accreting population III stars at 10−6–103 M⊙ yr−1

Context. The first stars formed over five orders of magnitude in mass by accretion in primordial dark matter halos. Aims. We study the evolution of massive, very massive and supermassive primordial (Pop III) stars over nine orders of magnitude in accretion rate. Methods. We use the stellar evolution code GENEC to evolve accreting Pop III stars from 10−6–103 M⊙ yr−1 and study how these rates determine final masses. The stars are evolved until either the end central Si burning or they encounter the general relativistic instability (GRI). We also examine how metallicity affects the evolution of the star at one accretion rate. Results. At rates below ∼2.5 × 10−5 M⊙ yr−1 the final mass of the star falls below that required for pair-instability supernovae. The minimum rate required to produce black holes with masses above 250 M⊙ is ∼5 × 10−5 M⊙ yr−1, well within the range of infall rates found in numerical simulations of halos that cool via H2, ≲10−3 M⊙ yr−1. At rates of 5 × 10−5 M⊙ yr−1 to 4 × 10−2 M⊙ yr−1, like those expected for halos cooling by both H2 and Lyα, the star collapses after Si burning. At higher accretion rates the GRI triggers the collapse of the star during central H burning. Stars that grow at above these rates are cool red hypergiants with effective temperatures log(Teff) = 3.8 and luminosities that can reach 1010.5 L⊙. At accretion rates of 100–1000 M⊙ yr−1 the gas encounters the general relativistic instability prior to the onset of central hydrogen burning and collapses to a black hole with a mass of ∼106 M⊙ without ever having become a star. Conclusions. Our models corroborate previous studies of Pop III stellar evolution with and without hydrodynamics over separate, smaller ranges in accretion rate. They also reveal for the first time the critical transition rate in accretion above which catastrophic baryon collapse, like that which can occur during galaxy collisions in the high-redshift Universe, produces supermassive black holes via dark collapse.

Shocking and Mass Loss of Compact Donor Stars in Type Ia Supernovae

Type Ia supernovae arise from thermonuclear explosions of white dwarfs accreting from a binary companion. Following the explosion, the surviving donor star leaves at roughly its orbital velocity. The discovery of the runaway helium subdwarf star US 708, and seven hypervelocity stars from Gaia data, all with spatial velocities ≳900 km s−1, strongly support a scenario in which the donor is a low-mass helium star or a white dwarf. Motivated by these discoveries, we perform three-dimensional hydrodynamical simulations with the Athena++ code, modeling the hydrodynamical interaction between a helium star or helium white dwarf and the supernova ejecta. We find that ≈0.01–0.02 M ⊙ of donor material is stripped, and explain the location of the stripped material within the expanding supernova ejecta. We continue the postexplosion evolution of the shocked donor stars with the MESA code. As a result of entropy deposition, they remain luminous and expanded for ≈105–106 yr. We show that the postexplosion properties of our helium white dwarf donor agree reasonably with one of the best-studied hypervelocity stars, D6-2.

Investigating episodic mass loss in evolved massive stars

Mass loss during the red supergiant (RSG) phase plays a crucial role in the evolution of an intermediate-mass star; however, the underlying mechanism remains unknown. We aim to increase the sample of well-characterized RSGs at subsolar metallicity by deriving the physical properties of 127 RSGs in nine nearby southern galaxies. For each RSG, we provide spectral types and used MARCS atmospheric models to measure stellar properties from their optical spectra, such as the effective temperature, extinction, and radial velocity. By fitting the spectral energy distribution, we obtained the stellar luminosity and radius for 92 RSGs, finding that ~50% of them have log(L/L⊙) ≥ 5.0 and six RSGs have R ≳ 1400 R⊙. We also find a correlation between the stellar luminosity and mid-IR excess of 33 dusty variable sources. Three of these dusty RSGs have luminosities exceeding the revised Humphreys-Davidson limit. We then derived a metallicity-dependent J − Ks color versus temperature relation from synthetic photometry and two new empirical J − Ks color versus temperature relations calibrated on literature TiO and J-band temperatures. To scale our derived cool TiO temperatures to values that are in agreement with the evolutionary tracks, we derived two linear scaling relations calibrated on J-band and i-band temperatures. We find that the TiO temperatures are more discrepant as a function of the mass-loss rate, and discuss future prospects of the TiO bands as a mass-loss probe. Finally, we speculate that three hot dusty RSGs may have experienced a recent mass ejection (12% of the K-type sample) and classify them as candidate Levesque-Massey variables.

Gaia DR3 reveals the complex dynamical evolution within star clusters

Context. Star clusters, composed of stars born from the same molecular cloud, serve as invaluable natural laboratories for understanding the fundamental processes governing stellar formation and evolution. Aims. This study aims to investigate correlations between the Mean Interdistance ($ \bar{D_{\mathrm{i}}} $), Mean Closest Interdistance ($ \bar{D_{\mathrm{c}}} $) and Median Weighted Central Interdistance ($ \bar{D_{\mathrm{cc}}} $) with the age of star clusters, examining their evolutionary trends and assessing the robustness of these quantities as possible age indicators. Methods. We selected a sample of open clusters in the solar region and with a representative number of members (e.g. well populated and without outliers). The interdistances are derived from the spatial distribution of member stars within a cluster. Their evolution over time allows us to use them as age indicators for star clusters. Results. Our investigation reveals a high-significant correlation between the interdistances and cluster age. Considering the full sample of clusters between 7 and 9 kpc, the relationship is very broad. This is due to uncertainties in parallax, which increase with increasing distance. In particular, we must limit the sample to a maximum distance from the Sun of about 200 pc to avoid artificial effects on cluster shape and on the spatial distribution of their stars along the line of sight. Conclusions. By conservatively restraining the distance to a maximum of ∼200 pc, we have established a relationship between the interdistances and the age of the clusters. In our sample, the relationship is mainly driven by the internal expansion of the clusters and is marginally affected by external perturbative effects. Such relation might enhance our comprehension of cluster dynamics and might be used to derive cluster dynamical ages.

The true nature of HE 0057-5959, the most metal-poor, Li-rich star

The Li-rich stars are a class of rare objects with a surface lithium abundance, A(Li), that exceeds that of other stars in the same evolutionary stage. The origin of these stars is still debated, and valuable routes to look at include the Cameron-Fowler mechanism, a mass-transfer process in a binary system, or the engulfment of rocky planets or brown dwarfs. Metal-poor ([Fe/H]<−1 dex) stars are only a small fraction of the entire population of Li-rich stars. We observed the metal-poor ([Fe/H]=−3.95±0.11 dex) giant star HE 0057–5959 with MIKE at the Magellan Telescope, deriving A(Li)NLTE=+2.09±0.07 dex. Such an Li abundance is significantly higher, by about 1 dex, than that of other stars at the same evolutionary stage. A previous analysis of the same target suggested that its high A(Li) reflects an ongoing first-dredge-up process. We revised the nature of HE 0057-5959 by comparing its stellar parameters and A(Li) with appropriate stellar evolution models describing Li depletion due to the deepening of the convective envelope. This comparison rules out that HE 0057-5959 is caught during its first dredge-up, the latter having already ended according to the parameters of this star. Its A(Li), remarkably higher than the typical lithium plateau drawn by similar giant stars, demonstrates that HE 0057-5959 joins the class of the rare metal-poor, Li-rich stars. HE 0057-5959 is the most metal-poor, Li-rich star discovered so far. We considered different scenarios to explain this star also comparing it with the other metal-poor, Li-rich stars. No internal mixing able to activate the Cameron-Fowler mechanism is known for metal-poor stars at this evolutionary stage. The engulfment of planets is also disfavoured because such metal-poor stars should not host planets. Finally, HE 0057-5959 is one of the most Na-rich among the Li-rich stars, and we found that a strong excess of Na abundance is common to all three Li-rich stars with [Fe/H]<–3 dex. This finding could support the scenario of mass transfer from a massive companion star (able to simultaneously produce large amounts of both elements) in a binary system, even if we found no evidence of radial velocity variations.