Search

In this paper, we studied the X-ray polarization of the jet in the Vela Pulsar Wind Nebula. After the background subtraction, we found a high degree of polarization in the helical jet region, close to the upper limit allowed by the synchrotron radiation mechanism. This result indicates the magnetic field is highly ordered, even though the jet's morphology is curved and varies rapidly, which suggests that the magnetic collimation process for the jet does not introduce significant turbulence. The polarization angle is perpendicular to the symmetry axis of the nebula, consistent with the results from the magnetohydrodynamic numerical simulations. This finding sheds light on the underlying mechanisms related to the dynamics and evolution of the jet structures. [ABSTRACT FROM AUTHOR]

Context. Cygnus X-1 (Cyg X-1) is a high-mass X-ray binary where accretion onto the black hole (BH) is mediated by the stellar wind from the blue supergiant companion star HDE 226868. Due to its inclination, the system is a perfect laboratory to study the not yet well-understood stellar wind structure. In fact, depending on the position of the BH along the orbit, X-ray observations can probe different layers of the stellar wind. Deeper wind layers can be investigated at superior conjunction (i.e. null orbital phases). Aims. We aim to characterise the stellar wind in the Cyg X-1/HDE 226868 system, analysing one passage at superior conjunction covered by XMM-Newton during the 'Cyg X-1 Hard state Observations of a Complete Binary Orbit in X-rays' (CHOCBOX) campaign. Methods. To analyse the properties of the stellar wind, we computed colour-colour diagrams. Since X-ray absorption is energy-dependent, colour indices provide information on the parameters of the stellar wind, such as the column density, NH, w, and the covering factor, fc. We fitted colour-colour diagrams with models that include both a continuum and a stellar wind component. We used the kernel density estimation method to infer the unknown probability distribution of the data points in the colour-colour diagram, and selected the model corresponding to the highest likelihood. In order to study the temporal evolution of the wind around superior conjunction, we extracted and fitted time-resolved colour-colour diagrams. Results. We found that the model that best describes the shape of the colour-colour diagram of Cyg X-1 at superior conjunction requires the wind to be partially ionised. The shape of the colour-colour diagram strongly varies during the analysed observation, due to concurrent changes of the mean NH, w and the fc of the wind. Our results suggest the existence of a linear scaling between the rapid variability amplitude of NH, w (on timescales between 10 s and 11 ks) and its long-term variations (on timescales > 11 ks). Using the inferred best-fit values, we estimated the stellar mass loss rate to be ∼7 × 10−6 M⊙ yr−1 and the clumps to have a characteristic mass of ∼1017 g. [ABSTRACT FROM AUTHOR]

Recent radio observations and coincident neutrino detections suggest that some tidal disruption events (TDEs) exhibit late-time activities, relative to the optical emission peak, and these may be due to delayed outflows launched from the central supermassive black hole. We investigate the possibility that jets launched with a time delay of days to months, interact with a debris that may expand outwards. We discuss the effects of the time delay and expansion velocity on the outcomes of jet breakout and collimation. We find that a jet with an isotropic-equivalent luminosity of |$\lesssim 5 \times 10^{45}\, {\rm erg\, s}^{-1}$| is likely to be choked for a delay time of |$\sim 3$| months. We also study the observational signatures of such delayed choked jets. The jet–debris interaction preceding the breakout would lead to particle acceleration and the resulting synchrotron emission can be detected by current and near-future radio, optical and X-ray telescopes, and the expanding jet-driven debris could explain late-time radio emission. We discuss high-energy neutrino production in delayed choked jets, and the time delay can significantly alleviate the difficulty of the hidden jet scenario in explaining neutrino coincidences. [ABSTRACT FROM AUTHOR]

The spectral variation in the optical region of a peculiar variable star V Sagittae was monitored in the recent three decades. The main aim of this article is to report several interesting phenomena discovered in the monitoring; we briefly discuss interpretations of them. Emission lines of H  i and He  ii were decomposed into narrow components and wide wings. The radial velocities of the wings were related to the orbital motion of the hot primary, but did not reflect it directly. Widely red-shifted sinusoidal variations according to the orbital period were detected in 1993–2005 (e.g. |$K = 195$|  km s |$^{-1}$| and |$\gamma = +252$|  km s |$^{-1}$| for He  ii 4686). Sporadic blue-shifts were also observed. The mode of the radial velocity variation drastically changed between 2005 and 2007. Red-shifted sinusoidal variations (⁠|$K = 224$|  km s |$^{-1}$| and |$\gamma = +260$|  km s |$^{-1}$| idem) or blue-shifted ones (⁠|$K =206$|  km s |$^{-1}$| and |$\gamma = -289$|  km s |$^{-1}$| idem) appeared alternately in 2007–2019. One red- or blue-shift state lasted several years. The radial velocities of O  iii 5592 and the doublet C  iv 5805 exhibited sinusoidal variations without large red- or blue-shifts (e.g. |$K =215$|  km s |$^{-1}$| and |$\gamma = -18$|  km s |$^{-1}$| for C  iv 5805). The widths of emission lines He  ii 5412 and C  iv in the doublet 5805 suggest that helium and carbon ions contained nearly the same kinetic energy, because the former line was wider than the latter by a factor of about |$\surd {3}$| on most spectra. It is probable that there was no large-scale mass motion in such a thermal-equilibrium-like condition. The C  iv lines and He  ii 5412 had roughly same widths, like the spectra of Wolf–Rayet stars, only three times in 2007–2019 on which there were likely large-scale mass outflows. [ABSTRACT FROM AUTHOR]

Context. The details of the mass-loss process in the late stages of low- and intermediate-mass stellar evolution are not well understood, in particular its dependence on stellar parameters. Mira variables are highly suitable targets for studying this mass-loss process. Aims. Here, we follow up on our earlier finding that a near-to-mid-infrared (NIR-MIR) colour versus pulsation period diagram shows two sequences of Miras that can be distinguished by the third dredge-up (3DUP) indicator technetium in those stars. While IR colours are good indicators of the dust mass-loss rate (MLR) from Miras, no corresponding sequences have been found using the gas MLR. However, investigations of the gas MLR have been hampered by data limitations. We aim to alleviate these limitations with new observational data. Methods. We present new optical spectra of a well-selected sample of Miras. We searched these spectra for absorption lines of Tc and other 3DUP indicators, and combine our findings with gas MLRs and expansion velocities from the literature. Furthermore, we extend the analysis of the MIR emission to WISE data and compare the broadband spectral energy distributions (SEDs) of Miras with and without Tc. Results. We find no systematic difference in gas MLRs between Miras with and without Tc. However, the gas envelopes of Tc-poor Miras appear to have a higher terminal expansion velocity than those of Miras with Tc. Furthermore, our analysis of the IR photometry strongly corroborates the earlier finding that Tc-poor Miras have a higher MIR emission than Tc-rich ones, by as much as a factor of two. We model the IR colours with DARWIN and stationary wind models and conclude that Miras with and without Tc have different dust content or dust properties. Conclusions. We discuss several hypotheses and interpretations of the observations and conclude that the reduction of free oxygen by 3DUP of carbon and iron-depleted dust grains in Tc-rich stars are the most convincing explanations for our observations. [ABSTRACT FROM AUTHOR]

We present the serendipitous discovery of a new radio-continuum ring-like object nicknamed Kýklos (J1802–3353), with MeerKAT UHF and L-band observations. The radio ring, which resembles the recently discovered odd radio circles (ORCs), has a diameter of ∼80″ and is located just ∼6° from the Galactic plane. However, Kýklos exhibits an atypical thermal radio-continuum spectrum (α = −0.1 ± 0.3), which led us to explore different possible formation scenarios. We concluded that a circumstellar shell around an evolved massive star, possibly a Wolf-Rayet, is the most convincing explanation with the present data. [ABSTRACT FROM AUTHOR]

Context. The detected population of "spiders", referring to millisecond pulsar binaries, has significantly grown in the past decade thanks to multiwavelength follow-up investigations of unidentified Fermi sources. These systems consist of low-mass stellar companions orbiting rotation-powered millisecond pulsars in short periods of a few hours up to day. Among them, a subset of intriguing objects called transitional millisecond pulsars (tMSPs) has been shown to exhibit a remarkable behavior, transitioning between pulsar-binary and faint low-mass X-ray binary states over a span of a few years. Aims. Our objective is to study the interaction of stellar winds in tMSPs in order to understand their observational properties. To this end we focus on the parameter range that places the system near Roche-lobe overflow. Methods. Employing the adaptative mesh refinement (AMR) AMRVAC 2.0 code, we performed 2D hydrodynamical (HD) simulations of the interaction between the flows from both stars, accounting for the effects of gravity and orbital motion. Results. By studying the mass loss and launch speed of the winds, we successfully recreated two phenomenologically distinct regimes: the accretion stream and the pulsar radio state. We also identified the tipping point that marks the sharp transition between these two states. In the accretion stream state, we discover a very strong variability induced by the pulsar wind. In the pulsar state, we reconstructed the corresponding X-ray light curves of the system that produces the characteristic double-peak pattern of these systems. The position of the peaks is shifted due to orbital motion and the leading peak is weaker due to eclipsing by the companion. Conclusions. This study highlights the importance of gravity and orbital motion in the interaction between the companion and pulsar winds. Our setup allows the study of the complex interaction between the pulsar wind and an accretion stream during mass transfer. We suggest that a smaller leading peak in X-rays is indicative of a nearly edge-on system. [ABSTRACT FROM AUTHOR]

The young stellar object [BHB2007]-1 has been extensively studied in the past at radio, millimeter, and infrared wavelengths. It has revealed a gap in the disk and previous observations have claimed possible emission from a sub-stellar object undergoing formation, in correspondence to the disk gap. In this work, we analyzed a set of eight Karl G. Jansky Very Large Array (VLA) observations at 15 GHz and spread out over a month. We inferred a slowly variable emission from the star, with a ∼15 − 20% circular polarization detected in two of the eight observations. The latter can be related to the magnetic fields in the system, while the unpolarized and moderately varying component can be indicative of free–free emission associated with jet induced shocks or interactions of the stellar wind, with dense surrounding material. We discarded any relevant short-flaring activities when sampling the radio light curves down to 10 seconds and found no clear evidence of emission from the sub-stellar object inferred from past observations, although deeper observations could shed further light on this. [ABSTRACT FROM AUTHOR]

Context. The hot inner flow in black-hole X-ray binaries is not just a static corona rotating around the black hole: it must be partially outflowing. It is therefore a mildly relativistic "outflowing corona". We have developed a model in which Comptonization takes place in this outflowing corona. In all of our previous work, we assumed a rather high outflow speed of 0.8c. Aims. Here, we investigate whether an outflow with a significantly lower speed can also reproduce the observations. Thus, in this work we consider an outflow speed of 0.1c or less. Methods. As in all of our previous work, we used a Monte Carlo code to compute not only the emergent X-ray spectra, but also the time lags that are introduced to the higher-energy photons with respect to the lower-energy ones via multiple scatterings. We also record the angle (with respect to the symmetry axis of the outflow) and the height at which photons escape. Results. Our results are very similar to those of our previous work, with some small quantitative differences that can be easily explained. We are again able to quantitatively reproduce five observed correlations: (a) the time lag as a function of Fourier frequency, (b) the time lag as a function of photon energy, (c) the time lag as a function of Γ, (d) the time lag as a function of the cutoff energy in the spectrum, and (e) the long-standing radio–X-ray correlation – and all of them with only two parameters, which vary in the same ranges for all the correlations. Conclusions. Our model does not require a compact, narrow relativistic jet, although its presence does not affect the results. The essential ingredient of our model is the parabolic shape of the Comptonizing corona. The outflow speed plays a minor role. Furthermore, the bottom of the outflow, in the hard state, looks like a "slab" to the incoming soft photons from the disk, and this can explain the observed X-ray polarization, which is along the outflow. In the hard-intermediate state, we predict that the polarization of GX 339−4 will be perpendicular to the outflow. [ABSTRACT FROM AUTHOR]

Context. The particle-in-cell approach has proven effective in modeling neutron-star and black-hole magnetospheres from first principles, but global simulations are plagued with an unrealistically small separation between the scales where microphysics operates and the system-size scales due to limited numerical resources. A legitimate concern is whether the scale separation achieved to date is large enough for results to be safely extrapolated to realistic scales. Aims. In this work, our aim is to explore the effect of scaling up physical parameters and to check whether salient features uncovered by pure kinetic models at smaller scales are still valid, with a special emphasis on particle acceleration and high-energy radiation emitted beyond the light cylinder. Methods. To reach this objective, we developed a new hybrid numerical scheme coupling the ideal force-free and the particle-in-cell methods to optimize the numerical cost of global models. We propose a domain decomposition of the magnetosphere based on the magnetic-field topology using the flux function. The force-free model is enforced along open field lines while the particle-in-cell model is restricted to the reconnecting field line region. Results. As a proof of concept, this new hybrid model is applied to simulate a weak millisecond pulsar magnetosphere with realistic scales using high-resolution axisymmetric simulations. Magnetospheric features reported by previous kinetic models are recovered, and strong synchrotron radiation above 100MeV consistent with the Fermi-LAT gamma-ray pulsar population is successfully reproduced. Conclusions.This work further consolidates the shining-reconnecting current sheet scenario as the origin of the gamma-ray emission in pulsars, as well as firmly establishing pulsar magnetospheres as at least teraelectronvolt particle accelerators. [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]

We performed numerical simulations of the common envelope (CE) interaction between thermally pulsing asymptotic giant branch (AGB) stars of 1.7 and 3.7 M |$_{\odot }$|⁠ , respectively, and a 0.6 M |$_{\odot }$| compact companion. We use tabulated equations of state to take into account recombination energy. For the first time, formation and growth of dust is calculated explicitly, using a carbon dust nucleation network with a C/O abundance ratio of 2.5 (by number). The first dust grains appear within |$\sim$| 1–3 yr after the onset of the CE, forming an optically thick shell at |$\sim$| 10–20 au, growing in thickness and radius to values of |$\sim$| 400–500 au over |$\sim$| 40 yr, with temperatures around 400 K. Most dust is formed in unbound material, having little effect on mass ejection or orbital evolution. By the end of the simulations, the total dust yield is |$\sim 8.4\times 10^{-3}$| and |$\sim 2.2\times 10^{-2}$| M |$_{\odot }$| for the CE with a 1.7 and a 3.7 M |$_{\odot }$| AGB star, respectively, corresponding to a nucleation efficiency close to 100 per cent, if no dust destruction mechanism is considered. Despite comparable dust yields to single AGB stars, in CE ejections the dust forms a thousand times faster, over tens of years as opposed to tens of thousands of years. This rapid dust formation may account for the shift in the infrared of the spectral energy distribution of some optical transients known as luminous red novae. Simulated dusty CEs support the idea that extreme carbon stars and 'water fountains' may be objects observed after a CE event. [ABSTRACT FROM AUTHOR]

We report on Chandra X-ray observations of SN 2016jae and SN 2018cqj, both low-luminosity Type Ia supernova (SN) that showed the presence of an H line in their early optical spectrum. No X-ray emission is detected at the location of either SN. Upper limits to the luminosity of up to 2 |$\times 10^{40}$|  erg s |$^{-1}$| are calculated for each SN, depending on the assumed spectral model, temperature, and column density. This luminosity is comparable to that of another low-luminosity Type Ia SN, SN 2018fhw, that was observed with Chandra. It is generally lower than upper limits calculated for Type Ia-CSM SNe observed in X-rays, and also below that of SN 2012ca, the only Type Ia-CSM SN to have been detected in X-rays. Comparisons are made to other Type Ia SN with an H line observed in X-rays. The observations suggest that while the density into which the SN is expanding may have been high at the time the H |$\alpha$| line was detected, it had decreased considerably by the time of X-ray observations. [ABSTRACT FROM AUTHOR]

The massive binary HDE 228766 is composed of an O type primary and an evolved secondary. However, previous qualitative analyses of the composite spectrum have led to a wide discussion about whether the secondary is an Of or a Wolf–Rayet star. We use new observations and our novel QER20 package to disentangle for the first time the spectra of the two stellar components and obtain artefact |$-$| free reconstructed spectra, yielding the more accurate and reliable spectral classifications of O7.5 V((f))z for the primary and O6 Iaf for the secondary. The emission features of the P |$-$| Cygni profiles of the H  |$\beta$| and He i 5876 Å lines, present in the reconstructed spectrum of the secondary, show that this star is at an initial phase of its transition to the WN evolutionary stage. A previously unobserved variable emission, composed of at least four independent features, is seen since 2014 superposed to the H  |$\gamma$| absorption line. Our analysis reveals that these emission features originate from a physically extended region. This could be explained by an episode of enhanced mass-loss in the scenario of a non |$-$| conservative evolution of the binary. [ABSTRACT FROM AUTHOR]

We perform a large-scale hydrodynamic simulation of a massive star cluster whose stellar population mimics that of the Cygnus OB2 association. The main-sequence stars are first simulated during 1.6 Myr, until a quasi-stationary state is reached. At this time, the three Wolf–Rayet stars observed in Cygnus OB2 are added to the simulation, which continues to 2 Myr. Using a high-resolution grid in the centre of the domain, we can resolve the most massive stars individually, which allows us to probe the kinetic structures at small (parsec) scales. We find that, although the cluster excavates a spherical 'superbubble' cavity, the stellar population is too loosely distributed to blow a large-scale cluster wind termination shock, and that collective effects from wind–wind interactions are much less efficient than usually assumed. This challenges our understanding of the ultra-high energy emission observed from the region. [ABSTRACT FROM AUTHOR]

Many high-state cataclysmic variables (CVs) exhibit blue-shifted absorption features in their ultraviolet (UV) spectra – a smoking-gun signature of outflows. However, the impact of these outflows on optical spectra remains much more uncertain. During its recent outburst, the eclipsing dwarf nova V455 And displayed strong optical emission lines whose cores were narrower than expected from a Keplerian disc. Here, we explore whether disc + wind models developed for matching UV observations of CVs can also account for these optical spectra. Importantly, V455 And was extremely bright at outburst maximum: the accretion rate implied by fitting the optical continuum with a standard disc model is |$\dot{M}_{\rm acc} \simeq 10^{-7}~{\rm M}_\odot ~{\rm yr^{-1}}$|⁠. Allowing for continuum reprocessing in the outflow helps to relax this constraint. A disc wind can also broadly reproduce the optical emission lines, but only if the wind is (i) highly mass-loaded, with a mass-loss rate reaching |$\dot{M}_{\rm wind} \simeq 0.4 \dot{M}_{\rm acc}$|⁠ , and/or (ii) clumpy, with a volume filling factor |$f_V \simeq 0.1$|⁠. The same models can describe the spectral evolution across the outburst, simply by lowering |$\dot{M}_{\rm acc}$| and |$\dot{M}_{\rm wind}$|⁠. Extending these models to lower inclinations and into the UV produces spectra consistent with those observed in face-on high-state CVs. We also find, for the first time in simulations of this type, P-Cygni-like absorption features in the Balmer series, as have been observed in both CVs and X-ray binaries. Overall, dense disc winds provide a promising framework for explaining multiple observational signatures seen in high-state CVs, but theoretical challenges persist. [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]

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]

Mixed-morphology supernova remnants (MMSNRs) are characterized by a shell-like morphology in the radio and centrally peaked thermal emission in the X-ray band. The nature of this peculiar class of supernova remnants (SNRs) remains a controversial issue. In this work, by pairing the predictions of stellar evolution theory with two-dimensional hydrodynamic simulations, we show that the mixed morphology properties of an SNR can arise by the interaction of the SNR with the circumstellar medium shaped by a red supergiant progenitor star, embedded in a dense environment. As a study case, we model the circumstellar medium formation and the subsequent interaction of the SNR with it of a |$15~\rm {\rm M}_{\odot }$| progenitor star. The reflected shock, formed by the collision of the SNR with the density walls of the surrounding circumstellar cavity, accumulates and re-shocks the supernova ejecta at the centre of the remnant, increasing its temperature so that the gas becomes X-ray bright. Such a formation mechanism may naturally explain the nature of MMSNRs resulted from Type II supernovae without the demand of additional physical mechanisms and/or ambient medium inhomogeneities. We discuss alternative evolutionary paths that potentially could be ascribed for the MMSNR formation within the framework of the reflected shock model. [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]

We investigate outflows and the physics of super‐Eddington versus sub‐Eddington regimes in black hole systems. Our focus is on prospective science using next‐generation high‐resolution soft x‐ray instruments. We highlight the properties of black hole ultraluminous x‐ray source (ULX) systems in particular. Owing to scale invariance in accreting black holes, ULX accretion properties, including their outflows, inform our understanding not only of the closely related population of (similar‐mass) x‐ray binary systems but also of tidal disruption events (TDEs) around supermassive black holes. A subsample of TDEs are likely to transcend super‐Eddington to sub‐Eddington regimes as they evolve, offering an important unifying analog to ULXs and sub‐Eddington x‐ray binaries. We demonstrate how next‐generation soft x‐ray observations with resolving power ≳1000$$ \gtrsim 1000 $$ and collecting area ≳1000cm2$$ \gtrsim 1000\kern0.5em {\mathrm{cm}}^2 $$ can simultaneously identify ultrafast and more typical wind components, distinguish between different wind mechanisms, and constrain changing wind properties over characteristic variability timescales. [ABSTRACT FROM AUTHOR]

Spectroscopic observations have shown for decades that the Wolf–Rayet (WR) phenomenon is ubiquitous among stars with different initial masses. Although much effort to understand the winds from massive WR stars has been presented in the literature, not much has been done for such type of stars in the low-mass range. Here we present an attempt to understand the winds from [WR]-type stars using results from spectral analyses with the full non-local thermodynamic equilibrium stellar atmosphere code PoWR. These results are put into context with the properties of massive WR stars. We found that WC + [WC] stars and WO + [WO] stars create independent sequences in the mass-loss rate (⁠|$\dot{M}$|⁠) and modified wind momentum (D mom) versus luminosity (L) diagrams. Our analysis indicates that even when the winds of WR and [WR] stars become optically thin, there is no breakdown of the general mass-loss trend, contrary to the observed 'weak wind phenomenon' in OB stars. We report that all WR-type stars studied here broadly define single sequences in the wind efficiency (η) versus transformed mass-loss rate (⁠|$\dot{M}_\mathrm{t}$|⁠), the |$\dot{M}_\mathrm{t}$| – T eff diagram, and the |$(L, T_\mathrm{eff}, \dot{M})$| space, which suggest these to be fundamental properties of the WR phenomenon (regardless of the mass range), at least for WR-type stars of the O and C sequences. Our analytical estimations could drive computations of future stellar evolution models for WR-type stars. [ABSTRACT FROM AUTHOR]

Stellar rotation at early ages plays a crucial role in the survival of primordial atmospheres around Earth-mass exoplanets. Earth-like planets orbiting fast-rotating stars may undergo complete photoevaporation within the first few hundred Myr driven by the enhanced stellar XUV [X-rays and extreme ultraviolet (EUV)] radiation, while planets orbiting slow-rotating stars are expected to experience difficulty in losing their primordial envelopes. Besides the action of stellar radiation, stellar winds induce additional erosion on these primordial atmospheres, altering their morphology, extent, and causing supplementary atmospheric losses. In this paper, we study the impact of activity-dependent stellar winds on primordial atmospheres to evaluate the extent to which the action of these winds can be significant in the whole planetary evolution at early evolutionary stages. We performed 3D magnetohydrodynamical (MHD) simulations of the interaction of photoevaporating atmospheres around unmagnetized Earth-mass planets in the time span between 50 and 500 Myr, analysing the joint evolution of stellar winds and atmospheres for both fast- and slow-rotating stars. Our results reveal substantial changes in the evolution of primordial atmospheres when influenced by fast-rotating stars, with a significant reduction in extent at early ages. In contrast, atmospheres embedded in the stellar winds from slow-rotating stars remain largely unaltered. The interaction of the magnetized stellar winds with the ionized upper atmospheres of these planets allows us to evaluate the formation and evolution of different MHD structures, such as double bow shocks and induced magnetospheres. This work will shed light on the first evolutionary stages of Earth-like exoplanets, which are of crucial relevance in terms of planet habitability. [ABSTRACT FROM AUTHOR]

We present a set of controlled hydrodynamical simulations to study the effects of strong galactic outflows on the density and temperature structures, and associated X-ray signatures, of extra-planar and circumgalactic gas. We consider three initial state models, isothermal, isentropic, and rotating cooling-flow, for the hot circumgalactic medium (CGM) into which the outflows are driven. The energy sources are either stellar winds and supernovae, or active galactic nuclei. We consider energy injection rates in the range |$10^{40} \lt \dot{E}_{\rm inj} \lt 10^{44.5}$| erg s−1, and compute the time-dependent soft X-ray (0.5–2 keV) surface brightness. For |$\dot{E}_{\rm inj} \gtrsim 10^{41} {\small --} 10^{42}$| erg s−1, with the exact threshold depending on the initial CGM state, the X-ray response is dominated by dense hot gas in the forward shock that eventually fades into the CGM as a sound wave. The shock surrounds an inner hot bubble leading to a radial flattening of the X-ray surface brightness. For lower energy injection rates, the X-ray surface brightness of the initial CGM state is almost unaffected. We present analytical approximations for the outflow shock propagation and the associated X-ray emissions. [ABSTRACT FROM AUTHOR]

We consider scenarios for non-thermal particle acceleration and re-acceleration in the central cores of compact massive star clusters, aided by insights from high-resolution hydrodynamic simulations. We show that (i) particles are unlikely to interact with many shocks during their lifetimes in the core; (ii) colliding flows do not produce hard spectra; and (iii) turbulent re-acceleration in the core is suppressed. Inefficient re-acceleration mechanisms are not expected to produce hard components nor to increase the maximum energy within the cores of massive star clusters. Models in which the observed ultra-high-energy gamma rays originate in the core of massive stellar clusters are thus disfavoured. [ABSTRACT FROM AUTHOR]

The remnant of the historical supernova 1181 is claimed to be associated with a white dwarf merger remnant J005311. The supernova remnant (SNR) shock, and a termination shock expected to be formed by the intense wind of J005311, are potential sites for radio emission via synchrotron emission from shock-accelerated electrons. In this paper, we estimate the radio emission from these two shocks, and find the peak radio flux to be 0.1–10 mJy (at 0.01–1 GHz) in the outer SNR shock and 0.01–0.1 mJy (at 1–10 GHz) in the inner termination shock. We also search for radio emission from this source in the archival data of the Karl G. Jansky Very Large Array (VLA) Sky Survey at 3 GHz, the NRAO VLA Sky Survey at 1.4 GHz and the Canadian Galactic Plane Survey at 408 MHz, finding no significant detection. While targeted observations with higher sensitivity are desired, we particularly encourage those at higher frequency and angular resolution to probe the inner termination shock and its evolution. [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. Nova 1670 is a historical transient bearing strong similarities to a recently recognized type of stellar eruption known as the red nova, which is thought to be powered by stellar mergers. The remnant of the transient, CK Vul, is observable today mainly through cool circumstellar gas and dust, and recombining plasma, but we have no direct view on the stellar object. Aims. Within the merger hypothesis, we aim to infer the most likely configuration of the progenitor system that resulted in Nova 1670. Methods. We collected the literature data on the physical properties of the outburst and the remnant (including their energetics), and on the chemical composition of the circumstellar material (including elemental and isotopic abundances). These data, which result from optical and submillimeter observations of the circumstellar gas of CK Vul, are summarized here. We performed simple simulations to analyze the form and the level of mixing within the material associated with the merger. We identified products of nuclear burning, among which we find ashes of hydrogen burning in the CNO cycles and in the MgAl chain, as well as signs of partial helium burning. Results. Based primarily on the luminosity and chemical composition of the remnant, we find that the progenitor primary had to be an evolutionarily advanced red giant branch star of a mass of 1–2 M⊙. The secondary was either a very similar giant, or –more likely– a helium white dwarf. While the eruption event was mainly powered by accretion, we estimate that about 12% of total energy is likely to have come from helium burning activated during the merger. The coalescence of a first-ascent giant with a helium white dwarf created a star with a rather unique internal structure and composition, which resemble those of early R-type carbon stars. Conclusions. Nova 1670 is the result of a merger between a helium white dwarf and a first-ascent red giant and is likely now evolving to become an early R-type carbon star. [ABSTRACT FROM AUTHOR]

Rest-frame far-ultraviolet spectra are fundamental to our understanding of star-forming galaxies, providing a unique window on massive star populations (MSs), chemical evolution, feedback processes, and reionization. JWST is ushering in a new era, pushing the FUV frontier beyond z=10. The success of such endeavors hinges on a comprehensive understanding of the MSs and gas conditions that power the observed spectra. The COS Legacy Archive Spectroscopic SurveY (CLASSY) is a powerful and promising solution providing high-quality, high-resolution FUV spectra of 45 nearby star-forming galaxies. The spectra contain a suite of features that simultaneously characterize the MSs that populate metal-poor galaxies, physical properties of large-scale outflows, and chemical abundance patterns. The CLASSY sample is consistent with the z 0 mass-metallicity relationship and spans 1.5 dex in metallicity. These unique properties make CLASSY the benchmark training set for studies of MSs in star-forming galaxies both across cosmic time and connecting metal-poor to metal-rich populations. [ABSTRACT FROM AUTHOR]

While we have growing numbers of massive star observations, it remains unclear how efficient the key physical processes such as mass loss, convection and rotation are, or indeed how they impact each other. We reconcile this with detailed stellar evolution models, yet these models have their own drawbacks with necessary assumptions for 3-dimensional processes like rotation which need to be adapted into 1-dimensional models. The implementation of empirical mass-loss prescriptions in stellar evolution codes can lead to the extrapolation of base rates to unconstrained evolutionary stages leading to a range of uncertain fates. In short, there remain many free parameters and physical processes which need to be calibrated in order to align our theory better with upcoming observations. We have tested various processes such as mass loss and internal mixing, including rotational mixing and convective overshooting, against multiple observational constraints such as using eclipsing binaries, the Humphreys-Davidson limit, and the final masses of Wolf-Rayet stars, across a range of metallicities. In fact, we developed a method of disentangling the effects of mixing and mass loss in the 'Mass-Luminosity Plane' allowing direct calibration of these processes. In all cases, it is important to note that a combined appreciation for both stellar winds and internal mixing are important to reproduce observations. [ABSTRACT FROM AUTHOR]

UV wind line variability in OB stars appears to be universal. We review the evidence that the variability is due to large, dense, optically thick structures rooted in or near the photosphere. Using repeated observations and a simple model we translate observed profile variations into optical depth variations and, consequently, variations in measured mass loss rates. Although global rates may be stable, measured rates vary. Consequently, profile variations infer how mass loss rates determined from UV wind lines vary. These variations quantify the intrinsic error inherent in any mass loss rate derived from a single observation. These derived rates can differ by factors of 3 or more. Our results also imply that rates from non-simultaneous observations (such as UV and ground based data) need not agree. Finally, we use our results to examine the nature of the structures responsible for the variability. [ABSTRACT FROM AUTHOR]

Adequate stellar atmosphere models are prerequisite to derive robust stellar parameters from spectroscopic analyses. I will briefly review recent results obtained with the Potdam Wolf-Rayet (PoWR) model atmosphere code, which is applicable to all types of hot stars. Using multi-wavelength observations including the UV, we analyzed large samples of massive stars at various metallicities, gaining important insights on their cosmic role and the feedback to their environment. A recent extension of PoWR allows to compose the model atmosphere from two zones. A rapidly rotating star, e.g., might possess a cooler equatorial region with a slow wind, and two polar cones with higher photospheric temperature and fast wind. For two examples of rapidly rotating O-type stars, we demonstrate that such model can reproduce wind-line profiles which otherwise would stay inconsistent. Fast rotation, which prevails in particular at low metallicities, thus might bias empirically derived parameters, having implications for feedback as well as for angular-momentum losses of SN and GRB progenitors. [ABSTRACT FROM AUTHOR]

Line-driven stellar winds are ubiquitous among hot massive stars. In some cases they can become so strong, that the whole star is cloaked by an optically thick wind. The strong outflow gives rise to large emission lines, defining the class of so-called Wolf-Rayet (WR) stars. While being major players in the evolution of massive stars, the formation of heavy black holes, and the distribution of elements, the occurrence and nature of WR winds is still quite enigmatic. A promising instrument towards a better theoretical understanding are stellar atmospheres allowing for a consistent inclusion of the hydrodynamics. By coupling stellar and wind parameters and the inclusion of a detailed non-LTE radiative transfer, they allow us to go beneath the observable layers and study the onset of WR-type winds. Establishing larger sets of models, we were able to make ground-breaking progress by identifying trends with mass and metallicity that deviate significantly from present empirical descriptions. Our modelling efforts reveal a complex picture for WR-type winds with strong, non-linear dependencies. Besides covering metallicity and mass, we further identify surface hydrogen as an important ingredient to retain WR-type mass loss at lower metallicity. Here, we present a summary of recent insights on the nature and onset of WR-type winds in massive stars including the consequences for stellar evolution, remaining open questions, and current efforts to overcome them. [ABSTRACT FROM AUTHOR]

Luminous blue variables (LBVs) and B[e] supergiants (B[e]SGs) are some of the most massive stars that display extreme and puzzling behavior. Their rarity indicates that they belong to short evolutionary phases or short-lived phenomena in the post-main sequence evolution of massive stars. However, their strong mass loss and episodic mass eruptions may be crucially impacting massive star evolution. LBVs are a group of evolved massive stars that exhibit irregular variability and eruptive mass loss. Various subtypes, including S Doradus variables, giant eruptions, and pre-supernova outbursts, exist. The physical cause of the LBV phenomenon remains heavily debated. B[e]SGs have strong forbidden line emission and infrared excess from dust that are thought to arise in a circumstellar disk or torus. The formation mechanism of their disk-like structures is yet to be settled. The evolutionary phases of LBVs and B[e]SGs and their connection to other evolved massive stars are important unanswered questions in massive star evolution. [ABSTRACT FROM AUTHOR]

We present photometry and spectroscopy of the slowly evolving superluminous Type IIn supernova (SN) 2015da. SN 2015da is extraordinary for its very high peak luminosity, and also for sustaining a high luminosity for several years. Even at 8 yr after explosion, SN 2015da remains as luminous as the peak of a normal SN II-P. The total radiated energy integrated over this time period (with no bolometric correction) is at least |$1.6 \times 10^{51}$| erg (or 1.6 FOE). Including a mild bolometric correction, adding kinetic energy of the expanding cold dense shell of swept-up circumstellar material (CSM), and accounting for asymmetry, the total explosion kinetic energy was likely 5–10 FOE. Powering the light curve with CSM interaction requires an energetic explosion and 20 M |$_{\odot }$| of H-rich CSM, which in turn implies a massive progenitor system |$\gt $| 30 M |$_{\odot }$|⁠. Narrow P Cyg features show steady CSM expansion at 90 km s |$^{-1}$|⁠ , requiring a high average mass-loss rate of |$\sim$| 0.1 M |$_{\odot }$|  yr |$^{-1}$| sustained for two centuries before explosion (although ramping up toward explosion time). No current theoretical model for single-star pre-SN mass-loss can account for this. The slow CSM, combined with broad wings of H  |$\alpha$| indicating H-rich material in the unshocked ejecta, disfavours a pulsational pair instability model for the pre-SN mass-loss. Instead, violent pre-SN binary interaction is a likely culprit. Finally, SN 2015da exhibits the characteristic asymmetric blueshift in its emission lines from shortly after peak until the present epoch, adding another well-studied superluminous SNe IIn with unambiguous evidence of post-shock dust formation. [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]

Dust populations in space plasmas are often described by a size distribution function, generally a power law distribution. In view of that, we include this feature in the kinetic description of a homogeneous magnetized dusty plasma with electrically charged immobile dust grains, in order to study its effects in the propagation and damping of Alfvén waves. The dispersion relation is numerically solved using parameters typically found in the dust-driven stellar winds of carbon-rich stars and in Earth's auroral acceleration region, two space systems with unalike plasma parameters and in which Alfvén waves are known to play important roles in the plasma acceleration and heating processes. We show that the characteristics of the normal modes, namely the ion cyclotron and whistler modes, will change when one considers a power law distribution of dust sizes in the theory, as compared to a mono-sized dust population; and that these differences will depend on the exponent p of the power law, which alters the plasma charge imbalance between electrons and ions. We also notice that power-law distribution functions will modify the waves' damping rate values. In particular, we show that in a stellar wind environment the ion cyclotron mode at very small wavenumber decreases with the reduction of p , while for higher wavenumber the damping of this mode increases with the reduction of p. For the Earth's magnetosphere, the results obtained show that the wave damping increases with the decrease of p for all wavenumbers, for the parameters considered in the analysis. [ABSTRACT FROM AUTHOR]

Growing supermassive black holes (Active Galactic Nuclei; AGN) release energy with the potential to alter their host galaxies and larger-scale environment; a process named "AGN feedback". Feedback is a required component of galaxy formation models and simulations to explain the observed properties of galaxy populations. We provide a broad overview of observational approaches that are designed to establish the physical processes that couple AGN energy to the multi-phase gas, or to find evidence that AGN impact upon galaxy evolution. The orders-of-magnitude range in spatial, temporal, and temperature scales, requires a diverse set of observational studies. For example, studying individual targets in detail sheds light on coupling mechanisms; however, evidence for the long-term impact of AGN is better established within galaxy populations that are not necessarily currently active. We emphasise how modern surveys have revealed the importance of radio emission for identifying and characterising feedback mechanisms. At the achieved sensitivities, the detected radio emission can trace a range of processes, including a shocked interstellar medium caused by AGN outflows (driven by various mechanisms including radiation pressure, accretion disc winds, and jets). We also describe how interpreting observations in the context of theoretical work can be challenging, in part, due to some of the adopted terminology. [ABSTRACT FROM AUTHOR]

We present time-series linear-polarization observations of the bright O4 supergiant ζ Puppis. The star is found to show polarization variation on time-scales of around an hour and longer. Many of the observations were obtained contemporaneously with Transiting Exoplanet Survey Satellite (TESS) photometry. We find that the polarization varies on similar time-scales to those seen in the TESS light curve. The previously reported 1.78-d photometric periodicity is seen in both the TESS and polarization data. The amplitude ratio of photometry to polarization is ∼9 for the periodic component and the polarization variation is oriented along position angle ∼70°–160°. Higher frequency stochastic variability is also seen in both data sets with an amplitude ratio of ∼19 and no preferred direction. We model the polarization expected for a rotating star with bright photospheric spots and find that models that fit the photometric variation produce too little polarization variation to explain the observations. We suggest that the variable polarization is more likely the result of scattering from the wind, with corotating interaction regions producing the periodic variation and a clumpy outflow producing the stochastic component. The Hα emission line strength was seen to increase by 10 per cent in 2021 with subsequent observations showing a return to the pre-2018 level. [ABSTRACT FROM AUTHOR]

We explore the impact of a magnetar giant flare (GF) on the neutron star (NS) crust, and the associated baryon mass ejection. We consider that sudden magnetic energy dissipation creates a thin high-pressure shell above a portion of the NS surface, which drives a relativistic shockwave into the crust, heating a fraction of these layers sufficiently to become unbound along directions unconfined by the magnetic field. We explore this process using spherically symmetric relativistic hydrodynamical simulations. For an initial shell pressure P GF we find the total unbound ejecta mass roughly obeys the relation |$M_{\rm {ej}}\sim 4\!-\!9\times 10^{24}\, \rm {g}\, (P_{\rm GF}/10^{30}\, \rm {erg}\, \rm {cm}^{-3})^{1.43}$|⁠. For |$P_{\rm {GF}}\sim 10^{30}\!-\!10^{31}\, \rm {erg}\, \rm {cm}^{-3}$| corresponding to the dissipation of a magnetic field of strength |$\sim 10^{15.5}\!-\!10^{16}\, \rm {G}$|⁠ , we find |$M_{\rm {ej}}\sim 10^{25}\!-\!10^{26}\, \rm {g}$| with asymptotic velocities v ej/ c ∼ 0.3–0.6 compatible with the ejecta properties inferred from the afterglow of the 2004 December GF from SGR 1806-20. Because the flare excavates crustal material to a depth characterized by an electron fraction Y e ≈ 0.40–0.46, and is ejected with high entropy and rapid expansion time-scale, the conditions are met for heavy element r -process nucleosynthesis via the alpha-rich freeze-out mechanism. Given an energetic GF rate of roughly once per century in the Milky Way, we find that magnetar GFs could be an appreciable heavy r -process source that tracks star formation. We predict that GFs are accompanied by short ∼minutes long, luminous |$\sim 10^{39}\, \rm {erg}\, \rm {s}^{-1}$| optical transients powered by r -process decay (nova brevis), akin to scaled-down kilonovae. Our findings also have implications for the synchrotron nebulae surrounding some repeating fast radio burst sources. [ABSTRACT FROM AUTHOR]

Solar-type stars form with a wide range of rotation rates Ω. A wide Ω range persists until a stellar age of t  ∼ 0.6 Gyr, after which solar-type stars exhibit Skumanich spin-down where Ω ∝  t −1/2. Rotational evolution models incorporating polytropic stellar winds struggle to simultaneously reproduce these two regimes, namely the initially wide Ω range and the Skumanich spin-down without imposing an a priori cap on the wind mass-loss rate. We show that a three-dimensional wind model driven by Alfvén waves and observational data yields wind torques that agree with the observed age distribution of Ω. In our models of the Sun and 27 open cluster stars aged from 0.04 to 0.6 Gyr that have observationally derived surface magnetic maps and rotation rates, we find evidence of exponential spin-down in young stars that are rapid rotators and Skumanich spin-down for slow rotators. The two spin-down regimes emerge naturally from our data-driven models. Our modelling suggests that the observed age distribution of stellar rotation rates Ω arises as a consequence of magnetic field strength saturation in rapid rotators. [ABSTRACT FROM AUTHOR]

Only 3–4 per cent of Galactic O stars are observed to display the emission features representative of the OBe phenomenon, compared to Galactic B stars, which display these characteristics in 25–35 per cent of B0 and B1 stars. We present new observations of the high-mass O star, VES 735, which confirms its classification as one of these rare emission-line stars. These are its first recorded observations that display strong spectroscopic variations in nearly 30 yr of monitoring, with the H α profile exhibiting a 10-fold increase in emission compared to observations taken between 1996 and 2014 and having variations which show episodes of inflowing and outflowing material. These observations, coupled with photometric variations in the visible and infrared, show behaviour that is consistent with the mass reservoir effect for viscous decretion discs. We propose that in 2015 VES 735 began an approximately 3 yr event in which mass was being injected into the circumstellar environment followed by re-accretion towards the star. We also find evidence that the re-accretion may have been interrupted with another, smaller, mass-injection event based on observations in 2022 and 2023. Observational cadences ranging from hours to months show no evidence that VES 735 is part of a binary system, making it an ideal candidate for future observations to further investigate the evolution of high-mass stars and the OBe phenomenon as it pertains to their circumstellar environment. [ABSTRACT FROM AUTHOR]

Shock boundary layers are regions bounded by a shock wave on the one side and tangential discontinuity on the other side. These boundary layers are commonly observed in astrophysics. For example, they exist in the regions of the interaction of the stellar winds with the surrounding interstellar medium. Additionally, the shock layers (SLs) are often penetrated by the flows of interstellar atoms, as, for instance, in the astrospheres of the stars embedded by the partially ionized interstellar clouds. This paper presents a simple toy model of an SL that aims to qualitatively describe the influence of charge exchange with interstellar hydrogen atoms on the plasma flow in astrospheric SLs. To clearly explore this effect, magnetic fields are neglected, and the geometry is kept as simple as possible. The model explains why the cooling of plasma due to charge exchange in the inner heliosheath leads to an increase in plasma density in front of the heliopause. It also demonstrates that the source of momentum causes changes in the pressure profile within the SL. The paper also discusses the decrease in plasma density near the astropause in the outer SL in the case of layer heating, which is particularly relevant in light of the Voyager measurements in the heliospheric SL. [ABSTRACT FROM AUTHOR]

A long-lived magnetar, potentially originating from a binary neutron star system, has been proposed to explain the extended emission observed in certain short-duration gamma-ray bursts (sGRBs), and is posited as a potential central engine to power the engine-fed kilonovae. Previously, the process by which energy is injected into the surrounding ejecta/jet was widely believed to be nearly isotropic. In this study, we employ special relativity magnetohydrodynamic (SRMHD) simulations to investigate the wind injection process from a magnetar central engine. We explore the dynamics and energy distribution within the system and found that the parameter α = u A/ u MWN can be used to indicate the collimation of the magnetar wind energy injection, where u A is the local Alfven four-speed and u MWN is the four-speed of the magnetar wind nebular (MWN) formed from wind-ejecta collision. A significant portion of the injected energy from the magnetar spin-down wind will be channeled to the jet axis due to collimation within the MWN. Achieving isotropic energy injection requires a significantly small α that necessitates either an ultra-relativistic expanding MWN or an extremely low magnetization MWN, both of which are challenging to attain in sGRBs. Consequently, a considerably reduced energy budget (i.e. energy per solid angle reduced by a factor of up to 10 with respect to the value under isotropic assumption) is anticipated to be injected into the ejecta for engine-fed kilonovae. Engine-fed kilonovae would appear fainter than originally anticipated. [ABSTRACT FROM AUTHOR]

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle-5 observations of HBC 494, as well as calculations of the kinematic and dynamic variables which represent the object's wide-angle bipolar outflows. HBC 494 is a binary FU Orionis type object located in the Orion A molecular cloud. We take advantage of combining the ALMA main array, Atacama Compact Array (ACA), and Total Power (TP) array in order to map HBC 494's outflows and thus, estimate their kinematic parameters with higher accuracy in comparison to prior publications. We use 12CO, 13CO, C18O, and SO observations to describe the object's outflows, envelope, and disc, as well as estimate the mass, momentum, and kinetic energy of the outflows. After correcting for optical opacity near systemic velocities, we estimate a mass of 3.0 × 10−2 M⊙ for the southern outflow and 2.8 × 10−2 M⊙ for northern outflow. We report the first detection of a secondary outflow cavity located approximately 15 arcsec north of the central binary system, which could be a remnant of a previous large-scale accretion outburst. Furthermore, we find CO spatial features in HBC 494's outflows corresponding to position angles of ∼35° and ∼145°. This suggests that HBC 494's outflows are most likely a composite of overlapping outflows from two different sources, i.e. HBC 494a and HBC 494b, the two objects in the binary system. [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]

B[e] supergiants (sgB[e]) are rare objects whose evolutionary stage remains uncertain. Observationally, they display strong Balmer emission lines, infrared excess, and intrinsic polarization, indicating a non-spherical circumstellar envelope. We present a study of the sgB[e] RMC 82, using new spectropolarimetric data complemented by photometry from the ultraviolet (UV) to the mid-infrared. Our two-component model comprises a slow, dense equatorial wind wherein dust grains form and a fast polar wind. We applied the hdust radiative transfer code and Bayesian statistics to infer the parameters from a grid of 3240 pre-computed models. The model accurately reproduces the spectral energy distribution and polarized spectrum, but struggles to match the H α emission. Our results suggest a large mass-loss rate of |$6.6 \times 10^{-6}\, \mathrm{{\rm M}_{\odot }\, yr^{-1}\, sr^{-1}}$|⁠. The dense wind is confined within an opening angle of 11°. The hottest dust grains are located at 277  R * with a temperature of 870 K. The dust grains are porous, with a density of 0.051  |$\rm {g\, cm^{-3}}$|⁠. The central star was found to be significantly hotter than previous estimates (T eff =  |$27\, 000$|  K). By comparing models with different components, we find that gas reprocesses a significant amount of UV radiation, shielding the dust. However, the dust also scatters UV photons back to the inner disc, increasing its temperature and H α emission. We conclude that self-consistent models, that account for the gas–dust interplay in the envelope, are essential for studying sgB[e] and similar objects. [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]