An overview is provided of the scientific goals of the Magellanic Cloud component of the STScI Directors Discretionary UV initiative ULLYSES, together with the complementary spectroscopic survey XShootU (VLT/Xshooter) and other ancillary datasets. Together, ULLYSES and XShootU permit the first comprehensive, homogeneous study of wind densities and velocities in metal-poor massive stars, plus UV/optical spectroscopic libraries for population synthesis models and a large number of interstellar sight-lines towards the Magellanic Clouds.

Expanding nebulae are produced by mass loss from stars, especially during late stages of evolution. We describe the algorithms and methods implemented in the radiation-magnetohydrodynamics (MHD) code PION for highly scalable simulations using static mesh-refinement. We present results from 3D MHD simulations of bow shocks around runaway massive stars, and of the expansion of a fast wind from a Wolf-Rayet star into the slow wind from a previous red supergiant phase of evolution. PION is free software that can be downloaded from https://www.pion.ie/

Massive stars drive strong winds that impact the surrounding interstellar medium, producing parsec-scale bubbles for isolated stars and superbubbles around young clusters. These bubbles can be observed across the electromagnetic spectrum, both the wind itself and the swept up interstellar gas. Runaway massive stars produce bow shocks that strongly compresses interstellar gas, producing bright infrared, optical and radio nebulae. With the detection of non-thermal radio emission from bow shocks, particle acceleration can now also be investigated. I review research on wind bubbles and bow shocks around massive stars, highlighting recent advances in infrared, radio and X-ray observations, and progress in multidimensional simulations of these nebulae. These advances enable quantitative comparisons between theory and observations and allow to test the importance of some physical processes such as thermal conduction and Kelvin-Helmholtz instability in shaping nebulae and in constraining the energetics of stellar-wind feedback to the interstellar medium.

We run numerical simulations of massive colliding wind binaries, and quantify the accretion onto the secondary under different conditions. We set 3D simulation of a LBV–WR system and vary the LBV mass loss rate to obtain different values of wind momentum ratio η. We show that the mean accretion rate for stationary systems fits a power law Macc∝ η–1.6 for a wide range of η, until for extremely small η saturation in the accretion is reached. We find that the stronger the primary wind, the smaller the opening angle of the colliding wind structure (CWS), and compare it with previous analytical estimates. We demonstrate the efficiency of clumpy wind in penetrating the CWS and inducing smaller scale clumps that can be accreted. We propose that simulations of colliding winds can reveal more relations as the ones we found, and can be used to constrain stellar parameters.

With the use of high-resolution ALMA observations, complex structures that resemble those observed in post-AGB stars and planetary nebulae are detected in the circumstellar envelopes of low-mass evolved stars. These deviations from spherical symmetry are believed to be caused primarily by the interaction with a companion star or planet. With the use of three-dimensional hydrodynamic simulations, we study the impact of a binary companion on the wind morphology and dynamics of an AGB outflow. We classifiy the wind structures and morphology that form in these simulations with the use of a classification parameter, constructed with characteristic parameters of the binary configuration. Finally we conclude that the companion alters the wind expansion velocity through the slingshot mechanism, if it is massive enough.

Low- and intermediate mass stars experience a significant mass loss during the last phases of their evolution, which obscures them in a vast, dusty envelope. Although it has long been thought this envelope is generally spherically symmetric in shape, recent high-resolution observations find that most of these stars exhibit complex and asymmetrical morphologies, most likely resulting from binary interaction. In order to improve our understanding about these systems, theoretical studies are needed in the form of numerical simulations. Currently, a handful of simulations exist, albeit they mainly focus on the hydrodynamics of the outflow. Hence, we here present the pathway to more detailed and accurate modelling of companion-perturbed outflows with, by discussing the missing but crucial physical and chemical processes. With these state-of-the-art simulations we aim to make a direct comparison with observations to unveil the true identity on the embedded systems.

The observations and simulations have revealed that large-scale magnetic field and outflows can exist in the inner regions of an advection-dominated accretion disc where the resistive diffusion may also be important. In the present paper, the roles of large-scale magnetic field and outflows in the structure of resistive advection-dominated accretion discs are explored by assuming that the accretion flow is radially self-similar. In the non-ideal magnetohydrodynamic (MHD) approximation, the results show that the angular velocity is always sub-Keplerian when both the outflow and the large-scale magnetic field are taken into account. A stronger toroidal field component leads to faster rotation, while the disc rotates with faster rate if the vertical field component is weaker. The increase of magnetic diffusivity causes the infall velocity to be close to Keplerian velocity. Although the previous studies in the ideal MHD approximation have shown that the disc temperature decreases due to the vertical field component, we find that the effect of vertical field component on the temperature of a resistive disc depends on the magnetic diffusivity. When the magnetic diffusivity is high, the more efficient mechanism for decreasing the disc temperature can be the outflows, and not the large-scale magnetic field. In such a limit of the magnetic diffusivity, the components of the large-scale magnetic field enhance the gas temperature. The increase of temperature can lead to heating and acceleration of the electrons and help us to explain the origin of phenomena such as the flares in Sgr A*. On the other hand, the infall velocity in such a limit rises as the temperature increases, and therefore the surface density falls to too low values. Any change in the density profile can alter the structure and the emitted spectrum of disc.

The hot Jupiter HD189733b is expected to be a source of strong radio emission, due to its close proximity to its magnetically active host star. Here, we model the stellar wind of its host star, based on reconstructed surface stellar magnetic field maps. We use the local stellar wind properties at the planetary orbit obtained from our models to compute the expected radio emission from the planet. Our findings show that the planet emits with a peak flux density within the detection capabilities of LOFAR. However, due to absorption by the stellar wind itself, this emission may be attenuated significantly. We show that the best time to observe the system is when the planet is near primary transit of the host star, as the attenuation from the stellar wind is lowest in this region.

The stellar magnetic field completely dominates the environment around late-type stars. It is responsible for driving the coronal high-energy radiation (e.g. EUV/X-rays), the development of stellar winds, and the generation transient events such as flares and coronal mass ejections (CMEs). While progress has been made for the first two processes, our understanding of the eruptive behavior in late-type stars is still very limited. One example of this is the fact that despite the frequent and highly energetic flaring observed in active stars, direct evidence for stellar CMEs is almost non-existent. Here we discuss realistic 3D simulations of stellar CMEs, analyzing their resulting properties in contrast with solar eruptions, and use them to provide a common framework to interpret the available stellar observations. Additionally, we present results from the first 3D CME simulations in M-dwarf stars, with emphasis on possible observable signatures imprinted in the stellar corona.

Ionized winds from late-type main-sequence stars are important for stellar spin-down and therefore the evolution of stellar activity; winds blow an “astrosphere” into the interstellar medium that absorbs a large part of galactic cosmic rays; and the winds play a key role in shaping planetary environments, in particular their upper atmospheres. These issues have been well studied for the solar wind but little is known about winds escaping from other solar-type stars. Several methods have been devised to either detect winds directly or to infer the presence of such winds from features that are shaped by the winds. This paper summarizes these methods and discusses exemplary findings. There is need for more studies using multiple methods for the same stars.

Stars and their exoplanets evolve together. Depending on the physical characteristics of these systems, such as age, orbital distance and activity of the host stars, certain types of star-exoplanet interactions can dominate during given phases of the evolution. Identifying observable signatures of such interactions can provide additional avenues for characterising exoplanetary systems. Here, I review some recent works on star-planet interactions and discuss their observability at different wavelengths across the electromagnetic spectrum.

The association of 6.7 GHz class II methanol (CH3OH) masers with ATLASGAL/ ALMA 0.9 mm massive dense cores is presented in this work from a statistical viewpoint. 42 of the 112 cores (37.5%) detected with the Atacama Compact Array (ACA) excite 6.7 GHz CH3OH masers. ACA cores have offsets 0\rlap.″17 to 4\rlap.″79 from the methanol multibeam survey (MMB), with a median of 2.″19. Approximately 90% of the MMB-associated cores are of masses > 40 M⊙. Because all the cores show evidence of outflow activity, and only a fraction of the cores excited CH3OH masers, we suggest that outflows precede the emergence of maser emission. This first ALMA survey of massive dense cores combined with the MMB survey along with other maser specie surveys is a promising tool to trace the evolutionary sequence of high-mass stars.

Massive star winds greatly influence the evolution of both their host star and local environment though their mass-loss rates, but current radiative line-driven wind models do not incorporate any magnetic effects. Recent surveys of O and B stars have found that about ten percent have large-scale, organized magnetic fields. These massive-star magnetic fields, which are thousands of times stronger than the Sun’s, affect the inherent properties of their own winds by changing the mass-loss rate. To quantify this, we present a simple surface mass-flux scaling over the stellar surface which can be easily integrated to get an estimate of the mass-loss rate for a magnetic massive star. The overall mass-loss rate is found to decrease by factors of 2-5 relative to the non-magnetic CAK mass-loss rate.

Gamma Cassiopeiae is an enigmatic Be star with unusually hard, strong X-ray emission compared with normal main-sequence B stars. The origin has been debated for decades between two theories: mass accretion onto a hidden compact companion and a magnetic dynamo driven by the star-Be disk differential rotation. There has been no decisive signature found that supports either theory, such as a pulse in X-ray emission or the presence of large-scale magnetic field. In a ~100 ksec duration observation of the star with the Suzaku X-ray observatory in 2011, we detected six rapid X-ray spectral hardening events called “softness dips”. All the softness dip events show symmetric softness ratio variations, and some of them have flat bottoms apparently due to saturation. The softness dip spectra are best described by either ~40% or ~70% partial covering absorption to kT ~12 keV plasma emission by matter with a neutral hydrogen column density of ~2 − 8 × 1021cm−2, while the spectrum outside of these dips is almost free of absorption. This result suggests that two distinct X-ray emitting spots in the γ Cas system, perhaps on a white dwarf companion with dipole mass accretion, are occulted by blobs in the Be stellar wind, the Be disk, or rotating around the white dwarf companion. The formation of a Be star and white dwarf binary system requires mass transfer between two stars; γ Cas may have experienced such activity in the past.

In the frame of radiation driven wind theory (Castor et al.1975), we present self-consistent hydrodynamical solutions to the line-force parameters (k, α, δ) under LTE conditions. Hydrodynamic models are provided by HydWind (Curé 2004). We evaluate these results with those ones previously found in literature, focusing in different regions of the optical depth to be used to perform the calculations. The values for mass-loss rate and terminal velocity obtained from our calculations are also presented.

We also examine the line-force parameters for the case when large changes in ionization throughout the wind occurs (δ-slow solutions, Curé et al.2011).

Recent work (Corradi et al. 2015; Jones et al. 2016) has shown that the phenomenon of extreme abundance discrepancies, where recombination line abundances exceed collisionally excited line abundances by factors of 10 or more, seem to be strongly associated with planetary nebulae with close binary central stars. To further investigate, we have obtained spectra of a sample of nebulae with known close binary central stars, using FORS2 on the VLT, and we have discovered several new extreme abundance discrepancy objects. We did not find any non-extreme discrepancies, suggesting that a very high fraction of nebulae with close binary central stars also have an extreme abundance discrepancy.

We derive an updated, analytic formula for the gravitational radius, rg, of an irradiated, thin disk. For such a disk, the region outside the gravitational radius will produce a thermal wind, while the region interior to the gravitational radius will be stable. We find that rg ≈ 1.4 (M*/M⊙)/(T0/104K) AU. The value of 1.4 AU is approximately one fifth of currently used values.

The analysis uses the adiabatic approximation. We argue that the same formula applies for the non-adiabatic, isothermal case.

We have performed large-scale nucleosynthesis calculations within the high-entropy-wind (HEW) scenario of Type II supernovae. The primary aim was to constrain the conditions for the production of the classical ‘p-only’ isotopes of the light trans-Fe elements. We find, however, that for electron fractions in the range 0.458 ≤ Ye ≤ 0.478, sizeable abundances of p-, s- and r-process nuclei between 64Zn and 98Ru are coproduced in the HEW at low entropies (S ≤ 100) by a primary charged-particle process after an α-rich freezeout. With the above Ye–S correlation, most of the predicted isotopic abundance ratios within a given element, e.g. 64Zn(p)/70Zn(r) or 92Mo(p)/94Mo(p), as well as of neighboring elements, e.g. 70Ge(s + p)/74Se(p) or 74Se(p)/78Kr(p) agree with the observed Solar-System ratios. Taking the Mo isotopic chain as a particularly challenging example, we show that our HEW model can account for the production of all 7 stable isotopes, from ‘p-only’ 92Mo, via ‘s-only’ 96Mo up to ‘r-only’ 100Mo. Furthermore, our model is able to reproduce the isotopic composition of Mo in presolar SiC X-grains.

Magnetohydrodynamic devices that can accelerate plasmas to speeds of the order of hundreds of kilometres per second have been designed and built for nearly forty years. Up to the time of writing, however, the theory for such devices has been exclusively non-relativistic. In this paper we derive the special relativistic magnetohydrodynamic (SRMHD) equations and use them to obtain the relativistic, magnetic nozzle equation which describes the production of jet flows with speeds approaching the speed of light.We obtain analytic solutions to this equation and show that, in principle, magnetic field gradients can accelerate a plasma to highly relativistic speeds. We also show that the exit kinetic energy, EK, of a particle is given by the equation EK = m0C2FR, where m0 is the rest mass of the particle and CFR is the fast magnetosonic speed at the start of the flow.

The relativistic nozzle differs in a number of ways from the non-relativistic case.A non-relativistic nozzle has a relatively symmetric converging/diverging shape, while a highly relativistic nozzle converges in the usual manner, but diverges, in an abrupt fashion, at the very end of the nozzle. The gentle divergence of non-relativistic nozzles causes the exit plasma densities and magnetic fields of the flow to have values that are small relative to their values at the start of the nozzle. The abrupt divergence of a highly relativistic nozzle implies that, for a less than perfect nozzle, the exit values of the mass density and the magnetic field strength are comparable to their initial values. This unexpected dichotomy in behaviour may have future application in understanding the ‘radio-loud’ and ‘radio-quiet’ relativistic jets that are produced from astrophysical sources.

Optical and X-ray spectroscopy indicate that the X-ray pulsar GX 1+4 is seen through a cloud of gravitationally bound matter. We discuss an unstable negative feedback mechanism (originally proposed by Kotani et al. 1999), based on X-ray heating of this matter which controls the accretion rate when the source is in a low X-ray luminosity state. A deep minimum lasting ∼6 hours occurred during observations with the RXTE satellite over 1996 July 19–21. The shape of the X-ray pulses changed remarkably from before to after the minimum. These changes may be related to the transition from neutron star spin-down to spin-up which occurred at about the same time. Smoothed particle hydrodynamic simulations of the effect of adding matter with opposite angular momentum to an existing disk, show that it is possible for a number of concentric rings with alternating senses of rotation to co-exist in a disk. This could provide an explanation for the step-like changes in Ṗ which are observed in GX 1+4. Changes at the inner boundary of the disk occur at the same timescale as that imposed at the outer boundary. Reversals of material torque on the neutron star occur at a minimum in LX.