Search

At metallicities lower than that of the Small Magellanic Cloud, it remains essentially unexplored how fossil magnetic fields, forming large-scale magnetospheres, could affect the evolution of massive stars, thereby impacting the fundamental building blocks of the early Universe. We extend our stellar evolution model grid with representative calculations of main-sequence, single-star models with initial masses of 20 and 60 M$_\odot$, including appropriate changes for low-metallicity environments ($Z = 10^{-3}-10^{-6}$). We scrutinise the magnetic, rotational, and chemical properties of the models. When lowering the metallicity, the rotational velocities can become higher and the tendency towards quasi-chemically homogeneous evolution increases. While magnetic fields aim to prevent the development of this evolutionary channel, the weakening stellar winds lead to less efficient magnetic braking in our models. Since the stellar radius is almost constant during a blueward evolution caused by efficient chemical mixing, the surface magnetic field strength remains unchanged in some models. We find core masses at the terminal-age main sequence between 22 and 52 M$_\odot$ for initially 60 M$_\odot$ models. This large difference is due to the vastly different chemical and rotational evolution. We conclude that in order to explain chemical species and, in particular, high nitrogen abundances in the early Universe, the adopted stellar models need to be under scrutiny. The assumptions regarding wind physics, chemical mixing, and magnetic fields will strongly impact the model predictions., Comment: Accepted for publication in MNRAS. A full Reproduction Package is available

CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quantitative spectroscopy of hot stars employed within the X-Shooting ULLYSES collaboration, from grid-based approaches to tailored fits. By performing a blind test, we gain an overview about the similarities and differences of the resulting parameters. Our study aims to provide an overview of the parameter spread caused by different approaches. METHODS: For three different stars from the sample (SMC O5 star AzV 377, LMC O7 star Sk -69 50, and LMC O9 star Sk -66 171), we employ different atmosphere codes (CMFGEN, Fastwind, PoWR) and strategies to determine their best-fitting model. For our analyses, UV and optical spectra are used to derive the properties with some methods relying purely on optical data for comparison. To determine the overall spectral energy distribution, we further employ additional photometry from the literature. RESULTS: Effective temperatures for each of three sample stars agree within 3 kK while the differences in log g can be up to 0.2 dex. Luminosity differences of up to 0.1 dex result from different reddening assumptions, which seem to be larger for the methods employing a genetic algorithm. All sample stars are nitrogen-enriched. CONCLUSIONS: We find a reasonable agreement between the different methods. Tailored fitting tends to be able to minimize discrepancies obtained with more course or automatized treatments. UV spectral data is essential for the determination of realistic wind parameters. For one target (Sk -69 50), we find clear indications of an evolved status., Comment: 19+15 pages, 21+4 figures, accepted version (A&A 689, A30) including language editing, condensed abstract

The Hubble Space Telescope has devoted 500 orbits to observing 250 massive stars with low metallicity in the ultraviolet (UV) range within the framework of the ULLYSES program. The X-Shooting ULLYSES (XShootU) project enhances the legacy value of this UV dataset by providing high-quality optical and near-infrared spectra, which are acquired using the wide-wavelength-coverage X-shooter spectrograph at ESO's Very Large Telescope. XShootU emphasises the importance of combining UV with optical spectra for the consistent determination of key stellar parameters such as effective temperature, surface gravity, luminosity, abundances, and wind characteristics including mass-loss rates as a function of metallicity. Since uncertainties in these parameters have implications across various branches of astrophysics, the data and modelling generated by the XShootU project are poised to significantly advance our understanding of massive stars at low metallicity. This is particularly crucial for confidently interpreting JWST data of the earliest stellar generations, making XShootU a unique resource for comprehending individual spectra of low-metallicity stars., Comment: 6 pages, 6 figures. ESO Large Programme Overview

In this paper we present the first set of 3D magnetohydrodynamic (MHD) simulations performed with the riemann geomesh code. We study the dynamics of the magnetically channeled winds of magnetic massive stars in full three dimensions using a code that is uniquely suited to spherical problems. Specifically, we perform isothermal simulations of a smooth wind on a rotating star with a tilted, initially dipolar field. We compare the mass-loss, angular momentum loss, and magnetospheric dynamics of a template star (with the properties that are reminiscent of the O4 supergiant {\zeta} Pup) over a range of rotation rates, magnetic field strengths, and magnetic tilt angles. The simulations are run up to a quasi-steady state and the results are observed to be consistent with the existing literature, showing the episodic centrifugal breakout events of the mass outflow, confined by the magnetic field loops that form the closed magnetosphere of the star. The catalogued results provide perspective on how angular-momentum loss varies for different configurations of rotation rate, magnetic field strength, and large magnetic tilt angles. In agreement with previous 2D MHD studies, we find that high magnetic confinement reduces the overall mass-loss rate, and higher rotation increases the mass-loss rate. This and future studies will be used to estimate the angular-momentum evolution, spin-down time, and mass-loss evolution of magnetic massive stars as a function of magnetic field strength, rotation rate, and dipole tilt.

Early-type B stars with strong magnetic fields and rapid rotation form centrifugal magnetospheres (CMs), as the relatively weak stellar wind becomes magnetically confined and centrifugally supported above the Kepler co-rotation radius. CM plasma is concentrated at and above the Kepler co-rotation radius at the intersection between the rotation and magnetic field axis. Stellar rotation can cause these clouds of material to intersect the viewer's line-of-sight, leading to photometric eclipses. However, for stars with strong ($\sim 10\,{\rm kG}$) magnetic fields and rapid rotation, CMs can become optically thick enough for emission to occur via electron scattering. Using high-precision space photometry from a sample of stars with strong H$\alpha$ emission, we apply simulated light curves from the Rigidly Rotating Magnetosphere model to directly infer magnetic and rotational properties of these stars. By comparing the values inferred from photometric modelling to those independently determined by spectropolarimetry, we find that magnetic obliquity angle $\beta$, viewer inclination $i$ and critical rotation fraction $W$ can be approximately recovered for 3 of the 4 stars studied here. However, there are large discrepancies between the optical depth at the Kepler radius $\tau_{\rm K}$ expected from magnetometry, and the values required to match the observations. We show that $\tau_{\rm K}$ of order unity is needed to reasonably match the light curve morphology of our sample stars., Comment: 15 pages, 8 figures. Accepted for publication in MNRAS

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver., Comment: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

We report detailed X-ray observations of the unique binary system $\epsilon$ Lupi, the only known short-period binary consisting of two magnetic early-type stars. The components have comparably strong, but anti-aligned magnetic fields. The orbital and magnetic properties of the system imply that the magnetospheres overlap at all orbital phases, suggesting the possibility of variable inter-star magnetospheric interaction due to the non-negligible eccentricity of the orbit. To investigate this effect, we observed the X-ray emission from $\epsilon$ Lupi both near and away from periastron passage, using the Neutron Star Interior Composition Explorer mission (NICER) X-ray Telescope. We find that the system produces excess X-ray emission at the periastron phase, suggesting the presence of variable inter-star magnetospheric interaction. We also discover that the enhancement at periastron is confined to a very narrow orbital phase range ($\approx 5\%$ of the orbital period), but the X-ray properties close to periastron phase are similar to those observed away from periastron. From these observations, we infer that the underlying cause is magnetic reconnection heating the stellar wind plasma, rather than shocks produced by wind-wind collision. Finally, by comparing the behavior of $\epsilon$ Lupi with that observed for cooler magnetic binary systems, we propose that elevated X-ray flux at periastron phase is likely a general characteristic of interacting magnetospheres irrespective of the spectral types of the constituent stars., Comment: Accepted for publication in MNRAS (20 pages, 17 figures)

The winds of massive stars have an impact on stellar evolution and on the surrounding medium. The maximum speed reached by these outflows, the terminal wind speed, is a global wind parameter and an essential input for models of stellar atmospheres and feedback. With the arrival of the ULLYSES programme, a legacy UV spectroscopic survey with HST, we have the opportunity to quantify the wind speeds of massive stars at sub-solar metallicity (in the Large and Small Magellanic Clouds, 0.5Z and 0.2Z) at an unprecedented scale. We empirically quantify the wind speeds of a large sample of OB stars, including supergiants, giants, and dwarfs at sub-solar metallicity. Using these measurements, we investigate trends of terminal wind speed with a number of fundamental stellar parameters, namely effective temperature, metallicity, and surface escape velocity. We empirically determined the terminal wind speed for a sample of 149 OB stars in the Magellanic Clouds either by directly measuring the maximum velocity shift of the absorption component of the Civ 1548-1550 line profile, or by fitting synthetic spectra produced using the Sobolev with exact integration method. Stellar parameters were either collected from the literature, obtained using spectral-type calibrations, or predicted from evolutionary models. We find strong trends of terminal wind speed with effective temperature and surface escape speed when the wind is strong enough to cause a saturated P Cygni profile in Civ 1548-1550. We find evidence for a metallicity dependence on the terminal wind speed proportional to Z^0.22+-0.03 when we compared our results to previous Galactic studies. Our results suggest that effective temperature rather than surface escape speed should be used as a straightforward empirical prediction of terminal wind speed and that the observed metallicity dependence is steeper than suggested by earlier works., Comment: 21 pages, 16 figures, 8 tables. Accepted in A&A

We analyse optical and X-ray spectroscopy of the Of?p star HD108, known for its strong dipolar magnetic field and its optical line profile variability with a timescale of $54 \pm 3$ yrs, interpreted as the stellar rotation period. Optical emission lines have now recovered from their minimum emission state reached in 2007 - 2008. The variations of the equivalent width of the H$\alpha$ emission provide constraints on the inclination of the rotation axis ($i$) and the obliquity of the magnetic axis ($\beta$). The best agreement between model and observations is found for ($i$, $\beta$) pairs with $i + \beta \simeq 85^{\circ}$ and $i \in [30^{\circ},55^{\circ}]$. The Balmer emission lines display stochastic variability at the $\sim 5$ % level on timescales of a few days. TESS photometry unveils transient modulations on similar timescales in addition to prominent red noise variations. A Chandra X-ray observation of December 2021, when the star was at a higher emission level, indicates a slight increase of the flux and a spectral hardening compared to the August 2002 XMM-Newton observation, taken near minimum emission state. Magnetohydrodynamic simulations are used to compute synthetic X-ray spectra. With our current best estimate of the $\dot{M}_{B=0}$ mass-loss rate, the simulated X-ray luminosity and spectral energy distribution agree very well with the observations. Finally, the radial velocities vary on a period of 8.5 years with a peak-to-peak amplitude of 10 - 11 km s$^{-1}$, suggesting orbital motion with an unseen companion of at least 4 M$_{\odot}$., Comment: Accepted for publication in MNRAS

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

Magnetism can greatly impact the evolution of stars. In some stars with OBA spectral types there is direct evidence via the Zeeman effect for stable, large-scale magnetospheres, which lead to the spin-down of the stellar surface and reduced mass loss. So far, a comprehensive grid of stellar structure and evolution models accounting for these effects was lacking. For this reason, we computed and studied models with two magnetic braking and two chemical mixing schemes in three metallicity environments with the MESA software instrument. We find notable differences between the subgrids, which affects the model predictions and thus the detailed characterisation of stars. We are able to quantify the impact of magnetic fields in terms of preventing quasi-chemically homogeneous evolution and producing slowly-rotating, nitrogen-enriched ("Group 2") stars. Our model grid is fully open access and open source., Comment: to appear in "Winds of Stars and Exoplanets" Proceedings IAUS 370, 2022, eds.: A.A. Vidotto, L. Fossati, J.S. Vink

Many stars across all classes possess strong enough magnetic fields to influence dynamical flow of material off the stellar surface. For the case of massive stars (O and B types), about 10\% of them harbour strong, globally ordered (mostly dipolar) magnetic fields. The trapping and channeling of their stellar winds in closed magnetic loops leads to {\it magnetically confined wind shocks} (MCWS), with pre-shock flow speeds that are some fraction of the wind terminal speed that can be a few thousand km s$^{-1}$. These shocks generate hot plasma, a source of X-rays. In the last decade, several developments took place, notably the determination of the hot plasma properties for a large sample of objects using \xmm\ and \ch, as well as fully self-consistent MHD modelling and the identification of shock retreat effects in weak winds. In addition, these objects are often sources of H$\alpha$ emission which is controlled by either sufficiently high mass loss rate or centrifugal breakout. Here we review the theoretical aspects of such magnetic massive star wind dynamics., Comment: Accepted for publication invited chapter of the Handbook of X-ray and Gamma-ray Astrophysics published by Nature Springer. arXiv admin note: text overlap with arXiv:1509.06482, arXiv:1605.04979

Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3-60 M$_\odot$), surface equatorial magnetic field strength (0-50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly-rotating, nitrogen-enriched ("Group 2") stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of i) revisiting the derived stellar parameters of known magnetic stars, and ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes., Comment: Accepted for publication in MNRAS. A full reproduction package is available on Zenodo at https://doi.org/10.5281/zenodo.7069766

Polstar is a proposed NASA MIDEX space telescope that will provide high-resolution, simultaneous full-Stokes spectropolarimetry in the far ultraviolet, together with low-resolution linear polarimetry in the near ultraviolet. This observatory offers unprecedented capabilities to obtain unique information on the magnetic and plasma properties of the magnetospheres of hot stars. We describe an observing program making use of the known population of magnetic hot stars to test the fundamental hypothesis that magnetospheres should act to rapidly drain angular momentum, thereby spinning the star down, whilst simultaneously reducing the net mass-loss rate. Both effects are expected to lead to dramatic differences in the evolution of magnetic vs. non-magnetic stars., Comment: 20 pages, 10 figures, accepted for publication in ApSS. arXiv admin note: substantial text overlap with arXiv:2111.06434

Several space missions and instruments for UV spectropolarimetry are in preparation, such as the proposed NASA MIDEX Polstar project, the proposed ESA M mission Arago, and the Pollux instrument on the future LUVOIR-like NASA flagship mission. In the frame of Polstar, we have studied the capabilities these observatories would offer to gain information on the magnetic and plasma properties of the magnetospheres of hot stars, helping us test the fundamental hypothesis that magnetospheres should act to rapidly drain angular momentum, thereby spinning the star down, whilst simultaneously reducing the net mass-loss rate. Both effects are expected to lead to dramatic differences in the evolution of magnetic vs. non-magnetic stars., Comment: Accepted for publication in Astrophysics and Space Science. arXiv admin note: substantial text overlap with arXiv:2111.06434

Magnetic B-stars often exhibit circularly polarized radio emission thought to arise from gyrosynchrotron emission by energetic electrons trapped in the circumstellar magnetosphere. Recent empirical analyses show that the onset and strength of the observed radio emission scale with both the magnetic field strength and the stellar rotation rate. This challenges the existing paradigm that the energetic electrons are accelerated in the current sheet between opposite-polarity field lines in the outer regions of magnetised stellar winds, which includes no role for stellar rotation. Building on recent success in explaining a similar rotation-field dependence of H$\alpha$ line emission in terms of a model in which magnetospheric density is regulated by centrifugal breakout (CBO), we examine here the potential role of the associated CBO-driven magnetic reconnection in accelerating the electrons that emit the observed gyrosynchrotron radio. We show in particular that the theoretical scalings for energy production by CBO reconnection match well the empirical trends for observed radio luminosity, with a suitably small, nearly constant conversion efficiency $\epsilon \approx 10^{-8}$. We summarize the distinct advantages of our CBO scalings over previous associations with an electromotive force, and discuss the potential implications of CBO processes for X-rays and other observed characteristics of rotating magnetic B-stars with centrifugal magnetospheres., Comment: 11 pages, 6 pages, accepted for publication in Monthly Notices of the Royal Astronomical Society

Strongly magnetic, rapidly rotating B-type stars with relatively weak winds form centrifugal magnetospheres (CMs), as the stellar wind becomes magnetically confined above the Kepler co-rotation radius. Approximating the magnetic field as a dipole tilted by an angle $\beta$ with respect to the rotation axis, the CM plasma is concentrated in clouds at and above the Kepler radius along the intersection of the rotational and magnetic equatorial planes. Stellar rotation can bring such clouds in front of the stellar disk, leading to absorption of order 0.1 magnitude ($\sim 10 \%$ of continuum flux). However some stars with prominent CMs, such as $\sigma$ Ori E, show an emission bump in addition to absorption dips, which has been so far unexplained. We show that emission can occur from electron scattering toward the observer when CM clouds are projected off the stellar limb. Using the Rigidly Rotating Magnetosphere model, modified with a centrifugal breakout density scaling, we present a model grid of photometric light curves spanning parameter space in observer inclination angle $i$, magnetic obliquity angle $\beta$, critical rotation fraction $W$, and optical depth at the Kepler radius $\tau_{\text{K}}$. We show that $\tau_{\text{K}}$ of order unity can produce emission bumps of the magnitude $\sim 0.05$ seen in $\sigma$ Ori E. We discuss the implications for modeling the light curves of CM stars, as well as future work for applying the radiative transfer model developed here to 3D MHD simulations of CMs., Comment: Accepted for publication in MNRAS, 11 pages, 9 figures

Numerous magnetic hot stars exhibit gyrosynchrotron radio emission. The source electrons were previously thought to be accelerated to relativistic velocities in the current sheet formed in the middle magnetosphere by the wind opening magnetic field lines. However, a lack of dependence of radio luminosity on the wind power, and a strong dependence on rotation, has recently challenged this paradigm. We have collected all radio measurements of magnetic early-type stars available in the literature. When constraints on the magnetic field and/or the rotational period are not available, we have determined these using previously unpublished spectropolarimetric and photometric data. The result is the largest sample of magnetic stars with radio observations that has yet been analyzed: 131 stars with rotational and magnetic constraints, of which 50 are radio-bright. We confirm an obvious dependence of gyrosynchrotron radiation on rotation, and furthermore find that accounting for rotation neatly separates stars with and without detected radio emission. There is a close correlation between H$\alpha$ emission strength and radio luminosity. These factors suggest that radio emission may be explained by the same mechanism responsible for H$\alpha$ emission from centrifugal magnetospheres, i.e. centrifugal breakout (CBO), however, whereas the H$\alpha$-emitting magnetosphere probes the cool plasma before breakout, radio emission is a consequence of electrons accelerated in centrifugally-driven magnetic reconnection., Comment: 47 pages (16 pages plus appendices), 43 figures, 9 tables, accepted for publication in MNRAS

The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, with consequences for ISM energization and chemical enrichment. Understanding these processes requires accurate observational constraints on the mass-loss rates of the most luminous stars, which can also be used to test theories of stellar wind generation. In the past, mass-loss rates have been characterized via collisional emission processes such as H$\alpha$ and free-free radio emission, but these so-called "density squared" diagnostics require correction in the presence of widespread clumping. Recent observational and theoretical evidence points to the likelihood of a ubiquitously high level of such clumping in hot-star winds, but quantifying its effects requires a deeper understanding of the complex dynamics of radiatively driven winds. Furthermore, large-scale structures arising from surface anisotropies and propagating throughout the wind can further complicate the picture by introducing further density enhancements, affecting mass-loss diagnostics. Time series spectroscopy of UV resonance lines with high resolution and high signal-to-noise are required to better understand this complex dynamics, and help correct "density squared" diagnostics of mass-loss rates. The proposed Polstar mission easily provides the necessary resolution at the sound-speed scale of 20 km s$^{-1}$, on three dozen bright targets with signal-to noise an order of magnitude above that of the celebrated IUE MEGA campaign, via continuous observations that track structures advecting through the wind in real time., Comment: white paper

Polstar is a proposed NASA MIDEX space telescope that will provide high-resolution, simultaneous full-Stokes spectropolarimetry in the far ultraviolet, together with low-resolution linear polarimetry in the near ultraviolet. In this white paper, we describe the unprecedented capabilities this observatory would offer in order to obtain unique information on the magnetic and plasma properties of the magnetospheres of hot stars. This would enable a test of the fundamental hypothesis that magnetospheres should act to rapidly drain angular momentum, thereby spinning the star down, whilst simultaneously reducing the net mass-loss rate. Both effects are expected to lead to dramatic differences in the evolution of magnetic vs. non-magnetic stars., Comment: White paper, 40 pages

Approximately 7% of massive stars host stable surface magnetic fields that are strong enough to alter stellar evolution through their effect on the stellar wind. It is therefore crucial to characterize the strength and structure of these large-scale fields in order to quantify their influence on massive star evolution. This is traditionally done by measuring the circular polarization caused by Zeeman splitting in optical photospheric lines, but we investigate here the possibility of detecting Stokes $V$ signatures in the wind-sensitive resonance lines formed in magnetically confined winds in the high opacity ultraviolet (UV) domain. This unique diagnostic would be accessible to high-sensitivity spaceborne UV spectropolarimeters such as POLSTAR., Comment: Proceedings article from IAUS 360 -- Astronomical Polarimetry 2020: New Era of Multiwavelength Polarimetry. 2 pages, 1 figure

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., Comment: 16 pages, 15 figures, accepted to MNRAS

Rapidly rotating early-type stars with strong magnetic fields frequently show H$\alpha$ emission originating in Centrifugal Magnetospheres (CMs), circumstellar structures in which centrifugal support due to magnetically enforced corotation of the magnetically confined plasma enables it to accumulate to high densities. It is not currently known whether the CM plasma escapes via Centrifugal Breakout (CB), or by an unidentified leakage mechanism. We have conducted the first comprehensive examination of the H$\alpha$ emission properties of all stars currently known to display CM-pattern emission. We find that the onset of emission is dependent primarily on the area of the CM, which can be predicted simply by the value $B_{\rm K}$ of the magnetic field at the Kepler corotation radius $R_{\rm K}$. Emission strength is strongly sensitive to both CM area and $B_{\rm K}$. Emission onset and strength are {\em not} dependent on effective temperature, luminosity, or mass-loss rate. These results all favour a CB scenario, however the lack of intrinsic variability in any CM diagnostics indicates that CB must be an essentially continuous process, i.e.\ it effectively acts as a leakage mechanism. We also show that the emission profile shapes are approximately scale-invariant, i.e.\ they are broadly similar across a wide range of emission strengths and stellar parameters. While the radius of maximum emission correlates closely as expected to $R_{\rm K}$, it is always larger, contradicting models that predict that emission should peak at $R_{\rm K}$., Comment: 33 pages, 26 figures, 2 tables, accepted for publication in MNRAS

Strongly magnetic B-type stars with moderately rapid rotation form `centrifugal magnetospheres' (CMs), from the magnetic trapping of stellar wind material in a region above the Kepler co-rotation radius. A longstanding question is whether the eventual loss of such trapped material occurs from gradual drift and/or diffusive leakage, or through sporadic `{\em centrifugal break out}' (CBO) events, wherein magnetic tension can no longer contain the built-up mass. We argue here that recent empirical results for Balmer-$\alpha$ emission from such B-star CMs strongly favor the CBO mechanism. Most notably, the fact that the onset of such emission depends mainly on the field strength at the Kepler radius, and is largely {\em independent} of the stellar luminosity, strongly disfavors any drift/diffusion process, for which the net mass balance would depend on the luminosity-dependent wind feeding rate. In contrast, we show that in a CBO model the {\em maximum confined mass} in the magnetosphere is independent of this wind feeding rate, and has a dependence on field strength and Kepler radius that naturally explains the empirical scalings for the onset of H$\alpha$ emission, its associated equivalent width, and even its line profile shapes. However, the general lack of observed Balmer emission in late-B and A-type stars could still be attributed to a residual level of diffusive or drift leakage that does not allow their much weaker winds to fill their CMs to the breakout level needed for such emission; alternatively this might result from a transition to a metal-ion wind that lacks the requisite Hydrogen., Comment: 13 pages, 11 figures, accepted for MNRAS

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries., Comment: 11 pages, 10 figures, accepted for publication in MNRAS

In this paper, we present results from the spectropolarimetric follow-up of photometrically selected candidate magnetic B stars from the MOBSTER project. Out of four observed targets, one (HD 38170) is found to host a detectable surface magnetic field, with a maximum longitudinal field measurement of 105$\pm$14 G. This star is chemically peculiar and classified as an $\alpha^2$ CVn variable. Its detection validates the use of TESS to perform a photometric selection of magnetic candidates. Furthermore, upper limits on the strength of a putative dipolar magnetic field are derived for the remaining three stars, and we report the discovery of a previously unknown spectroscopic binary system, HD 25709. Finally, we use our non-detections as case studies to further inform the criteria to be used for the selection of a larger sample of stars to be followed up using high-resolution spectropolarimetry., Comment: 9 pages, 6 figures, 4 tables, accepted by MNRAS (2021 March 23)

The time evolution of angular momentum and surface rotation of massive stars is strongly influenced by fossil magnetic fields via magnetic braking. We present a new module containing a simple, comprehensive implementation of such a field at the surface of a massive star within the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. We test two limiting scenarios for magnetic braking: distributing the angular momentum loss throughout the star in the first case, and restricting the angular momentum loss to a surface reservoir in the second case. We perform a systematic investigation of the rotational evolution using a grid of OB star models with surface magnetic fields ($M_\star=5-60$ M$_\odot$, $\Omega/\Omega_{\rm crit} =0.2-1.0$, $B_{\rm p} =1-20$ kG). We then employ a representative grid of B-type star models ($M_\star=5, 10, 15$ M$_\odot$, $\Omega/\Omega_{\rm crit} =0.2 , 0.5, 0.8$, $B_{\rm p} = 1, 3 ,10, 30$ kG) to compare to the results of a recent self-consistent analysis of the sample of known magnetic B-type stars. We infer that magnetic massive stars arrive at the zero age main sequence with a range of rotation rates, rather than with one common value. In particular, some stars are required to have close-to-critical rotation at the ZAMS. However, magnetic braking yields surface rotation rates converging to a common low value, making it difficult to infer the initial rotation rates of evolved, slowly-rotating stars., Comment: accepted for publication in MNRAS, a zenodo record is available at: https://zenodo.org/record/3250412

Magnetic massive and intermediate-mass stars constitute a separate population whose properties are still not fully understood. Increasing the sample of known objects of this type would help answer fundamental questions regarding the origins and characteristics of their magnetic fields. The MOBSTER Collaboration seeks to identify candidate magnetic A, B and O stars and explore the incidence and origins of photometric rotational modulation using high-precision photometry from the Transiting Exoplanet Survey Satellite (\textit{TESS}) mission. In this contribution, we present an overview of our methods and planned targeted spectropolarimetric follow-up surveys., Comment: 7 pages, 1 figure, Conference proceedings: "Stellar Magnetic Fields - A workshop in honour of the career and contributions of John Landstreet", London, Canada, 2019 (to be published in the Proceedings of the Polish Astronomical Society)

In this contribution, we present the MOBSTER Collaboration, a large community effort to leverage high-precision photometry from the Transiting Exoplanet Survey Satellite (\textit{TESS}) in order to characterize the variability of magnetic massive and intermediate-mass stars. These data can be used to probe the varying column density of magnetospheric plasma along the line of sight for OB stars, thus improving our understanding of the interaction between surface magnetic fields and massive star winds. They can also be used to map out the brightness inhomogeneities present on the surfaces of Ap/Bp stars, informing present models of atomic diffusion in their atmospheres. Finally, we review our current and ongoing studies, which lead to new insights on this topic., Comment: 4 pages, Conference proceedings: "Stars and their variability observed from space", Vienna, Austria, 2019

The investigation of the emergence of life is a major endeavour of science. Astronomy is contributing to it in three fundamental manners: (1) by measuring the chemical enrichment of the Universe, (2) by investigating planet formation and searching for exoplanets with signatures of life and, (3) by determining the abundance of aminoacids and the chemical routes to aminoacid and protein growth in astronomical bodies. This proposal deals with the first two. In the Voyage to 2050, the world-wide scientific community is getting equipped with large facilities for the investigation of the emergence of life in the Universe (i.e. VLT, JWST, ELT, GMT, TMT, ALMA, FAST, VLA, ATHENA, SKA) including the ESA's CHEOPS, PLATO and ARIEL missions. This white paper is a community effort to call for the development of a large ultraviolet optical observatory to gather fundamental data for this investigation that will not be accessible through other ranges of the electromagnetic spectrum. A versatile space observatory with UV sensitivity a factor of 50-100 greater than existing facilities will revolutionize our understanding of the pathway to life in the Universe., Comment: White Paper submitted to the European Space Agency call for Voyage 2050. This WP proposes the development of a European Ultraviolet Visible Observatory (EUVO) to study the origins of life and it is largely based in arXiv 1306.3358 (ApSS,354,229)

Magnetic confinement of stellar winds leads to the formation of magnetospheres, which can be sculpted into Centrifugal Magnetospheres (CMs) by rotational support of the corotating plasma. The conditions required for the CMs of magnetic early B-type stars to yield detectable emission in H$\alpha$ -- the principal diagnostic of these structures -- are poorly constrained. A key reason is that no detailed study of the magnetic and rotational evolution of this population has yet been performed. Using newly determined rotational periods, modern magnetic measurements, and atmospheric parameters determined via spectroscopic modelling, we have derived fundamental parameters, dipolar oblique rotator models, and magnetospheric parameters for 56 early B-type stars. Comparison to magnetic A- and O-type stars shows that the range of surface magnetic field strength is essentially constant with stellar mass, but that the unsigned surface magnetic flux increases with mass. Both the surface magnetic dipole strength and the total magnetic flux decrease with stellar age, with the rate of flux decay apparently increasing with stellar mass. We find tentative evidence that multipolar magnetic fields may decay more rapidly than dipoles. Rotational periods increase with stellar age, as expected for a magnetic braking scenario. Without exception, all stars with H$\alpha$ emission originating in a CM are 1) rapid rotators, 2) strongly magnetic, and 3) young, with the latter property consistent with the observation that magnetic fields and rotation both decrease over time., Comment: 34 pages, 21 figures, 5 tables, accepted for publication in MNRAS

We have performed 3D isothermal MHD simulation of a magnetic rotating massive star with a non-zero dipole obliquity and predicted the radio/sub-mm observable lightcurves and continuum spectra for a frequency range compatible with ALMA. From these results we also compare the model input mass-loss to that calculated from the synthetic thermal emission. Spherical and cylindrical symmetry is broken due to the obliquity of the stellar magnetic dipole resulting in an inclination and phase dependence of both the spectral flux and inferred mass-loss rate, providing testable predictions of variability for oblique rotator. Both quantities vary by factors between 2 and 3 over a full rotational period of the star, demonstrating that the role of rotation as critical in understanding the emission. This illustrates the divergence from a symmetric wind, resulting in a two armed spiral structure indicative of a oblique magnetic rotator. We show that a constant spectral index, $\alpha$, model agrees well with our numerical prediction for a spherical wind for $\nu~<~10^{3} \ \mathrm{GHz}$, however it is unable to capture the behavior of emission at $\nu~>~10^{3} \ \mathrm{GHz}$. As such we caution the use of such constant $\alpha$ models for predicting emission from non-spherical winds such as those which form around magnetic massive stars., Comment: 19 pages, 14 figures

In this paper we introduce the MOBSTER collaboration and lay out its scientific goals. We present first results based on the analysis of nineteen previously known magnetic O, B and A stars observed in 2-minute cadence in sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS) mission. We derive precise rotational periods from the newly obtained light curves and compare them to previously published values. We also discuss the overall photometric phenomenology of the known magnetic massive and intermediate-mass stars and propose an observational strategy to augment this population by taking advantage of the high-quality observations produced by TESS., Comment: 15 pages, 3 figures, accepted for publication by MNRAS (2019 April 22)

Surface magnetic fields have been detected in 5 to 10% of isolated massive stars, hosting outer radiative envelopes. They are often thought to have a fossil origin, resulting from the stellar formation phase. Yet, magnetic massive stars are scarcer in (close) short-period binaries, as reported by the BinaMIcS (Binarity and Magnetic Interaction in various classes of Stars) collaboration. Different physical conditions in the molecular clouds giving birth to isolated stars and binaries are commonly invoked. In addition, we propose that the observed lower magnetic incidence in close binaries may be due to nonlinear tides. Indeed, close binaries are probably prone to tidal instability, a fluid instability growing upon the equilibrium tidal flow via nonlinear effects. Yet, stratified effects have hitherto been largely overlooked. We theoretically and numerically investigate tidal instability in rapidly rotating, stably stratified fluids permeated by magnetic fields. We use the short-wavelength stability method to propose a comprehensive (local) theory of tidal instability at the linear onset, discussing damping effects. Then, we propose a mixing-length theory for the mixing generated by tidal instability in the nonlinear regime. We successfully assess our theoretical predictions against proof-of-concept, direct numerical simulations. Finally, we compare our predictions with the observations of short-period, double-lined spectroscopic binary systems. Using new analytical results, cross-validated by a direct integration of the stability equations, we show that tidal instability can be generated by nonlinear couplings of inertia-gravity waves with the equilibrium tidal flow in short-period massive binaries, even against the Joule diffusion. In the nonlinear regime, a fossil magnetic field can be dissipated by the turbulent magnetic diffusion induced by the saturated tidal flows. We predict that the turbulent Joule diffusion of fossil fields would occur in a few million years for several short-period massive binaries. Therefore, turbulent tidal flows could explain the observed dearth of some short-period magnetic binaries.

Classification is an important supervised machine learning method, which is necessary and challenging issue for ecological research. It offers a way to classify a dataset into subsets that share common patterns. Notably, there are many classification algorithms to choose from, each making certain assumptions about the data and about how classification should be formed. In this paper, we applied eight machine learning classification algorithms such as Decision Trees, Random Forest, Artificial Neural Network, Support Vector Machine, Linear Discriminant Analysis, k-nearest neighbors, Logistic Regression and Naive Bayes on ecological data. The goal of this study is to compare different machine learning classification algorithms in ecological dataset. In this analysis we have checked the accuracy test among the algorithms. In our study we conclude that Linear Discriminant Analysis and k-nearest neighbors are the best methods among all other methods, Comment: 18 pages, 20 figures

In this paper, we present high-resolution HST/COS observations of the extreme magnetic O star NGC 1624-2. These represent the first ultraviolet spectra of this archetypal object. We examine the variability of its wind-sensitive resonance lines, comparing it to that of other known magnetic O stars. In particular, the observed variations in the profiles of the CIV and SiIV doublets between low state and high state are the largest observed in any magnetic O-type star, consistent with the expected properties of NGC 1624-2's magnetosphere. We also observe a redshifted absorption component in the low state, a feature not seen in most stars. We present preliminary modelling efforts based on the Analytic Dynamical Magnetosphere (ADM) formalism, demonstrating the necessity of using non-spherically symmetric models to determine wind/magnetospheric properties of magnetic O stars., Comment: 13 pages, 6 figures, accepted for publication by MNRAS (23 November 2018)

Spectropolarimetric surveys reveal that 8-10\% of OBA stars harbor large-scale magnetic fields, but thus far no such fields have been detected in any classical Be stars. Motivated by this, we present here MHD simulations for how a pre-existing Keplerian disc -- like that inferred to form from decretion of material from rapidly rotating Be stars -- can be disrupted by a rotation-aligned stellar dipole field. For characteristic stellar and disc parameters of a near-critically rotating B2e star, we find that a polar surface field strength of just 10 G can significantly disrupt the disc, while a field of 100 G, near the observational upper limit inferred for most Be stars, completely destroys the disc over just a few days. Our parameter study shows that the efficacy of this magnetic disruption of a disc scales with the characteristic plasma beta {(defined as the ratio between thermal and magnetic pressure)} in the disc, but is surprisingly insensitive to other variations, e.g. in stellar rotation speed, or the mass loss rate of the star's radiatively driven wind. The disc disruption seen here for even a modest field strength suggests that the presumed formation of such Be discs by decretion of material from the star would likely be strongly inhibited by such fields; this provides an attractive explanation for why no large-scale fields are detected from such Be stars., Comment: MNRAS, accepted for publication; 7 pp

Most radiatively-driven massive star winds can be modelled with m-CAK theory resulting in so called fast solution. However, those most rapidly rotating among them, especially when the stellar rotational speed is higher than $\sim 75\%$ of the critical rotational speed, can adopt a different solution called $\Omega$-slow solution characterized by a dense and slow wind. Here, in this work we study the transition region of the solutions where the fast solution changes to the $\Omega$-slow. Using both time-steady and time-dependent numerical codes, we study this transition region for different equatorial models of B-type stars. In all the cases, at certain range of rotational speeds, we found a region where the fast and $\Omega$-slow solution can co-exist. We find that the type of solution obtained in this co-existence region depends heavily on the initial conditions of our models. We also test the stability of the solutions within the co-existence region by performing base density perturbations in the wind. We find that under certain conditions, the fast solution can switch to a $\Omega$-slow solution, or vice versa. Such solution switching may be a possible contributor of material injected into the circumstellar environment of Be stars, without requiring rotational speeds near critical values., Comment: 12 pages, 14 figures, accepted for publication in MNRAS

Although initially thought to be promising for production of the r-process nuclei, standard models of neutrino-heated winds from proto-neutron stars (PNSs) do not reach the requisite neutron-to-seed ratio for production of the lanthanides and actinides. However, the abundance distribution created by the r-, rp-, or $\nu p$-processes in PNS winds depends sensitively on the entropy and dynamical expansion timescale of the flow, which may be strongly affected by high magnetic fields. Here, we present results from magnetohydrodynamic simulations of non-rotating neutrino-heated PNS winds with strong dipole magnetic fields from $10^{14}-10^{16}$ G, and assess their role in altering the conditions for nucleosynthesis. The strong field forms a closed zone and helmet streamer configuration at the equator, with episodic dynamical mass ejections in toroidal plasmoids. We find dramatically enhanced entropy in these regions and conditions favorable for third-peak r-process nucleosynthesis if the wind is neutron-rich. If instead the wind is proton-rich, the conditions will affect the abundances from the $\nu p$-process. We quantify the distribution of ejected matter in entropy and dynamical expansion timescale, and the critical magnetic field strength required to affect the entropy. For $B\gtrsim10^{15}$ G, we find that $\gtrsim10^{-6}$ M$_{\odot}$ and up to $\sim10^{-5}$ M$_{\odot}$ of high entropy material is ejected per highly-magnetized neutron star birth in the wind phase, providing a mechanism for prompt heavy element enrichment of the universe. Former binary companions identified within (magnetar-hosting) supernova remnants, the remnants themselves, and runaway stars may exhibit overabundances. We provide a comparison with a semi-analytic model of plasmoid eruption and discuss implications and extensions., Comment: accepted to MNRAS

The groundbreaking detection of gravitational waves produced by the inspiralling and coalescence of the black hole (BH) binary GW150914 confirms the existence of "heavy" stellar-mass BHs with masses >25 Msun. Initial modelling of the system by Abbott et al. (2016a) supposes that the formation of black holes with such large masses from the evolution of single massive stars is only feasible if the wind mass-loss rates of the progenitors were greatly reduced relative to the mass-loss rates of massive stars in the Galaxy, concluding that heavy BHs must form in low-metallicity (Z < 0.25-0.5 Zsun) environments. However, strong surface magnetic fields also provide a powerful mechanism for modifying mass loss and rotation of massive stars, independent of environmental metallicity (ud-Doula & Owocki 2002; ud-Doula et al. 2008). In this paper we explore the hypothesis that some heavy BHs, with masses >25 Msun such as those inferred to compose GW150914, could be the natural end-point of evolution of magnetic massive stars in a solar-metallicity environment. Using the MESA code, we developed a new grid of single, non-rotating, solar metallicity evolutionary models for initial ZAMS masses from 40-80 Msun that include, for the first time, the quenching of the mass loss due to a realistic dipolar surface magnetic field. The new models predict TAMS masses that are significantly greater than those from equivalent non-magnetic models, reducing the total mass lost by a strongly magnetized 80 Msun star during its main sequence evolution by 20 Msun. This corresponds approximately to the mass loss reduction expected from an environment with metallicity Z = 1/30 Zsun., Comment: 10 pages, 7 figures, 1 table. Accepted for publication in MNRAS

We conducted a survey of seven magnetic O and eleven B-type stars with masses above $8M_{\odot}$ using the Very Large Array in the 1cm, 3cm and 13cm bands. The survey resulted in a detection of two O and two B-type stars. While the detected O-type stars - HD 37742 and HD 47129 - are in binary systems, the detected B-type stars, HD 156424 and ALS 9522, are not known to be in binaries. All four stars were detected at 3cm, whereas three were detected at 1cm and only one star was detected at 13cm. The detected B-type stars are significantly more radio luminous than the non-detected ones, which is not the case for O-type stars. The non-detections at 13cm are interpreted as due to thermal free-free absorption. Mass-loss rates were estimated using 3cm flux densities and were compared with theoretical mass-loss rates, which assume free-free emission. For HD 37742, the two values of the mass-loss rates were in good agreement, possibly suggesting that the radio emission for this star is mainly thermal. For the other three stars, the estimated mass-loss rates from radio observations were much higher than those expected from theory, suggesting either a possible contribution from non- thermal emission from the magnetic star or thermal or non-thermal emission due to interacting winds of the binary system, especially for HD 47129. All the detected stars are predicted to host centrifugal magnetospheres except HD 37742, which is likely to host a dynamical magnetosphere. This suggests that non-thermal radio emission is favoured in stars with centrifugal magnetospheres., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society

High-resolution spectra of the magnetic star HD191612 were acquired using the Chandra X-ray observatory at both maximum and minimum emission phases. We confirm the flux and hardness variations previously reported with XMM-Newton, demonstrating the great repeatability of the behavior of HD191612 over a decade. The line profiles appear typical for magnetic massive stars: no significant line shift, relatively narrow lines for high-Z elements, and formation radius at about 2R*. Line ratios confirm the softening of the X-ray spectrum at the minimum emission phase. Shift or width variations appear of limited amplitude at most (slightly lower velocity and slightly increased broadening at minimum emission phase, but within 1--2 sigma of values at maximum). In addition, a fully self-consistent 3D magnetohydrodynamic (MHD) simulation of the confined wind in HD191612 was performed. The simulation results were directly fitted to the data leading to a remarkable agreement overall between them., Comment: accepted for publication in ApJ