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.

Detailed knowledge about stellar winds and evolution at different metallicities is crucial for understanding stellar populations and feedback in the Local Group of galaxies and beyond. Despite efforts in the literature, we still lack a comprehensive, empirical view of the dependence of wind properties on metallicity ($Z$). Here, we investigate the winds of O and B stars in the Milky Way (MW) and Small Magellanic Cloud (SMC). We gathered a sample of 96 stars analyzed by means of the NLTE code CMFGEN. We explored their wind strengths and terminal velocities to address the $Z$ dependence, over a large luminosity range. The empirical wind-luminosity relation (WLR) obtained updates and extends previous results in the literature. It reveals a luminosity and $Z$ dependence, in agreement with the radiatively driven wind theory. For bright objects ($\log L/L_\odot \gtrsim 5.4$), we infer that $\dot{M} \sim Z^{0.5-0.8}$. However, this dependence seems to get weaker or vanish at lower luminosities. The analysis of the terminal velocities suggests a shallow $Z^n$ dependence, with $n \sim 0.1-0.2$, but it should be confirmed with a larger sample and more accurate $V_{\infty}$ determinations. Recent results on SMC stars based on the PoWR code support our inferred WLR. On the other hand, recent bow-shocks measurements stand mostly above our derived WLR. Theoretical calculations of the WLR are not precise, specially at low $L$, where the results scatter. Deviations between our results and recent predictions are identified to be due to the weak wind problem and the extreme terminal velocities predicted by the models. The Z dependence suggested by our analysis deserves further investigations, given its astrophysical implications., 17 Pages, 10 Figures and 3 Tables, accepted for publication in MNRAS

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

V906 Carinae was one of the best observed novae of recent times. It was a prolific dust producer and harboured shocks in the early evolving ejecta outflow. Here, we take a close look at the consequences of these early interactions through study of high-resolution Ultraviolet and Visual Echelle spectrograph spectroscopy of the nebular stage and extrapolate backwards to investigate how the final structure may have formed. A study of ejecta geometry and shaping history of the structure of the shell is undertaken following a spectral line $\rm {\small SHAPE}$ model fit. A search for spectral tracers of shocks in the nova ejecta is undertaken and an analysis of the ionized environment. Temperature, density, and abundance analyses of the evolving nova shell are presented.

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.

Context. In accelerating and supersonic media, understanding the interaction of photons with spectral lines can be of utmost importance, especially in an accelerating flow. However, fully accounting for such line forces is computationally expensive and challenging, as it involves complicated solutions of the radiative transfer problem for millions of contributing lines. This currently can only be done by specialised codes in 1D steady-state flows. More general cases and higher dimensions require alternative approaches. Aims. We present a comprehensive and fast method for computing the radiation line force using tables of spectral-line-strength distribution parameters, which can be applied in arbitrary (multi-D, time-dependent) simulations, including those that account for the line-deshadowing instability, to compute the appropriate opacities. Methods. We assume local thermodynamic equilibrium to compute a flux-weighted line opacity from ~4 million spectral lines. We fit the opacity computed from the line list with an analytic result derived for an assumed distribution of the spectral line strength and found the corresponding line-distribution parameters, which we tabulate here for a range of assumed input densities ρ ∈ [10−20, 10−10] g cm−3 and temperatures T ∊ [104, 1047] K. Results. We find that the variation in the line-distribution parameters plays an essential role in setting the wind dynamics in our models. In our benchmark study, we also find a good overall agreement between the O-star mass-loss rates of our models and those derived from steady-state studies that use a more detailed radiative transfer. Conclusions. Our models reinforce the idea that self-consistent variation in the line-distribution parameters is important for the dynamics of line-driven flows. Within a well-calibrated O-star regime, our results support the proposed methodology. In practice, utilising the provided tables, yielded a factor >100 speed-up in computational time compared to specialised 1D model-atmosphere codes of line-driven winds, which constitutes an important step towards efficient multi-dimensional simulations. We conclude that our method and tables are ready to be exploited in various radiation-hydrodynamic simulations where the line force is important.

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

Disk winds are an important mechanism for accretion and disk evolution around young stars. The accreting intermediate-mass T-Tauri star RY Tau has an active jet and a previously known disk wind. Archival optical and new near-infrared observations of the RY Tau system show two horn-like components stretching out as a cone from RY Tau. Scattered light from the disk around RY Tau is visible in near-infrared but not seen at optical wavelengths. In the near-infrared, dark wedges that separates the horns from the disk, indicating we may see the scattered light from a disk wind. We use archived ALMA and SPHERE/ZIMPOL I-band observations combined with newly acquired SPEHRE/IRDIS H-band observations and available literature to build a simple geometric model of the RY Tau disk and disk wind. We use Monte Carlo radiative transfer modelling \textit{MCMax3D} to create comparable synthetic observations that test the effect of a dusty wind on the optical effect in the observations. We constrain the grain size and dust mass needed in the disk wind to reproduce the effect from the observations. A model geometrically reminiscent of a dusty disk wind with small micron to sub-micron size grains elevated above the disk can reproduce the optical effect seen in the observations. The mass in the obscuring component of the wind has been constrained to $1\times10^{-9} M_{\odot} \leq M \leq 5\times10^{-8} M_{\odot}$ which corresponds to a lower limit mass loss rate in the wind of about $\sim 1\times10^{-8}M_{\odot}\mathrm{yr}^{-1}$. While an illuminate dust cavity cannot be ruled out without measurements of the gas velocity, we argue that a magnetically launched disk wind is the most likely scenario., 11 pages, 5 figures

Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decreasing metal content mimicking the environment at the earliest epochs. Ultimately, a physical model that could be extrapolated to zero metallicity would enable tackling long-standing questions such as "What did the First, very massive stars of the Universe look like?" or "What was their role in the re-ionization of the Universe?". Yet, most our knowledge of metal-poor massive stars is drawn from one single point in metallicity. Massive stars in the Small Magellanic Cloud (SMC, $\sim 1/5 Z_{\odot}$) currently serve as templates for low-metallicity objects in the early Universe, even though significant differences with respect to massive stars with poorer metal content have been reported. This White Paper summarizes the current knowledge on extremely (sub-SMC) metal poor massive stars, highlighting the most outstanding open questions and the need to supersede the SMC as standard. A new paradigm can be built from nearby extremely metal-poor galaxies that make a new metallicity ladder, but massive stars in these galaxies are out of reach to current observational facilities. Such task would require an L-size mission, consisting of a 10m-class space telescope operating in the optical and the ultraviolet ranges. Alternatively, we propose that ESA unites efforts with NASA to make the LUVOIR mission concept a reality, thus continuing the successful partnership that made Hubble Space Telescope one of the greatest observatories of all time., Comment: A white paper submitted for the Voyage 2050 long-term plan in the ESA Science Programme. 21 pages, 1 table, 6 figures. arXiv admin note: substantial text overlap with arXiv:1903.05235

The properties and the evolution of asymptotic giant branch (AGB) stars are strongly influenced by their mass loss through a stellar wind. This is believed to be caused by radiation pressure due to the absorption and scattering of the stellar radiation by the dust grains formed in the atmosphere. The optical properties of dust are often estimated using the small particle limit (SPL) approximation, and it has been used frequently in modelling AGB stellar winds when performing radiation-hydrodynamics (RHD) simulations. We aim to investigate the effects of replacing the SPL approximation by detailed Mie calculations of the size-dependent opacities for grains of amorphous carbon forming in C-rich AGB star atmospheres. We performed RHD simulations for a large grid of carbon star atmosphere+wind models with different effective temperatures, luminosities, stellar masses, carbon excesses, and pulsation properties. Also, a posteriori radiative transfer calculations for many radial structures (snapshots) of these models were done, resulting in spectra and filter magnitudes. We find, when giving up the SPL approximation, the wind models become more strongly variable and more dominated by gusts, although the average mass-loss rates and outflow speeds do not change significantly; the increased radiative pressure on the dust throughout its formation zone does, however, result in smaller grains and lower condensation fractions (and thus higher gas-to-dust ratios). The photometric K magnitudes are generally brighter, but at V the effects of using size-dependent dust opacities are more complex: brighter for low mass-loss rates and dimmer for massive stellar winds. Given the large effects on spectra and photometric properties, it is necessary to use the detailed dust optical data instead of the simple SPL approximation in stellar atmosphere+wind modelling where dust is formed., Comment: 14 pages, 24 figures. Accepted for publication in Astronomy \& Astrophysics

The temporal and spectral evolution of gamma-ray burst (GRB) afterglows can be used to infer the density and density profile of the medium through which the shock is propagating. In long-duration (core-collapse) GRBs, the circumstellar medium (CSM) is expected to resemble a wind-blown bubble, with a termination shock, separating the stellar wind and the interstellar medium (ISM). A long standing problem is that flat density profiles, indicative of the ISM, are often found at lower radii than expected for a massive star progenitor. Furthermore, the presence of both wind-like environments at high radii and ISM-like environments at low radii remains a mystery. In this paper, we perform a ‘CSM population synthesis’ with long GRB progenitor stellar evolution models. Analytic results for the evolution of wind blown bubbles are adjusted through comparison with a grid of 2D hydrodynamical simulations. Predictions for the emission radii, ratio of ISM to wind-like environments, wind, and ISM densities are compared with the largest sample of afterglow derived parameters yet compiled, which we make available for the community. We find that high ISM densities of n ∼ 1000 cm−3 best reproduce observations. If long GRBs instead occur in typical ISM densities of n ∼ 1 cm−3, then the discrepancy between theory and observations is shown to persist at a population level. We discuss possible explanations for the origin of variety in long GRB afterglows, and for the overall trend of CSM modelling to over-predict the termination shock radius. © 2022 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society., BPG and MN are supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 948381). DAK acknowledges support from Spanish National Research Project RTI2018-098104-J-I00 (GRBPhot). AJvM is supported by the ANR-19-CE31-0014GAMALO project. AJL has received funding from the European Research Council (ERC) under the European Union’s Seventh Framework Programme (FP7-2007-2013) (Grant agreement No. 725246). ERS has been supported by STFC consolidated grant ST/P000495/1. PJG is supported by NRF SARChI Grant 111692. We gratefully acknowledge the use of GOTOHEAD, the computing cluster of the Gravitational-wave Optical Transient Observer (GOTO), as well as support from Joe Lyman and Krzysztof Ulaczyk., With funding from the Spanish government through the Severo Ochoa Centre of Excellence accreditation SEV-2017-0709.

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

Close binary evolution is widely invoked to explain the formation of axisymmetric planetary nebulae, after a brief common envelope phase. The evolution of the primary would be interrupted abruptly, its still quite massive envelope being fully ejected to form the PN, which should be more massive than a planetary nebula coming from the same star, were it single. We test this hypothesis by investigating the ionised and molecular masses of a sample consisting of 21 post-common-envelope planetary nebulae, roughly one fifth of their known total population, and comparing them to a large sample of regular planetary nebulae (not known to host close-binaries). We find that post-common-envelope planetary nebulae arising from single-degenerate systems are, on average, neither more nor less massive than regular planetary nebulae, whereas post-common-envelope planetary nebulae arising from double-degenerate systems are considerably more massive, and show substantially larger linear momenta and kinetic energy than the rest. Reconstruction of the common envelope of four objects further suggests that the mass of single-degenerate nebulae actually amounts to a very small fraction of the envelope of their progenitor stars. This leads to the uncomfortable question of where the rest of the envelope is, raising serious doubts on our understanding of these intriguing objects.

We present spectroscopic observations of the Be/X-ray binary X Per obtained during the period December 2017 - January 2020 (MJD~58095 - MJD~58865). In December 2017 the $H\alpha$, $H\beta$, and HeI 6678 emission lines were symmetric with violet-to-red peak ratio $V/R \approx 1$. During the first part of the period (December 2017 - August 2018) the V/R-ratio decreased to 0.5 and the asymmetry developed simultaneously in all three lines. In September 2018, a third component with velocity $\approx 250$~km~s$^{-1}$ appeared on the red side of the HeI line profile. Later this component emerged in $H\beta$, accompanied by the appearance of a red shoulder in $H\alpha$. Assuming that it is due to an eccentric wave in the circumstellar disc, we find that the eccentric wave appeared first in the innermost part of the disc, it spreads out with outflowing velocity $v_{wave} \approx 1.1 \pm 0.2 $~km~s$^{-1}$, and the eccentricity of the eccentric wave is $e_{wave} \approx 0.29 \pm 0.07$. A detailed understanding of the origin of such eccentricities would have applications to a wide range of systems from planetary rings to AGNs., Comment: MNRAS (accepted)

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

Recent solar physics missions have shown the definite role of waves and magnetic fields deep in the inner corona, at the chromosphere-corona interface, where dramatic and physically dominant changes occur. HiRISE (High Resolution Imaging and Spectroscopy Explorer), the ambitious new generation ultra-high resolution, interferometric, and coronagraphic, solar physics mission, proposed in response to the ESA Voyage 2050 Call, would address these issues and provide the best-ever and most complete solar observatory, capable of ultra-high spatial, spectral, and temporal resolution observations of the solar atmosphere, from the photosphere to the corona, and of new insights of the solar interior from the core to the photosphere. HiRISE, at the L1 Lagrangian point, would provide meter class FUV imaging and spectro-imaging, EUV and XUV imaging and spectroscopy, magnetic fields measurements, and ambitious and comprehensive coronagraphy by a remote external occulter (two satellites formation flying 375 m apart, with a coronagraph on a chaser satellite). This major and state-of-the-art payload would allow us to characterize temperatures, densities, and velocities in the solar upper chromosphere, transition zone, and inner corona with, in particular, 2D very high resolution multi-spectral imaging-spectroscopy, and, direct coronal magnetic field measurement, thus providing a unique set of tools to understand the structure and onset of coronal heating. HiRISE’s objectives are natural complements to the Parker Solar Probe and Solar Orbiter-type missions. We present the science case for HiRISE which will address: i) the fine structure of the chromosphere-corona interface by 2D spectroscopy in FUV at very high resolution; ii) coronal heating roots in the inner corona by ambitious externally-occulted coronagraphy; iii) resolved and global helioseismology thanks to continuity and stability of observing at the L1 Lagrange point; and iv) solar variability and space climate with, in addition, a global comprehensive view of UV variability.

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.

We present 850 $��$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~��$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities., 18 pages, accepted for publication in ApJ

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

Low- and intermediate-mass stars evolve into asymptotic giant branch (AGB) stars near the end of their lives, losing mass through slow and massive winds. The ejected material creates a chemically-rich expanding envelope around the star, namely the circumstellar envelope (CSE). Investigating the anisotropy of the mass-loss phenomenon on the AGB is crucial in gaining a better understanding of the shaping of the CSE during the transition from AGB star to planetary nebula (PN). We investigate possible signs of deviation from spherical symmetry in the CO-emitting CSEs of 70 AGB stars by analysing their emission maps in CO J=2−1 and 3−2 observed with the Atacama Compact Array, as part of the DEATHSTAR project. We find that about one third of the sources are likely aspherical, as they exhibit large-scale asymmetries that are unlikely to have been created by a smooth wind. Further high-resolution observations would be necessary to investigate the nature of, and the physical processes behind, these asymmetrical structures.

High-mass microquasar jets, produced in an accreting compact object in orbit around a massive star, must cross a region filled with stellar wind. The combined effects of the wind and orbital motion can strongly affect the jet properties on binary scales and beyond. The study of such effects can shed light on how high-mass microquasar jets propagate and terminate in the interstellar medium. We study for the first time, using relativistic hydrodynamical simulations, the combined impact of the stellar wind and orbital motion on the properties of high-mass microquasar jets on binary scales and beyond. We have performed 3-dimensional relativistic hydrodynamic simulations, using the PLUTO code, of a microquasar scenario in which a strong weakly relativistic wind from a star interacts with a relativistic jet under the effect of the binary orbital motion. The parameters of the orbit are chosen such that the results can provide insight on the jet-wind interaction in compact systems like for instance Cyg X-1 or Cyg X-3. The wind and jet momentum rates are set to values that may be realistic for these sources and lead to moderate jet bending, which together with the close orbit and jet instabilities could trigger significant jet precession and disruption. For high-mass microquasars with orbit size a ∼ 0.1 AU, and (relativistic) jet power Lj∼1037(M˙w/10−6M⊙yr−1) erg s−1, where M˙w is the stellar wind mass rate, the combined effects of the stellar wind and orbital motion can induce relativistic jet disruption on scales ∼1 AU., VB-R acknowledges the support by the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant PID2019-105510GB-C31 and through the ‘Center of Excellence María de Maeztu 2020–2023’ award to the ICCUB (CEX2019-000918-M), and by the Catalan DEC grant 2017 SGR 643. VB-R is Correspondent Researcher of CONICET, Argentina, at the IAR.

Magnetic interactions between stars and close-in planets may lead to a detectable signal on the stellar disk. HD 189733 is one of the key exosystems thought to harbor magnetic interactions, which may have been detected in August 2013. We present a set of twelve wind models at that period, covering the possible coronal states and coronal topologies of HD 189733 at that time. We assess the power available for the magnetic interaction and predict its temporal modulation. By comparing the predicted signal with the observed signal, we find that some models could be compatible with an interpretation based on star-planet magnetic interactions. We also find that the observed signal can be explained only with a stretch-and-break interaction mechanism, while that the Alfv\'en wings scenario cannot deliver enough power. We finally demonstrate that the past observational cadence of HD 189733 leads to a detection rate of only between 12 to 23%, which could explain why star-planet interactions have been hard to detect in past campaigns. We conclude that the firm confirmation of their detection will require dedicated spectroscopic observations covering densely the orbital and rotation period, combined with scarcer spectropolarimetric observations to assess the concomitant large-scale magnetic topology of the star., Comment: 16 pages, 8 figures, 4 tables, accepted for publication in MNRAS

[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 SOFIA+ALMA continuum and spectral-line polarisation data on the massive molecular cloud BYF 73, revealing important details about the magnetic field morphology, gas structures, and energetics in this unusual massive star formation laboratory. The 154$\mu$m HAWC+ polarisation map finds a highly organised magnetic field in the densest, inner 0.55$\times$0.40 pc portion of the cloud, compared to an unremarkable morphology in the cloud's outer layers. The 3mm continuum ALMA polarisation data reveal several more structures in the inner domain, including a pc-long, $\sim$500 M$_{\odot}$ "Streamer" around the central massive protostellar object MIR 2, with magnetic fields mostly parallel to the east-west Streamer but oriented north-south across MIR 2. The magnetic field orientation changes from mostly parallel to the column density structures to mostly perpendicular, at thresholds $N_{\rm crit}$ = 6.6$\times$10$^{26}$ m$^{-2}$, $n_{\rm crit}$ = 2.5$\times$10$^{11}$ m$^{-3}$, and $B_{\rm crit}$ = 42$\pm$7 nT. ALMA also mapped Goldreich-Kylafis polarisation in $^{12}$CO across the cloud, which traces in both total intensity and polarised flux, a powerful bipolar outflow from MIR 2 that interacts strongly with the Streamer. The magnetic field is also strongly aligned along the outflow direction; energetically, it may dominate the outflow near MIR 2, comprising rare evidence for a magnetocentrifugal origin to such outflows. A portion of the Streamer may be in Keplerian rotation around MIR 2, implying a gravitating mass 1350$\pm$50 M$_{\odot}$ for the protostar+disk+envelope; alternatively, these kinematics can be explained by gas in free fall towards a 950$\pm$35 M$_{\odot}$ object. The high accretion rate onto MIR 2 apparently occurs through the Streamer/disk, and could account for $\sim$33% of MIR 2's total luminosity via gravitational energy release., Comment: 33 pages, 32 figures, accepted by ApJ. Line-Integral Convolution (LIC) images and movie versions of Figures 3b, 7, and 29 are available at https://gemelli.spacescience.org/~pbarnes/research/champ/papers/

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

We investigate the impact of optically thick clumping on stellar wind diagnostics in O supergiants and constrain wind parameters associated with porosity in velocity space. This is the first time the effects of optically thick clumping have been investigated for a sample of massive hot stars, using models including a full optically thick clumping description. We re-analyse spectroscopic observations of a sample of eight O supergiants. Using a genetic algorithm wrapper around the NLTE atmosphere code FASTWIND we obtain simultaneous fits to optical and UV spectra and determine photospheric and wind properties and surface abundances. We provide empirical constraints on a number of wind parameters including the clumping factors, mass-loss rates and terminal wind velocities. Additionally, we establish the first systematic empirical constraints on velocity filling factors and interclump densities. These parameters describe clump distribution in velocity-space and density of the interclump medium in physical-space, respectively. We observe a mass-loss rate reduction of a factor of 3.6 compared to theoretical predictions from Vink et al. (2000), and mass-loss rates within a factor 1.4 of predictions from Bj\"orklund et al. (2021). We confirm that including optically thick clumping allows simultaneous fitting of recombination lines and resonance lines, including the unsaturated UV phosphorus lines (Pv 1118-1128), without reducing the phosphorus abundance. We find that, on average, half of the wind velocity field is covered by dense clumps. We also find that these clumps are 25 times denser than the average wind, and that the interclump medium is 3-10 times less dense than the mean wind. The former result agrees well with theoretical predictions, the latter suggests that lateral filling-in of radially compressed gas might be critical for setting the scale of the rarefied interclump matter., Comment: Accepted for publication in A&A; 35 pages, 20 figures

Context. Debris disks have commonly been studied around intermediate-mass stars. Their intense radiation fields are believed to efficiently remove the small dust grains that are constantly replenished by collisions. For lower-mass stars, in particular M-stars, the dust removal mechanism needs to be further investigated given the much weaker radiation field produced by these objects. Aims. We present new polarimetric observations of the nearly edge-on disk around the pre-main sequence M-type star GSC 07396-00759, taken with VLT/SPHERE IRDIS, with the aim to better understand the morphology of the disk, its dust properties, and the star-disk interaction via the stellar mass-loss rate. Methods. We model our observations to characterize the location and properties of the dust grains using the Henyey-Greenstein approximation of the polarized phase function and evaluate the strength of the stellar winds. Results. We find that the observations are best described by an extended and highly inclined disk ($i\approx 84.3\,^{\circ}\pm0.3$) with a dust distribution centered at a radius $r_{0}\approx107\pm2$ au. The polarized phase function $S_{12}$ is best reproduced by an anisotropic scattering factor $g\approx0.6$ and small micron-sized dust grains with sizes $s>0.3\,\mathrm��$m. We furthermore discuss some of the caveats of the approach and a degeneracy between the grain size and the porosity. Conclusions. Even though the radius of the disk may be over-estimated, our results suggest that using a given scattering theory might not be sufficient to fully explain key aspects such as the shape of the phase function, or the dust grain size. With the caveats in mind, we find that the average mass-loss rate of GSC 07396-00759 can be up to 500 times stronger than that of the Sun, supporting the idea that stellar winds from low-mass stars can evacuate small dust grains from the disk., 22 pages, 17 figures, 4 tables

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]: The Vela X-1 system is one of the best-studied X-ray binaries because it was detected early, has persistent X-ray emission, and a rich phenomenology at many wavelengths. The system is frequently quoted as the archetype of wind-accreting high-mass X-ray binaries, and its parameters are referred to as typical examples. Specific values for these parameters have frequently been used in subsequent studies, however, without full consideration of alternatives in the literature, even more so when results from one field of astronomy (e.g., stellar wind parameters) are used in another (e.g., X-ray astronomy). The issues and considerations discussed here for this specific, very well-known example will apply to various other X-ray binaries and to the study of their physics., [Aims]: We provide a robust compilation and synthesis of the accumulated knowledge about Vela X-1 as a solid baseline for future studies, adding new information where available. Because this overview is targeted at a broader readership, we include more background information on the physics of the system and on methods than is usually done. We also attempt to identify specific avenues of future research that could help to clarify open questions or determine certain parameters better than is currently possible., [Methods]: We explore the vast literature for Vela X-1 and on modeling efforts based on this system or close analogs. We describe the evolution of our knowledge of the system over the decades and provide overview information on the essential parameters. We also add information derived from public data or catalogs to the data taken from the literature, especially data from the Gaia EDR3 release., [Results]: We derive an updated distance to Vela X-1 and update the spectral classification for HD 77518. At least around periastron, the supergiant star may be very close to filling its Roche lobe. Constraints on the clumpiness of the stellar wind from the supergiant star have improved, but discrepancies persist. The orbit is in general very well determined, but a slight difference exists between the latest ephemerides. The orbital inclination remains the least certain factor and contributes significantly to the uncertainty in the neutron star mass. Estimates for the stellar wind terminal velocity and acceleration law have evolved strongly toward lower velocities over the years. Recent results with wind velocities at the orbital distance in the range of or lower than the orbital velocity of the neutron star support the idea of transient wind-captured disks around the neutron star magnetosphere, for which observational and theoretical indications have emerged. Hydrodynamic models and observations are consistent with an accretion wake trailing the neutron star., [Conclusions]: With its extremely rich multiwavelength observational data and wealth of related theoretical studies, Vela X-1 is an excellent laboratory for exploring the physics of accreting X-ray binaries, especially in high-mass systems. Nevertheless, much room remains to improve the accumulated knowledge. On the observational side, well-coordinated multiwavelength observations and observing campaigns addressing the intrinsic variability are required. New opportunities will arise through new instrumentation, from optical and near-infrared interferometry to the upcoming X-ray calorimeters and X-ray polarimeters. Improved models of the stellar wind and flow of matter should account for the non-negligible effect of the orbital eccentricity and the nonspherical shape of HD 77581. There is a need for realistic multidimensional models of radiative transfer in the UV and X-rays in order to better understand the wind acceleration and effect of ionization, but these models remain very challenging. Improved magnetohydrodynamic models covering a wide range of scales are required to improve our understanding of the plasma-magnetosphere coupling, and they are thus a key factor for understanding the variability of the X-ray flux and the torques applied to the neutron star. A full characterization of the X-ray emission from the accretion column remains another so far unsolved challenge., IEM has received funding from the Research Foundation Flanders (FWO), from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement No 665501 and from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 863412). SMN acknowledges funding by the Spanish Ministry MCIU under project RTI2018-096686-B-C21 (MCIU/AEI/FEDER, UE), co-funded by FEDER funds and by the Unidad de Excelencia María de Maeztu, ref. MDM-2017-0765. VG is supported through the Margarete von Wrangell fellowship by the ESF and the Ministry of Science, Research and the Arts Baden-Württemberg. IEM and PK acknowledge support from the ESA/ESAC Faculty Visiting Scientist Programme. JMA acknowledges support from the Spanish Government Ministerio de Ciencia through grant PGC2018-095049-B-C22. FJE acknowledges support from ESCAPE – The European Science Cluster of Astronomy & Particle Physics ESFRI Research Infrastructures, which received funding from the European Union’s Horizon 2020 research and innovation programme under Grant Agreement no. 824064. JvdE is supported by a Lee Hysan Junior Research Fellowship awarded by St. Hilda’s College.

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.

The properties of carbon stars in the Magellanic Clouds (MCs) and their total dust production rates are predicted by fitting their spectral energy distributions (SED) over pre-computed grids of spectra reprocessed by dust. The grids are calculated as a function of the stellar parameters by consistently following the growth for several dust species in their circumstellar envelopes, coupled with a stationary wind. Dust radiative transfer is computed taking as input the results of the dust growth calculations. The optical constants for amorphous carbon are selected in order to reproduce different observations in the infrared and optical bands of \textit{Gaia} Data Release 2. We find a tail of extreme mass-losing carbon stars in the Large Magellanic Cloud (LMC) with low gas-to-dust ratios that is not present in the Small Magellanic Cloud (SMC). Typical gas-to-dust ratios are around $700$ for the extreme stars, but they can be down to $\sim160$--$200$ and $\sim100$ for a few sources in the SMC and in the LMC, respectively. The total dust production rate for the carbon star population is $\sim 1.77\pm 0.45\times10^{-5}$~M$_\odot$~yr$^{-1}$, for the LMC, and $\sim 2.52\pm 0.96 \times 10^{-6}$~M$_\odot$~yr$^{-1}$, for the SMC. The extreme carbon stars observed with the Atacama Large Millimeter Array and their wind speed are studied in detail. For the most dust-obscured star in this sample the estimated mass-loss rate is $\sim 6.3 \times 10^{-5}$~M$_\odot$~yr$^{-1}$. The grids of spectra are available at: https://ambrananni085.wixsite.com/ambrananni/online-data-1 and included in the SED-fitting python package for fitting evolved stars https://github.com/s-goldman/Dusty-Evolved-Star-Kit ., 22 pages, 16 figures, accepted for publication in MNRAS

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).

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.

GRS 1915+105 has been in a bright flux state for more than 2 decades, but in 2018 a significant drop in flux was observed, partly due to changes in the central engine along with increased X-ray absorption. The aim of this work is to explore how X-ray spectro-polarimetry can be used to derive the basic geometrical properties of the absorbing and reflecting matter. In particular, the expected polarisation of the radiation reflected off the disc and the putative outflow is calculated. We use \textit{NuSTAR} data collected after the flux drop to derive the parameters of the system from hard X-ray spectroscopy. The spectroscopic parameters are then used to derive the expected polarimetric signal, using results from a MonteCarlo radiative transfer code both in the case of neutral and fully ionised matter. From the spectral analysis, we find that the continuum emission becomes softer with increasing flux, and that in all flux levels the obscuring matter is highly ionised. This analysis, on the other hand, confirms that spectroscopy alone is unable to put constraints on the geometry of the reflectors. Simulations show that X-ray polarimetric observations, like those that will be provided soon by the Imaging X-ray Polarimetry Explorer (IXPE), will help to determine the geometrical parameters which are left unconstrained by the spectroscopic analysis., 14 pages, 13 figures, Accepted for publication in Astronomy & Astrophysics

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.

Filaments play an important role in star formation, but the formation process of filaments themselves is still unclear. The high-mass star forming clump G286.21+0.17 (G286 for short) that contains an "L" type filament was thought to undergo global collapse. Our high resolution ALMA band 3 observations resolve the gas kinematics of G286 and reveal two sub-clumps with very different velocities inside it. We find that the "blue profile" (an indicator of gas infall) of HCO+ lines in single dish observations of G286 is actually caused by gas emission from the two sub-clumps rather than gas infall. We advise great caution in interpreting gas kinematics (e.g., infall) from line profiles toward distant massive clumps in single dish observations. Energetic outflows are identified in G286 but the outflows are not strong enough to drive expansion of the two sub-clumps. The two parts of the "L" type filament ("NW-SE" and "NE-SW" filaments) show prominent velocity gradients perpendicular to their major axes, indicating that they are likely formed due to large-scale compression flows. We argue that the large-scale compression flows could be induced by the expansion of nearby giant HII regions. The "NW-SE" and "NE-SW" filaments seem to be in collision, and a large amount of gas has been accumulated in the junction region where the most massive core G286c1 forms., Comment: Accepted by MNRAS

We investigate the wind of lambda And, a solar-mass star that has evolved off the main sequence becoming a sub-giant. We present spectropolarimetric observations and use them to reconstruct the surface magnetic field of lambda And. Although much older than our Sun, this star exhibits a stronger (reaching up to 83 G) large-scale magnetic field, which is dominated by the poloidal component. To investigate the wind of lambda And, we use the derived magnetic map to simulate two stellar wind scenarios, namely a polytropic wind (thermally-driven) and an Alfven-wave driven wind with turbulent dissipation. From our 3D magnetohydrodynamics simulations, we calculate the wind thermal emission and compare it to previously published radio observations and more recent VLA observations, which we present here. These observations show a basal sub-mJy quiescent flux level at ~5 GHz and, at epochs, a much larger flux density (>37 mJy), likely due to radio flares. By comparing our model results with the radio observations of lambda And, we can constrain its mass-loss rate Mdot. There are two possible conclusions. 1) Assuming the quiescent radio emission originates from the stellar wind, we conclude that lambda And has Mdot ~ 3e-9 Msun/yr, which agrees with the evolving mass-loss rate trend for evolved solar-mass stars. 2) Alternatively, if the quiescent emission does not originate from the wind, our models can only place an upper limit on mass-loss rates, indicating that Mdot, Comment: Accepted for publication in MNRAS

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.

Nearby pulsar wind nebulae exhibit complex morphological features: jets, torus, arcs and knots. These structures are well captured and understood in the scope of global magnetohydrodynamic models. However, the origin of knots in the inner radius of the Crab Nebula remains elusive. In this work, we investigate the dynamics of the shock front and downstream flow with a special emphasis on the reconnecting equatorial current sheet. We examine whether giant plasmoids produced in the reconnection process could be good candidates for the knots. To this end, we perform large semi-global three-dimensional particle-in-cell simulations in a spherical geometry. The hierarchical merging plasmoid model is used to extrapolate numerical results to pulsar wind nebula scales. The shocked material collapses into the midplane, forming and feeding a large-scale, but thin, ring-like current layer. The sheet breaks up into a dynamical chain of merging plasmoids reminiscent of three-dimensional reconnection. Plasmoids grow to a macroscopic size. The final number of plasmoids predicted is solely governed by the inverse of the dimensionless reconnection rate. The formation of giant plasmoids is a robust feature of pulsar wind termination shocks. They provide a natural explanation for the inner-ring knots in the Crab Nebula, provided that the nebula is highly magnetized., Comment: 8 pages, 7 figures, Accepted for publication in Astronomy & Astrophysics, revision includes A&A language editing

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.

The discovery of GRB 170817A, the first unambiguous off-axis short gamma-ray burst arising from a neutron star merger, has challenged our understanding of the angular structure of relativistic jets. Studies of the jet propagation usually assume that the jet is ejected from the central engine with a top-hat structure and its final structure, which determines the observed light curve and spectra, is primarily regulated by the interaction with the nearby environment. However, jets are expected to be produced with a structure that is more complex than a simple top-hat, as shown by global accretion simulations. We present numerical simulations of short GRBs launched with a wide range of initial structures, durations and luminosities. We follow the jet interaction with the merger remnant wind and compute its final structure at distances $\gtrsim 10^{11}$~cm from the central engine. We show that the final jet structure, as well as the resulting afterglow emission, depend strongly on the initial structure of the jet, its luminosity and duration. While the initial structure at the jet is preserved for long-lasting SGRBs, it is strongly modified for jets barely making their way through the wind. This illustrates the importance of combining the results of global simulations with propagation studies in order to better predict the expected afterglow signatures from neutron star mergers. Structured jets provide a reasonable description of the GRB 170817A afterglow emission with an off-axis angle $\theta_{\rm obs} \approx 22.5^\circ$., Comment: 9 pages, 7 Figures; Added discussion; Accepted for publication in MNRAS

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

Context. NGC 2264 is a young cluster whose accretion properties can be investigated in detail by taking advantage of the FLAMES data in the context of the Gaia-ESO Survey. In fact, the analysis of the Hα emission line profile can provide us with information about the accretion and ejection activity of young stars. However, a strong nebular emission that contributes to the Hα emission can alter the profiles, with consequences for their physical interpretation. Aims. Our study is aimed at investigating the accretion and ejection properties of NGC 2264 by applying a proper treatment of the sky contribution to the Hα and forbidden emission lines (FELs; [SII] and [NII] doublets). Methods. We developed a tool, the OHαNA-method, to handle the strong nebular contribution and spectra with spurious profiles of the Hα and FELs, namely altered Hα profiles or absorption features artificially created where emission lines (FELs) are expected. We derived the quantitative measurements of relevant parameters to describe the accretion and ejection processes in young members of NGC 2264, focusing on reliable quantities derived from the width of the lines, which is relatively unaffected by the nebular emission, unlike the intensity peak, which can be altered significantly. Results. We derive the quantitative measurements related to the Hα emission line and discuss the comparison between the original and sky-subtracted spectra. We thus reveal possible profile alterations with consequences for their physical interpretation. Furthermore, we show the analysis of the variability for multi-epoch observations, also deriving the velocity of the infalling and outflowing plasma from the wings of the broad Hα emission line (in accreting stars). We also explore the mass accretion rate versus full width at zero intensity of the Hα line, namely Ṁ versus FWZI(Hα), a correlation based on the width of the emission line, which is expected to be more robust with respect to any measurement derived from the peak (e.g., Hα10%) and possibly altered by the nebular contribution. Conclusions. We are able to ascertain that more than 20% of the confirmed accretors, which have already been identified in NGC 2264, are affected by the alteration of their line profiles due to the contribution of the nebular emission. Therefore, this is an important issue to consider when investigating accretion and ejection processes in young stellar clusters. While a small fraction of spectra can be unequivocally classified as either unaffected by nebular emission or dominated by nebular emission, the majority (> 90%) represent intermediate cases whose spectral features have to be investigated in detail to derive reliable measurements of the relevant parameters and their physical implications.