Spherical Wind Research Articles

The symbiotic recurrent nova V745 Sco at radio wavelengths

Abstract V745 Sco is a Galactic symbiotic recurrent nova with nova eruptions in 1937, 1989 and 2014. We study the behaviour of V745 Sco at radio wavelengths (0.6–37 GHz), covering both its 1989 and 2014 eruptions and informed by optical, X-ray, and γ-ray data. The radio light curves are synchrotron-dominated. Surprisingly, compared to expectations for synchrotron emission from explosive transients such as radio supernovae, the light curves spanning 0.6–37 GHz all peak around the same time (∼18–26 days after eruption) and with similar flux densities (5–9 mJy). We model the synchrotron light curves as interaction of the nova ejecta with the red giant wind, but find that simple spherically symmetric models with wind-like circumstellar material (CSM) cannot explain the radio light curve. Instead, we conclude that the shock suddenly breaks out of a dense CSM absorbing screen around 20 days after eruption, and then expands into a relatively low density wind ($\dot{M}_{out} \approx 10^{-9}-10^{-8}$ M⊙ yr−1 for vw = 10 km s−1) out to ∼1 year post-eruption. The dense, close-in CSM may be an equatorial density enhancement or a more spherical red giant wind with $\dot{M}_{in} \approx [5-10] \times 10^{-7}$ M⊙ yr−1, truncated beyond several × 1014 cm. The outer lower-density CSM would not be visible in typical radio observations of Type Ia supernovae: V745 Sco cannot be ruled out as a Type Ia progenitor based on CSM constraints alone. Complementary constraints from the free–free radio optical depth and the synchrotron luminosity imply the shock is efficient at accelerating relativistic electrons and amplifying magnetic fields, with εe and εB ≈ 0.01 − 0.1.

The Optically Thick Rotating Magnetic Wind from a Massive White Dwarf Merger Product. II. Axisymmetric Magnetohydrodynamic Simulations

We numerically construct a series of axisymmetric rotating magnetic wind solutions, aiming at exploring the observation properties of massive white dwarf (WD) merger remnants with a strong magnetic field, a fast spin, and an intense mass loss, as inferred for WD J005311. We investigate the magnetospheric structure and the resultant spin-down torque exerted to the merger remnant with respect to the surface magnetic flux Φ*, spin angular frequency Ω* and the mass-loss rate Ṁ . We confirm that the wind properties for σ≡Φ*2Ω*2/Ṁvesc3≳1 significantly deviate from those of the spherical Parker wind, where v esc is the escape velocity at stellar surface. For such a rotating magnetic wind sequence, we find (i) a quasiperiodic mass eruption triggered by magnetic reconnection along with the equatorial plane and (ii) a scaling relation for the spin-down torque T≈(1/2)×ṀΩ*R*2σ1/4 . We apply our results to discuss the spin-down evolution and wind anisotropy of massive WD merger remnants, the latter of which could be probed by a successive observation of WD J005311 using Chandra.

A Gaussian-processes approach to fitting for time-variable spherical solar wind in pulsar timing data

ABSTRACT Propagation effects are one of the main sources of noise in high-precision pulsar timing. For pulsars below an ecliptic latitude of 5°, the ionized plasma in the solar wind can introduce dispersive delays of order $100\, \mu \mathrm{s}$ around solar conjunction at an observing frequency of 300 MHz. A common approach to mitigate this assumes a spherical solar wind with a time-constant amplitude. However, this has been shown to be insufficient to describe the solar wind. We present a linear, Gaussian-process piecewise Bayesian approach to fit a spherical solar wind of time-variable amplitude, which has been implemented in the pulsar software run_enterprise. Through simulations, we find that the current EPTA+InPTA data combination is not sensitive to such variations; however, solar wind variations will become important in the near future with the addition of new InPTA data and data collected with the low-frequency LOFAR telescope. We also compare our results for different high-precision timing data sets (EPTA+InPTA, PPTA, and LOFAR) of 3 ms pulsars (J0030+0451, J1022+1001, J2145−0450), and find that the solar-wind amplitudes are generally consistent for any individual pulsar, but they can vary from pulsar to pulsar. Finally, we compare our results with those of an independent method on the same LOFAR data of the three millisecond pulsars. We find that differences between the results of the two methods can be mainly attributed to the modelling of dispersion variations in the interstellar medium, rather than the solar wind modelling.

Damage mechanism analysis of concrete-filled steel tube (CFST) spherical joint wind turbine plane tower

In order to investigate the failure modes and damage mechanism of the concrete filled steel tubular latticed planar wind turbine tower with split spherical joints under earthquakes, the low cyclic loading tests were carried out on the two plane towers with web wall thickness as parameters, and their hysteresis curves and skeleton curves were analyzed. At the same time, it is verified via ABAQUS finite element simulation, and based on this analysis, the parameter expansion analysis of the web diameter ratio was carried out. The results indicate that the failure characteristics of the tower include three types: weld tear failure, web member buckling failure and web member high-strength bolt pull-out failure; to increase the thickness of the web member can make the hysteresis curve of the tower fuller and the energy dissipation capacity stronger effectively. The ultimate bearing capacity and ductility coefficient of the tower are positively correlated with the diameter ratio of web members. When the diameter ratio of web members was greater than 0.13, the increase in ultimate bearing capacity and ductility coefficient reduced obviously. In order to maximize the overall mechanical performance of the tower and improve the economy, in practical engineering applications, the diameter ratio of the web member of such towers should be taken as the range of 0.11~0.13.

Modelling variable linear polarization produced by Co-rotating Interaction Regions (CIRs) across optical recombination lines of Wolf–Rayet stars

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

The energy and dynamics of trapped radiative feedback with stellar winds

ABSTRACT In this paper, we explore the significant, non-linear impact that stellar winds have on H ii regions. We perform a parameter study using three-dimensional radiative magnetohydrodynamic simulations of wind and ultraviolet radiation feedback from a 35 M⊙ star formed self-consistently in a turbulent, self-gravitating cloud, similar to the Orion Nebula (M42) and its main ionizing source θ1 Ori C. Stellar winds suppress early radiative feedback by trapping ionizing radiation in the shell around the wind bubble. Rapid breakouts of warm photoionized gas (‘champagne flows’) still occur if the star forms close to the edge of the cloud. The impact of wind bubbles can be enhanced if we detect and remove numerical overcooling caused by shocks crossing grid cells. However, the majority of the energy in the wind bubble is still lost to turbulent mixing between the wind bubble and the gas around it. These results begin to converge if the spatial resolution at the wind bubble interface is increased by refining the grid on pressure gradients. Wind bubbles form a thin chimney close to the star, which then expands outwards as an extended plume once the wind bubble breaks out of the dense core the star formed in, allowing them to expand faster than a spherical wind bubble. We also find wind bubbles mixing completely with the photoionized gas when the H ii region breaks out of the cloud as a champagne flow, a process we term ‘hot champagne’.

Variations in the intermediate wind region of the blue supergiant 55 Cygni

Context. The quantitative near-infrared (NIR) spectroscopic synthesis is an important technique for determining wind properties of massive stars. The Brα line is an excellent mass-loss tracer and provides valuable information on the physical conditions of intermediate-wind regions. The knowledge of the wind properties gained by studying the NIR lines could provide extra ingredients to the theory of line-driven winds, mainly because the standard theory does not predict observed properties of blue supergiants, such as high values for the β parameter (β> 2), low terminal velocities, and mass-loss variability. Aims. We seek to enhance our understanding of the wind properties of B supergiants. To this end, we propose analysing their NIR spectra over different epochs to study wind variability and its connection with phenomena arising from regions close to the photosphere. Methods. We present the first sets of multi-epoch high-resolution K- and L-band spectra of 55 Cyg acquired with the Gemini Near-InfraRed Spectrograph (GNIRS). We measured line equivalent widths and modelled the Brα line to derive (unclumped) mass-loss rates. Synthetic line profiles were computed for a homogeneous spherical wind by solving the radiative transfer equations in the co-moving frame for a multi-level atom in non-local thermodynamic equilibrium (NLTE). Results. We observe variations in the spectral lines originating in the upper photosphere and the wind. The perturbations, on average, have periods of ~13 and ~23 days; the latter is similar to that found previously from optical data (22.5 days). The NIR lines observed in 2013 are described with the same wind structure used to model a quasi-simultaneous observation in Hα. By contrast, from observations taken in 2015, we derived a higher mean mass-loss rate. Variations in the mass-loss rate are also detected within a few weeks. Interestingly, we find that the profile shape of the Hu14 line sets constraints on the mass loss. Moreover, we find the Mg II doublet in emission, which suggests a tenuous circumstellar gas ring or shell. Conclusions. The variability detected in the NIR H emission lines of 55 Cyg is related to changes in the mass-loss rate, which doubled its value between 2013 and 2015. Furthermore, the short-term variability (within three weeks) in the spectral lines and mass loss supports the hypothesis of strange-mode oscillations. This pilot project demonstrates the importance of comprehensive monitoring of blue supergiants’ variability to deeply understand the physical properties of their stellar winds and the role of pulsations in recurrently enhancing mass loss.

The dusty circumstellar environment of Betelgeuse during the Great Dimming as seen by VLTI/MATISSE

Context. The ‘Great Dimming’ of the prototypical red supergiant Betelgeuse, which occurred between December 2019 and April 2020, gives us unprecedented insight into the processes occurring on the stellar surface and in the inner wind of this type of star. In particular it may bring further understanding of their dust nucleation and mass-loss processes. Aims. Here, we present and analyse VLTI/MATISSE observations in the N band (8–13 µm) taken near the brightness minimum in order to assess the status of the dusty circumstellar environment. Methods. We explored the compatibility of a dust clump obscuring the star with our mid-infrared interferometric observations using continuum 3D radiative transfer modelling, and probed the effect of adding multiple clumps close to the star on the observables. We also tested the viability of a large cool spot on the stellar surface without dust present in the ambient medium. Results. Using the visibility data, we derived a uniform disk diameter of 59.02 ± 0.64 mas in the spectral range 8–8.75 µm. We find that both the dust clump and the cool spot models are compatible with the data. Further to this, we note that the extinction and emission of our localised dust clump in the line of sight of the star directly compensate for each other, making the clump undetectable in the spectral energy distribution and visibilities. The lack of infrared brightening during the Great Dimming therefore does not exclude extinction due to a dust clump as one of the possible mechanisms. The visibilities can be reproduced by a spherical wind with dust condensing at 13 stellar radii and a dust mass-loss rate of (2.1–4.9) × 10−10 M⊙ yr−1; however, in order to reproduce the complexity of the observed closure phases, additional surface features or dust clumps would be needed.

Faint but Not Forgotten. I. First Results from a Search for Astrospheres around AGB Stars in the Far-ultraviolet

Using the GALEX archive, we have discovered extended structures around ten asymptotic giant branch (AGB) stars (out of a total 92 searched) emitting in the far-ultraviolet (FUV) band. In all but one, we find the typical morphology expected for a spherical wind moving relative to, and interacting with, the interstellar medium (ISM) to produce an astrosphere. The exception is V Hya whose mass ejection is known to be highly aspherical, where we find evidence of its large parabolic outflows interacting with the ISM, and its collimated, extreme velocity outflows interacting with the circumstellar medium. For eight objects with relatively large proper motions, we find (as expected) that the termination-shock region lies in a hemisphere that contains the proper motion vector. Radial intensity cuts for each source have been used to locate the termination shock and the astropause’s outer edge. In a few objects, the cuts also reveal faint emission just outside the astropause that likely arises in shocked ISM material. We have used these data, together with published mass-loss rates and wind expansion velocities, to determine the total mass lost and duration for each source—we find that the duration of and total mass in the shocked wind are significantly larger than their corresponding values for the unshocked wind. The combination of FUV and far-IR data on AGB astrospheres provides a unique database for theoretical studies (numerical simulations) of wind–ISM interactions. We show that a Cyclical Spatial Heterodyne Spectrometer on a small space-based telescope can provide high-resolution spectra of astrospheres to confirm the emission mechanism.

Probing the tidal disruption event iPTF16axa with cloudy and disc-wind models

ABSTRACT We present both a disc-wind model on the optical/ultraviolet (UV) emission continuum and cloudy modelling on the spectral lines of the tidal disruption event (TDE) intermediate Palomar Transient Factory (iPTF)16axa to understand the disc-wind emission and the properties of the atmosphere that impacts the line luminosity of the TDE. Assuming the optical/UV emission from the wind due to the disc super-Eddington phase, we use the steady structured disc-wind model with a spherical wind with constant velocity to fit the observations on multiple days. The extracted parameters are stellar-mass M⋆ = 6.20 ± 1.19M⊙, disc radiative efficiency log10(η) = −1.22 ± 1.327, wind inner radius rl = (2.013 ± 0.551) × 1014 cm, and velocity vw = 18999.4 ± 1785.1 km s−1. The photosphere temperature for wind emission is ∼2 × 104 K and the disc single blackbody temperature is ∼0.995 × 105 K. We also perform cloudy modelling to explain the observed He and H line luminosities that estimate a wind inner radius rl = 7.07 × 1014 cm and velocity vw = 1.3 × 104 km s−1. The independent analyses of iPTF16axa using cloudy and disc-wind models show comparable results that agree with observations. The cloudy modelling finds that both the super-solar abundance of He and a smaller He II line optical depth is responsible for the enhancement of He II line luminosity over the Hα line luminosity. The super-solar abundance of He II agrees with a relatively large stellar mass and suggests that the disrupted star might have been a red giant.

Three scenarios for shaping ‘ears’ in planetary nebulae

We present three scenarios for shaping planetary nebulae (PNe) with ears : (i) Ejection of early and short-duration jets into a spherical wind. (ii) Fast wind ejection into a slow dense wind and then jet ejection into the hot bubble. (iii) Interaction of fast wind with two opposite punctured regions. We find that all three scenarios can lead to the formation of ears. Ears are two opposite polar protrusions from the main nebula that are smaller than the main nebula and with a cross section that decreases monotonically from the base of each ear at the shell to its far end. Ears are generally indicative of the operation of jets, but in a rich spectrum of flow structures, jet properties and time and duration of activity.

Radio emission from simulated tidal disruption events

ABSTRACT Several tidal disruption events such as ASASSN-14li and XMMSL1 J0740-85 have recently been observed in the radio. While the radio emission of some tidal disruption events are attributed to a relativistic jet, a few others are associated with a non-relativistic outflow. This outflow can either be due to a spherical wind or unbound tidal debris. We explore this latter hypothesis in this paper. We show that the maximum velocity of the unbound debris is a function of the impact parameter, such that smaller impact parameters (closer approaches) produce larger maximum velocities. We then model this outflow which expands and shocks the local interstellar medium and compute the peak radio flux and frequency as functions of the impact parameter. Moreover, multiple epochs of observations can put additional constraints on the profile of the local interstellar medium. We apply this analysis to four tidal disruption events whose radio emission is attributed to a non-relativistic outflow and show that the velocities of the unbound material are consistent with our simulated events. We also place constraints on the density profile of three of the four tidal disruption events with multiple epochs of observations.

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars

This overview paper presentsATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter/submillimeter Array (ALMA). The goal ofATOMIUMis to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. Seventeen oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (∼0″​​.025–0″​​.050), medium (∼0″​​.13–0″​​.24), and low (∼1″) angular resolution. The sensitivity per ∼1.3 km s−1channel was 1.5–5 mJy beam−1, and the line-free channels were used to image the millimetre wave continuum. Our primary molecules for studying the gas dynamics and dust formation are CO, SiO, AlO, AlOH, TiO, TiO2, and HCN; secondary molecules include SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described. In addition, we highlight one scientific result – the wind kinematics of theATOMIUMsources. Our analysis suggests that theATOMIUMsources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion. Moreover, the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In fifteen sources, some molecular transitions other than12COv = 0J = 2 − 1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. We propose that a binary interaction with a (sub)stellar companion may (partly) explain the non-monotonic behaviour of the projected velocity field. TheATOMIUMdata hence provide a crucial benchmark for the wind dynamics of evolved stars in single and binary star models.

Stellar Wind Confinement of Evaporating Exoplanet Atmospheres and Its Signatures in 1083 nm Observations

Abstract Atmospheric escape from close-in exoplanets is thought to be crucial in shaping observed planetary populations. Recently, significant progress has been made in observing this process in action through excess absorption in-transit spectra and narrowband light curves. We model the escape of initially homogeneous planetary winds interacting with a stellar wind. The ram pressure balance of the two winds governs this interaction. When the impingement of the stellar wind on the planetary outflow is mild or moderate, the planetary outflow expands nearly spherically through its sonic surface before forming a shocked boundary layer. When the confinement is strong, the planetary outflow is redirected into a cometary tail before it expands to its sonic radius. The resultant transmission spectra at the He 1083 nm line are accurately represented by a 1D spherical wind solution in cases of mild to moderate stellar wind interaction. In cases of strong stellar wind interaction, the degree of absorption is enhanced and the cometary tail leads to an extended egress from transit. The crucial features of the wind–wind interaction are, therefore, encapsulated in the light curve of He 1083 nm equivalent width as a function of time. The possibility of extended He 1083 nm absorption well beyond the optical transit carries important implications for planning out-of-transit observations that serve as a baseline for in-transit data.

Multiple components in the molecular outflow of the red supergiant NML Cyg

ABSTRACT Despite their large impact on stellar and galactic evolution, the properties of outflows from red supergiants are not well characterized. We used the Onsala 20m telescope to perform a spectral survey at 3 and 4 mm (68–116 GHz) of the red supergiant NML Cyg, alongside the yellow hypergiant IRC + 10420. Our observations of NML Cyg were combined with complementary archival data to enable a search for signatures of morphological complexity in the circumstellar environment, using emission lines from 15 molecular species. The recovered parameters imply the presence of three distinct, coherent, and persistent components, comprised of blue-shifted and red-shifted components, in addition to an underlying outflow centred at the stellar systemic velocity. Furthermore, to reproduce 12CO emission with 3D radiative transfer models required a spherical outflow with three superposed conical outflows, one towards and one away from the observer, and one in the plane of the sky. These components are higher in density than the spherical outflow by up to an order of magnitude. We hence propose that NML Cyg’s circumstellar environment consists of a small number of high-density large-scale coherent outflows embedded in a spherical wind. This would make the mass-loss history similar to that of VY CMa, and distinct from μ Cep, where the outflow contains many randomly distributed smaller clumps. A possible correlation between stellar properties, outflow structures, and content is critical in understanding the evolution of massive stars and their environmental impact.

Radio modelling of the brightest and most luminous non-thermal colliding-wind binary Apep

ABSTRACT Apep is the brightest and most luminous non-thermal colliding-wind binary by over an order of magnitude. It has been suggested from infrared observations that one of the Wolf–Rayet stars in Apep is launching an anisotropic wind. Here we present radio observations of Apep from 0.2 to 20 GHz taken over 33 yr. The spectrum reveals an extremely steep turnover in the flux density at low frequencies, where the flux density decreases by two orders of magnitude over only 325 MHz of bandwidth. This exponential decline is best described by free–free absorption, with a turnover frequency at 0.54 ± 0.01 GHz. Above the turnover, the spectrum is well described by a power-law and a high-frequency cut-off likely caused by inverse-Compton cooling. The light curve of Apep shows significant variation over the observing period, with Apep brightening by over 50 mJy in a span of 25 yr at 1.4 GHz. Models that assume spherical winds do not replicate all of the structure evident in the radio light curve. We derived a model that allows one of the winds in the system to be anisotropic. This anisotropic model recovers most of the structure of the light curve and is a significantly better statistical fit to the data than the spherical wind model. We suggest such a result is independent support that one of the Wolf–Rayet stars in Apep is launching an anisotropic wind. If the anisotropic wind model is correct, we predict a ∼25 per cent decrease of the 1.4 GHz flux density of Apep over the next 5 yr.

3D MHD simulation of a pulsationally driven MRI decretion disc

ABSTRACT We explore the pulsationally driven orbital mass ejection mechanism for Be star disc formation using isothermal, 3D magnetohydrodynamic (MHD) and hydrodynamic simulations. Non-radial pulsations are added to a star rotating at 95 per cent of critical as an inner boundary condition that feeds gas into the domain. In MHD, the initial magnetic field within the star is weak. The hydrodynamics simulation has limited angular momentum transport, resulting in repeating cycles of mass accumulation into a rotationally supported disc at small radii followed by fall-back on to the star. The MHD simulation, conversely, has efficient (Maxwell αM ∼ 0.04) angular momentum transport provided by both turbulent and coherent magnetic fields: a slowly decreting midplane driven by the magnetorotational instability and a supersonic wind on the surface of the disc driven by global magnetic torques. The angle and time-averaged properties near the midplane agree reasonably well with a 1D viscous decretion disc model with a modified $\tilde{\alpha }=0.5$, in which the gas transitions from a subsonic thin disc to a supersonic spherical wind at the critical point. 1D models, however, cannot capture the multiphase decretion/angular structure seen in our simulations. Our results demonstrate that, at least under certain conditions, non-radial pulsations on the surface of a rapidly rotating, weakly magnetised star can drive a Keplerian disc with the basic properties of the viscous decretion disc paradigm, albeit coupled to a laminar wind away from the midplane. Future modelling of Be star discs should consider the possible existence of such a surface wind.

Nucleosynthesis signatures of neutrino-driven winds from proto-neutron stars: a perspective from chemical evolution models

ABSTRACT We test the hypothesis that the observed first-peak (Sr, Y, Zr) and second-peak (Ba) s-process elemental abundances in low-metallicity Milky Way stars, and the abundances of the elements Mo and Ru, can be explained by a pervasive r-process contribution originating in neutrino-driven winds from highly magnetic and rapidly rotating proto-neutron stars (proto-NSs). We construct chemical evolution models that incorporate recent calculations of proto-NS yields in addition to contributions from asymptotic giant branch stars, Type Ia supernovae, and two alternative sets of yields for massive star winds and core-collapse supernovae. For non-rotating massive star yields from either set, models without proto-NS winds underpredict the observed s-process peak abundances by 0.3–$1\, \text{dex}$ at low metallicity, and they severely underpredict Mo and Ru at all metallicities. Models incorporating wind yields from proto-NSs with spin periods P ∼ 2–$5\, \text{ms}$ fit the observed trends for all these elements well. Alternatively, models omitting proto-NS winds but adopting yields of rapidly rotating massive stars, with vrot between 150 and $300\, \text{km}\, \text{s}^{-1}$, can explain the observed abundance levels reasonably well for [Fe/H] < −2. These models overpredict [Sr/Fe] and [Mo/Fe] at higher metallicities, but with a tuned dependence of vrot on stellar metallicity they might achieve an acceptable fit at all [Fe/H]. If many proto-NSs are born with strong magnetic fields and short spin periods, then their neutrino-driven winds provide a natural source for Sr, Y, Zr, Mo, Ru, and Ba in low-metallicity stellar populations. Conversely, spherical winds from unmagnetized proto-NSs overproduce the observed Sr, Y, and Zr abundances by a large factor.

Molecules and Outflows in NML Cygni: New Insights from a 1 mm Spectral Line Survey

Abstract A sensitive (1σ rms ≤ 3 mK) 1 mm spectral survey (214.5–285.5 GHz) of the oxygen-rich circumstellar envelope of the red hypergiant NML Cygni (NML Cyg) has been conducted using the Sub-millimeter Telescope (SMT) of the Arizona Radio Observatory (ARO). Over 100 spectral lines were detected, arising from 17 different molecules, including the carbon-bearing species CO, HCN, HCO+, CN, and HNC; sulfur- and silicon-containing compounds H2S, SO, SO2, SiO, and SiS; and more exotic NaCl and AlO. The 1 mm spectrum of NML Cyg closely resembles that of VY Canis Majoris (VY CMa) suggesting that the chemistries of hypergiant stars are similar. The line profiles in NML Cyg consist of multiple velocity features, particularly evident in SO2 and SO. In addition to a spherical wind at the star’s systemic velocity, the spectra suggest an asymmetric, blueshifted component near V LSR = −21 ± 3 km s−1 and a collimated, redshifted component near 15 ± 3 km s−1, positioned ∼34° and ∼12°, respectively, from the line of sight. The red- and blueshifted flows appear to be randomly oriented, and likely trace sporadic mass loss events. Their LSR velocities align closely with those of 22 GHz water masers, suggesting an NE–SE orientation. The winds may also be associated with the asymmetric nebula in F555W HST images but extending to 5″ (∼600R *). NML Cyg appears to be another example of rare, massive stars with collimated, episodic ejections, analogous to Betelgeuse and VY CMa, lending support for an important new mass loss mechanism—surface activity.

Moving-mesh radiation-hydrodynamic simulations of wind-reprocessed transients

ABSTRACT Motivated by recent theoretical work on tidal disruption events and other peculiar transients, we present moving-mesh radiation-hydrodynamic simulations of radiative luminosity emitted by a central source being reprocessed by a wind-like outflow. We couple the moving-mesh hydrodynamic code jet with our newly developed radiation module based on mixed-frame grey flux-limited diffusion with implicit timestep update. This allows us to study the self-consistent multidimensional radiation-hydrodynamic evolution over more than 10 orders of magnitude in both space and time in a single run. We simulate an optically thick spherical wind with constant or evolving mass-loss rate, which is irradiated by a central isotropic or angularly dependent radiation source. Our spherically symmetric simulations confirm previous analytic results by identifying different stages of radiation reprocessing: radiation trapped in the wind, diffusing out through the wind, and reaching constant maximum attenuation. We find that confining the central radiation source in a cone with moderate opening angles decrease up to one order of magnitude the early flux along sightlines oriented away from the direction of radiation injection but that the reprocessed radiation becomes isotropic roughly after one lateral diffusion time through the ejecta. We discuss further applications and guidelines for the use of our novel radiation-hydrodynamics tool in the context of transient modelling.