Search

Asteroseismic masses of globular cluster (GC) stars are invaluable to investigate stellar evolution. Previously, only two GCs have been seismically studied. We present new detections of solar-like oscillations in the clusters M9 and M19, focusing on two key areas: stellar mass-loss and GC multiple populations. Using K2 photometry, we detect solar-like oscillations in stars on the red giant branch and early asymptotic giant branch. We measure an integrated mass-loss for M9 of |$0.16\pm 0.02$| (rand) |$^{+0.03}_{-0.03}$| (sys) |$~\mathrm{M}_{\odot }$| and M19 of |$0.33\pm 0.03$| (rand) |$^{+0.09}_{-0.07}$| (sys) |$~\mathrm{M}_{\odot }$|⁠. Comparing these to the mass-loss estimates from previous seismically studied clusters, we derive a preliminary relationship between stellar mass-loss and metallicity for Type I GCs. We find that the mass-loss for M19 – a Type II GC – is significantly larger, suggesting Type II clusters follow a different mass-loss-metallicity trend. We also examine the mass distributions in each evolutionary phase for evidence of a bimodality that could indicate mass differences between sub-populations. While no clear bimodality is observed, there is tentative evidence suggesting the presence of two mass populations. Classification through spectroscopic abundances into the sub-populations is needed to verify these findings. This study reinforces that asteroseismology of GC stars provides an excellent testbed for studying stellar evolution. However, to advance the field we need high-quality photometry of more GCs, a goal that could be realised with the upcoming Roman Telescope. [ABSTRACT FROM AUTHOR]

While stellar expansion after core-hydrogen exhaustion related to thermal imbalance has been documented for decades, the physical phenomenon of stellar inflation that occurs close to the Eddington limit has only come to the fore in recent years. We aim to elucidate the differences between these physical mechanisms for stellar radius enlargement, especially given that additional terms such as 'bloated' and 'puffed-up' stars have been introduced in the recent massive star literature. We employ single and binary star MESA structure and evolution models for constant mass, as well as models allowing the mass to change due to winds or binary interaction. We find cases that were previously attributed to stellar inflation in fact to be due to stellar expansion. We also highlight that while the opposite effect of expansion is contraction, the removal of an inflated zone should not be referred to as contraction but rather deflation, as the star is still in thermal balance. [ABSTRACT FROM AUTHOR]

With the ever-increasing number of light-curve solutions of contact binary systems, an increasing number of potential bright red nova progenitors are being reported. There remains, however, only one confirmed event. In this study, we undertake a comprehensive review of the orbital stability of contact binary systems, considering the effects of the stellar internal composition (metallicity) and age on the evolution of the gyration radius and its effect on the instability mass ratio of contact binaries. We find that both metallicity and age have an independent effect on orbital stability, with metal-poor and older systems being more stable. The combined effects of age and metallicity are quite profound, such that for most systems with primaries of solar mass or greater, which are halfway or more through the main-sequence lifespans have instability mass ratios at levels where the secondary component would be below the hydrogen fusion mass limit. We find that from the currently available solutions we cannot confidently assign any system as unstable. Although we identify eight potential red nova progenitors, all have methodological or astrophysical concerns, which lowers our confidence in designating any of them as potential merger candidates. [ABSTRACT FROM AUTHOR]

Recently inaugurated telescopes, such as the MeerKAT radio telescope and the upcoming Rubin Observatory Legacy Survey of Space and Time, will be able to detect millions of transient events in the night sky. Stellar collisions between white dwarfs (WDs) and main-sequence (MS) stars may be detectable among such transients. Simulations will play a key role in characterizing these events and selecting targets for follow-up. We present 3D smoothed particle hydrodynamics models of dynamical interactions between a |$0.6\ {\rm M}_{\odot }$|  WD and 0.3, 0.6, and |$1.2\ {\rm M}_{\odot }$|  MS stars within globular cluster environments. Utilizing a 34-isotope nuclear network, we investigate the energetics, gas morphologies, and mass-loss properties of these collisions for different stellar mass ratios. Our models predict an overabundance of |$^{13}$| C, |$^{15}$| N, and |$^{17}$| O isotopes relative to solar abundances. Moreover, we find that the time-scale of the collisions is too short and maximum temperatures too low for any significant hydrogen burning or triple-alpha reactions to occur. This combined with a negligible production of elements heavier than neon may be key signatures in distinguishing these events from other transient events with similar peak bolometric luminosities (⁠|$\sim\!\! 10^{38}\!\!-\!\!10^{41}$|  erg s |$^{-1}$|⁠). [ABSTRACT FROM AUTHOR]

We present an analysis of the Be star HD 212044 that reveals intriguing correlations between photometric variations, circumstellar disk emission, and frequencies. Our findings, based on data from the Transiting Exoplanet Survey Satellite (TESS) mission, show an unexpectedly strong negative correlation between photometric brightness and Hα equivalent width, challenging the established understanding of the behavior of Be star seen at low inclination angles. Notably, beating episodes precede brightening events. This study suggests that the variability of HD 212044 may be due to additional mechanisms. [ABSTRACT FROM AUTHOR]

Aims. We present contemporaneous high-angular-resolution millimeter imaging and visible polarimetric imaging of the nearby asymptotic giant branch (AGB) star W Hya to better understand the dynamics and dust formation within a few stellar radii. Methods. The star W Hya was observed in two vibrationally excited H2O lines at 268 and 251 GHz with Atacama Large Millimeter/submillimeter Array (ALMA) at a spatial resolution of 16 × 20 mas and at 748 and 820 nm at a resolution of 26 × 27 mas with the Very Large Telescope (VLT)/Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE)-Zurich Imaging Polarimeter (ZIMPOL). Results. ALMA's high spatial resolution allowed us to image strong emission of the vibrationally excited H2O line at 268 GHz (v2 = 2, JKa, Kc = 65, 2–74, 3) over the stellar surface instead of absorption against the continuum, which is expected for thermal excitation. Strong, spotty emission was also detected along and just outside the stellar disk limb at an angular distance of ∼40 mas (∼1.9 R⋆), extending to ∼60 mas (∼2.9 R⋆). Another H2O line (v2 = 2, JKa, Kc = 92, 8–83, 5) at 251 GHz with a similar upper-level energy was tentatively identified, which shows absorption over the stellar surface. This suggests that the emission over the surface seen in the 268 GHz H2O line is suprathermal or even maser emission. The estimated gas temperature and H2O density are consistent with the radiatively pumped masers. The 268 GHz H2O line reveals global infall at up to ∼15 km s−1 within 2–3 R⋆, but outflows at up to ∼8 km s−1 are also present. The polarized intensity maps obtained in the visible reveal clumpy dust clouds forming within ∼40 mas (∼1.9 R⋆) with a particularly prominent cloud in the SW quadrant and a weaker cloud in the east. The 268 GHz H2O emission overlaps very well with the visible polarized intensity maps, which suggests that both the nonthermal and likely maser H2O emission and the dust originate from dense, cool pockets in the inhomogeneous atmosphere within ∼2–3 R⋆. [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]

The viscous decretion disc (VDD) scenario has been used to model the observables of the Be stars. Its capability to predict individual observables has been confirmed for several Be stars. Here, we simultaneously analyse the spectroscopic and |$BVI$| -band photometric data for the Be star MT91-213 with the Monte Carlo radiative-transfer code hdust to determine the stellar parameters, geometry, and physical conditions for its circumstellar disc. MT91-213 is the primary component of a binary system whose companion is the pulsar PSR J2032+4127. We find that the VDD model can simultaneously reproduce the multiple observables qualitatively, but not quantitatively. We determine the mass of the primary star to be 13.1  |${\rm M}_{\odot }$| which is smaller than reported in the literature. We present a dynamical scenario for the evolving disc density from a diffuse to a dense phase. Also, we determine that the inclination of the disc is about 40 |$^{\circ }$| which means it is 20 |$^{\circ }$| tilted from the orbit of the secondary star. Our results indicate that the mass loss rate for MT91-213 is |$\sim 10^{-7}$| to |$10^{-6} {\rm M}_{\odot }\,\mathrm{yr}^{-1}$| which is in agreement with the suggested values in the literature, required to explain the observed X-ray synchrotron luminosity, |$L_\mathrm{x}$|⁠ , for PSR J2032+4127. [ABSTRACT FROM AUTHOR]

At metallicities lower than that of the Small Magellanic Cloud, it remains essentially unexplored how fossil magnetic fields, forming large-scale magnetospheres, could affect the evolution of massive stars, thereby impacting the fundamental building blocks of the early Universe. We extend our stellar evolution model grid with representative calculations of main-sequence, single-star models with initial masses of 20 and 60 M |$_\odot$|⁠ , including appropriate changes for low-metallicity environments (⁠|$Z = 10^{-3}$| – |$10^{-6}$|⁠). We scrutinize the magnetic, rotational, and chemical properties of the models. When lowering the metallicity, the rotational velocities can become higher and the tendency towards quasi-chemically homogeneous evolution increases. While magnetic fields aim to prevent the development of this evolutionary channel, the weakening stellar winds lead to less efficient magnetic braking in our models. Since the stellar radius is almost constant during a blueward evolution caused by efficient chemical mixing, the surface magnetic field strength remains unchanged in some models. We find core masses at the terminal-age main sequence between 22 and 52 M |$_\odot$| for initially 60 M |$_\odot$| models. This large difference is due to the vastly different chemical and rotational evolution. We conclude that in order to explain chemical species and, in particular, high nitrogen abundances in the early Universe, the adopted stellar models need to be under scrutiny. The assumptions regarding wind physics, chemical mixing, and magnetic fields will strongly impact the model predictions. [ABSTRACT FROM AUTHOR]

Context. Despite mounting evidence that extreme red supergiants (RSGs) undergo episodic mass-loss events, their role in RSG evolution remains uncertain. Critical questions remain unanswered, such as whether or not these events can strip the star, and their timescale and frequency. Aims. This study delves into [W60] B90, one of the most luminous and extreme RSGs in the Large Magellanic Cloud (LMC), with our aim being to search for evidence of episodic mass loss. Our discovery of a bar-like nebular structure at 1 pc, which is reminiscent of the bar around Betelgeuse, raised the question of whether [W60] B90 also has a bow shock, motivating the present study. Methods. We collected and analyzed proper motion data from Gaia, as well as new multi-epoch spectroscopic and imaging data, and archival time-series photometry in the optical and mid-infrared (MIR). We used MARCS models to derive the physical properties of the star from the spectra. Results. We find [W60] B90 to be a walkaway star, with a supersonic peculiar velocity in the direction of the bar. We detect shocked emission between the bar and the star, based on the [S II]/Hα > 0.4 criterion, providing strong evidence for a bow shock. The 30 yr optical light curve reveals semi-regular variability, showing three similar dimming events with ΔV ~ 1 mag, a recurrence of ~12 yr, and a rise time of 400 days. We find the MIR light curve to vary by 0.51 mag and 0.37 mag in the WISE1 and WISE2 bands, respectively, and by 0.42 mag and 0.25 mag during the last dimming event. During this event, optical spectroscopy reveals spectral variability (M3 I to M4 I), a correlation between the Teff and the brightness, increased extinction, and, after the minimum, spectral features incompatible with the models. We also find a difference of >300 K between the Teff measured from the TiO bands in the optical and the atomic lines from our J-band spectroscopy. Conclusions. [W60] B90 is a more massive analog of Betelgeuse in the LMC and therefore the first single extragalactic RSG with a suspected bow shock. Its high luminosity of log(L/L⊙) = 5.32 dex, mass-loss rate, and MIR variability compared to other RSGs in the LMC indicate that it is in an unstable evolutionary state, undergoing episodes of mass loss. Investigating other luminous and extreme RSGs in low-metallicity environments using both archival photometry and spectroscopy is crucial to understanding the mechanism driving episodic mass loss in extreme RSGs in light of the Humphreys-Davidson limit and the "RSG problem". [ABSTRACT FROM AUTHOR]

Context. The outcome of the formation of massive stars is an important anchor point in the formation and evolution process of these stars. It provides insight into the physics of the assembly process, and sets the conditions for stellar evolution. For massive stars, the outcome of formation is rarely observed because the processes involved unfold deep down in highly extincted molecular clouds. Aims. We characterize a population of highly reddened stars in the very young massive star-forming region M17. The group of 18 O4.5 to B9 stars constitutes one of the best samples of almost zero-age main-sequence (ZAMS) high- and intermediate-mass stars. Their properties allow us to identify the empirical location of the ZAMS of massive stars, and the rotation and mass-loss rate of stars close to or at the onset of core-hydrogen burning. Methods. We performed quantitative spectroscopic modeling of a uniform set of over 100 spectral features in optical VLT/X-shooter spectra using the nonlocal thermal equilibrium stellar atmosphere code FASTWIND and a fitting approach based on a genetic algorithm, KIWI-GA. The spectral energy distributions of photometric observations were used to determine the line-of-sight extinction. From a comparison of their positions in the Hertzsprung-Russell diagram with MIST evolutionary tracks, we inferred the stellar masses and ages. Results. We find an age of 0.4−0.2+0.6 Myr for our sample, however we also identify a strong relation between the age and the mass of the stars. All sources are highly reddened, with AV ranging from 3.6 to 10.6 mag. The sample can be subdivided into two groups. Stars more massive than 10 M⊙ have reached the ZAMS. Their projected ZAMS spin rate distribution extends to 0.3 of the critical velocity; their mass-loss rates agree with those of other main-sequence O and early-B stars. Stars with a mass in the range 3 < M/M⊙ < 7 are still on the pre-main sequence (PMS), and some of them have circumstellar disks. Evolving their υ sin i to the ZAMS assuming angular momentum conservation yields values up to ~0.6 υcrit. For PMS stars without disks, we find tentative mass-loss rates up to 10−8.5 M⊙ yr−1. The total-to-selective extinction RV is higher for PMS stars with disks than for the remainder of the sample. Conclusions. We constrain the empirical location of the ZAMS for massive (10 < M/M⊙ < 50) stars and find it to agree with its location in MIST evolutionary tracks. The ZAMS rotation rates for intermediate-mass stars are twice as high as for massive stars, suggesting that the angular momentum gain processes differ between the two groups. The relation between the age and mass of the stars suggests a lag in the formation of more massive stars relative to lower mass stars. Taking the derived mass-loss rates at face value, stellar winds are already initiated in the PMS phase. The PMS-star winds are found to be substantially more powerful than indicated by predictions for line-driven outflows. [ABSTRACT FROM AUTHOR]

Context. Asymptotic giant branch (AGB) stars are significant contributors to the metal enrichment of the interstellar medium. They have strong dust-driven winds that have their origin in regions close to the AGB star's surface, where dense dust clouds form. Aims. In this methods paper, we adapted models from advanced radiation-hydrodynamical (RHD) simulations as input for radiative transfer software to create synthetic observables. A major goal is to describe an AGB star's non-sphericity and to simulate its effects on the surrounding dusty envelope. Methods. We developed tools in Python to translate models of an AGB star and its dust-driven wind from 3D RHD simulations with CO5BOLD into the format used for radiative transfer with RADMC-3D. We preserved the asymmetric shape of the AGB star by including the star as a 'dust species' and by using temperature data computed in CO5BOLD. The circumstellar silicate dust from the 3D RHD simulation is included using Mg2SiO4 opacity data in RADMC-3D with spatially dependent grain sizes. We compared images and spectral energy distributions (SEDs) created with RADMC-3D of a model snapshot with similar output made with a spherically symmetric stellar atmosphere from the 1D program DARWIN and with a point source star in RADMC-3D. Results. Our CO5BOLD model features substantial and clumpy dust formation just above 3.4 au from the grid centre (∼1 R⋆ above the star), and large-scale structures due to giant convection cells are visible on the stellar surface. With the properties of VLTI as a basis, we have created simple synthetic observables where the dust clouds close to the star and features on the stellar surface are resolved. The flux density and the contrast to the star are high enough that optical interferometers, such as the VLTI, should be able to detect these dust clouds. We find that it is important to include asymmetric stellar models since their irregular shapes, radiation fields, and their dusty envelopes even put their marks on spatially unresolved observables and affect the flux levels and shapes of the SEDs. The effects on flux levels can mostly be linked to the clumpiness of the circumstellar dust. In contrast, the angle-dependent illumination resulting from temperature variations on the stellar surface causes shifts in the wavelengths of the flux maximum, as shown by replacing the asymmetric star with a spherical one. Conclusions. The methods presented here are an important step towards producing realistic synthetic observables and testing predictions of advanced 3D RHD models. With the model used here, we find that optical interferometers should be able to resolve thermal emission from dense clouds in the dust-formation zone close to an AGB star. Taking the angle-dependence of SEDs as a proxy for temporal variations in unresolved data, we conclude that not all variability observed in AGB stars should be interpreted as global changes in the sense of spherical models. [ABSTRACT FROM AUTHOR]

Context. The nature of mass-loss in massive stars is one of the most important and difficult to constrain processes in the evolution of massive stars. The largest observational uncertainties are related to the influence of metallicity and wind structure with optically thick clumps. Aims. We aim to constrain the wind parameters of sample of 18 O-type stars in the LMC, through analysis with stellar atmosphere and wind models including the effects of optically thick clumping. This will allow us to determine the most accurate spectroscopic mass-loss and wind structure properties of massive stars at sub-solar metallicity to date. This will allow us to gain insight into the impact of metallicity on massive stellar winds. Methods. Combining high signal to noise (S/N) ratio observations in the ultraviolet and optical wavelength ranges gives us access to diagnostics of multiple different ongoing physical processes in the stellar wind. We produce synthetic spectra using the stellar atmosphere modelling code FASTWIND, and reproduce the observed spectra using a genetic algorithm based fitting technique to optimise the input parameters. Results. We empirically constrain 15 physical parameters associated with the stellar and wind properties of O-type stars from the dwarf, giant and supergiant luminosity classes. These include temperature, surface gravity, surface abundances, rotation, macroturbulence and wind parameters. Conclusions. We find, on average, mass-loss rates a factor of 4–5 lower than those from theoretical predictions commonly used in stellar-evolution calculations, but in good agreement with more recent theoretical predictions. In the 'weak-wind' regime we find massloss rates orders of magnitude below any theoretical predictions. We find a positive correlation of clumping factors with effective temperature with an average fcl = 14 ± 8 for the full sample. It is clear that there is a difference in the porosity of the wind in velocity space, and interclump density, above and below a temperature of roughly 38 kK. Above 38 kK an average 46 ± 24% of the wind velocity span is covered by clumps and the interclump density is 10–30% of the mean wind. Below an effective temperature of roughly 38 kK there must be additional light leakage for supergiants. For dwarf stars at low temperatures there is a statistical preference for very low clump velocity spans, however it is unclear if this can be physically motivated as there are no clearly observable wind signatures in UV diagnostics. [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]

In our Galactic Centre, about |$10\,000$| to |$100\,000$| stars are estimated to have survived tidal disruption events, resulting in partially disrupted remnants. These events occur when a supermassive black hole (SMBH) tidally interacts with a star, but not enough to completely disrupt the star. We use the 1D stellar evolution code Kepler and the 3D smoothed particle hydrodynamics code Phantom to model the tidal disruption of 1, 3, and |$10\, \mathrm{M}_\odot$| stars at zero-age main sequence (ZAMS), middle-age main sequence (MAMS), and terminal-age main sequence (TAMS). We map the disruption remnants into Kepler in order to understand their post-distribution evolution. We find distinct characteristics in the remnants, including increased radius, rapid core rotation, and differential rotation in the envelope. The remnants undergo composition mixing that affects their stellar evolution. Although the remnants formed by disruption of ZAMS models evolve similarly to unperturbed models of the same mass, for MAMS and TAMS stars, the remnants have higher luminosity and effective temperature. Potential observational signatures include peculiarities in nitrogen and carbon abundances, higher luminosity, rapid rotation, faster evolution, and unique tracks in the Hertzsprung–Russell diagram. [ABSTRACT FROM AUTHOR]

The properties of blue supergiants are key for constraining the end of the main sequence (MS) of massive stars. Whether the observed drop in the relative number of fast-rotating stars below ≈21 kK is due to enhanced mass-loss rates at the location of the bistability jump, or the result of the end of the MS is still debated. Here, we combine newly derived estimates of photospheric and wind parameters with Gaia distances and wind terminal velocities from the literature to obtain upper limits on the mass-loss rates for a sample of 116 Galactic luminous blue supergiants. The parameter space covered by the sample ranges between 35–15 kK in Teff and 4.8–5.8 dex in log(L/L⊙). Our results show no increase in the mass-loss rates over the bistability jump. Therefore, we argue that the drop in rotational velocities cannot be explained by enhanced mass loss. Since a large jump in the mass-loss rates is commonly included in evolutionary models, we suggest an urgent revision of the default prescriptions currently in use. [ABSTRACT FROM AUTHOR]

Context. Radio surveys of early-type stars have revealed a number of non-thermal emitters. Most of these have been shown to be binaries, where the collision between the two stellar winds is responsible for the non-thermal emission. Aims. HD 168112 is a non-thermal radio emitter, whose binary nature has only recently been confirmed spectroscopically. We obtained independent spectroscopic observations to determine its orbit, in addition to radio observations to see if the thermal or non-thermal nature of the emission changes during the periastron passage. Methods. We monitored HD 168112 spectroscopically for a 13 yr time span. From these data, we determined the orbital parameters, which we compared to the previous results in the literature. The stellar parameters of both components were determined by comparing the spectra to TLUSTY models. From the spectral index of the radio observations, we found how the nature of the emission changes as the system goes through periastron. Combining our results with other literature data allowed us to further constrain the orbital and stellar parameters. Results. We find HD 168112 to have an orbital period of P = 512.17−0.11+0.41 days, an eccentricity of e = 0.7533−0.0124+0.0053, and a mass ratio close to one. From our spectroscopic modelling, we derived the stellar parameters, but we had difficulty arriving at a spectroscopic mass ratio of one. The radio observations around periastron show only thermal emission, suggesting that most of the synchrotron photons are absorbed in the two stellar winds at that phase. Combining our data with the optical interferometry detection, we could constrain the inclination angle to i ~ 63°, and the mass of each component to ~26 M⊙. Conclusions. We have provided an independent spectroscopic confirmation of the binary nature of HD 168112. Although detected as a non-thermal radio emitter, near periastron the radio emission of this highly eccentric system is thermal and is mainly formed in the colliding-wind region. This effect will also occur in other colliding-wind binaries. [ABSTRACT FROM AUTHOR]

Aims. We aim to measure the impact of our mass-loss recipe in the evolution of massive stars at the metallicity of the Small Magellanic Cloud (SMC). Methods. We used the Geneva-evolution code (GENEC) to run evolutionary tracks for stellar masses ranging from 20 to 85 M⊙ at SMC metallicity (ZSMC = 0.002). We upgraded the recipe for stellar winds by replacing Vink's formula with our self-consistent m-CAK prescription, which reduces the value of the mass-loss rate, Ṁ, by a factor of between two and six depending on the mass range. Results. The impact of our new [weaker] winds is wide, and it can be divided between direct and indirect impact. For the most massive models (60 and 85 M⊙) with Ṁ ≳ 2 × 10−7M⊙ yr−1, the impact is direct because lower mass loss make stars remove less envelope, and therefore they remain more massive and less chemically enriched at their surface at the end of their main sequence (MS) phase. For the less massive models (20 and 25 M⊙) with Ṁ ≲ 2 × 10−8M⊙ yr−1, the impact is indirect because lower mass loss means the stars keep high rotational velocities for a longer period of time, thus extending the H-core burning lifetime and subsequently reaching the end of the MS with higher surface enrichment. In either case, given that the conditions at the end of the H-core burning change, the stars will lose more mass during their He-core burning stages anyway. For the case of Mzams = 20–40 M⊙, our models predict stars will evolve through the Hertzsprung gap, from O-type supergiants to blue supergiants (BSGs), and finally red supergiants (RSGs), with larger mass fractions of helium compared to old evolution models. New models also sets the minimal initial mass required for a single star to become a Wolf-Rayet (WR) at metallicity Z = 0.002 at Mzams = 85 M⊙. Conclusions. These results reinforce the importance of upgrading mass-loss prescriptions in evolution models, in particular for the earlier stages of stellar lifetime, even for Z ≪ Z⊙. New values for Ṁ need to be complemented with upgrades in additional features such as convective-core overshooting and distribution of rotational velocities, besides more detailed spectroscopical observations from projects such as XShootU, in order to provide a robust framework for the study of massive stars at low-metallicity environments. [ABSTRACT FROM AUTHOR]

Context. The chemical evolution of asymptotic giant branch (AGB) stars is driven by repeated thermal pulses (TPs). The duration of a TP is only a few hundred years, whereas an inter-pulse period lasts 104 − 105 yr. Direct observations of TPs are hence unlikely. However, the detached shells seen in CO line emission that are formed as a result of a TP provide indirect constraints on the changes experienced by the star during the pulse. Aims. We aim to resolve the spatial and kinematic sub-structures in five detached-shell sources to provide detailed constraints for hydrodynamic models that describe the formation and evolution of the shells. Methods. We used observations of the 12CO (1 − 0) emission towards five carbon-AGB stars with ALMA (Atacama Large Millimeter/submillimeter Array), including previously published observations of the carbon AGB star U Ant. The data have angular resolutions of 0″​​.3 to 1″ and a velocity resolution of 0.3 km s−1. This enabled us to quantify spatial and kinematic structures in the shells. Combining the ALMA data with single-dish observations of the 12CO (1 − 0) to 12CO (4 − 3) emission towards the sources, we used radiative transfer models to compare the observed structures with previous estimates of the shell masses and temperatures. Results. The observed emission is separated into two distinct components: a more coherent, bright outer shell and a more filamentary, fainter inner shell. The kinematic information shows that the inner sub-shells move at a higher velocity relative to the outer sub-shells. The observed sub-structures reveal a negative velocity gradient outwards across the detached shells, confirming the predictions from hydrodynamical models. However, the models do not predict a double-shell structure, and the CO emission likely only traces the inner and outer edges of the shell, implying a lack of CO in the middle layers of the detached shell. Previous estimates of the masses and temperatures are consistent with originating mainly from the brighter subshell, but the total shell masses are likely lower limits. Also, additional structures in the form of partial shells outside the detached shell around V644 Sco, arcs within the shell of R Scl, and a partially filled shell for DR Ser indicate a more complicated evolution of the shells and mass-loss process throughout the TP cycle than previously assumed. Conclusions. The observed spatial and kinematical splittings of the shells appear consistent with results from the hydrodynamical models, provided the CO emission does not trace the H2 density distribution in the shell but rather traces the edges of the shells. The hydrodynamical models predict very different density profiles depending on the evolution of the shells and the different physical processes involved in the wind-wind interaction (e.g. heating and cooling processes). It is therefore not possible to constrain the total shell mass based on the CO observations alone. Additional features outside and inside the shells complicate the interpretation of the data. Complementary observations of, for example, CI as a dissociation product of CO would be necessary to understand the distribution of CO compared to H2, in addition to new detailed hydrodynamical models of the pre-pulse, pulse, and post-pulse wind. Only a comprehensive combination of observations and models will allow us to constrain the evolution of the shells and the changes in the star during the thermal-pulse cycle. [ABSTRACT FROM AUTHOR]

Massive stars in the red supergiant (RSG) phase are known to undergo strong mass-loss through winds and observations indicate that a substantial part of this mass-loss could be driven by localized and episodic outflows. Various mechanisms have been considered to explain this type of mass-loss in RSGs, but these models often focus on single-star evolution. However, massive stars commonly evolve in binary systems, potentially interacting with their companions. Motivated by observations of the highly asymmetric circumstellar ejecta around the RSG VY CMa, we investigate a scenario where a companion on an eccentric orbit grazes the surface of an RSG at periastron. The companion ejects part of the outer RSG envelope, which radiatively cools, reaching the proper conditions for dust condensation and eventually giving rise to dust-driven winds. Using simple treatments for radiative cooling and dust-driven winds, we perform three-dimensional smoothed particle hydrodynamic simulations of this scenario with a |$20\, {\rm M}_\odot$| RSG and a |$2\, {\rm M}_\odot$| companion. We follow the evolution of the binary throughout a total of 14 orbits and observe that the orbit tightens after each interaction, in turn enhancing the mass-loss of subsequent interactions. We show that one such grazing interaction yields outflows of |$3\times 10^{-4}\, {\rm M}_\odot$|⁠ , which later results in wide asymmetric dusty ejecta, carrying a total mass of |$0.185\, {\rm M}_\odot$| by the end of simulations. We discuss the implications for the evolution of the binary, potential observational signatures, as well as future improvements of the model required to provide sensible predictions for the evolution of massive binaries. [ABSTRACT FROM AUTHOR]

Be stars are rapidly rotating B-type stars surrounded by circumstellar disks formed from self-ejected material. Understanding the mechanisms driving mass ejection and disk formation, known as the Be phenomenon, requires a detailed investigation of their variability and underlying physical processes. In this study, we analyze the photometric, spectroscopic, and seismic characteristics of three Be stars—HD 212044, 28 Cyg, and HD 174237—using high-cadence data from the TESS mission and spectral data from the BeSS database. Photometric variability was analyzed through iterative prewhitening and wavelet techniques, revealing distinct frequency groups associated with non-radial pulsations (NRPs). Spectral data provided equivalent width measurements of the Hα line, which correlated with photometric changes, reflecting dynamic interactions between the stars and their disks. Seismic analysis identified core rotation rates and buoyancy travel times for HD 212044 and 28 Cyg, offering insights into internal stellar processes and angular momentum distribution. HD 212044 exhibits a strong negative correlation between photometric brightness and Hα equivalent width, whereas this correlation is weaker in the case of 28 Cyg. The findings for these two stars highlight the interplay between NRPs, rapid rotation, and circumstellar disk dynamics. In contrast, the case of HD 174237 presents an example of how a binary system with mass transfer and a B-type component is revealed when observed simultaneously with space-based photometry and ground-based spectroscopy, demonstrating the importance of distinguishing classical Be stars from interacting binaries.

We investigate the secular evolution of a four-body planetary system, where two planets move around a binary star configuration on quasi-elliptic orbits. It is assumed that the masses of all bodies can change isotropically at different rates and the bodies attract each other according to Newton's law of universal gravitation. The purpose of this work is to investigate an influence of variability of the masses of the binary stars and their planets on the dynamical evolution of the planetary system. We consider the case of small eccentricities and inclinations of the bodies orbits and assume that their orbits do not intersect during evolution. Differential equations of the perturbed motion in the osculating analogues of canonical Poincaré elements were obtained in a general case of the many-body problem with variable masses in our previous work. Here we solve these equations numerically and investigate the secular evolution of some fictitious circumbinary 2-planet system under assumption that the two stars of the binary system lose their masses independently at different rates. In order to demonstrate the dynamical features of the equations we use the known parameters of the TOI-1338 system. Comparing the results of calculations in the cases of constant and variable masses, we show that the isotropic variability of the masses of bodies can influence substantially the secular perturbation of the orbital elements. [ABSTRACT FROM AUTHOR]

Context. Betelgeuse, a red supergiant star of semi-regular variability, reached a historical minimum brightness in February 2020, known as the Great Dimming. Even though the brightness has returned to the values prior to the Great Dimming now, it continues to exhibit highly unusual behavior. Aims. Understanding the long-term atmospheric motions of Betelgeuse and its variability could be a clue to the nature of the Great Dimming and the mass-loss process in red supergiants. Our goal is to study long-term dynamics of the photosphere, including during the Great Dimming. Methods. We applied the tomographic method, which allows different layers in the stellar atmosphere to be probed in order to reconstruct depth-dependent velocity fields. The method is based on the construction of spectral masks by grouping spectral lines from specific optical depths. These masks are cross-correlated with the observed spectra to recover the velocity field inside each atmospheric layer. Results. We obtained about 2800 spectra over the past 15 yr that were observed with the STELLA robotic telescope in Tenerife. We analyzed the variability of five different layers of Betelgeuse's photosphere. We found phase shift between the layers, as well as between the variability of velocity and photometry. The time variations of the widths of the cross-correlation function reveal propagation of two shockwaves during the Great Dimming. For about 2 yr after the dimming, the timescale of variability was different between the inner and outer photospheric layers. By 2022, all the layers seemingly started to follow a similar behavior as before the dimming, but pulsating with higher frequency corresponding with the first overtone. Conclusions. The combination of the extensive high-resolution spectroscopic data set with the tomographic method revealed the variable velocity fields in the photosphere of Betelgeuse, for the first time in such detail. We were also able to find new insights related to the Great Dimming event and its aftermath, namely the discovery of another shockwave and the subsequent rearrangement of the photosphere. Our results demonstrate that powerful shocks are the triggering mechanism for episodic mass-loss events, which may be the missing component to explain the mass-loss process in red supergiants. [ABSTRACT FROM AUTHOR]

Context. The mass loss experienced on the asymptotic giant branch (AGB) at the end of the lives of low- and intermediate-mass stars is widely accepted to rely on radiation pressure acting on newly formed dust grains. Dust formation happens in the extended atmospheres of these stars, where the density, velocity, and temperature distributions are strongly affected by convection, stellar pulsation, and heating and cooling processes. The interaction between these processes and how that affects dust formation and growth is complex. Hence, characterising the extended atmospheres empirically is paramount to advance our understanding of the dust formation and wind-driving processes. Aims. We aim to determine the density, temperature, and velocity distributions of the gas in the extended atmosphere of the AGB star R Dor. Methods. We acquired observations using ALMA towards R Dor to study the gas through molecular line absorption and emission. We modelled the observed 12CO v = 0,  J = 2 − 1, v = 1,  J = 2 − 1, and 3 − 2 and 13CO v = 0,  J = 3 − 2 lines using the 3D radiative transfer code LIME to determine the density, temperature, and velocity distributions up to a distance of four times the radius of the star at sub-millimetre wavelengths. Results. The high angular resolution of the sub-millimetre maps allows for even the stellar photosphere to be spatially resolved. By analysing the absorption against the star, we infer that the innermost layer in the near-side hemisphere is mostly falling towards the star, while gas in the layer above that seems to be mostly outflowing. Interestingly, the high angular resolution of the ALMA Band 7 observations reveal that the velocity field of the gas seen against the star is not homogenous across the stellar disc. The gas temperature and density distributions have to be very steep close to the star to fit the observed emission and absorption, but they become shallower for radii larger than ∼1.6 times the stellar sub-millimetre radius. While the gas mass in the innermost region is hundreds of times larger than the mass lost on average by R Dor per pulsation cycle, the gas densities just above this region are consistent with those expected based on the mass-loss rate and expansion velocity of the large-scale outflow. Our fits to the line profiles require the velocity distribution on the far side of the envelope to be mirrored, on average, with respect to that on the near side. Using a sharp absorption feature seen in the CO v = 0,  J = 2 − 1 line, we constrained the standard deviation of the stochastic velocity distribution in the large-scale outflow to be ≲0.4 km s−1. We characterised two blobs detected in the CO v = 0,  J = 2 − 1 line and found densities substantially larger than those of the surrounding gas. The two blobs also display expansion velocities that are high relative to that of the large-scale outflow. Monitoring the evolution of these blobs will lead to a better understanding of the role of these structures in the mass-loss process of R Dor. [ABSTRACT FROM AUTHOR]

Context. The features in the light curves and spectra of many Type I and Type II supernovae (SNe) can be understood by assuming an interaction of the SN ejecta with circumstellar matter (CSM) surrounding the progenitor star. This suggests that many massive stars may undergo various degrees of envelope stripping shortly before exploding, and may therefore produce a considerable diversity in their pre-explosion CSM properties. Aims. We explore a generic set of about 100 detailed massive binary evolution models in order to characterize the amount of envelope stripping and the expected CSM configurations. Methods. Our binary models were computed with the MESA stellar evolution code, considering an initial primary star mass of 12.6 M⊙ and secondaries with initial masses of between ∼12 M⊙ and ∼1.3 M⊙, and focus on initial orbital periods above ∼500 d. We compute these models up to the time of iron core collapse in the primary. Results. Our models exhibit varying degrees of stripping due to mass transfer, resulting in SN progenitor models ranging from fully stripped helium stars to stars that have not been stripped at all. We find that Roche lobe overflow often leads to incomplete stripping of the mass donor, resulting in a large variety of pre-SN envelope masses. In many of our models, the red supergiant (RSG) donor stars undergo core collapse during Roche lobe overflow, with mass transfer and therefore system mass-loss rates of up to 0.01 M⊙ yr−1 at that time. The corresponding CSM densities are similar to those inferred for Type IIn SNe, such as SN 1998S. In other cases, the mass transfer becomes unstable, leading to a common-envelope phase at such late time that the mass donor explodes before the common envelope is fully ejected or the system has merged. We argue that this may cause significant pre-SN variability, as witnessed for example in SN 2020tlf. Other models suggest a common-envelope ejection just centuries before core collapse, which may lead to the strongest interactions, as observed in superluminous Type IIn SNe, such as SN 1994W and SN 2006gy. Conclusions. Wide massive binaries exhibit properties that may not only explain the diverse envelope stripping inferred in Type Ib, IIb, IIL, and IIP SNe, but also offer a natural framework to understand a broad range of hydrogen-rich interacting SNe. On the other hand, the flash features observed in many Type IIP SNe, such as SN 2013fs, may indicate that RSG atmospheres are more extended than currently assumed; this could enhance the parameter space for wide binary interaction. [ABSTRACT FROM AUTHOR]

This contribution presents new results on two members of the class of post-Red Supergiants, IRAS 17163-3907, the central star of the Fried Egg nebula and IRC +10420. New optical spectra in the blue spectral range confirm their spectral type to be of A-supergiant class. Our VLTI/GRAVITY K -band interferometry reveals that the neutral Na i 2.2 μm line emitting region is smaller than that of the hydrogen Br γ emission. This can be explained with the hydrogen emission the result of collisional excitation populating the higher levels in a neutral region instead them being populated through recombination in an ionised environment as mostly inferred in stellar winds. Finally, the central star of the Fried Egg nebula, has undergone 3 distinct mass loss episodes over the last hundreds of years. As it is likely that at least the last mass loss event occurred when the star was already a Yellow Hypergiant and not a Red Supergiant, we put forward the bi-stability mechanism as explanation for the mass loss. [ABSTRACT FROM AUTHOR]

Supernova properties in radio strongly depend on their circumstellar environment and they are an important probe to investigate the mass loss of supernova progenitors. Recently, core-collapse supernova observations in radio have been assembled and the rise time and peak luminosity distribution of core-collapse supernovae in radio has been obtained. In this talk, we will discuss the constraints on the mass-loss prescriptions of red supergiants obtained from the assembled radio properties of Type II supernovae. We take a couple of mass-loss prescriptions for red supergiants, calculate the rise time and peak luminosity distribution based on them, and compare the results with the observed distribution. We found that the widely spread radio rise time and peak luminosity distribution of Type II supernovae can only be explained by mass-loss prescriptions having strong dependence on the luminosity. Red supergiant mass-loss prescriptions should have steep luminosity dependence in the supernova progenitor range. [ABSTRACT FROM AUTHOR]

We fit the high-resolution Chandra X-ray spectra of the O supergiant ζ Puppis using the variable boundary condition (VBC) line model to test the stability of its mass-loss rate between two epochs of observation: 2000 March and 2018 July – 2019 August. At issue is whether the observed variations are induced by global changes in the cool (unshocked) wind itself or are isolated to the local pockets of hot gas (i.e. changes in the frequency and location of the shocks). Evidence in the literature favoured the possibility of a 40 per cent increase in the mass flux of the entire stellar wind, based on X-ray reabsorption from a line-deshadowing-instability-inspired parametrization, whereas our fit parameters are consistent with a constant mass flux with a change in the velocity variations that determine the locations where shocks form. Our results suggest the shocks in the more recent data are formed at somewhat larger radii, mimicking the enhanced blueshifts and increased line fluxes interpreted in the previous analysis as being due to increases in both the X-ray generation and reabsorption from an overall stronger wind. [ABSTRACT FROM AUTHOR]

We conducted four-night multiwavelength observations of an active M-dwarf star EV Lacertae on 2022 October 24–27 with simultaneous coverage of soft X-rays (NICER; 0.2–12 keV, Swift XRT; 0.2–10 keV), near-ultraviolet (Swift UVOT/UVW2; 1600–3500 Å), optical photometry (TESS; 6000–10000 Å), and optical spectroscopy (Nayuta/MALLS; 6350–6800 Å). During the campaign, we detected a flare starting at 12:28 UTC on October 25 with a white-light bolometric energy of 3.4 × 1032 erg. At about 1 h after this flare peak, our Hα spectrum showed a blueshifted excess component at a corresponding velocity of ∼100 km s−1. This may indicate that the prominence erupted with a 1 h delay of the flare peak. Furthermore, the simultaneous 20 s cadence near-ultraviolet (NUV) and white-light curves show gradual and rapid brightening behaviors during the rising phase at this flare. The ratio of flux in NUV to white light at the gradual brightening was ∼0.49, which may suggest that the temperature of the blackbody is low (<9000 K) or the maximum energy flux of a non-thermal electron beam is less than 5 × 1011 erg cm−2 s−1. Our simultaneous observations of the NUV and white-light flare raise the issue of a simple estimation of UV flux from optical continuum data by using a blackbody model. [ABSTRACT FROM AUTHOR]

Just 10 recurrent novae (RNe) – which erupt repeatedly on time-scales shorter than one century – are known in our Galaxy. The most extreme RN known (located in the Andromeda galaxy), M31N 2008-12a, undergoes a nova eruption every year, and is surrounded by a vast nova 'super-remnant', 134 pc in extent. Simulations predict that all RNe should be surrounded by similar vast shells, but previous searches have failed to detect them. KT Eri has recently been suggested to be a RN, and we have used the Condor Array Telescope to image its environs through multiple narrow-band filters. We report the existence of a large (∼50-pc diameter), H  |$\, \alpha$| -bright shell centred on KT Eri, exactly as predicted. This strongly supports the claim that KT Eri is the 11th Galactic recurrent nova, and only the second nova known to be surrounded by a super-remnant. SALT spectra of the super-remnant demonstrate that its velocity width is consistent with that of M31-2008-12a. [ABSTRACT FROM AUTHOR]

We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a |$2\, \mathrm{M_\odot }$| red giant branch primary and a |$1\, \mathrm{M_\odot }$| companion modelled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, ∼15 per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of |${\sim}2\, {\rm yr}$|⁠ , after which the envelope would be unbound. [ABSTRACT FROM AUTHOR]

Planetary nebulae (PNe) have long been known to consist of gas, molecules, and dust. The vast majority of the work in which properties of the dust of these objects are inferred were carried out through infrared observations. We use a technique that gives information about the dust, without nebular contamination, and regardless of its temperature. The primary objective is to present an atlas of images with the distribution of dust extinction in PNe. We obtained images in narrow-band filter centred on the nebular emission of the Balmer lines H α and H β to derive extinction images. By combining these images we have been able to achieve the distribution of dust in PNe. The atlas contains 29 PNe, representing the most comprehensive and extensive surveys of dust distribution in PNe. [ABSTRACT FROM AUTHOR]

Accurate determination of mass-loss rates from massive stars is important to understand stellar and galactic evolution and enrichment of the interstellar medium. Large-scale structure and variability in stellar winds have significant effects on mass-loss rates. Time-series observations provide direct quantification of such variability. Observations of this nature are available for some Galactic early supergiant stars but not yet for stars in lower metallicity environments such as the Magellanic Clouds. We utilize ultraviolet spectra from the Hubble Space Telescope ULLYSES program to demonstrate that the presence of structure in stellar winds of supergiant stars at low metallicities may be discerned from single-epoch spectra. We find evidence that, for given stellar luminosities and mean stellar wind optical depths, structure is more prevalent at higher metallicities. We confirm, at Large Magellanic Cloud (0.5 Z ⊙), Small Magellanic Cloud (0.2 Z ⊙), and lower (0.14–0.1 Z ⊙) metallicities, earlier Galactic results that there does not appear to be correlation between the degree of structure in stellar winds of massive stars and stellar effective temperature. Similar lack of correlation is found with regard to terminal velocity of stellar winds. Additional and revised values for radial velocities of stars and terminal velocities of stellar winds are presented. Direct evidence of temporal variability, on time-scales of several days, in stellar wind at low metallicity is found. We illustrate that narrow absorption components in wind-formed profiles of Galactic OB stellar spectra remain common in early B supergiant spectra at low metallicities, providing means for better constraining hot, massive star mass-loss rates. [ABSTRACT FROM AUTHOR]

Hot subdwarf B (sdB) stars are helium core burning stars that have lost almost their entire hydrogen envelope due to binary interaction. Their assumed canonical mass of |${M_{\mathrm{sdB}}}\sim 0.47$|  M⊙ has recently been debated given a broad range found both from observations as well as from the simulations. Here, we revise and refine the mass range for sdBs derived two decades ago with the Eggleton code, using the stellar evolution code mesa , and discuss the effects of metallicity and the inclusion of core overshooting during the main sequence. We find an excellent agreement for low-mass progenitors, up to ∼2.0 M⊙. For stars more massive than ∼2.5 M⊙, we obtain a wider range of sdB masses compared to the simulations from the literature. Our mesa  models for the lower metallicity predict, on average, slightly more massive sdBs. Finally, we show the results for the sdB lifetime as a function of sdB mass and discuss the effect this might have in the comparison between simulations and observational samples. This study paves the way for reproducing the observed Galactic mass distribution of sdB binaries. [ABSTRACT FROM AUTHOR]

Asteroseismology provides a new avenue for accurately measuring the masses of evolved globular cluster (GC) stars. We present the first detections of solar-like oscillations in 47 red giant branch (RGB) and early asymptotic giant branch (EAGB) stars in the metal-poor GC M80; only the second with measured seismic masses. We investigate two areas of stellar evolution and GC science: multiple populations and stellar mass-loss. We detect a distinct bimodality in the EAGB mass distribution. We suggest that this could be due to sub-population membership. If confirmed in future work with spectroscopy, it would be the first direct measurement of a mass difference between sub-populations. A mass difference was not detected between the sub-populations in our RGB sample. We instead measured an average RGB mass of |$0.782\pm 0.009~\mathrm{M}_{\odot }$|⁠ , which we interpret as the average of the sub-populations. Differing mass-loss rates on the RGB have been proposed as the second parameter that could explain the horizontal branch morphology variations between GCs. We calculated an integrated RGB mass-loss separately for each sub-population: |$0.12\pm 0.02~\mathrm{M}_{\odot }$| (SP1) and |$0.25\pm 0.02~\mathrm{M}_{\odot }$| (SP2). Thus, SP2 stars appear to have enhanced mass-loss on the RGB. Mass-loss is thought to scale with metallicity, which we confirm by comparing our results to a higher metallicity GC, M4. Finally, our study shows the robustness of the Δ ν -independent mass scaling relation in the low-metallicity (and low surface gravity) regime. [ABSTRACT FROM AUTHOR]

Asymptotic giant branch (AGB) stars experience heavy episodes of mass-loss through a slow stellar wind during the thermal pulse phase that form large, pc-scale structures around them. As the AGB stellar wind interacts with the interstellar medium (ISM), the otherwise isotropic ejecta gets distorted, resulting in asymmetric shapes, bow-shock structures and, in the case of fast motion relative to the surrounding ISM, extended wakes and tails as unexpectedly detected in GALEX ultraviolet (UV) images of o  Cet, also known as Mira. Since that discovery, another fourteen AGB stars have been reported to exhibit extended UV emission around them. Here, we present the discovery of extended UV emission around another three AGB stars, namely R For, R Hor, and DM Tuc. The analysis of the overall properties of the regions of extended UV emission indicates that these are preferentially detected in closer AGB stars far away from the Galactic plane, i.e. sources less extincted. Faster AGB stars tend to have more elongated shapes, with bow-shocks and wakes or tails of material, while AGB stars farther away from the Galactic plane tend to be surrounded by larger regions of UV emission. [ABSTRACT FROM AUTHOR]

Context. V1309 Sco is an example of a red nova, a product of the merger between non-compact stars. V1309 Sco is particularly important within the class of red novae due to the abundance of photometric data of the progenitor binary before the merger. Aims. We aim to investigate the spatio-kinematic and chemical properties of the circumstellar environment, including deriving the physical conditions and establishing the origins of the different circumstellar components. Methods. We used radiative transfer modelling of molecular emission in sub-millimetre spectra to examine the properties of the molecular gas, and we used forbidden line diagnostics from optical spectra to constrain the electron density and temperature using forbidden line diagnostics. We compared line intensities from shock models to observations to look for and constrain shocks. Results. We derived a new kinematical distance of 5.6 kpc to the source. The detection of ro-vibrational H2 and sub-millimetre HCO+ emission in 2016 and 2019, respectively, indicate active shock interactions within the circumstellar environment. The velocity profiles of both H2 and HCO+, as well as the moment-1 maps of sub-millimetre CO and 29SiO, indicate a bipolar structure that may be asymmetric. The sub-millimetre and optical molecular emission exhibits temperatures of 35–113 and 200 K, respectively, whilst the atomic gas is much hotter, with temperatures of 5–15 kK, which may be due to shock heating. Conclusions. The detection of a bipolar structure in V1309 Sco indicates further similarities with the structure of another Galactic red nova, V4332 Sgr. It provides evidence that bipolar structures may be common in red novae. All collected data are consistent with V1309 Sco being a kinematically and chemically complex system. [ABSTRACT FROM AUTHOR]

Context. U Equ is an unusual maser-hosting infrared source discovered in the 1990s. It was tentatively classified as a post-asymptotic giant branch (post-AGB) star with a unique optical spectrum displaying rare emission and absorption features from molecular gas at a temperature of about 500 K. In 2022, we serendipitously discovered that its optical spectrum had drastically changed since the last observations in the 1990s. Aims. We aim to characterize the drastic change in the spectrum and analyze the photometric behavior of the object since 1989. Methods. Optical high-resolution spectra of U Equ from the Southern African Large Telescope were supplemented by archival data and near-infrared photometry from the Nordic Optical Telescope. New spectral line observations with the Effelsberg 100 m radio telescope and Atacam Large Millimeter Array are presented. Radiative transfer modeling of multiple epoch spectral energy distributions was performed. Results. No circumstellar molecular features are present in the contemporary optical spectra of U Equ. Nonphotospheric absorption and emission from neutral and ionized species dominate the current spectrum. Some of the observed features indicate an outflow with a projected terminal velocity of 215 km s−1. Broad H&K lines of [Ca II] indicate a photosphere of spectral type F or similar. For the first time, we find SiO J = 1−0 υ=1 maser emission in U Equ. Our collected photometric measurements show that the source has monotonically increased its optical and near-infrared fluxes since about the beginning of this century and continues to do so. The current rise in the optical regime is about 1 mag. Spectral energy distributions at different epochs show dusty circumstellar material that is very likely arranged in a highly inclined disk. Adopting a distance of 4 kpc, informed by the Gaia parallax of U Equ, we find that the source luminosity is about 104L⊙. This luminosity has likely increased by a factor of a few in the last decades, which is most probably related to the drastic change in the optical circumstellar spectrum of the object. Conclusions. The object has changed considerably in the past three decades, either due to geometrical reconfiguration of the circumstellar medium, evolutionary changes in the central star, or owing to an accretion event that started very recently in the system. Observationally, U Equ appears to resemble category 0 of disk-hosting post-AGB stars reported previously, especially the post-common-envelope binary HD 101584. It is uncertain whether the drastic spectral change and the associated optical and mid-infrared rise in brightness witnessed in U Equ are common in post-AGB stars, but this radical change may be related to the real-time onset of the evolution of the system into a planetary nebula. We find that the post-AGB star V576 Car has undergone a similar transformation as U Equ in the past few decades, which means that the phenomenon is not extremely rare. [ABSTRACT FROM AUTHOR]

The orbital stability of contact binary systems has been receiving considerable attention recently. Theoretical studies indicate that merger is likely to occur at very low-mass ratios, but the actual mass ratio at which merger may take place is likely to be variable and dependent on the mass of the primary. We consider the effects of metal content on the orbital stability of contact binary systems by modelling the gyration radius of a rotating and tidally distorted primary component at various values of |$\rm [Fe/H]$| in the range −1.25 to +0.5. We determine the instability mass ratio range for contact binary systems with a low-mass primary 0.6 M⊙ ≤ M 1 ≤ 1.4 M⊙ at various metallicity levels and show that systems with low metallicity have an instability mass ratio lower than those with higher metal content and therefore are likely to be more stable. We illustrate the effect through light-curve analysis of two otherwise very similar contact binary systems, except for different metallicity. While both would be considered unstable if metallicity was not taken into consideration, only one remains in that category after appropriate adjustments based on metallicity have been made. [ABSTRACT FROM AUTHOR]

The theory of line-driven winds can explain many observed spectral features in early-type stars, though our understanding the winds of B supergiants remains incomplete. The hydrodynamic equations for slowly rotating stellar winds predict two regimes based on the line-force parameter δ: the fast and the δ- slow solution. In this paper, we aim to explore the capability of the latter to explain the observed properties of B supergiant winds. We calculate Hα line profiles, the most sensitive wind diagnostics in the optical, for both fast and δ- slow wind models. We fit them to observed data from a well-studied sample of B supergiants, by adapting the line-force parameters (k , α, and δ) of the hydrodynamic model. Unexpectedly, the observed Hα spectra can be reproduced by both hydrodynamic wind regimes with similar precision. We argue that this similarity results from the similar shape of the normalized velocity law produced by both regimes in the lower, Hα-forming wind region. Our findings raise a dichotomy, because mass-loss rates and terminal velocities (v ∞) for each solution are quite different. The δ- slow solution predicts maximum values for v ∞ that are systematically lower than those measured in the ultraviolet, whereas the v ∞ values of the fast solution are closer, and probably more appropriate. However, our results also indicate that the δ- slow solution might better describe the dense winds of B hypergiants. Multiwavelength analyses and a larger sample of stars are needed to reach a definitive conclusion. [ABSTRACT FROM AUTHOR]

η Carinae is an extremely luminous and energetic colliding-wind binary. The combination of its orbit and orientation, with respect to our line of sight, enables direct investigation of the conditions and geometry of the colliding winds. We analyse optical He i 5876 and 7065 Å line profiles from the Global Jet Watch observatories covering the last 1.3 orbital periods. The sustained coverage throughout apastron reveals the distinct dynamics of the emitting versus absorbing components: the emission lines follow orbital velocities, while one of the absorption lines is detected only around apastron (0.08 < ϕ < 0.95) and exhibits velocities that deviate substantially from the orbital motion. To interpret these deviations, we conjecture that this He i absorption component is formed in the post-shock primary wind, and is only detected when our line of sight intersects with the shock cone formed by the collision of the two winds. We formulate a geometrical model for the colliding winds in terms of a hyperboloid in which the opening angle and location of its apex are parametrized in terms of the ratio of the wind momentum of the primary star to that of companion. We fit this geometrical model to the absorption velocities, finding results that are concordant with the panchromatic observations and simulations of η Carinae. The model presented here is an extremely sensitive probe of the exact geometry of the wind momentum balance of binary stars, and can be extended to probe the latitudinal dependence of mass-loss. [ABSTRACT FROM AUTHOR]

We present the analysis of spatially resolved spectroscopic observations of the born-again planetary nebula (PN) A 78 that are used to investigate the chemistry and physical properties of its three main morphological components, namely the inner knots, its eye-like structure, and the low surface-brightness outer nebula. The H-poor chemical abundances of the inner knots confirm the born-again nature of A 78, with an N/O abundances ratio consistent with the predictions of very late thermal pulses (VLTP). On the other hand, the high Ne/O is not expected in VLTP events, which prompts a possible different evolutionary path may be involving a binary system. The low N/O ratio and He/H abundances of the outer shell are indicative of a low-mass progenitor, whereas the chemical abundances of the eye-like structure, which results from the interaction between the H-poor born-again material and the outer nebula, evidence their mixture. Unlike previous works, the extinction is found to be inhomogeneous, being much higher towards the H-poor inner knots, where the presence of large amounts of C-rich dust has been previously reported. Dust-rich material seems to diffuse into outer nebular regions, resulting in zones of enhanced extinction. [ABSTRACT FROM AUTHOR]

Massive star winds are structured both stochastically ('clumps') and often coherently (Co-rotation Interaction Regions, or CIRs). Evidence for CIRs threading the winds of Wolf–Rayet (WR) stars arises from multiple diagnostics including linear polarimetry. Some observations indicate changes in polarization position angle across optical recombination emission lines from a WR star wind but limited to blueshifted Doppler velocities. We explore a model involving a spherical wind with a single conical CIR stemming from a rotating star as qualitative proof-of-concept. To obtain a realistic distribution of limb polarization and limb darkening across the pseudo-photosphere formed in the optically thick wind of a WR star, we used Monte Carlo radiative transfer (MCRT). Results are shown for a parameter study. For line properties similar to WR 6 (EZ CMa; HD 50896), combining the MCRT results, a simple model for the CIR, and the Sobolev approximation for the line formation, we were able to reproduce variations in both polarization amplitude and position angle commensurate with observations. Characterizing CIRs in WR winds has added importance for providing stellar rotation periods since the v sin i values are unobtainable because the pseudo-photosphere forms in the wind itself. [ABSTRACT FROM AUTHOR]

The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMSs, M>100  |$\rm {\rm M}_{\odot }$|⁠) are relatively rare compared to O stars, they lose disproportionately large amounts of mass already from the onset of core H-burning. VMS have optically thick winds with elevated mass-loss rates in comparison to optically thin standard O-star winds. We compute wind yields and ejected masses on the main sequence, and we compare enhanced mass-loss rates to standard ones. We calculate solar metallicity wind yields from MESA stellar evolution models in the range 50–500  |$\rm {\rm M}_{\odot }$|⁠ , including a large nuclear network of 92 isotopes, investigating not only the CNO-cycle, but also the Ne–Na and Mg–Al cycles. VMS with enhanced winds eject 5–10 times more H-processed elements (N, Ne, Na, Al) on the main sequence in comparison to standard winds, with possible consequences for observed anticorrelations, such as C–N and Na–O, in globular clusters. We find that for VMS 95 per cent of the total wind yields is produced on the main sequence, while only ∼ 5 per cent is supplied by the post-main sequence. This implies that VMS with enhanced winds are the primary source of 26Al, contrasting previous works where classical Wolf–Rayet winds had been suggested to be responsible for galactic 26Al enrichment. Finally, 200  |$\rm {\rm M}_{\odot }$| stars eject 100 times more of each heavy element in their winds than 50  |$\rm {\rm M}_{\odot }$| stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50  |$\rm {\rm M}_{\odot }$| stars. [ABSTRACT FROM AUTHOR]

Neutron stars accreting from OB supergiants are often divided between persistently and transiently accreting systems, called supergiant X-ray binaries (SgXBs) and supergiant fast X-ray transients (SFXTs). This dichotomy in accretion behaviour is typically attributed to systematic differences in the massive stellar wind, binary orbit, or magnetic field configuration, but direct observational evidence for these hypotheses remains sparse. To investigate their stellar winds, we present the results of pilot 100-GHz observations of one SFXT and one SgXB with the Northern Extended Millimetre Array. The SFXT, IGR J18410-0535, is detected as a point source at 63.4 ± 9.6 μJy, while the SgXB, IGR J18410-0535 remains undetected. Radio observations of IGR J18410-0535 imply a flat or inverted low-frequency spectrum, arguing for wind emission and against non-thermal flaring. Due to the uncertain SFXT distance, however, the observations do not necessarily imply a difference between the wind properties of the SFXT and SgXB. We compare the mm constraints with other HMXBs and isolated OB supergiants, before considering how future mm campaigns can constrain HMXB wind properties by including X-ray measurements. Specifically, we discuss caveats and future steps to successfully measure wind mass-loss rates and velocities in HMXBs with coordinated mm, radio, and X-ray campaigns. [ABSTRACT FROM AUTHOR]

Some interaction-powered supernovae have long rise times of more than 100 d. We show that such long rise times are naturally expected if cirumstellar matters (CSM) have a flat density structure (s ≲ 1.5, where ρCSM∝ r − s ). In such cases, bolometric luminosities from the CSM interaction keep increasing as long as the CSM interacts with the outer layers of the SN ejecta. Thus, the rise time is determined by the dynamical time-scale in which the reverse shock propagates the outer layers of the SN ejecta, not by the time-scales in which photons diffuse in the CSM as often considered. Interaction-powered supernovae with very long rise times can be an important probe of extensive non-steady mass-loss in massive stars. [ABSTRACT FROM AUTHOR]

As a cool star evolves, it loses mass and angular momentum due to magnetized stellar winds that affect its rotational evolution. This change has consequences that range from the alteration of its activity to influences over the atmosphere of any orbiting planet. Despite their importance, observations constraining the properties of stellar winds in cool stars are extremely limited. Therefore, numerical simulations provide a valuable way to understand the structure and properties of these winds. In this work, we simulate the magnetized winds of 21 cool main-sequence stars (F-type to M-dwarfs) using a state-of-the-art 3D MHD code driven by observed large-scale magnetic field distributions. We perform a qualitative and quantitative characterization of our solutions, analysing the dependencies between the driving conditions (e.g. spectral type, rotation, and magnetic field strength) and the resulting stellar wind parameters (e.g. Alfvén surface size, mass-loss rate, angular momentum loss rate, and stellar wind speeds). We compare our models with the current observational knowledge on stellar winds in cool stars and explore the behaviour of the mass-loss rate as a function of the Rossby number. Furthermore, our 3D models encompass the entire classical Habitable Zones (HZ) of all the stars in our sample. This allows us to provide the stellar wind dynamic pressure at both edges of the HZ and analyse the variations of this parameter across spectral type and orbital inclination. The results here presented could serve to inform future studies of stellar wind-magnetosphere interactions and stellar wind erosion of planetary atmospheres via ion escape processes. [ABSTRACT FROM AUTHOR]

Mass-transfer interactions in binary stars can lead to accretion disc formation, mass-loss from the system, and spin-up of the accretor. To determine the trajectory of the mass-transfer stream, and whether it directly impacts the accretor, or forms an accretion disc, requires numerical simulations. The mass-transfer stream is approximately ballistic, and analytic approximations based on such trajectories are used in many binary population synthesis codes as well as in detailed stellar evolution codes. We use binary population synthesis to explore the conditions under which mass transfer takes place. We then solve the reduced three-body equations to compute the trajectory of a particle in the stream for systems with varying system mass ratio, donor synchronicity, and initial stream velocity. Our results show that, on average, both more mass and more time are spent during mass transfer from a sub-synchronous donor than from a synchronous donor. Moreover, we find that at low initial stream velocity, the asynchronous rotation of the donor leads to self-accretion over a large range of mass ratios, especially for supersynchronous donors. The stream (self-)intersects in a narrow region of parameter space where it transitions between accreting on to the donor or the accretor. Increasing the initial stream velocity leads to larger areas of the parameter space where the stream accretes on to the accretor, but also more (self-)intersection. The radii of closest approach generally increase, but the range of specific angular momenta that these trajectories carry at the radius of closest approach gets broader. Our results are made publicly available. [ABSTRACT FROM AUTHOR]

Very massive stars (VMS) up to 200–300 M⊙ have been found in the Local Universe. If they would lose little mass, they produce intermediate-mass black holes or pair-instability supernovae (PISNe). Until now, VMS modellers have extrapolated mass-loss versus metallicity (Z) exponents from optically thin winds, resulting in a range of PISN thresholds that might be unrealistically high in Z , as VMS develop optically thick winds. We utilize the transition mass-loss rate of Vink and Gräfener (2012) that accurately predicts mass-loss rates of Of/WNh ('slash') stars that characterize the morphological transition from absorption-dominated O-type spectra to emission-dominated WNh spectra. We develop a wind efficiency framework, where optically thin winds transition to enhanced winds, enabling us to study VMS evolution at high redshift where individual stars cannot be resolved. We present a MESA grid covering Z ⊙/2 to Z ⊙/100. VMS above the transition evolve towards lower luminosity, skipping the cool supergiant phase but directly forming pure He stars at the end of hydrogen burning. Below the transition, VMS evolve as cooler luminous blue variables (LBVs) or yellow hypergiants (YHGs), naturally approaching the Eddington limit. Strong winds in this YHG/LBV regime – combined with a degeneracy in luminosity – result in a mass-loss runaway, where a decrease in mass increases wind mass loss. Our models indicate an order-of-magnitude lower metallicity threshold for PISN than usually assumed, at Z ⊙/20 due to our mass-loss runaway. While future work on LBV mass loss could affect the PISN threshold, our framework will be critical for establishing definitive answers on the PISN threshold and galactic chemical evolution modelling. [ABSTRACT FROM AUTHOR]

Calculation of the mass transfer (MT) rate |$\dot{M}_\text{d}$| of a Roche lobe overflowing star is a fundamental task in binary star evolution theory. Most of the existing MT prescriptions are based on a common set of assumptions that combine optically thick and optically thin regimes with different flow geometries. In this work, we develop a new model of MT based on the assumption that the Roche potential sets up a nozzle converging on the inner Lagrangian point and that the gas flows mostly along the axis connecting both stars. We derive a set of 1D hydrodynamic equations governing the gas flow with |$\dot{M}_\text{d}$| determined as the eigenvalue of the system. The inner boundary condition directly relates our model to the structure of the donor obtained from 1D stellar evolution codes. We obtain an algebraic solution for the polytropic equation-of-state (EOS). This gives |$\dot{M}_\text{d}$| within a factor of 0.9 to 1.0 of existing optically thick prescriptions and reduces to the existing optically thin prescription for isothermal gas. For a realistic EOS, we find that |$\dot{M}_\text{d}$| differs by up to a factor of 4 from existing models. We illustrate the effects of our new MT model on a |$30\, M_\odot$| low-metallicity star undergoing intensive thermal time-scale MT and find that it is more likely to become unstable to L2 overflow and common-envelope evolution than expected according to MT prescriptions. Our model provides a framework to include additional physics such as radiation or magnetic fields. [ABSTRACT FROM AUTHOR]