Search

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]

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]

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]

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]

We present high-angular-resolution observations of Sakurai's object using the Atacama Large Millimeter Array, shedding new light on its morpho-kinematical structure. The millimetre continuum emission observed at an angular resolution of 20 milliarcsec (corresponding to 70 AU) reveals a bright compact central component whose spectral index indicates that it is composed of amorphous carbon dust. Based on these findings, we conclude that this emission traces the previously suggested dust disc observed in mid-infrared observations, and therefore our observations provide the first direct imaging of this disc. The H12CN(J = 4 → 3) line emission observed at an angular resolution of 300 milliarcsec (corresponding to 1000 AU) displays a bipolar structure with a north–south velocity gradient. From the position–velocity diagram of this emission, we identify the presence of an expanding disc and a bipolar molecular outflow. The inclination of the disc is determined to be i = 72°. The derived values for the de-projected expansion velocity and the radius of the disc are vexp = 53 km s−1 and R = 277 AU, respectively. On the other hand, the de-projected expansion velocity of the bipolar outflow detected in the H12CN(J = 4 → 3) emission is of approximately 1000 km s−1. We propose that the molecular outflow has an hourglass morphology with an opening angle of around 60°. Our observations unambiguously show that an equatorial disc and bipolar outflows formed in Sakurai's object during the 30 years following the occurrence of the born-again event, providing important constraints for future modelling efforts of this phenomenon. [ABSTRACT FROM AUTHOR]

Using observations with e-MERLIN and the Karl G. Jansky Very Large Array (VLA), together with archival data from the Atacama Large Millimeter/submillimeter Array (ALMA), we obtain high-resolution radio images of two binary YSOs: L1551 IRS 5 and L1551 NE, covering a wide range of frequencies from 5 to 336 GHz, and resolving emission from the radio jet on scales of only ∼15 au. By comparing these observations to those from a previous epoch, it is shown that there is a high degree of variability in the free-free emission from the jets of these sources. In particular, the northern component of L1551 IRS 5 shows a remarkable decline in flux density of a factor of ∼5, suggesting that the free-free emission of this source has almost disappeared. By fitting the spectra of the sources, the ionised mass-loss rates of the jets were derived and it was shown that there is significant variability of up to a factor of ∼6 on timescales of ∼20 years. Using radiative transfer modelling, we also obtained a model image for the jet of the southern component of L1551 IRS 5 to help study the inner region of the ionised high-density jet. The findings favour the X-wind model launched from a very small innermost region. [ABSTRACT FROM AUTHOR]

Context. B supergiants (BSGs) are evolved stars with effective temperatures between ~10 to ~30 kK. Knowing the properties of these objects is important to understand massive star evolution. Located on the cool end of the line-driven wind regime, the study of their atmospheres can help us to understand the physics of their winds and phenomena such as the bi-stability jump. Aims. Despite being well-studied stars, key UV features of their spectra have so far not been reproduced by atmosphere models for spectral types later than B1. In this study, we aim to remedy this situation by performing quantitative spectral analyzes that account for the effects of X-rays and clumping in the wind. In addition, we also briefly investigate the evolutionary status of our sample stars based on the stellar parameters we obtained. Methods. We determined photospheric and wind parameters via quantitative spectroscopy using atmosphere models computed with CMFGEN and PoWR. These models were compared to high-resolution UV and optical spectra of four BSGs: HD206165, HD198478, HD53138, and HD164353. We further employed GENEC and MESA tracks to study the evolutionary status of our sample. Results. When including both clumping and X-rays, we obtained a good agreement between synthetic and observed spectra for our sample stars. For the first time, we reproduced key wind lines in the UV, where previous studies were unsuccessful. To model the UV spectra, we require a moderately clumped wind (fV∞ ≳ 0.5). We also infer a relative X-ray luminosity of about 10−7.5 to 10−8, which is lower than the typical ratio of 10−7. Moreover, we find a possible mismatch between evolutionary mass predictions and the derived spectroscopic masses, which deserves deeper investigation as this might relate to the mass-discrepancy problem present in other types of OB stars. Conclusions. Our results provide direct spectroscopic evidence that both X-rays and clumping need to be taken into account to describe the winds of cool BSGs. However, their winds seem to be much less structured than in earlier OB-type stars. Our findings are in line with observational X-rays and clumping constraints as well as recent hydrodynamical simulations. The evolutionary status of BSGs seems to be diverse with some objects potentially being post-red supergiants or merger products. The obtained wind parameters provide evidence for a moderate increase of the mass-loss rate around the bi-stability jump. [ABSTRACT FROM AUTHOR]

Stellar winds of cool and pulsating asymptotic giant branch (AGB) stars enrich the interstellar medium with large amounts of processed elements and various types of dust. We present the first study on the influence of gas-to-dust drift on ab initio simulations of stellar winds of M-type stars driven by radiation pressure on forsterite particles. Our study is based on our radiation hydrodynamic model code T-800 that includes frequency-dependent radiative transfer, dust extinction based on Mie scattering, grain growth and ablation, gas-to-dust drift using one mean grain size, a piston that simulates stellar pulsations, and an accurate high spatial resolution numerical scheme. To enable this study, we calculated new gas opacities based on the EXOMOL database, and we extended the model code to handle the formation of minerals that may form in M-type stars. We determine the effects of drift by comparing drift models to our new and extant non-drift models. Three out of four new drift models show high drift velocities, 87–310 km s−1. Our new drift model mass-loss rates are 1.7–13 per cent of the corresponding values of our non-drift models, but compared to the results of two extant non-drift models that use the same stellar parameters, these same values are 0.33–1.5 per cent. Meanwhile, a comparison of other properties such as the expansion velocity and grain size show similar values. Our results, which are based on single-component forsterite particles, show that the inclusion of gas-to-drift is of fundamental importance in stellar wind models driven by such transparent grains. Assuming that the drift velocity is insignificant, properties such as the mass-loss rate may be off from more realistic values by a factor of 50 or more. [ABSTRACT FROM AUTHOR]

Context. Young stellar systems actively accrete from their circumstellar disk and simultaneously launch outflows. The physical link between accretion and ejection processes remains to be fully understood. Aims. We investigate the structure and dynamics of magnetospheric accretion and associated outflows on a scale smaller than 0.1 au around the young transitional disk system GM Aur. Methods. We devised a coordinated observing campaign to monitor the variability of the system on timescales ranging from days to months, including partly simultaneous high-resolution optical and near-infrared spectroscopy, multiwavelength photometry, and lowresolution near-infrared spectroscopy, over a total duration of six months, covering 30 rotational cycles. We analyzed the photometric and line profile variability to characterize the accretion and ejection processes. Results. The optical and near-infrared light curves indicate that the luminosity of the system is modulated by surface spots at the stellar rotation period of 6.04 ± 0.15 days. Part of the Balmer, Paschen, and Brackett hydrogen line profiles as well as the HeI 5876Å and HeI 10830Å line profiles are modulated on the same period. The Pa line flux correlates with the photometric excess in the u0 band, which suggests that most of the line emission originates from the accretion process. High-velocity redshifted absorptions reaching below the continuum periodically appear in the near-infrared line profiles at the rotational phase in which the veiling and line fluxes are the largest. These are signatures of a stable accretion funnel flow and associated accretion shock at the stellar surface. This large-scale magnetospheric accretion structure appears fairly stable over at least 15 and possibly up to 30 rotational periods. In contrast, outflow signatures randomly appear as blueshifted absorption components in the Balmer and HeI 10830Å line profiles. They are not rotationally modulated and disappear on a timescale of a few days. The coexistence of a stable, large-scale accretion pattern and episodic outflows supports magnetospheric ejections as the main process occurring at the star-disk interface. Conclusions. Long-term monitoring of the variability of the GM Aur transitional disk system provides clues to the accretion and ejection structure and dynamics close to the star. Stable magnetospheric accretion and episodic outflows appear to be physically linked on a scale of a few stellar radii in this system. [ABSTRACT FROM AUTHOR]

Stripped-envelope stars can be observed as Wolf-Rayet (WR) stars or as less luminous hydrogen-poor stars with low mass-loss rates and transparent winds. Both types are potential progenitors of Type I core-collapse supernovae (SNe). We used grids of core-collapse models obtained from single helium stars at different metallicities to study the effects of metallicity on the transients and remnants these stars produce. We characterised the surface and core properties of our core-collapse models and investigated their 'explodability' using three criteria. In the cases where explosions are predicted, we estimated the ejecta mass, explosion energy, nickel mass, and neutron star (NS) mass. Otherwise, we predicted the mass of the resulting black hole (BH). We constructed a simplified population model and find that the properties of SNe and compact objects depend strongly on metallicity. The ejecta masses and explosion energies for Type Ic SNe are best reproduced by models with Z = 0.04 that exhibit strong winds during core helium burning. This implies that either their mass-loss rates are underestimated or that Type Ic SN progenitors experience mass loss through other mechanisms before exploding. The distributions of ejecta masses, explosion energies, and nickel mass for Type Ib SNe are not well reproduced by progenitor models with WR mass loss, but are better reproduced if we assume no mass loss in progenitors with luminosities below the minimum WR star luminosity. We find that Type Ic SNe become more common as metallicity increases, and that the vast majority of progenitors of Type Ib SNe must be transparent-wind stripped-envelope stars. We find that several models with pre-collapse CO masses of up to ∼30 M⊙ may form ∼3 M⊙ BHs in fallback SNe. This may have important consequences for our understanding of SNe, binary BH and NS systems, X-ray binary systems, and gravitational wave transients. [ABSTRACT FROM AUTHOR]

Context. Classical T Tauri stars (CTTs) magnetically interact with their surrounding disks, a process that is thought to regulate their rotational evolution. Aims. We compute torques acting on the stellar surface of CTTs that arise from different accreting (accretion funnels) and ejecting (stellar winds and magnetospheric ejections) flow components. Furthermore, we compare the magnetic braking due to stellar winds in two different systems: isolated (i.e., weak-line T Tauri and main-sequence) and accreting (i.e., classical T Tauri) stars. Methods. We use 2.5D magnetohydrodynamic, time-dependent, axisymmetric simulations that were computed with the PLUTO code. For both systems, the stellar wind is thermally driven. In the star-disk-interaction (SDI) simulations, the accretion disk is Keplerian, viscous, and resistive, and is modeled with an alpha prescription. Two series of simulations are presented, one for each system (i.e., isolated and accreting stars). Results. In classical T Tauri systems, the presence of magnetospheric ejections confines the stellar-wind expansion, resulting in an hourglass-shaped geometry of the outflow, and the formation of the accretion columns modifies the amount of open magnetic flux exploited by the stellar wind. These effects have a strong impact on the stellar-wind properties, and we show that the stellar-wind braking is more efficient in the SDI systems than in the isolated ones. We further derive torque scalings over a wide range of magnetic field strengths for each flow component in an SDI system (i.e., magnetospheric accretion and ejections, and stellar winds), which directly applies a torque on the stellar surface. Conclusions. In all the performed SDI simulations, the stellar wind extracts less than 2% of the mass accretion rate and the disk is truncated by up to 66% of the corotation radius. All simulations show a net spin-up torque. We conclude that in order to achieve a stellar-spin equilibrium, we need either more massive stellar winds or disks that are truncated closer to the corotation radius, which increases the torque efficiency of the magnetospheric ejections. [ABSTRACT FROM AUTHOR]

Context. Pulsar wind nebulae are efficient particle accelerators, and yet the processes at work remain elusive. Self-generated, microturbulence is too weak in relativistic magnetized shocks to accelerate particles over a wide energy range, suggesting that the global dynamics of the nebula may be involved in the acceleration process instead. Aims. In this work, we study the role played by the large-scale anisotropy of the transverse magnetic field profile on the shock dynamics. Methods. We performed large two-dimensional particle-in-cell simulations for a wide range of upstream plasma magnetizations, from weakly magnetized to strongly magnetized pulsar winds. Results. The magnetic field anisotropy leads to a dramatically different structure of the shock front and downstream flow. A large-scale velocity shear and current sheets form in the equatorial regions and at the poles, where they drive strong plasma turbulence via Kelvin-Helmholtz vortices and kinks. The mixing of current sheets in the downstream flow leads to efficient nonthermal particle acceleration. The power-law spectrum hardens with increasing magnetization, akin to those found in relativistic reconnection and kinetic turbulence studies. The high end of the spectrum is composed of particles surfing on the wake produced by elongated spearhead-shaped cavities forming at the shock front and piercing through the upstream flow. These particles are efficiently accelerated via the shear-flow acceleration mechanism near the Bohm limit. Conclusions. Magnetized relativistic shocks are very efficient particle accelerators. Capturing the global dynamics of the downstream flow is crucial to understanding them, and therefore local plane parallel studies may not be appropriate for pulsar wind nebulae and possibly other astrophysical relativistic magnetized shocks. A natural outcome of such shocks is a variable and Doppler-boosted synchrotron emission at the high end of the spectrum originating from the shock-front cavities, reminiscent of the mysterious Crab Nebula gamma-ray flares. [ABSTRACT FROM AUTHOR]

Context. Hot subdwarfs are helium burning objects in late stages of their evolution. These subluminous stars can develop winds driven by light absorption in the lines of heavier elements. The wind strength depends on chemical composition which can significantly vary from star to star. Aims. We aim to understand the influence of metallicity on the strength of the winds of the hot hydrogen-rich subdwarfs HD 49798 and BD+18° 2647. Methods. We used high-resolution UV and optical spectra to derive stellar parameters and abundances using the TLUSTY and SYNSPEC codes. For derived stellar parameters, we predicted wind structure (including mass-loss rates and terminal velocities) with our METUJE code. Results. We derived effective temperature Teff = 45 900 K and mass M = 1.46 M⊙ for HD 49798 and Teff = 73 000 K and M = 0.38 M⊙ for BD+18° 2647. The derived surface abundances can be interpreted as a result of interplay between stellar evolution and diffusion. The subdwarf HD 49798 has a strong wind that does not allow for chemical separation and consequently the star shows solar chemical composition modified by hydrogen burning. On the other hand, we did not find any wind in BD+18° 2647 and its abundances are therefore most likely affected by radiative diffusion. Accurate abundances do not lead to a significant modification of wind mass-loss rate for HD 49798, because the increase of the contribution of iron and nickel to the radiative force is compensated by the decrease of the radiative force due to other elements. The resulting wind mass-loss rate Ṁ = 2.1 × 10−9M⊙ yr−1 predicts an X-ray light curve during the eclipse which closely agrees with observations. On the other hand, the absence of the wind in BD+18° 2647 for accurate abundances is a result of its peculiar chemical composition. Conclusions. Wind models with accurate abundances provide more reliable wind parameters, but the influence of abundances on the wind parameters is limited in many cases. [ABSTRACT FROM AUTHOR]

η Car is one of the most luminous and massive stars in our Galaxy and is the brightest mid-IR source in the sky outside our solar system. Since the late 1990s, the central source has dramatically brightened at UV and optical wavelengths. This might be explained by a decrease in circumstellar dust extinction. We aim to establish the mid-IR flux evolution and further our understanding of the star's UV and optical brightening. Mid-IR images from 8−20 μm were obtained in 2018 with VISIR at the Very Large Telescope. Archival data from 2003 and 2005 were retrieved from the ESO Science Archive Facility, and historical records were collected from publications. We present mid-IR images of η Car with the highest angular resolution to date at the corresponding wavelengths (≥0.22″). We reconstruct the mid-IR evolution of the spectral energy distribution of the spatially integrated Homunculus nebula from 1968 to 2018 and find no long-term changes. The bolometric luminosity of η Car has been stable over the past five decades. We do not observe a long-term decrease in the mid-IR flux densities that could be associated with the brightening at UV and optical wavelengths, but circumstellar dust must be declining in our line of sight alone. Short-term flux variations within about 25% of the mean levels could be present. [ABSTRACT FROM AUTHOR]

Context. High-mass microquasars (HMMQs) are systems from which relativistic jets are launched. At the scales of several times the binary system size, the jets are expected to follow a helical path caused by the interaction with a strong stellar wind and orbital motion. Such a trajectory has its influence on the non-thermal emission of the jets, which also depends strongly on the observing angle due to Doppler boosting effects. Aims. We explore how the expected non-thermal emission of HMMQ jets at small scales is affected by the impact of the stellar wind and the orbital motion on the jet propagation. Methods. We studied the broadband non-thermal emission, from radio to gamma rays, produced in HMMQ jets up to a distance of several orbital separations, taking into account a realistic jet trajectory, different model parameters, and orbital modulation. The jet trajectory is computed by considering momentum transfer with the stellar wind. Electrons are injected at the position where a recollimation shock in the jets is expected due to the wind impact. Their distribution along the jet path is obtained assuming local acceleration at the recollimation shock, and cooling via adiabatic, synchrotron, and inverse Compton processes. The synchrotron and inverse Compton emission is calculated taking into account synchrotron self-absorption within the jet, free-free absorption with the stellar wind, and absorption by stellar photons via pair production. Results. The spectrum is totally dominated by the jet over the counter-jet due to Doppler boosting. Broadband emission from microwaves to gamma rays is predicted, with radio emission being totally absorbed. This emission is rather concentrated in the regions close to the binary system and features strong orbital modulation at high energies. Asymmetric light curves are obtained owing to the helical trajectory of the jets. Conclusions. The presence of helical shaped jets could be inferred from asymmetries in the light curves, which become noticeable only for large jet Lorentz factors and low magnetic fields. Model parameters could be constrained if accurate phase-resolved light curves from GeV to TeV energies were available. The predictions for the synchrotron and the inverse Compton radiation are quite sensitive of the parameters determining the wind-jet interaction structure. [ABSTRACT FROM AUTHOR]

Carbon-enhanced metal-poor stars, CH stars, barium stars, and extrinsic S stars, among other classes of chemically peculiar stars, are thought to be the products of the interaction of low- and intermediate-mass binaries, which occurred when the most evolved star was in the asymptotic giant branch (AGB) phase. Binary evolution models predict that because of the large sizes of AGB stars, if the initial orbital periods of such systems are shorter than a few thousand days, their orbits should have circularised due to tidal effects. However, observations of the progeny of AGB binary stars show that many of these objects have substantial eccentricities, up to e ≈ 0.9. In this work we explore the impact of wind mass transfer on the orbital parameters of AGB binary stars by performing numerical simulations in which the AGB wind is modelled using a hydrodynamical code and the dynamics of the stars is evolved using an N-body code. We find that in most models the effect of wind mass transfer contributes to the circularisation of the orbit, but on longer timescales than tidal circularisation if e ≲ 0.4. For relatively low initial wind velocities and pseudo-synchronisation of the donor star, we find a structure resembling wind Roche-lobe overflow as the stars approach periastron. In this case, the interaction between the gas and the star is stronger than when the initial wind velocity is high and the orbit shrinks while the eccentricity decreases. In one of our models wind interaction is found to pump the eccentricity of the orbit on a similar timescale as tidal circularisation. However, since the orbit of this model is shrinking tidal effects will become stronger during the evolution of the system. Although our study is based on a small sample of models, it offers some insight into the orbital evolution of eccentric binary stars interacting via winds. A larger grid of numerical models for different binary parameters is needed to test if a regime exists where hydrodynamical eccentricity pumping can effectively counteract tidal circularisation, and if this can explain the puzzling eccentricities of the descendants of AGB binaries. [ABSTRACT FROM AUTHOR]

We present the highest angular resolution (∼20 × 15 mas–44 × 33 au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations that are currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16 km probe scales < 50 au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission that is especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O 55, 0−64, 3 ν2 = 1 (Eu = 3461.9 K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45 ± 10 M⊙. The H30α (231.9 GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the ∼10σ level. The estimated disc mass is 0.6 − 2.6 M⊙ under the optically thin assumption. Analysis of the Toomre Q parameter in the optically thin regime indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures > 150 K; this does not preclude that the substructures formed earlier through disc fragmentation. [ABSTRACT FROM AUTHOR]

One of the first massive stars detected in X-rays, Schulte 12 has remained a puzzle in several aspects. In particular, its extreme brightness both in the visible and X-ray ranges is intriguing. Thanks to Swift and XMM-Newton observations covering ∼5000 d, we report the discovery of a regular 108 d modulation in X-ray flux of unknown origin. The minimum in the high-energy flux appears due to a combination of increased absorption and decreased intrinsic emission. We examined in parallel the data from a dedicated spectroscopic and photometric monitoring in the visible and near-IR domains, complemented by archives. While a similar variation timescale is found in those data, they do not exhibit the strict regular clock found at high energies. Changes in line profiles cannot be related to binarity but rather correspond to non-radial pulsations. Considering the substantial revision of the distance of Schulte 12 from the second Gaia data release, the presence of such oscillations agrees well with the evolutionary status of Schulte 12, as it lies in an instability region of the HR diagram. [ABSTRACT FROM AUTHOR]

Context. X-ray emission from massive stars was first reported four decades ago, but the precise physics governing its formation as a function of stellar properties and binarity remains not fully understood. With the recent suggestion that such objects may be important sites of cosmic ray production, a better understanding of their high-energy properties is particularly timely. Aims. The young massive cluster Westerlund 1 provides an ideal testbed for understanding this emission, with over 50 cluster members detected in historical X-ray observations. In the decade since these data were obtained, significant new multi-epoch observations of the cluster have been made, allowing a fundamental reappraisal of the nature of both X-ray bright and dark stars. Methods. Optical spectroscopy permits accurate classification of cluster members, while multi-epoch observations of a sub-set allow identification and characterisation of the binary population. Results. A total of 45 X-ray sources within Wd1 now have precise spectral classifications. Of these, 16 have been identified as candidate or confirmed massive binaries. X-ray emission is confined to O9-B0.5 supergiants, Wolf-Rayets and a small group of highly luminous interacting/post-interaction OB+OB binaries. Despite their presence in large numbers, no emission is seen from earlier, less evolved O stars or later, cooler B super-/hypergiants. A total of 22 stars have X-ray properties that are suggestive of a contribution from emission originating in a wind collision zone. Conclusions. We suppose that the lack of X-ray emission from O giants is due to their comparatively low intrinsic bolometric luminosity if, as expected, they follow the canonical LX/Lbol relation for hot stars. The transition away from X-ray emission for OB supergiants occurs at the location of the bistability jump; we speculate that below this limit, stellar wind velocities are insufficient for internal, X-ray emitting shocks to form. Our results are consistent with recent findings that massive binaries are not uniformly brighter than single stars of comparable luminosity or spectral type, although it is noteworthy that the brightest and hardest stellar X-ray sources within Wd1 are all either confirmed or candidate massive, interacting/post-interaction binaries. [ABSTRACT FROM AUTHOR]

The understanding of the rotational evolution of early-type stars is deeply related to that of anisotropic mass and angular momentum loss. In this paper, we aim to clarify the rotational evolution of rapidly rotating early-type stars along the main sequence (MS). We have used the 2D ESTER code to compute and evolve isolated rapidly rotating early-type stellar models along the MS, with and without anisotropic mass loss. We show that stars with Z = 0.02 and masses between 5 and 7 M⊙ reach criticality during the main sequence provided their initial angular velocity is larger than 50% of the Keplerian one. More massive stars are subject to radiation-driven winds and to an associated loss of mass and angular momentum. We find that this angular momentum extraction from the outer layers can prevent massive stars from reaching critical rotation and greatly reduce the degree of criticality at the end of the MS. Our model includes the so-called bi-stability jump of the Ṁ − Teff relation of 1D-models. This discontinuity now shows up in the latitude variations of the mass-flux surface density, endowing rotating massive stars with either a single-wind regime (no discontinuity) or a two-wind regime (a discontinuity). In the two-wind regime, mass loss and angular momentum loss are strongly increased at low latitudes inducing a faster slow-down of the rotation. However, predicting the rotational fate of a massive star is difficult, mainly because of the non-linearity of the phenomena involved and their strong dependence on uncertain prescriptions. Moreover, the very existence of the bi-stability jump in mass-loss rate remains to be substantiated by observations. [ABSTRACT FROM AUTHOR]

Context. The microquasar SS 433 is the only known steady supercritical accretor in the Galaxy. It is well-known for its relativistic baryonic jets, but the system also drives equatorial outflows. These have been routinely detected in radio images, and components associated with a circumbinary disk have also been suggested in optical emission lines. Aims. We aim to spatially resolve the regions producing the stationary emission lines of SS 433 to shed light on its circumbinary structure and outflows. With an estimated binary orbit size of ≲0.1 mas, this requires optical interferometry. Methods. We use the optical interferometer VLTI+GRAVITY to spatially resolve SS 433 in the near-infrared K band at high spectral resolution (R ≈ 4000) on three nights in July 2017. This is the second such observation, after the first one in July 2016. Results. The stationary Brγ line in the 2017 observation is clearly dominated by an extended ∼1 mas ∼ 5 AU circumbinary structure perpendicular to the jets with a strong rotation component. The rotation direction is retrograde relative to the jet precession, in accordance with the slaved disk precession model. The structure has a very high specific angular momentum and is too extended to be a stable circumbinary disk in Keplerian rotation; interpreting it as such leads to a very high enclosed mass M ≳ 400 M⊙. We instead interpret it as the centrifugal ejection of the circumbinary disk, with the implication that there must be an efficient transfer of specific angular momentum from the binary to the disk. We suggest that the equatorial outflows sometimes seen in radio images result from similar episodes of circumbinary disk centrifugal ejection. In addition to the equatorial structure, we find a very extended ∼6 mas ∼ 30 AU spherical wind component to the Brγ line: the entire binary is engulfed in an optically thin spherical line emission envelope. [ABSTRACT FROM AUTHOR]

The study of massive stars in different metallicity environments is a central topic of current stellar research. The spectral analysis of massive stars requires adequate model atmospheres. The computation of such models is difficult and time-consuming. Therefore, spectral analyses are greatly facilitated if they can refer to existing grids of models. Here we provide grids of model atmospheres for OB-type stars at metallicities corresponding to the Small and Large Magellanic Clouds, as well as to solar metallicity. In total, the grids comprise 785 individual models. The models were calculated using the state-of-the-art Potsdam Wolf-Rayet (PoWR) model atmosphere code. The parameter domain of the grids was set up using stellar evolution tracks. For all these models, we provide normalized and flux-calibrated spectra, spectral energy distributions, feedback parameters such as ionizing photons, Zanstra temperatures, and photometric magnitudes. The atmospheric structures (the density and temperature stratification) are available as well. All these data are publicly accessible through the PoWR website. [ABSTRACT FROM AUTHOR]

Observational surveys are now able to detect an increasing number of transients, such as core-collapse supernovae (SN) and powerful non-terminal outbursts (SN impostors). Dedicated spectroscopic facilities can follow up these events shortly after detection. Here we investigate the properties of these explosions at early times. We use the radiative transfer code CMFGEN to build an extensive library of spectra simulating the interaction of supernovae and their progenitor's wind or circumstellar medium (CSM). We have considered a range of progenitor mass-loss rates (Ṁ = 5 × 10−4−10−2M⊙ yr−1), abundances (solar, CNO-processed, and He-rich), and SN luminosities (L = 1.9 × 108 − 2.5 × 1010 L⊙). The models simulate events approximately one day after explosion, and we assume a fixed location of the shock front as Rin = 8.6 × 1013 cm. We show that the large range of massive star properties at the pre-SN stage causes a diversity of early-time interacting SN and impostors. We identify three main classes of early-time spectra consisting of relatively high-ionisation (e.g. He II and O VI), medium-ionisation (e.g. C III and N III), and low-ionisation lines (e.g. He I and Fe II/III). They are regulated by L and the CSM density. Given a progenitor wind velocity υ∞, our models also place a lower limit of Ṁ ≳ 5 × 10−4 (υ∞/150 km s−1) M⊙ yr−1 for detection of CSM interaction signatures in observed spectra. Early-time SN spectra should provide clear constraints on progenitors by measuring H, He, and CNO abundances if the progenitors come from single stars. The connections are less clear considering the effects of binary evolution. Nevertheless, our models provide a clear path for linking the final stages of massive stars to their post-explosion spectra at early times, and guiding future observational follow-up of transients with facilities such as the Zwicky Transient Facility. [ABSTRACT FROM AUTHOR]

Aims. We aim to use the high spatial resolution of the Atacama Large Millimeter/submillimeter Array (ALMA) to map the flow pattern of molecular gas near DG Tauri and its disk, a young stellar object driving a jet and a molecular outflow. Methods. We use observations from ALMA in the J = 2 − 1 transition of 12CO, 13CO, and C18O to study the Keplerian disk of DG Tauri and outflows that may be related to the disk and the jet. Results. We find a new wind component flowing radially at a steep angle (≈25° from the vertical) above the disk with a velocity of ≈3.1 km s−1. It continues the trend of decreasing velocity for increasing distance from the jet axis (“onion-like velocity structure”). Conclusions. The new component is located close to the protostellar disk surface and may be related to photoevaporative winds. [ABSTRACT FROM AUTHOR]

Fraction of hot stars posses strong magnetic fields that channel their radiatively driven outflows. We study the influence of line splitting in the magnetic field (Zeeman effect) on the wind properties. We use our own global wind code with radiative transfer in the comoving frame to understand the influence of the Zeeman splitting on the line force. We show that the Zeeman splitting has a negligible influence on the line force for magnetic fields that are weaker than about 100 kG. This means that the wind mass-loss rates and terminal velocities are not affected by the magnetic line splitting for magnetic fields as are typically found on the surface of nondegenerate stars. Neither have we found any strong flux variability that would be due to the magnetically split line blanketing. [ABSTRACT FROM AUTHOR]

Context. Massive stars live short but intense lives. While less numerous than low-mass stars, they enormously impact their surroundings by several feedback mechanisms. They form in opaque and far-away regions of the galaxy, such that one of these feedback mechanisms also becomes a record of their evolution: their bright large-scale jets and outflows. Aims. In a comprehensive convergence study, we investigate the computational conditions necessary to resolve (pseudo-) disk formation and jet-launching processes, and analyze possible caveats. We explore the magneto-hydrodynamic (MHD) processes of the collapse of massive prestellar cores in detail, including an analysis of the forces involved and their temporal evolution for up to two free-fall times. Methods. We conduct MHD simulations using the state-of-the-art code PLUTO, combining nonideal MHD, self-gravity, and very high resolutions as they have never been achieved before. Our setup includes a 100 M⊙ cloud core that collapses under its own self-gravity to self-consistently form a dense disk structure and launch tightly collimated magneto-centrifugal jets and wide-angle tower flows. Results. We show a comprehensive evolutionary picture of the collapse of a massive prestellar core with a detailed analysis of the physical processes involved and our high-resolution simulations can resolve a magneto-centrifugal jet and a magnetic pressure-driven outflow, separately. The nature of the outflows depends critically on spatial resolution. Only high-resolution simulations are able to differentiate a magneto-centrifugally launched, highly collimated jet from a slow wide-angle magnetic-pressure-driven tower flow. Of these two outflow components, the tower flow dominates angular-momentum transport. The mass outflow rate is dominated by the entrained material from the interaction of the jet with the stellar environment and only part of the ejected medium is directly launched from the accretion disk. A tower flow can only develop to its full extent when much of the original envelope has already dispersed. Taking into account both the mass launched from the surface of the disk and the entrained material from the envelope, we find an ejection-to-accretion efficiency of 10%. Nonideal MHD is required to form centrifugally supported accretion disks and the disk size is strongly dependent on spatial resolution. A converged result for disk and both outflow components requires a spatial resolution of Δx ≤ 0.17 au at 1 au and sink-cell sizes ≤3.1 au. Conclusions. Massive stars not only possess slow wide-angle tower flows, but also produce magneto-centrifugal jets, just as their low-mass counterparts. The actual difference between low-mass and high-mass star formation lies in the "embeddedness" of the high-mass star which implies that the jet and tower flow interact with the infalling large-scale stellar environment, potentially resulting in entrainment. [ABSTRACT FROM AUTHOR]

A small fraction of the radiative flux emitted by hot stars is absorbed by their winds and redistributed towards longer wavelengths. This effect, which leads also to the heating of the stellar photosphere, is termed wind blanketing. For stars with variable winds, the effect of wind blanketing may lead to the photometric variability. We have studied the consequences of line driven wind instability and wind blanketing for the light variability of O stars. We combined the results of wind hydrodynamic simulations and of global wind models to predict the light variability of hot stars due to the wind blanketing and instability. The wind instability causes stochastic light variability with amplitude of the order of tens of millimagnitudes and a typical timescale of the order of hours for spatially coherent wind structure. The amplitude is of the order of millimagnitudes when assuming that the wind consists of large number of independent concentric cones. The variability with such amplitude is observable using present space borne photometers. We show that the simulated light curve is similar to the light curves of O stars obtained using BRITE and CoRoT satellites. [ABSTRACT FROM AUTHOR]

We performed a spectroscopic and photometric analysis to study new eruptions in two luminous blue variables (LBVs) in the Magellanic Clouds. We detected a strong new eruption in the LBV R40 that reached V ~ 9.2 in 2016, which is around 1.3 mag brighter than the minimum registered in 1985. During this new eruption, the star changed from an A-type to a late F-type spectrum. Based on photometric and spectroscopic empirical calibrations and synthetic spectral modeling, we determine that R 40 reached Teff = 5800-6300 K during this new eruption. This object is thereby probably one of the coolest identified LBVs. We could also identify an enrichment of nitrogen and r- and s-process elements. We detected a weak eruption in the LBV R 110 with a maximum of V ~ 9.9 mag in 2011, that is, around 1.0 mag brighter than in the quiescent phase. On the other hand, this new eruption is about 0.2 mag fainter than the first eruption detected in 1990, but the temperature did not decrease below 8500 K. Spitzer spectra show indications of cool dust in the circumstellar environment of both stars, but no hot or warm dust was present, except by the probable presence of PAHs in R 110. We also discuss a possible post-red supergiant nature for both stars. [ABSTRACT FROM AUTHOR]

Context. Evolution of massive stars is affected by a significant loss of mass either via (nearly) spherically symmetric stellar winds or by aspherical mass-loss mechanisms, namely the outflowing equatorial disks. However, the scenario that leads to the formation of a disk or rings of gas and dust around massive stars is still under debate. It is also unclear how various forming physical mechanisms of the circumstellar environment affect its shape and density, as well as its kinematic and thermal structure. Aims. We study the hydrodynamic and thermal structure of optically thick, dense parts of outflowing circumstellar disks that may be formed around various types of critically rotating massive stars, for example, Be stars, B[e] supergiant (sgB[e]) stars or Pop III stars. We calculate self-consistent time-dependent models of temperature and density structure in the disk's inner dense region that is strongly affected by irradiation from a rotationally oblate central star and by viscous heating. Methods. Using the method of short characteristics, we specify the optical depth of the disk along the line-of-sight from stellar poles. Within the optically thick dense region with an optical depth of τ > 2=3 we calculate the vertical disk thermal structure using the diffusion approximation while for the optically thin outer layers we assume a local thermodynamic equilibrium with the impinging stellar irradiation. For time-dependent hydrodynamic modeling, we use two of our own types of hydrodynamic codes: two-dimensional operator-split numerical code based on an explicit Eulerian finite volume scheme on a staggered grid, and unsplit code based on the Roe's method, both including full second-order Navier-Stokes shear viscosity. Results. Our models show the geometric distribution and contribution of viscous heating that begins to dominate in the central part of the disk for mass-loss rates higher than Ṁ& ≳ 10-10 M☉ yr-1. In the models of dense viscous disks with Ṁ > 10-8 M☉ yr-1, the viscosity increases the central temperature up to several tens of thousands of Kelvins, however the temperature rapidly drops with radius and with distance from the disk midplane. The high mass-loss rates and high viscosity lead to instabilities with significant waves or bumps in density and temperature in the very inner disk region. Conclusions. The two-dimensional radial-vertical models of dense outflowing disks including the full Navier-Stokes viscosity terms show very high temperatures that are however limited to only the central disk cores inside the optically thick area, while near the edge of the optically thick region the temperature may be low enough for the existence of neutral hydrogen, for example. [ABSTRACT FROM AUTHOR]

Context. Clumping in the radiation-driven winds of hot, massive stars arises naturally due to the strong, intrinsic instability of linedriving (the line-deshadowing instability, hereafter LDI). But LDI wind models have so far mostly been limited to 1D, mainly because of the severe computational challenges regarding calculation of the multi-dimensional radiation force. Aims. In this paper we simulate and examine the dynamics and multi-dimensional nature of wind structure resulting from the LDI. Methods. We introduce a pseudo-planar, box-in-a-wind method that allows us to efficiently compute the line force in the radial and lateral directions and then use this approach to carry out 2D radiation-hydrodynamical simulations of the time-dependent wind. Results. Our 2D simulations show that the LDI first manifests itself by mimicking the typical shell structure seen in 1D models, but that these shells quickly break up into complex 2D density and velocity structures, characterized by small-scale density "clumps" embedded in larger regions of fast and rarefied gas. Key results of the simulations are that density variations in the well-developed wind are statistically quite isotropic and that characteristic length scales are small; a typical clump size is ℓcl=R* ∼ 0:01 at 2R*, thus also resulting in rather low typical clump masses mcl ∼ 1017 g. Overall, our results agree well with the theoretical expectation that the characteristic scale for LDI generated wind-structure is on the order of the Sobolev length ℓSob.We further confirm some earlier results that lateral "filling in" of radially compressed gas leads to somewhat lower clumping factors in 2D simulations than in comparable 1D models.We conclude by discussing an extension of our method toward rotating LDI wind models that exhibit an intriguing combination of large and small-scale structures extending down to the wind base. [ABSTRACT FROM AUTHOR]

Context. Massive stars severely influence their environment by their strong ionizing radiation and by the momentum and kinetic energy input provided by their stellar winds and supernovae. Quantitative analyses of massive stars are required to understand how their feedback creates and shapes large scale structures of the interstellar medium. The giant Hii region N206 in the Large Magellanic Cloud contains an OB association that powers a superbubble filled with hot X-ray emitting gas, serving as an ideal laboratory in this context. Aims. We aim to estimate stellar and wind parameters of all OB stars in N206 by means of quantitative spectroscopic analyses. In this first paper, we focus on the nine Of-type stars located in this region. We determine their ionizing flux and wind mechanical energy. The analysis of nitrogen abundances in our sample probes rotational mixing. Methods. We obtained optical spectra with the multi-object spectrograph FLAMES at the ESO-VLT. When possible, the optical spectroscopy was complemented by UV spectra from the HST, IUE, and FUSE archives. Detailed spectral classifications are presented for our sample Of-type stars. For the quantitative spectroscopic analysis we used the Potsdam Wolf-Rayet model atmosphere code. We determined the physical parameters and nitrogen abundances of our sample stars by fitting synthetic spectra to the observations. Results. The stellar and wind parameters of nine Of-type stars, which are largely derived from spectral analysis are used to construct wind momentum - luminosity relationship. We find that our sample follows a relation close to the theoretical prediction, assuming clumped winds. The most massive star in the N206 association is an Of supergiant that has a very high mass-loss rate. Two objects in our sample reveal composite spectra, showing that the Of primaries have companions of late O subtype. All stars in our sample have an evolutionary age of less than 4 million yr, with the O2-type star being the youngest. All these stars show a systematic discrepancy between evolutionary and spectroscopic masses. All stars in our sample are nitrogen enriched. Nitrogen enrichment shows a clear correlation with increasing projected rotational velocities. Conclusions. The mechanical energy input from the Of stars alone is comparable to the energy stored in the N206 superbubble as measured from the observed X-ray and Hα emission. [ABSTRACT FROM AUTHOR]

Context. Through spectrally unresolved observations of high-J CO transitions, Herschel Photodetector Array Camera and Spectrometer (PACS) has revealed large reservoirs of warm (300 K) and hot (700 K) molecular gas around low-mass protostars. The excitation and physical origin of this gas is still not understood. Aims. We aim to shed light on the excitation and origin of the CO ladder observed toward protostars, and on the water abundance in different physical components within protostellar systems using spectrally resolved Herschel-HIFI data. Methods. Observations are presented of the highly excited CO line J = 16-15 (Eup/kB = 750 K) with the Herschel Heterodyne Instrument for the Far Infrared (HIFI) toward a sample of 24 low-mass protostellar objects. The sources were selected from the Herschel "Water in Star-forming regions with Herschel" (WISH) and "Dust, Ice, and Gas in Time" (DIGIT) key programs. Results. The spectrally resolved line profiles typically show two distinct velocity components: a broad Gaussian component with an average FWHM of 20 km s-1 containing the bulk of the flux, and a narrower Gaussian component with a FWHM of 5 km s-1 that is often offset from the source velocity. Some sources show other velocity components such as extremely-high-velocity features or "bullets". All these velocity components were first detected in H2O line profiles. The average rotational temperature over the entire profile, as measured from comparison between CO J = 16-15 and 10-9 emission, is ~300 K. A radiative-transfer analysis shows that the average H2O/CO column-density ratio is ~0.02, suggesting a total H2O abundance of ~2 × 10-6, independent of velocity. Conclusions. Two distinct velocity profiles observed in the HIFI line profiles suggest that the high-J CO ladder observed with PACS consists of two excitation components. The warm PACS component (300 K) is associated with the broad HIFI component, and the hot PACS component (700 K) is associated with the offset HIFI component. The former originates in either outflow cavity shocks or the disk wind, and the latter in irradiated shocks. The low water abundance can be explained by photodissociation. The ubiquity of the warm and hot CO components suggest that fundamental mechanisms govern the excitation of these components; we hypothesize that the warm component arises when H2 stops being the dominant coolant. In this scenario, the hot component arises in cooling molecular H2-poor gas just prior to the onset of H2 formation. High spectral resolution observations of highly excited CO transitions uniquely shed light on the origin of warm and hot gas in low-mass protostellar objects. [ABSTRACT FROM AUTHOR]

Radio observations of the region surrounding PSR J1301-6305 at 5.5 GHz and 7.5 GHz were conducted with ATCA on September 5, 2013. The observations were dedicated to the search of the radio counterpart of the evolved pulsar wind nebula (PWN) HESS J1303-631, which has been detected in X-rays and GeV-TeV γ-rays. The collected data do not reveal any significant extended emission associated with PSR J1301-6305. In addition, archival 1.384 GHz and 2.368 GHz data do not show any evidence for a radio counterpart of HESS J1303-631. Archival 1.384 GHz observations reveal the detection of an extended structure centred at an angular distance of 19' from the pulsar. This extended structure might be a supernova remnant (SNR) and a potential birth place of PSR J1301-6305. The implications of the lack of a radio counterpart of HESS J1303-631 on the understanding of the nature of the PWN are discussed. [ABSTRACT FROM AUTHOR]