Search

The nearby M dwarf WX UMa has recently been detected at radio wavelengths with LOFAR. The combination of its observed brightness temperature and circular polarisation fraction suggests that the emission is generated via the electron-cyclotron maser instability. Two distinct mechanisms have been proposed to power such emission from low-mass stars: either a sub-Alfv\'enic interaction between the stellar magnetic field and an orbiting planet, or reconnection at the edge of the stellar magnetosphere. In this paper, we investigate the feasibility of both mechanisms, utilising the information about the star's surrounding plasma environment obtained from modelling its stellar wind. Using this information, we show that a Neptune-sized exoplanet with a magnetic field strength of 10-100 G orbiting at ~0.034 au can accurately reproduce the observed radio emission from the star, with corresponding orbital periods of 7.4 days. Due to the stellar inclination, a planet in an equatorial orbit is unlikely to transit the star. While such a planet could induce radial velocity semi-amplitudes from 7 to 396 m s$^{-1}$, it is unlikely that this signal could be detected with current techniques due to the activity of the host star. The application of our planet-induced radio emission model here illustrates its exciting potential as a new tool for identifying planet-hosting candidates from long-term radio monitoring. We also develop a model to investigate the reconnection-powered emission scenario. While this approach produces less favourable results than the planet-induced scenario, it nevertheless serves as a potential alternative emission mechanism which is worth exploring further., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

The clumpy stellar wind from the companion star in high mass X-ray binaries causes variable, partial absorption of the emission from the X-ray source. We studied XMM-Newton observations from the 7.22 d-long "Cyg X-1 Hard state Observations of a Complete Binary Orbit in X-rays" (CHOCBOX) monitoring campaign, in order to constrain the effects of the stellar wind on the short-timescale X-ray spectral-timing properties of the source. We find these properties to change significantly in the presence of the wind. In particular, the longest sampled timescales (corresponding to temporal frequencies of $\nu\sim$ 0.1-1 Hz) reveal an enhancement of the fractional variability power, while on the shortest sampled timescales ($\nu\sim$ 1-10 Hz) the variability is suppressed. In addition, we observe a reduction (by up to a factor of $\sim$ 1.8) of the otherwise high coherence between soft and hard band light curves, as well as of the amplitude of the hard X-ray lags intrinsic to the X-ray continuum. The observed increase of low frequency variability power can be explained in terms of variations of the wind column density as a consequence of motions of the intervening clumps. In this scenario (and assuming a terminal velocity of $v_{\infty}=2400\ {\rm km\ s^{-1}}$), we obtain an estimate of $l \sim$ 0.5-1.5 $\times 10^{-4} R_{\ast}$ for the average radial size of a clump. On the other hand, we suggest the behaviour at high frequencies to be due to scattering in an optically thicker medium, possibly formed by collision of the stellar wind with the edge of the disc., Comment: 16 pages, 13 figures

We present the analysis of Very Large Telescope Multi Unit Spectroscopic Explorer (MUSE) observations of the planetary nebula (PN) IC 4406. MUSE images in key emission lines are used to unveil the presence of at least five ring-like structures north and south of the main nebula of IC 4406. MUSE spectra are extracted from the rings to unambiguously assess for the first time in a PN their physical conditions, electron density (ne), and temperature (Te). The rings are found to have similar Te as the rim of the main nebula, but smaller ne. Ratios between different ionic species suggest that the rings of IC 4406 have a lower ionization state than the main cavity, in contrast to what was suggested for the rings in NGC 6543, the Cat’s Eye Nebula. © 2022 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., GR-L acknowledges support from CONACyT grant 263373 and PRODEP (Mexico). VMAGG acknowledges support from the Programa de Becas posdoctorales of the Dirección General de Asuntos del Personal Académico (DGAPA) of the Universidad Nacional Autónoma de México (UNAM, Mexico). VMAGG and JAT acknowledge funding by DGAPA UNAM PAPIIT project IA100720. JAT also acknowledges support from the Marcos Moshinsky Fundation (Mexico). JBR-G thanks CONACyT for a student scholarship. MAG acknowledges support of the Spanish Ministerio de Ciencia, Innovación y Universidades grant PGC2018-102184-B-I00, co-funded by FEDER funds. Based on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA), and the Canadian Astronomy Data Centre (CADC/NRC/CSA). Based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere under ESO programme 60.A-9100. Based on observations obtained at the Canada–France–Hawaii Telescope (CFHT). This work has made extensive use of the NASA’s Astrophysics Data System.

The central star of the planetary nebula (PN) HuBi 1 has been recently proposed to have experienced a very late thermal pulse (VLTP), but the dilution of the emission of the recent ejecta by that of the surrounding H-rich old outer shell has so far hindered confirming its suspected H-poor nature. We present here an analysis of the optical properties of the ejecta in the innermost regions of HuBi 1 using MEGARA high-dispersion integral field and OSIRIS intermediate-dispersion long-slit spectroscopic observations obtained with the 10.4-m Gran Telescopio de Canarias. The unprecedented tomographic capability of MEGARA to resolve structures in velocity space allowed us to disentangle for the first time the Hα and Hβ emission of the recent ejecta from that of the outer shell. The recent ejecta is found to have much higher extinction than the outer shell, implying the presence of large amounts of dust. The spatial distribution of the emission from the ejecta and the locus of key line ratios in diagnostic diagrams probe the shock excitation of the inner ejecta in HuBi 1, in stark contrast with the photoionization nature of the H-rich outer shell. The abundances of the recent ejecta have been computed using the MAPPINGS V code under a shock scenario. They are found to be consistent with a born-again ejection scenario experienced by the progenitor star, which is thus firmly confirmed as a new ‘born-again’ star. © 2022 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., BM-M and MAG are funded by the Spanish Ministerio de Ciencia, Innovación y Universidades (MCIU) grant PGC2018-102184-B-I00, co-funded by FEDER funds. BM-M, BP, MAG, and SC acknowledge support from the State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). JAT is funded by UNAM DGAPA PAPIIT project IA100720 and the Marcos Moshinsky Fundation (Mexico). This work has made used of observations made with the Gran Telescopio Canarias (GTC), installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, in the island of La Palma. This work made use of IRAF, distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy under cooperative agreement with the National Science Foundation. This work has made extensive use of NASA’s Astrophysics Data System. Based on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA).

We present results from new self-consistent 3D MHD simulations of the magnetospheres from massive stars with a dipole magnetic axis that has a non-zero obliquity angle ($\beta$) to the star's rotation axis. As an initial direct application, we compare the global structure of co-rotating disks for nearly aligned ($\beta=5^o$) versus half-oblique ($\beta=45^o$) models, both with moderately rapid rotation ($\sim$ 0.5 critical). We find that accumulation surfaces broadly resemble the forms predicted by the analytic Rigidly Rotating Magnetosphere (RRM) model, but the mass buildup to near the critical level for centrifugal breakout against magnetic confinement distorts the field from the imposed initial dipole. This leads to an associated warping of the accumulation surface toward the rotational equator, with the highest density concentrated in {\em wings} centered on the intersection between the magnetic and rotational equators. These MHD models can be used to synthesize rotational modulation of photometric absorption and H$\alpha$ emission for a direct comparison with observations., Comment: Accepted for publication in MNRAS

Luminous Blue Variables (LBVs) are massive stars that are likely to be a transitionary phase between O stars and hydrogen-free classical Wolf–Rayet stars. The variability of these stars has been an area of study for both professional and amateur astronomers for more than a century. In this paper, we present 5 yr of precision photometry of the classical LBV P Cygni taken with the BRITE-Constellation nanosatellites. We have analyzed these data with Fourier analysis to search for periodicities that could elucidate the drivers of variability for these stars. These data show some long-time-scale variability over the course of all six calendar years of observations, but the frequencies needed to reproduce the individual light curves are not consistent from 1 yr to the next. These results likely show that there is no periodic phenomenon present for P Cygni, meaning that the variability is largely stochastic. We interpret the data as being caused by internal gravity waves similar to those seen in other massive stars, with P Cygni exhibiting a larger amplitude and lower characteristic frequency than the main-sequence or blue supergiant stars previously studied. These results show evidence that LBVs may be an extrapolation of the blue supergiants, which have previously been shown to be an extension of main-sequence stars in the context of the stochastic low-frequency photometric variability.

M33 X-7 is the only known eclipsing black hole high mass X-ray binary. The system is reported to contain a very massive O supergiant donor and a massive black hole in a short orbit. The high X-ray luminosity and its location in the metal-poor galaxy M33 make it a unique laboratory for studying the winds of metal-poor donor stars with black hole companions and it helps us to understand the potential progenitors of black hole mergers. Using phase-resolved simultaneous HST- and XMM-Newton-observations, we traced the interaction of the stellar wind with the black hole. Our comprehensive spectroscopic investigation of the donor star (X-ray+UV+optical) yields new stellar and wind parameters for the system that differs significantly from previous estimates. In particular, the masses of the components are considerably reduced to 38 for the O-star donor and 11.4 for the black hole. The O giant is overfilling its Roche lobe and shows surface He enrichment. The donor shows a densely clumped wind with a mass-loss rate that matches theoretical predictions. We investigated the wind-driving contributions from different ions and the changes in the ionization structure due to X-ray illumination. Toward the black hole, the wind is strongly quenched due to strong X-ray illumination. For this system, the standard wind-fed accretion scenario alone cannot explain the observed X-ray luminosity, pointing toward an additional mass overflow, which is in line with our acceleration calculations. The X-ray photoionization creates an He II emission region emitting $10^{47}$ ph/s. We computed binary evolutionary tracks for the system using MESA. Currently, the system is transitioning toward an unstable mass transfer phase, resulting in a common envelope of the black hole and donor. Since the mass ratio is q~3.3 and the period is short, the system is unlikely to survive the common envelope, but will rather merge., Accepted for publication in Astronomy & Astrophysics

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $\mu$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust., Comment: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

We present the results of a four-month, spectroscopic campaign of the Wolf-Rayet dust-making binary, WR137. We detect only small-amplitude, random variability in the CIII5696 emission line and its integrated quantities (radial velocity, equivalent width, skewness, kurtosis) that can be explained by stochastic clumps in the wind of the WC star. We find no evidence of large-scale, periodic variations often associated with Corotating Interaction Regions that could have explained the observed intrinsic continuum polarization of this star. Our moderately high-resolution and high signal-to-noise average Keck spectrum shows narrow double-peak emission profiles in the Halpha, Hbeta, Hgamma, HeII6678 and HeII5876 lines. These peaks have a stable blue-to-red intensity ratio with a mean of 0.997 and a root-mean-square of 0.004, commensurate with the noise level; no variability is found during the entire observing period. We suggest that these profiles arise in a decretion disk around the O9 companion, which is thus an O9e star. The characteristics of the profiles are compatible with those of other Be/Oe stars. The presence of this disk can explain the constant component of the continuum polarization of this system, for which the angle is perpendicular to the plane of the orbit, implying that the rotation axis of the O9e star is aligned with that of the orbit. It remains to be explained why the disk is so stable within the strong ultraviolet radiation field of the O star. We present a binary evolutionary scenario that is compatible with the current stellar and system parameters., 12 pages, 5 Figures

The interpretation of recent observations of bow shocks around O-stars and the creation of corresponding models require a detailed understanding of the associated (magneto-)hydrodynamic structures. We base our study on three-dimensional numerical magneto-)hydrodynamical models, which are analyzed using the dynamically relevant parameters, in particular, the (magneto)sonic Mach numbers. The analytic Rankine-Hugoniot relation for HD and MHD are compared with those obtained by the numerical model. In that context we also show that the only distance which can be approximately determined is that of the termination shock, if it is a hydrodynamical shock. For MHD shocks the stagnation point does not, in general, lie on the inflow line, which is the line parallel to the inflow vector and passing through the star. Thus an estimate via the Bernoulli equation as in the HD case is, in general, not possible. We also show that in O-star astrospheres, distinct regions exist in which the fast, slow, Alfv\'enic, and sonic Mach numbers become lower than one, implying sub-slow magnetosonic as well as sub-fast and sub-sonic flows. Nevertheless, the analytic MHD Rankine Hugoniot relations can be used for further studies of turbulence and cosmic ray modulation., Comment: 14 pages, 9 figures

2MASS J20395358+4222505 is an obscured early B supergiant near the massive OB star association Cygnus OB2. Despite its bright infrared magnitude (Ks = 5.82) it has remained largely ignored because of its dim optical magnitude (B = 16.63, V = 13.68). In a previous paper, we classified it as a highly reddened, potentially extremely luminous, early B-type supergiant. We obtained its spectrum in the U, B and R spectral bands during commissioning observations with the instrument MEGARA at the Gran Telescopio CANARIAS. It displays a particularly strong Hα emission for its spectral type, B1 Ia. The star seems to be in an intermediate phase between supergiant and hypergiant, a group that it will probably join in the near (astronomical) future. We observe a radial velocity difference between individual observations and determine the stellar parameters, obtaining Teff = 24 000 K and log gc = 2.88 ± 0.15. The rotational velocity found is large for a B supergiant, v sin i = 110 ± 25 kms−1⁠. The abundance pattern is consistent with solar, with a mild C underabundance (based on a single line). Assuming that J20395358+4222505 is at the distance of Cyg OB2, we derive the radius from infrared photometry, finding R = 41.2 ± 4.0 R⊙, log(L/L⊙) = 5.71 ± 0.04 and a spectroscopic mass of 46.5 ± 15.0 M⊙. The clumped mass-loss rate (clumping factor 10) is very high for the spectral type, M˙ = 2.4 × 10−6 M⊙ a−1. The high rotational velocity and mass-loss rate place the star at the hot side of the bi-stability jump. Together with the nearly solar CNO abundance pattern, they may also point to evolution in a binary system, J20395358+4222505 being the initial secondary. © 2022 The Author(s)., SS-D and AH acknowledge support from the Spanish Government Ministerio de Ciencia e Innovación through grants PGC-2018-091 3741-B-C22 and CEX2019-000920-S and from the Canarian Agency for Research, Innovation and Information Society (ACIISI), of the Canary Islands Government, and the European Regional Development Fund (ERDF), under grant with reference ProID2020010016. MG and FN acknowledge financial support through Spanish grant PID2019-105552RB-C41 (MINECO/MCIU/AEI/FEDER) and from the Spanish State Research Agency (AEI) through the Unidad de Excelencia ‘María de Maeztu’-Centro de Astrobiología (CSIC-INTA) project No. MDM-2017-0737. SRB acknowledges support by the Spanish Government under grants AYA2015-68012-C2-2-P and PGC2018-093741-B-C21/C22 (MICIU/AEI/FEDER, UE). SRA acknowledges funding support from the FONDECYT Iniciación project 11171025 and the FONDECYT Regular project 1201490. JIP acknowledges finantial support from projects Estallidos6 AYA2016-79724-C4 (Spanish Ministerio de Economia y Competitividad), Estallidos7 PID2019-107408GB-C44 (Spanish Ministerio de Ciencia e Innovacion), grant P18-FR-2664 (Junta de Andalucía), and grant SEV-2017-0709 ‘Centro de Excelencia Severo Ochoa Program’ (Spanish Science Ministry). AGP, SP, AG-M, JG and NC acknowledge support from the Spanish MCI through project RTI2018-096188-B-I00.

Context. There is a population of runaway stars that move at extremely high speeds with respect to their surroundings. The fast motion and the stellar wind of these stars, plus the wind-medium interaction, can lead to particle acceleration and non-thermal radiation. Aims. We characterise the interaction between the winds of fast runaway stars and their environment, in particular to establish their potential as cosmic-ray accelerators and non-thermal emitters. Methods. We model the hydrodynamics of the interaction between the stellar wind and the surrounding material. We self-consistently calculate the injection and transport of relativistic particles in the bow shock using a multi-zone code, and compute their broadband emission from radio to $\gamma$-rays. Results. Both the forward and reverse shocks are favourable sites for particle acceleration, although the radiative efficiency of particles is low and therefore the expected fluxes are in general rather faint. Conclusions. We show that high-sensitivity observations in the radio band can be used to detect the non-thermal radiation associated with bow shocks from hypervelocity and semi-relativistic stars. Hypervelocity stars are expected to be modest sources of sub-TeV cosmic rays, accounting perhaps for a $\sim 0.1$% of that of galactic cosmic rays., Comment: 11 pages, 8 figures, accepted for publication section 2. Astrophysical processes of Astronomy and Astrophysics

[Context] Relativistic jets are ubiquitous in the Universe. In microquasars, especially in high-mass X-ray binaries, the interaction of jets with the strong winds driven by the massive and hot companion star in the vicinity of the compact object is fundamental for understanding the jet dynamics, nonthermal emission, and long-term stability. However, the role of the jet magnetic field in this process is unclear. In particular, it is still debated whether the magnetic field favors jet collimation or triggers more instabilities that can jeopardize the jet evolution outside the binary. [Aims] We study the dynamical role of weak and moderate to strong toroidal magnetic fields during the first several hundred seconds of jet propagation through the stellar wind, focusing on the magnetized flow dynamics and the mechanisms of energy conversion. [Methods] We developed the code Lóstrego v1.0, a new 3D relativistic magnetohydrodynamics code to simulate astrophysical plasmas in Cartesian coordinates. Using this tool, we performed the first 3D relativistic magnetohydrodynamics numerical simulations of relativistic magnetized jets propagating through the clumpy stellar wind in a high-mass X-ray binary. To highlight the effect of the magnetic field in the jet dynamics, we compared the results of our analysis with those of previous hydrodynamical simulations. [Result] The overall morphology and dynamics of weakly magnetized jet models is similar to previous hydrodynamical simulations, where the jet head generates a strong shock in the ambient medium and the initial overpressure with respect to the stellar wind drives one or more recollimation shocks. On the timescales of our simulations (i.e., ta, The project that gave rise to these results received the support of a fellowship from “la Caixa” Foundation (ID 100010434). The fellowship code is LCF/BQ/DR19/11740030. J.L.M. acknowledges additional support from the Spanish Ministerio de Ciencia through grant PID2019-105510GB-C31. M.P. and J.M-M. acknowledge support by the Spanish Ministerio de Ciencia through grants PID2019-107427GB-C33 and PGC2018-095984-B-I00, and from the Generalitat Valenciana through grant PROMETEU/2019/071. M.P. acknowledges additional support from the Spanish Ministerio de Ciencia through grant PID2019-105510GB-C31. V.B.-R. acknowledges financial support from the State Agency for Research of the Spanish Ministry of Science and Innovation under grant PID2019-105510GB-C31 and through the ”Unit of Excellence María de Maeztu 2020-2023” award to the Institute of Cosmos Sciences (CEX2019-000918-M), and by the Catalan DEC grant 2017 SGR 643. V.B-R. is Correspondent Researcher of CONICET, Argentina, at the IAR.

We present the first instalment of a deep imaging catalogue containing 58 True, Likely, and Possible extended PNe detected with the Isaac Newton Telescope Photometric H α Survey (IPHAS). The three narrow-band filters in the emission lines of H α, [N ii] λ6584 Å, and [O iii] λ5007 Å used for this purpose allowed us to improve our description of the morphology and dimensions of the nebulae. In some cases even the nature of the source has been reassessed. We were then able to unveil new macro- and micro-structures, which will without a doubt contribute to a more accurate analysis of these PNe. It has been also possible to perform a primary classification of the targets based on their ionization level. A Deep Learning classification tool has also been tested. We expect that all the PNe from the IPHAS catalogue of new extended planetary nebulae will ultimately be part of this deep H α, [N ii], and [O iii] imaging catalogue. © 2021 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., LS acknowledges support from PAPIIT grant IN-101819 (Mexico). MAG acknowledges support from grant AYA PGC2018-102184-B-I00 co-funded with FEDER funds. GR-L acknowledges support from CONACyT (grant 263373) and PRODEP (Mexico). JAT acknowledges funding by Dirección General de Asuntos del Personal Academico of the Universidad Nacional Autónoma de México (DGAPA, UNAM) project IA100720 and the Marcos Moshinsky Fundation (Mexico). AAZ acknowledges funding from STFC under grant number ST/T000414/1 and Newton grant ST/R006768/1, and from a Hung Hing Ying visiting professorship at the University of Hong Kong. DNFAI acknowledges funding under ‘Deep Learning for Classification of Astronomical Archives’ from the Newton-Ungku Omar Fund: F08/STFC/1792/2018., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

The recently discovered massive binary system Apep is the most powerful synchrotron emitter among the known Galactic colliding-wind binaries. This makes this particular system of great interest to investigate stellar winds and the non-thermal processes associated with their shocks. This source was detected at various radio bands, and in addition the wind-collision region was resolved by means of very-long baseline interferometric observations. We use a non-thermal emission model for colliding-wind binaries to derive physical properties of this system. The observed morphology in the resolved maps allows us to estimate the system projection angle on the sky to be $\psi \approx 85^\circ$. The observed radio flux densities also allow us to characterise both the intrinsic synchrotron spectrum of the source and its modifications due to free--free absorption in the stellar winds at low frequencies; from this we derive mass-loss rates of the stars of $\dot{M}_\mathrm{WN} \approx 4\times10^{-5}$ $\mathrm{M}_\odot$ yr$^{-1}$ and $\dot{M}_\mathrm{WC} \approx 2.9\times10^{-5}$ $\mathrm{M}_\odot$ yr$^{-1}$. Finally, the broadband spectral energy distribution is calculated for different combinations of the remaining free parameters, namely the intensity of the magnetic field and the injected power in non-thermal particles. We show that the degeneracy of these two parameters can be solved with observations in the high-energy domain, most likely in the hard X-rays but also possibly in $\gamma$-rays under favourable conditions., Comment: 10 pages, 5 figures, 2 tables. This version has been accepted for publication in PASA after improvements and updates in the discussion of the results

Non-ideal magnetohydrodynamic (MHD) processes -- namely Ohmic resistivity, ambipolar diffusion and the Hall effect -- modify the early stages of the star formation process and the surrounding environment. Collectively, they have been shown to promote disc formation and promote or hinder outflows. But which non-ideal process has the greatest impact? Using three-dimensional smoothed particle radiation non-ideal MHD simulations, we model the gravitational collapse of a rotating, magnetised cloud through the first hydrostatic core phase to shortly after the formation of the stellar core. We investigate the impact of each process individually and collectively. Including any non-ideal process decreases the maximum magnetic field strength by at least an order of magnitude during the first core phase compared to using ideal MHD, and promotes the formation of a magnetic wall. When the magnetic field and rotation vectors are anti-aligned and the Hall effect is included, rotationally supported discs of $r \gtrsim 20$~au form; when only the Hall effect is included and the vectors are aligned, a counter-rotating pseudo-disc forms that is not rotationally supported. Rotationally supported discs of $r \lesssim 4$~au form if only Ohmic resistivity or ambipolar diffusion are included. The Hall effect suppresses first core outflows when the vectors are anti-aligned and suppresses stellar core outflows independent of alignment. Ohmic resistivity and ambipolar diffusion each promote first core outflows and delay the launching of stellar core outflows. Although each non-ideal process influences star formation, these results suggest that the Hall effect has the greatest influence., 19 pages, 13 figures. Accepted for publication in MNRAS

Context. Eta Carinae is a highly eccentric, massive binary system (semimajor axis ~15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, η Car A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims. Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods. Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results. We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar windin several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R ~ 960). Our observations were performed close to periastron passage in February 2020 (orbital phase ~ 14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (~14 au). The radius of the faintest outer wind regions is ~26 mas (~60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84 ± 0.06 mas (6.54 ± 0.14 au). We comparethe derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

Recent large-scale spectropolarimetric surveys have established that a small but significant percentage of massive stars host stable, surface dipolar magnetic fields with strengths on the order of kG. These fields channel the dense, radiatively driven stellar wind into circumstellar magnetospheres, whose density and velocity structure can be probed using ultraviolet (UV) spectroscopy of wind-sensitive resonance lines. Coupled with appropriate magnetosphere models, UV spectroscopy provides a valuable way to investigate the wind-field interaction, and can yield quantitative estimates of the wind parameters of magnetic massive stars. We report a systematic investigation of the formation of UV resonance lines in slowly rotating magnetic massive stars with dynamical magnetospheres. We pair the Analytic Dynamical Magnetosphere (ADM) formalism with a simplified radiative transfer technique to produce synthetic UV line profiles. Using a grid of models, we examine the effect of magnetosphere size, the line strength parameter, and the cooling parameter on the structure and modulation of the line profile. We find that magnetic massive stars uniquely exhibit redshifted absorption at most viewing angles and magnetosphere sizes, and that significant changes to the shape and variation of the line profile with varying line strengths can be explained by examining the individual wind components described in the ADM formalism. Finally, we show that the cooling parameter has a negligible effect on the line profiles., 16 pages, 15 figures, accepted to MNRAS

We have modelled the 3D distribution of OH/IR stars in the Galactic plane, traced by 1612 MHz OH maser sources with classic double horned spectral profiles. We statistically analysed over 700 maser sources detected by the HI/OH/Recombination line survey of the Milky Way (THOR) and the Australia Telescope Compact Array interferometric follow-up observations of the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH). With a simple model constructed from a classical density distribution of stars and luminosity functions of OH maser sources in the Galaxy, we estimate the scale height, or the half thickness of the OH/IR star distribution along the Galactic disc to be 90–290 pc. The simple model also implies that there are ∼4000 OH/IR stars hosting 1612 MHz OH masers along the Galactic Plane. Therefore, next generation telescopes such as the Australian Square Kilometre Array Pathfinder (ASKAP) and SKA Phase 1 will detect about 80 per cent of such OH/IR stars in the Galaxy at a 10 mJy detection limit. Comparing the data of previously detected circumstellar 1612 MHz OH maser sources with those of THOR and SPLASH, the maser source lifetime is estimated to be ∼300 yr. This is likely a lower limit, since non-detections of masers in some cases could be affected by the flux variation of the maser source. © 2021 The Author(s)., For deriving a kinematic distance, we used a python library provided by Trey Wenger (https://github.com/tvwenger/kd). This research has made use of: NASA’s Astrophysics Data System Abstract Service; NASA/IPAC Infrared Science Archive, which is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology; and Database of Circumstellar Masers. HI was supported by the Bilateral Collaboration Program by the Japan Society for Promotion of Science (JSPS) and KAKENHI programs 25610043 and 16H02167 by the Ministry of Education, Culture, Sports, Science and Technology (MEXT). JFG and HI were supported by the Invitation Program for Foreign Researcher by JSPS (S14128), by Amanogawa Galaxy Astronomy Research Center (AGARC), and the State Agency for Research of the Spanish MCIU through the ’Center of Excellence Severo Ochoa’ award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). JFG also acknowledges support by MCIU-AEI (Spain) grant AYA2017-84390-C2-R (co-funded by FEDER). JRD acknowledges the support of an Australian Research Council (ARC) DECRA Fellowship (project number DE170101086). H-HQ is partially supported by the Special Funding for Advanced Users, budgeted and administrated by Center for Astronomical Mega-Science, Chinese Academy of Sciences (CAMS-CAS), CAS ‘Light of West China’ Program and the National Natural Science Foundation of China (Grant No. 11903038).

Binary systems of Population III can evolve to microquasars when one of the stars collapses into a black hole. When the compact object accretes matter at a rate greater than the Eddington rate, powerful jets and winds driven by strong radiation pressure should form. We investigate the structure of the jet-wind system for a model of Population III microquasar on scales beyond the jet-wind formation region. Using relativistic hydrodynamic simulations we find that the ratio of kinetic power between the jet and the disk wind determines the configuration of the system. When the power is dominated by the wind, the jet fills a narrow channel, collimated by the dense outflow. When the jet dominates the power of the system, part of its energy is diverted turning the wind into a quasi-equatorial flow, while the jet widens. From the results of our simulations, we implement semi-analytical calculations of the impact of the quasi-equatorial wind on scales of the order of the size of the binary system. Our results indicate that Population III microquasars might inject gamma rays and relativistic particles into the early intergalactic medium, contributing to its reionization at large distances from the binary system., PSC is Fellow of CONICET and PhD student at Universidad Nacional de La Plata (UNLP). PSC thanks the UNLP for the training received during his undergraduate studies, and the public education system in Argentina. GER is very grateful to the ICCUB where part of this research was done. This work was supported by the Argentine agency CONICET (PIP 2014-00338) and the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant PID2019-105510GBC31 and though the “Center of Excellence María de Maeztu 2020-2023” award to the ICCUB (CEX2019-000918-M). V.B-R. is Correspondent Researcher of CONICET, Argentina, at the IAR.

Pastorello, A., et al., We present the results of our monitoring campaigns of the luminous red novae (LRNe) AT 2020hat in NGC 5068 and AT 2020kog in NGC 6106. The two objects were imaged (and detected) before their discovery by routine survey operations. They show a general trend of slow luminosity rise, lasting at least a few months. The subsequent major LRN outbursts were extensively followed in photometry and spectroscopy. The light curves present an initial short-duration peak, followed by a redder plateau phase. AT 2020kog is a moderately luminous event peaking at ∼7 × 1040 erg s-1, while AT 2020hat is almost one order of magnitude fainter than AT 2020kog, although it is still more luminous than V838 Mon. In analogy with other LRNe, the spectra of AT 2020kog change significantly with time. They resemble those of type IIn supernovae at early phases, then they become similar to those of K-type stars during the plateau, and to M-type stars at very late phases. In contrast, AT 2020hat already shows a redder continuum at early epochs, and its spectrum shows the late appearance of molecular bands. A moderate-resolution spectrum of AT 2020hat taken at +37 d after maximum shows a forest of narrow P Cygni lines of metals with velocities of 180 km s-1, along with an Hα emission with a full-width at half-maximum velocity of 250 km s-1. For AT 2020hat, a robust constraint on its quiescent progenitor is provided by archival images of the Hubble Space Telescope. The progenitor is clearly detected as a mid-K type star, with an absolute magnitude of MF606W = -3.33 ± 0.09 mag and a colour of F606W - F814W = 1.14 ± 0.05 mag, which are inconsistent with the expectations from a massive star that could later produce a core-collapse supernova. Although quite peculiar, the two objects nicely match the progenitor versus light curve absolute magnitude correlations discussed in the literature., MF gratefully acknowledges the support of a Royal Society – Science Foundation Ireland University Research Fellowship. Research by SV and YD is supported by NSF grants AST–1813176 and AST-2008108. Time domain research by DJS is supported by NSF grants AST-1813466, 1908972, and by the Heising-Simons Foundation under grant #2020-1864. YZC is supported in part by National Natural Science Foundation of China (NSFC grants 12033003, 11633002, 11325313, and 11761141001), National Program on Key Research and Development Project (grant no. 2016YFA0400803), and Ma Huateng Foun- dation. DJ acknowledges support from the State Research Agency (AEI) of the Spanish Ministry of Science, Innovation and Universities (MCIU) and the European Regional Development Fund (FEDER) under grant AYA2017-83383-P. DJ also acknowledges support under grant P/308614 financed by funds transferred from the Spanish Ministry of Science, Innovation and Universities, charged to the General State Budgets and with funds transferred from the General Budgets of the Autonomous Community of the Canary Islands by the Ministry of Economy, Industry, Trade and Knowledge. SJS and KWS acknowledge funding from STFC Grants ST/P000312/1, ST/T000198/1 and ST/S006109/1. EC and AF are partially supported by the PRIN-INAF 2017 with the project Towards the SKA and CTA era: discovery, localisation, and physics of transients sources (P.I. M. Giroletti). MS is supported by generous grants from Villum FONDEN (13261,28021) and by a project grant (8021-00170B) from the Independent Research Fund Denmark. This work is based on observations made with the Nordic Optical Telescope, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias; the Gran Telescopio Canarias (GTC), installed in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias, in the Island of La Palma; the Liverpool Telescope operated on the island of La Palma by Liverpool John Moores University at the Spanish Observatorio del Roque de los Muchachos of the Insti-tuto de Astrofísica de Canarias with financial support from the UK Science and Technology Facilities Council; the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundación Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias.

Binary systems that host a massive star and a non-accreting pulsar can be powerful non-thermal emitters. The relativistic pulsar wind and the non-relativistic stellar outflows interact along the orbit, producing ultrarelativistic particles that radiate from radio to gamma rays. To properly characterize the physics of these sources, and better understand their emission and impact on the environment, careful modelling of the outflow interactions, spanning a broad range of spatial and temporal scales, is needed. Full 3-dimensional approaches are very computationally expensive, but simpler approximate approaches, while still realistic at the semi-quantitative level, are available. We present here the results of calculations done with a quasi 3-dimensional scheme to compute the evolution of the interacting flows in a region spanning in size up to a thousand times the size of the binary. In particular, we analyze for the first time the role of different eccentricities in the large scale evolution of the shocked flows. We find that the higher the eccentricity, the closer the flows behave like a one-side outflow, which becomes rather collimated for eccentricity values $\gtrsim 0.75$. The simulations also unveil that the pulsar and the stellar winds become fully mixed within the grid for low eccentricity systems, presenting a more stochastic behavior at large scales than in the highly eccentric systems., 18 pages, 11 figures

We present the follow-up campaign of the luminous red nova (LRN) AT 2019zhd, the third event of this class observed in M 31. The object was followed by several sky surveys for about five months before the outburst, during which it showed a slow luminosity rise. In this phase, the absolute magnitude ranged from Mr =-2.8 ± 0.2 mag to Mr =-5.6 ± 0.1 mag. Then, over a four to five day period, AT 2019zhd experienced a major brightening, reaching a peak of Mr =-9.61 ± 0.08 mag and an optical luminosity of 1.4 × 1039 erg s-1. After a fast decline, the light curve settled onto a short-duration plateau in the red bands. Although less pronounced, this feature is reminiscent of the second red maximum observed in other LRNe. This phase was followed by a rapid linear decline in all bands. At maximum, the spectra show a blue continuum with prominent Balmer emission lines. The post-maximum spectra show a much redder continuum, resembling that of an intermediate-type star. In this phase, Hα becomes very weak, Hβ is no longer detectable, and a forest of narrow absorption metal lines now dominate the spectrum. The latest spectra, obtained during the post-plateau decline, show a very red continuum (Teff ≈ 3000 K) with broad molecular bands of TiO, similar to those of M-type stars. The long-lasting, slow photometric rise observed before the peak resembles that of LRN V1309 Sco, which was interpreted as the signature of the common-envelope ejection. The subsequent outburst is likely due to the gas outflow following a stellar merging event. The inspection of archival HST images taken 22 years before the LRN discovery reveals a faint red source (MF555W = 0.21 ± 0.14 mag, with F555W-F814W = 2.96 ± 0.12 mag) at the position of AT 2019zhd, which is the most likely quiescent precursor. The source is consistent with expectations for a binary system including a predominant M5-type star.

A correction has been published: Monthly Notices of the Royal Astronomical Society, Volume 509, Issue 2, January 2022, Page 2896, https://doi.org/10.1093/mnras/stab3247, We present a comprehensive analysis of the compact planetary nebula M 2-31 investigating its spectral properties, spatio-kinematical structure, and chemical composition using Gran Telescopio Canarias (GTC) Multi-Espectrógrafo en GTC de Alta Resolución para Astronomía (MEGARA) integral field spectroscopic observations and Nordic Optical Telescope(NOT) Alhambra Faint Object Spectrograph and Camera (ALFOSC) medium-resolution spectra and narrow-band images. The GTC MEGARA high-dispersion observations have remarkable tomographic capabilities, producing an unprecedented view of the morphology and kinematics of M 2-31 that discloses a fast spectroscopic bipolar outflow along position angles 50° and 230°, an extended shell, and a toroidal structure or waist surrounding the central star perpendicularly aligned with the fast outflows. These observations also show that the C ii emission is confined in the central region and enclosed by the [N ii] emission. This is the first time that the spatial segregation revealed by a two-dimensional map of the C ii line implies the presence of multiple plasma components. The deep NOT ALFOSC observations allowed us to detect broad Wolf-Rayet (WR) features from the central star of M 2-31, including previously undetected broad O vi lines that suggest a reclassification as a [WO4]-type star. © 2021 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., JSRG and VMAGG acknowledge support from the Programas de Beca Posdoctorales funded by Dirección General de Asuntos del Personal Académico of the Universidad Nacional Autónoma de México (DGAPA, UNAM). JAT acknowledges funding by DGAPA, UNAM PAPIIT project IA100720 and the Marcos Moshinsky Foundation (Mexico). SC and MAG acknowledge financial support from State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). MAG acknowledges support of the Spanish Ministerio de Ciencia, Innovación y Universidades (MCIU) grant PGC2018-102184-B-100. LS acknowledges support from PAPIIT-UNAM grant IN101819. GR-L acknowledges support from CONACYT (grant 263373) and PRODEP (Mexico). The GTC Science Operations group is recognized for scheduling the GTC MEGARA observations beneath the exacting conditions requested by this program.

Context LS 5039 is a powerful high-mass gamma-ray binary that probably hosts a young non-accreting pulsar. However, despite the wealth of data available, the means by which the non-thermal emitter is powered are still unknown. Aims. We use a dynamical-radiative numerical model, and multiwavelength data, to constrain the properties of a hypothetical pulsar wind that would power the non-thermal emitter in LS 5039. Methods. We ran simulations of an ultrarelativistic (weakly magnetized) cold e±-wind that Compton scatters stellar photons and that dynamically interacts with the stellar wind. The effects of energy losses on the unshocked e±-wind dynamics, and the geometry of the two-wind contact discontinuity, are computed for different wind models. The predicted unshocked e±-wind radiation at periastron, when expected to be the highest, is compared to LS 5039 data. Results. The minimum possible radiation from an isotropic cold e±-wind overpredicts the X-ray to gamma-ray fluxes at periastron by a factor of ∼3. In the anisotropic (axisymmetric) wind case X-ray and ≳ 100 MeV data are not violated by wind radiation if the wind axis is at ≲ 20-40° from the line of sight (chance probability of ≲ 6-24%), depending on the anisotropic wind model, or if the wind Lorentz factor ∈ 102 - 103, in which case the wind power can be higher, but it requires e±-multiplicities of ∼106 and 109 for a 10-2 s and 10 s pulsar period, respectively. Conclusions. The studied model predicts that a weakly magnetized cold pulsar e±-wind in LS 5039 should be strongly anisotropic, with either a wind Lorentz factor ∈ 102 - 103 and very high multiplicities or with a fine-tuned wind orientation. A weakly magnetized, cold baryon-dominated wind would be a possible alternative, but then the multiplicities should be rather low, while the baryon-to-e± energy transfer should be very efficient at wind termination. A strongly magnetized cold wind seems to be the most favorable option as it is consistent with recent research on pulsar winds and does not require fine-tuning of the pulsar wind orientation, and the wind multiplicity and Lorentz factor are less constrained., Acknowledgements. We thank the anonymous referee for constructive and useful comments that helped to improve the manuscript. We are grateful to Dmitry Khangulyan for insightful comments on this work. V.B-R. acknowledges financial support by the Spanish Ministerio de Economía, Industria y Competitivi-dad (MINEICO/FEDER, UE) under grant AYA2016-76012-C3-1-P, from the State Agency for Research of the Spanish Ministry of Science and Innovation under grant PID2019-105510GB-C31 and through the ”Unit of Excellence María de Maeztu 2020-2023” award to the Institute of Cosmos Sciences (CEX2019-000918-M), and by the Catalan DEC grant 2017 SGR 643. V.B-R. is Correspondent Researcher of CONICET, Argentina, at the IAR.

We present a multiwavelength characterization of the planetary nebula (PN) NGC 6905 and its [Wolf–Rayet]-type ([WR]) central star (CSPN) HD 193949. Our Nordic Optical Telescope Alhambra Faint Object Spectrograph and Camera (ALFOSC) spectra and images unveil in unprecedented detail the high-ionization structure of NGC 6905. The high-quality spectra of HD 193949 allowed us to detect more than 20 WR features including the characteristic O-bump, blue bump, and red bump, which suggests a spectral type no later than a [WO2]-subtype. Moreover we detect the Ne VII and Ne VIII broad emission lines, rendering HD 193949 yet another CSPN with Teff ≲ 150 kK exhibiting such stellar emission lines. We studied the physical properties (Te and ne) and chemical abundances of different regions within NGC 6905 including its low-ionization clumps; abundances are found to be homogeneous. We used the PoWR stellar atmosphere code to model the spectrum of HD 193949, which is afterwards used in a photoionization model performed with CLOUDY that reproduces the nebular and dust properties for a total mass in the 0.31–0.47 M⊙ range and a mass of C-rich dust of ∼2 × 10−3 M⊙. Adopting a current stellar mass of 0.6 M⊙, our model suggests an initial mass ∼1 M⊙ for HD 193949, consistent with the observations. © 2021 The Author(s)., VMAGG acknowledges support from the Programa de Becas posdoctorales of the Dirección General de Asuntos del Personal Académico of the Universidad Nacional Autónoma de México (DGAPA, UNAM, Mexico). VMAGG and JAT acknowledge funding by DGAPA UNAM PAPIIT project IA100720. JAT also acknowledges support from the Marcos Moshinsky Fundation (Mexico). GR acknowledge support from Consejo Nacional de Ciencia y Tecnología (CONACYT) for student scholarship. MAG acknowledges support of the Spanish Ministerio de Ciencia, Innovación y Universidades grant PGC2018-102184-B-I00, cofunded by FEDER funds. LS acknowledges funding by DGAPA UNAM PAPIIT project IN-101819. GR-L acknowledges support from CONACYT grant 263373 and PRODEP (Mexico). This work is based on observations made with the Nordic Optical Telescope, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos in La Palma, Spain, of the Instituto de Astrofisica de Canarias. This work uses public data from the IR telescope Spitzer Space Telescope through the NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory at the California Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA). WISE is a joint project of the University of California (Los Angeles, USA) and the JPL/Caltech. The Infrared Astronomical Satellite (IRAS) was a joint project of the US, UK, and the Netherlands. This research is based on observations with AKARI, a JAXA project with the participation of ESA. This work has make extensive use of the NASA’s Astrophysics Data System., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

Aims.The mechanism behind the long secondary period (LSP) observed in pulsating red giants remains unknown. In this work, I investigate the connection between the red giant branch and the asymptotic giant branch evolution and the appearance of the LSP – a phenomenon that has been observed in a significant proportion of red giants.Methods.For the purposes of this work, I used the OGLE-III sample of the OSARG variables in the Large Magellanic Cloud. I constructed density maps in the period-luminosity as well as color-magnitude planes for the stars showing an LSP and compared them to the remaining giants. I also fit the spectral energy distribution to test whether an additional source of reddening is present in the LSP stars.Results.Here, I posit the hypothesis that the LSP phenomenon may be related to a transition between the different pulsation period-luminosity sequences. I also show that an overabundance of the stars showing LSP can be observed around the tip of the red giant branch, and much more prominently, at the upper part of the asymptotic giant branch. The main over-density region appears to be slightly fainter and redder than the bulk of the asymptotic giant branch. It also seems to correspond to the area of the Hertzsprung-Russell diagram where stable winds and high mass loss are present.Conclusions.The LSP can possibly be shown to be a recurring phenomenon that appears and disappears at various points of the red giant evolution. The LSP stars appear to be more reddened than other giants, which suggests the intrinsic nature of the reddening is likely to be related to large dust emission. The analysis appears to confirm the hypothesis that there is a relation between the mass loss due to the presence of strong stellar wind and the appearance of an LSP.

[Context]: Gamma-ray binaries are systems that emit nonthermal radiation peaking at energies above 1 MeV. One proposed scenario to explain their emission consists of a pulsar orbiting a massive star, with particle acceleration taking place in shocks produced by the interaction of the stellar and pulsar winds. [Aims]: We develop a semi-Analytical model of the nonthermal emission of the colliding-wind structure, which includes the dynamical effects of orbital motion. We apply the model to a general case and to LS 5039. [Methods]: The model consists of a one-dimensional emitter, the geometry of which is affected by Coriolis forces owing to orbital motion. Two particle accelerators are considered: one at the two-wind standoff location and the other one at the turnover produced by the Coriolis force. Synchrotron and inverse Compton emission is studied taking into account Doppler boosting and absorption processes associated to the massive star. [Results]: If both accelerators are provided with the same energy budget, most of the radiation comes from the region of the Coriolis turnover and beyond, up to a few orbital separations from the binary system. Significant orbital changes of the nonthermal emission are predicted in all energy bands. The model allows us to reproduce some of the LS 5039 emission features, but not all of them. In particular, the MeV radiation is probably too high to be explained by our model alone, the GeV flux is recovered but not its modulation, and the radio emission beyond the Coriolis turnover is too low. The predicted system inclination is consistent with the presence of a pulsar in the binary. [Conclusions]: The model is quite successful in reproducing the overall nonthermal behavior of LS 5039. Some improvements are suggested to better explain the phenomenology observed in this source, such as accounting for particle reacceleration beyond the Coriolis turnover, unshocked pulsar wind emission, and the three-dimensional extension of the emitter., We acknowledge support by the Spanish Ministerio de Economía y Competitividad (MINECO/FEDER, UE) under grant AYA2016-76012-C3-1-P, with partial support by the European Regional Development Fund (ERDF/FEDER), and from the Catalan DEC grant 2017 SGR 643. EM acknowledges support from MINECO through grant BES-2016-076342.

28 pags., 10 figs., 7 tabs., Context. This is the first publication from the DEATHSTAR project. The overall goal of the project is to reduce the uncertainties of the observational estimates of mass-loss rates from evolved stars on the Asymptotic Giant Branch (AGB). Aim. The aim in this first publication is to constrain the sizes of the 12CO emitting region from the circumstellar envelopes around 42 mostly southern AGB stars, of which 21 are M-Type and 21 are C-Type, using the Atacama Compact Array (ACA) at the Atacama Large Millimeter/submillimeter Array. The symmetry of the outflows is also investigated. Methods. Line emission from 12CO J = 2→1 and 3→2 from all of the sources were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. A detailed radiative transfer analysis will be presented in a future publication. The major and minor axis of the best-fit Gaussian at the line center velocity of the 12CO J = 2→1 emission gives a first indication of the size of the emitting region. Furthermore, the fitting results, such as the Gaussian major and minor axis, center position, and the goodness of fit across both lines, constrain the symmetry of the emission distribution. For a subsample of sources, the measured emission distribution is compared to predictions from previous best-fit radiative transfer modeling results. Results. We find that the CO envelope sizes are, in general, larger for C-Type than for M-Type AGB stars, which is as expected if the CO/H2 ratio is larger in C-Type stars. Furthermore, the measurements show a relation between the measured (Gaussian) 12CO J = 2→1 size and circumstellar density that, while in broad agreement with photodissociation calculations, reveals large scatter and some systematic differences between the different stellar types. For lower mass-loss-rate irregular and semi-regular variables of both M-and C-Type AGB stars, the 12CO J = 2→1 size appears to be independent of the ratio of the mass-loss rate to outflow velocity, which is a measure of circumstellar density. For the higher mass-loss-rate Mira stars, the 12CO J = 2→1 size clearly increases with circumstellar density, with larger sizes for the higher CO-Abundance C-Type stars. The M-Type stars appear to be consistently smaller than predicted from photodissociation theory. The majority of the sources have CO envelope sizes that are consistent with a spherically symmetric, smooth outflow, at least on larger scales. For about a third of the sources, indications of strong asymmetries are detected. This is consistent with what was found in previous interferometric investigations of northern sources. Smaller scale asymmetries are found in a larger fraction of sources. Conclusions. These results for CO envelope radii and shapes can be used to constrain detailed radiative transfer modeling of the same stars so as to determine mass-loss rates that are independent of photodissociation models. For a large fraction of the sources, observations at higher spatial resolution will be necessary to deduce the nature and origin of the complex circumstellar dynamics revealed by our ACA observations., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.2.00025.S and #2017.1.00595.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www. cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/ consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

Winds from young massive stars contribute a large amount of energy to their host molecular clouds. This has consequences for the dynamics and observable structure of star-forming clouds. In this paper, we present radiative magnetohydrodynamic simulations of turbulent molecular clouds that form individual stars of 30, 60 and 120 solar masses emitting winds and ultraviolet radiation following realistic stellar evolution tracks. We find that winds contribute to the total radial momentum carried by the expanding nebula around the star at 10 % of the level of photoionisation feedback, and have only a small effect on the radial expansion of the nebula. Radiation pressure is largely negligible in the systems studied here. The 3D geometry and evolution of wind bubbles is highly aspherical and chaotic, characterised by fast-moving "chimneys" and thermally-driven "plumes". These plumes can sometimes become disconnected from the stellar source due to dense gas flows in the cloud. Our results compare favourably with the findings of relevant simulations, analytic models and observations in the literature while demonstrating the need for full 3D simulations including stellar winds. However, more targeted simulations are needed to better understand results from observational studies., Comment: 18 pages, 11 figures, accepted to MNRAS

The recurrent fast radio burst FRB 180916 was recently shown to exhibit a 16 day period (with possible aliasing) in its bursting activity. Given magnetars as widely considered FRB sources, this period has been attributed to precession of the magnetar spin axis or the orbit of a binary companion. Here, we make the simpler connection to a {\it rotational period}, an idea observationally motivated by the 6.7 hour period of the Galactic magnetar candidate, 1E 161348--5055. We explore three physical mechanisms that could lead to the creation of ultra long period magnetars: (i) enhanced spin-down due to episodic mass-loaded charged particle winds (e.g. as may accompany giant flares), (ii) angular momentum kicks from giant flares and (iii) fallback leading to long lasting accretion disks. We show that particle winds and fallback accretion can potentially lead to a sub-set of the magnetar population with ultra long periods, sufficiently long to accommodate FRB 180916 or 1E 161348--5055. If confirmed, such periods implicate magnetars in relatively mature states (ages $1-10$ kyr) and which possessed large internal magnetic fields at birth $B_{\rm int}\gtrsim 10^{16}$ G. In the low-twist magnetar model for FRBs, such long period magnetars may dominate FRB production for repeaters at lower isotropic-equivalent energies and broaden the energy distribution beyond that expected for a canonical population of magnetars which terminate their magnetic activity at shorter periods $P \lesssim 10$ s., Comment: 12 pages, 5 figures. Accepted for publication in MNRAS

[Context] Photospheric radius expansion during X-ray bursts can be used to measure neutron star radii and help constrain the equation of state of neutron star matter. Understanding the stellar wind dynamics is important for interpreting observations, and therefore stellar wind models, though studied in past decades, have regained interest and need to be revisited with updated data and methods., [Aims] Here, we study the radiative wind model in the context of X-ray bursts with modern techniques and physics input. We focus on characterization of the solutions and the study of observable magnitudes as a function of free model parameters., [Methods] We implemented a spherically symmetric nonrelativistic wind model in a stationary regime, with updated opacity tables and modern numerical techniques. Total mass and energy outflows (Ṁ, Ė) were treated as free parameters., [Results] A high-resolution parameter-space exploration was performed to allow better characterization of observable magnitudes. High correlation was found between different photospheric magnitudes and free parameters. For instance, the photospheric ratio of gravitational energy outflow to radiative luminosity is directly proportional to the photospheric wind velocity., [Conclusions] The correlations found here could help determine the physical conditions of the inner layers, where nuclear reactions take place, by means of observable photospheric values. Further studies are needed to determine the range of physical conditions in which the correlations are valid., The authors would like to thank Arman Aryaeipour and Manuel Linares for fruitful discussions. This work has been partially supported by the Spanish MINECO grant AYA2017–86274–P, by the E.U. FEDER funds, and by the AGAUR/Generalitat de Catalunya grant SGR-661/2017. This article benefited from discussions within the “ChETEC” COST Action (CA16117).

MSc (Astrophysical Sciences), North-West University, Potchefstroom Campus In the 11 years since its launch, the Fermi Large Area Telescope has yielded a large number of new detections. Among these new detections are about 30 black widow and redback systems. These sources are compact binaries in which a millisecond pulsar heats and may even ablate its low-mass companion by its intense wind of relativistic particles and radiation. In such systems, an intrabinary shock can form as a site of particle acceleration and associated non-thermal emission. We model the X-ray and gamma-ray synchrotron and inverse-Compton spectral components for select spider binaries, including diffusion, convection, and radiative energy losses in an axially-symmetric, steady-state approach. Our new multi-zone code simultaneously yields energy-dependent light curves and orbital-phase-resolved spectra. Using parameter studies and matching the observed X-ray spectra and light curves, and Fermi Large Area Telescope spectra where available, with a synchrotron component, we can constrain certain model parameters. For PSR J1723–2837 these are notably the magnetic field and bulk flow speed of plasma moving along the shock tangent, the shock acceleration efficiency, and the multiplicity and spectrum of pairs accelerated by the pulsar. This affords a more robust prediction of the expected high-energy and very-high-energy gamma-ray flux. We find that nearby pulsars with hot or flaring companions may be promising targets for the future Cherenkov Telescope Array. Moreover, many spiders are likely to be of significant interest to future MeV-band missions such as AMEGO and e-ASTROGAM. Masters

We present a detailed comparative study of the arcs and fragmented ring-like features in the haloes of the planetary nebulae (PNe) NGC 6543, NGC 7009, and NGC 7027 and the spiral pattern around the carbon star AFGL 3068 using high-quality multi-epoch HST images. This comparison allows us to investigate the connection and possible evolution between the regular patterns surrounding AGB stars and the irregular concentric patterns around PNe. The radial proper motion of these features, ≃15 km s-1, are found to be consistent with the AGB wind and their linear sizes and inter-lapse times (500 - 1900 yr) also agree with those found around AGB stars, suggesting a common origin. We find evidence using radiative-hydrodynamic simulations that regular patterns produced at the end of the AGB phase become highly distorted by their interactions with the expanding PN and the anisotropic illumination and ionization patterns caused by shadow instabilities. These processes will disrupt the regular (mostly spiral) patterns around AGB stars, plausibly becoming the arcs and fragmented rings observed in the haloes of PNe.© 2020 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, The authors acknowledge support from grants AYA 2014-57280- P and PGC2018-102184-B-I00, co-funded with FEDER funds. JAT and MAG are funded by DGAPA, UNAM PAPIIT project IA100318, and LS by DGAPA, UNAM PAPIIT project IN101819. GRL acknowledges support from Fundacion Marcos Moshinsky, ´ Consejo Nacional de Ciencia y Tecnolog´ıa (CONACyT), and Programa para el Desarrollo Profesional Docente (PRODEP) (Mexico). BB acknowledges support from program AR-14563 that was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. This paper is based on observations made with the NASA/ESAHubble Space Telescope, and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA), and the Canadian Astronomy Data Centre (CADC/NRC/CSA). The large scale image of NGC 7009 is based on observations collected at the European Southern Observatory by the MUSE science verification team. This work has made extensive use of the NASA’s Astrophysics Data System.

High-mass X-ray binaries (HMXBs) are exceptional astrophysical laboratories that offer a rare glimpse into the physical processes that govern accretion on compact objects, massive-star winds, and stellar evolution. In a subset of the HMXBs, the compact objects accrete matter solely from winds of massive donor stars. These so-called wind-fed HMXBs are divided in persistent HMXBs and supergiant fast X-ray transients (SFXTs) according to their X-ray properties. While it has been suggested that this dichotomy depends on the characteristics of stellar winds, they have been poorly studied. With this investigation, we aim to remedy this situation by systematically analyzing donor stars of wind-fed HMXBs that are observable in the UV, concentrating on those with neutron star (NS) companions. We obtained Swift X-ray data, HST UV spectra, and additional optical spectra for all our targets. Our multi-wavelength approach allows us to provide stellar and wind parameters for six donor stars (four wind-fed systems and two OBe X-ray binaries). The wind properties are in line with the predictions of the line-driven wind theory. Three of the donor stars are in an advanced evolutionary stage, while for some of the stars, the abundance pattern indicates that processed material might have been accreted. When passing by the NS in its tight orbit, the donor star wind has not yet reached its terminal velocity but it is still significantly slower; its speed is comparable with the orbital velocity of the NS companion. There are no systematic differences between the two types of wind-fed HMXBs (persistent versus transients) with respect to the donor stars. For the SFXTs in our sample, the orbital eccentricity is decisive for their transient X-ray nature. Based on the orbital parameters and the further evolution of the donor stars, the investigated HMXBs will presumably form Thorne-\.Zytkow objects in the future., Comment: 44 pages, 13 figures, 15 tables; accepted for publication in A&A

We present results from radiation non-ideal magnetohydrodynamics (MHD) calculations that follow the collapse of rotating, magnetised, molecular cloud cores to stellar densities. These are the first such calculations to include all three non-ideal effects: ambipolar diffusion, Ohmic resistivity and the Hall effect. We employ an ionisation model in which cosmic ray ionisation dominates at low temperatures and thermal ionisation takes over at high temperatures. We explore the effects of varying the cosmic ray ionisation rate from $\zeta_\text{cr}= 10^{-10}$ to $10^{-16}$ s$^{-1}$. Models with ionisation rates $\gtrsim 10^{-12}$ s$^{-1}$ produce results that are indistinguishable from ideal MHD. Decreasing the cosmic ray ionisation rate extends the lifetime of the first hydrostatic core up to a factor of two, but the lifetimes are still substantially shorter than those obtained without magnetic fields. Outflows from the first hydrostatic core phase are launched in all models, but the outflows become broader and slower as the ionisation rate is reduced. The outflow morphology following stellar core formation is complex and strongly dependent on the cosmic ray ionisation rate. Calculations with high ionisation rates quickly produce a fast (~14km s$^{-1}$) bipolar outflow that is distinct from the first core outflow, but with the lowest ionisation rate a slower (~3-4 km s$^{-1}$) conical outflow develops gradually and seamlessly merges into the first core outflow., Comment: 22 pages, 24 figures. Accepted for publication in MNRAS

Aims: AT 2018hso is a new transient showing transitional properties between those of LRNe and the class of intermediate luminosity red transients (ILRTs) similar to SN 2008S. Through the detailed analysis of the observed parameters, our study support that it actually belongs to the LRN class, and was likely produced by the coalescence of two massive stars. Methods: We obtained ten months of optical and near-infrared photometric monitoring, and eleven epochs of low-resolution optical spectroscopy of AT~2018hso. We compared its observed properties with those of other ILRTs and LRNe. We also inspected the archival Hubble Space Telescope (HST) images obtained about 15 years ago to constrain the progenitor's properties. Results: The light curves of AT 2018hso show a first sharp peak (Mr = -13.93 mag), followed by a broader and shallower second peak, that resembles a plateau in the optical bands. The spectra dramatically change with time. Early time spectra show prominent Balmer emission lines and a weak Ca II] doublet, which is usually observed in ILRTs. However, the major decrease in the continuum temperature, the appearance of narrow metal absorption lines, the major change in the H$\alpha$ strength and profile, and the emergence of molecular bands support a LRN classification. The possible detection of an I ~ -8 mag source at the position of AT 2018hso in HST archive images is consistent with expectations for a pre-merger massive binary, similar to the precursor of the 2015 LRN in M101. Conclusions: We provide reasonable arguments to support a LRN classification for AT~2018hso. This study reveals growing heterogeneity in the observables of LRNe than thought in the past, making sometimes tricky the discrimination between LRNe and ILRTs. This suggests the need of monitoring the entire evolution of gap transients to avoid misclassifications., Comment: 10 pages, 6 figures, 3 tables. Submitted to A&A Letter

著者人数: 52名, Accepted: 2019-06-05, 資料番号: SA1190104000

The Small Magellanic Cloud multiple system HD 5980 contains a luminous blue variable (LBV) that underwent a major eruption in 1994, and whose current spectrum is that of a hydrogen-rich Wolf–Rayet (WR) star. Since the eruption, the wind mass-loss rate has been declining while wind speeds have been steadily increasing. Observations obtained in 2014 when Star A (the LBV) eclipses Star B indicate that the fitted mass-loss rate and luminosity have reached the lowest values ever determined for such spectra: Ṁ = 4.5 × 10⁻⁵ M⊙ yr⁻¹, L = 1.7 × 10⁶ L⊙. In addition, the radius of the LBV’s continuum-emitting region is similar to that derived from the eclipse light curves of the late 1970s. Hence, it appears to have attained a similar ‘low’ state to that of the late 1970s. While a good fit to the emission spectrum is obtained using a cmfgen model, there are discrepancies in the UV. In particular, the extent of the observed absorption profiles is ∼1000 km s⁻¹ greater than predicted by the emission-line intensities. Further, HST UV observations obtained in 2016, when Star A is eclipsed by Star B, show unusual P Cygni profiles that are not easily explained. Surprisingly the 2016 emission-line spectrum is similar to that at the opposite eclipse obtained in 2014. The complex UV profiles are likely to arise as a consequence of the dynamics of the wind–wind collision and radiative braking, both of which will cause significant departures from spherical symmetry, and have a strong orbital phase dependence. However, other scenarios, such as intrinsically aspherical winds, cannot be ruled out., Facultad de Ciencias Astronómicas y Geofísicas

We present results from the MOSFIRE Deep Evolution Field survey on broad flux from the nebular emission lines Hα, [N ii], [O iii], Hβ, and [S ii]. The sample consists of 127 star-forming galaxies at 1.37

We present the analysis of infrared (IR) observations of the planetary nebula NGC 40 together with spectral analysis of its [WC]-type central star HD 826. Spitzer IRS observations were used to produce spectral maps centred at polycyclic aromatic hydrocarbons (PAH) bands and ionic transitions to compare their spatial distribution. The ionic lines show a clumpy distribution of material around the main cavity of NGC 40, with the emission from [Ar II] being the most extended, whilst the PAHs show a rather smooth spatial distribution. Analysis of ratio maps shows the presence of a toroidal structure mainly seen in PAH emission, but also detected in a Herschel PACS 70 μm image. We argue that the toroidal structure absorbs the UV flux from HD 826, preventing the nebula to exhibit lines of high-excitation levels as suggested by previous authors. We discuss the origin of this structure and the results from the spectral analysis of HD 826 under the scenario of a late thermal pulse.© 2019 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, JAT, MAG and HT are funded by UNAM DGAPA PAPIIT project IA100318. GRL acknowledges support from Fundacion Marcos Moshinsky, CONACyT and PRODEP (Mexico). MAG acknowledges support of the grant AYA 2014-57280-P, cofunded with FEDER funds.

The interacting binary Eta Carinae remains one of the most enigmatic massive stars in our Galaxy despite over four centuries of observations. In this work, its light curve from the ultraviolet to the near-infrared is analysed using spatially resolved HST observations and intense monitoring at the La Plata Observatory, combined with previously published photometry. We have developed a method to separate the central stellar object in the ground-based images using HST photometry and applying it to the more numerous ground-based data, which supports the hypothesis that the central source is brightening faster than the almost-constant Homunculus. After detrending from long-term brightening, the light curve shows periodic orbital modulation ($\Delta V$ $\sim$ 0.6 mag) attributed to the wind-wind collision cavity as it sweeps around the primary star and it shows variable projected area to our line-of-sight. Two quasi-periodic components with time scales of 2-3 and 8-10 yr and low amplitude, $\Delta V$ $, Comment: 26 pages, 19 figures. Accepted to appear on MNRAS

Classical T Tauri stars with ages of less than 10 Myr possess accretion discs. Magnetohydro-dynamic processes at the boundary between the disc and the stellar magnetosphere control the accretion and ejections gas flows. We carried out a long series of simultaneous spectroscopic and photometric observations of the classical T Tauri stars, RY Tauri and SU Aurigae, with the aim to quantify the accretion and outflow dynamics at time-scales from days to years. It is shown that dust in the disc wind is the main source of photometric variability of these stars. In RY Tau, we observed a new effect: during events of enhanced outflow, the circumstellar extinction becomes lower. The characteristic time of changes in outflow velocity and stellar brightness indicates that the obscuring dust is near the star. The outflow activity in both stars is changing on a time-scale of years. Periods of quiescence in the variability of the Hα profile were observed during the 2015-2016 period in RY Tau and during the 2016-2017 period in SU Aur. We interpret these findings in the framework of the magnetospheric accretion model, and we discuss how the global stellar magnetic field can influence the long-term variations of the outflow activity.© 2018 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society, This research was supported by the RFBR grant 16-02-00140. It was based partly on observations made with the Nordic Optical Telescope, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias, using ALFOSC, which is provided by the Instituto de Astrofisica de Andalucia (IAA) under a joint agreement with the University of Copenhagen and NOTSA. We acknowledge support from a Visitors Programme at the Department of Astronomy at Stockholm University. JFG and RMGA were supported by Fundacao para a Ciencia e a Tecnologia (FCT) through national funds (UID/FIS/04434/2013) and by FEDER through COMPETE2020 (POCI-01-0145-FEDER-007672). RMGA is supported by the fellowship PD/BD/113745/2015, under the FCT PD Program PhD::SPACE, funded by FCT (Portugal) and POPH/FSE (EC) and by CRUP through a cooperation programme (PAULIF: TC-16/17). DEM acknowledges his work as part of the research activity supported by the National Astronomical Research Institute of Thailand (NARIT), Ministry of Science and Technology of Thailand. SYuG was supported in part by the Ministry of Education and Science of Russia (the basic part of the Stateassignment, RK no. AAAA-A17-117030310283-7) and by the Act no. 211 of the Government of the Russian Federation, agreement No 02.A03.21.0006.

Despite the first detection of X-rays from massive stars occurring four decades ago, the physical dependence of the emission mechanism(s) on the underlying stellar and binary properties of the emitters remains uncertain. The young massive cluster Westerlund 1 provides an ideal testbed for understanding this phenomenon, 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. A total of 45 X-ray sources within Wd1 now have precise spectral classifications. These observations identify 16 candidate and confirmed massive binaries; by comparison 22 stars have X-ray properties that imply a contribution from a wind collision zone. X-ray emission appears confined to O9-B0.5 supergiants, Wolf-Rayets and a small group of highly luminous interacting/post-interaction binaries. Despite their presence in large numbers, no emission is seen from earlier, less evolved O stars or later, cooler B super-/hypergiants. We suppose that the lack of X-ray emission from O giants is due to their comparatively low bolometric luminosities if, as expected, they follow the canonical L_X/L_bol 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 and 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 (Abridged)., Comment: 10 pages, 2 figures, accepted for publication in Astronomy and Astrophysics

We obtained high-resolution spectra and multicolor photometry of V1082 Sgr to study the donor star in this 20.8 hr orbital period binary, which is assumed to be a detached system. We measured the rotational velocity (v sin i = 26.5+/-2.0 km/s), which, coupled with the constraints on the white dwarf mass from the X-ray spectroscopy, leads to the conclusion that the donor star barely fills 70% of its corresponding Roche lobe radius. It appears to be a slightly evolved K2-type star. This conclusion was further supported by a recently published distance to the binary system measured by the Gaia mission. At the same time, it becomes difficult to explain a very high (> 10e-9 Msun/yr) mass transfer and mass accretion rate in a detached binary via stellar wind and magnetic coupling., 16 pages, 10 figures, accepted for publication in Ap.J

Reliable predictions of mass-loss rates are important for massive-star evolution computations. We aim to provide predictions for mass-loss rates and wind-momentum rates of O-type stars, carefully studying the behaviour of these winds as functions of stellar parameters like luminosity and metallicity. We use newly developed steady-state models of radiation-driven winds to compute the global properties of a grid of O-stars. The self-consistent models are calculated by means of an iterative solution to the equation of motion using full NLTE radiative transfer in the co-moving frame to compute the radiative acceleration. In order to study winds in different galactic environments, the grid covers main-sequence stars, giants and supergiants in the Galaxy and both Magellanic Clouds. We find a strong dependence of mass-loss on both luminosity and metallicity. Mean values across the grid are $\dot{M}\sim L_{\ast}^{2.2}$ and $\dot{M}\sim Z_{\ast}^{0.95}$, however we also find a somewhat stronger dependence on metallicity for lower luminosities. Similarly, the mass loss-luminosity relation is somewhat steeper for the SMC than for the Galaxy. In addition, the computed rates are systematically lower (by a factor 2 and more) than those commonly used in stellar-evolution calculations. Overall, our results agree well with observations in the Galaxy that account properly for wind-clumping, with empirical $\dot{M}$ vs. $Z_\ast$ scaling relations, and with observations of O-dwarfs in the SMC. Our results provide simple fit relations for mass-loss rates and wind momenta of massive O-stars stars as functions of luminosity and metallicity, valid in the range $T_{\rm eff} = 28000 - 45000$\,K. Due to the systematically lower $\dot{M}$, our new models suggest that new rates might be needed in evolution simulations of massive stars., Accepted for publication in A&A. 16 pages, 13 figures

We present the results of our monitoring campaigns of the luminous red novae (LRNe) AT 2020hat in NGC 5068 and AT 2020kog in NGC 6106. The two objects were imaged (and detected) before their discovery by routine survey operations. They show a general trend of slow luminosity rise, lasting at least a few months. The subsequent major LRN outbursts were extensively followed in photometry and spectroscopy. The light curves present an initial short-duration peak, followed by a redder plateau phase. AT 2020kog is a moderately luminous event peaking at ∼7 × 1040 erg s−1, while AT 2020hat is almost one order of magnitude fainter than AT 2020kog, although it is still more luminous than V838 Mon. In analogy with other LRNe, the spectra of AT 2020kog change significantly with time. They resemble those of type IIn supernovae at early phases, then they become similar to those of K-type stars during the plateau, and to M-type stars at very late phases. In contrast, AT 2020hat already shows a redder continuum at early epochs, and its spectrum shows the late appearance of molecular bands. A moderate-resolution spectrum of AT 2020hat taken at +37 d after maximum shows a forest of narrow P Cygni lines of metals with velocities of 180 km s−1, along with an Hα emission with a full-width at half-maximum velocity of 250 km s−1. For AT 2020hat, a robust constraint on its quiescent progenitor is provided by archival images of the Hubble Space Telescope. The progenitor is clearly detected as a mid-K type star, with an absolute magnitude of MF606W = −3.33 ± 0.09 mag and a colour of F606W − F814W = 1.14 ± 0.05 mag, which are inconsistent with the expectations from a massive star that could later produce a core-collapse supernova. Although quite peculiar, the two objects nicely match the progenitor versus light curve absolute magnitude correlations discussed in the literature.

Context. Red supergiants are observed to undergo vigorous mass-loss. However, to date, no theoretical model has succeeded in explaining the origins of these objects' winds. This strongly limits our understanding of red supergiant evolution and Type II-P and II-L supernova progenitor properties. Aims. We examine the role that vigorous atmospheric turbulence may play in initiating and determining the mass-loss rates of red supergiant stars. Methods. We analytically and numerically solve the equations of conservation of mass and momentum, which we later couple to an atmospheric temperature structure, to obtain theoretically motivated mass-loss rates. We then compare these to state-of-the-art empirical mass-loss rate scaling formulae as well as observationally inferred mass-loss rates of red supergiants. Results. We find that the pressure due to the characteristic turbulent velocities inferred for red supergiants is sufficient to explain the mass-loss rates of these objects in the absence of the normally employed opacity from circumstellar dust. Motivated by this initial success, we provide a first theoretical and fully analytic mass-loss rate prescription for red supergiants. We conclude by highlighting some intriguing possible implications of these rates for future studies of stellar evolution, especially in light of the lack of a direct dependence on metallicity., Comment: 14 pages, 9 figures, 2 tables

Context: The observed orbital characteristics of post-asymptotic giant branch (post-AGB) and post-red giant branch (post-RGB) binaries are not understood. We suspect that the missing ingredients to explain them probably lie in the continuous interaction of the central binary with its circumstellar environment. Aims: We aim at studying the circumbinary material in these complex systems by investigating the connection between the innermost and large-scale structures. Methods: We perform high-angular resolution observations in the near-infrared continuum of HD101584, which has a complex structure as seen at millimeter wavelengths with a disk-like morphology and a bipolar outflow due to an episode of strong binary interaction. To account for the complexity of the target we first perform an image reconstruction and use this result to fit a geometrical model to extract the morphological and thermal features of the environment. Results: The image reveals an unexpected double-ring structure. We interpret the inner ring to be produced by emission from dust located in the plane of the disk and the outer ring to be produced by emission from dust that is located 1.6[D/1kpc] au above the disk plane. The inner ring diameter (3.94[D/1kpc] au), and temperature (T=1540$\pm$10K) are compatible with the dust sublimation front of the disk. The origin of the out-of-plane ring (with a diameter of 7.39[D/1kpc] au and a temperature of 1014$\pm10$K) could be due to episodic ejection or a dust condensation front in the outflow. Conclusion: The observed outer ring is possibly linked with the blue-shifted side of the large scale outflow seen by ALMA and is tracing its launching location to the central star. Such observations give morphological constraints on the ejection mechanism. Additional observations are needed to constrain the origin of the out-of-plane structure., Comment: Accepted to A&A. 14 pages, 13 figures

Astronomy and astrophysics 635, A144 (2020). doi:10.1051/0004-6361/201937031, The binary system η Carinae is a unique laboratory that facilitates the study of particle acceleration to high energies under a wide range of conditions, including extremely high densities around periastron. To date, no consensus has emerged as to the origin of the gigaelectronvolt γ-ray emission in this important system. With a re-analysis of the full Fermi-LAT data set for η Carinae, we show that the spectrum is consistent with a pion decay origin. A single population leptonic model connecting X-ray to γ-ray emission can be ruled out. We revisit our physical model from 2015, based on two acceleration zones associated with the termination shocks in the winds of both stars. We conclude that inverse Compton emission from in-situ accelerated electrons dominates the hard X-ray emission detected with NuSTAR at all phases away from periastron and that pion-decay from shock accelerated protons is the source of γ-ray emission. Very close to periastron there is a pronounced dip in hard X-ray emission, concomitant with the repeated disappearance of the thermal X-ray emission, which we interpret as due to the suppression of significant electron acceleration in the system. Within our model, the residual emission seen by NuSTAR at this phase can be accounted for with secondary electrons produced in interactions of accelerated protons, which agrees with the variation in pion-decay γ-ray emission. Future observations with H.E.S.S., CTA, and NuSTAR should confirm or refute this scenario.Key words: radiation mechanisms: non-thermal / acceleration of particles / gamma-rays: stars / stars: individual: η Car / stars: winds / outflows / X-rays: binaries, Published by EDP Sciences, Les Ulis