Observational and theoretical work has now established that the fossil fields of magnetic massive stars are surviving remnants from an earlier event, or an earlier evolutionary phase. However, many important questions remain regarding the effects of these fields on the late stages of stellar evolution, as well as their impact on the core-collapse mechanism and the formation of exotic compact objects such as magnetars and gravitational wave progenitors. There is currently a critical need to incorporate the impact of fossil fields in models of the structure and evolution of magnetic stars, and to determine the evolutionary history of magnetic massive stars. We present a preliminary population study of a cluster of co-evolving stars based on MESA evolutionary tracks that account for the effect of magnetic mass-loss quenching.

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.

It was recently claimed that the magnetic O-type star HD 54879 exhibits important radial velocity variability indicative of its presence in a spectroscopic binary. More remarkably, it was furthermore reported that the star underwent a short, sudden variation in spectral type and magnetic field. In this Letter we examine new Narval and ESPaDOnS data of this star in addition to the previously-published FORS2 data and conclude that both the reported velocity variations and the sudden spectral and magnetic changes are spurious.

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.

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.

Surface magnetic fields have a strong impact on stellar mass loss and rotation and, as a consequence, on the evolution of massive stars. In this work we study the influence of an evolving dipolar surface fossil magnetic field with an initial field strength of 4 kG on the characteristics of 15 M$_{\odot}$ solar metallicity models using the Geneva stellar evolution code. Non-rotating and rotating models considering two different scenarios for internal angular momentum transport are computed, including magnetic field evolution, mass-loss quenching, and magnetic braking. Magnetic field evolution results in weakening the initially strong magnetic field, however, in our models an observable magnetic field is still maintained as the star evolves towards the red supergiant phase. At the given initial mass of the models, mass-loss quenching is modest. Magnetic braking greatly enhances chemical element mixing if radial differential rotation is allowed for, on the other hand, the inclusion of surface magnetic fields yields a lower surface enrichment in the case of near solid-body rotation. Models including surface magnetic fields show notably different trends on the Hunter diagram (plotting nitrogen abundance vs $v \sin i$) compared to those that do not. The magnetic models agree qualitatively with the anomalous `Group 2 stars', showing slow surface rotation and high surface nitrogen enhancement on the main sequence.

We present the initial results of the Large Impact of magnetic Fields on the Evolution of hot stars (LIFE) project. The focus of this project is the search for magnetic fields in evolved OBA giants and supergiants with visual magnitudes between 4 and 8, with the aim to investigate how the magnetic fields observed in upper main sequence (MS) stars evolve from the MS until the late post-MS stages. In this paper, we present spectropolarimetric observations of 15 stars observed using the ESPaDOnS instrument of the CFHT. For each star, we have determined the fundamental parameters and have used stellar evolution models to calculate their mass, age and radius. Using the LSD technique, we have produced averaged line profiles for each star. From these profiles, we have measured the longitudinal magnetic field strength and have calculated the detection probability. We report the detection of magnetic fields in two stars of our sample: a weak field of $B_l=1.0\pm0.2\,$G is detected in the post-MS A5 star 19$\,$Aur and a stronger field of $B_l=-230\,\pm\,10\,$G is detected in the MS/post-MS B8/9 star HR$\,$3042.

Magnetic confinement of the winds of hot, massive stars has far-reaching consequences on timescales ranging from hours to Myr. Understanding the long-term effects of this interplay has already led to the identification of two new evolutionary pathways to form `heavy' stellar mass black holes and pair-instability supernova even at galactic metallicity. We are performing 1D stellar evolution model calculations that, for the first time, account for the surface effects and the time evolution of fossil magnetic fields. These models will be thoroughly confronted with observations and will potentially lead to a significant revision of the derived parameters of observed magnetic massive stars.

Over the last decade, tremendous strides have been achieved in our understanding of magnetism in main sequence hot stars. In particular, the statistical occurrence of their surface magnetism has been established (~10%) and the field origin is now understood to be fossil. However, fundamental questions remain: how do these fossil fields evolve during the post-main sequence phases, and how do they influence the evolution of hot stars from the main sequence to their ultimate demise? Filling the void of known magnetic evolved hot (OBA) stars, studying the evolution of their fossil magnetic fields along stellar evolution, and understanding the impact of these fields on the angular momentum, rotation, mass loss, and evolution of the star itself, is crucial to answering these questions, with far reaching consequences, in particular for the properties of the precursors of supernovae explosions and stellar remnants. In the framework of the BRITE spectropolarimetric survey and LIFE project, we have discovered the first few magnetic hot supergiants. Their longitudinal surface magnetic field is very weak but their configuration resembles those of main sequence hot stars. We present these first observational results and propose to interpret them at first order in the context of magnetic flux conservation as the radius of the star expands with evolution. We then also consider the possible impact of stellar structure changes along evolution.

Large-scale dipolar surface magnetic fields have been detected in a fraction of OB stars, however only few stellar evolution models of massive stars have considered the impact of these fossil fields. We are performing 1D hydrodynamical model calculations taking into account evolutionary consequences of the magnetospheric-wind interactions in a simplified parametric way. Two effects are considered: i) the global mass-loss rates are reduced due to mass-loss quenching, and ii) the surface angular momentum loss is enhanced due to magnetic braking. As a result of the magnetic mass-loss quenching, the mass of magnetic massive stars remains close to their initial masses. Thus magnetic massive stars - even at Galactic metallicity - have the potential to be progenitors of `heavy' stellar mass black holes. Similarly, at Galactic metallicity, the formation of pair instability supernovae is plausible with a magnetic progenitor.

Near-solar metallicity (and low-redshift) Pair-Instability Supernova (PISN) candidates challenge stellar evolution models. Indeed, at such a metallicity, even an initially very massive star generally loses so much mass by stellar winds that it will avoid the electron-positron pair-creation instability. We use recent results showing that a magnetic field at the surface of a massive star can significantly reduce its effective mass-loss rate to compute magnetic models of very massive stars (VMSs) at solar metallicity and explore the possibility that such stars end as PISNe. We implement the quenching of the mass loss produced by a surface dipolar magnetic field into the Geneva stellar evolution code and compute new stellar models with an initial mass of $200\,M_\odot$ at solar metallicity, with and without rotation. It considerably reduces the total amount of mass lost by the star during its life. For the non-rotating model, the total (CO-core) mass of the models is $72.8\,M_\odot$ ($70.1\,M_\odot$) at the onset of the electron-positron pair-creation instability. For the rotating model, we obtain $65.6\,M_\odot$ ($62.4\,M\odot$). In both cases, a significant fraction of the internal mass lies in the region where pair instability occurs in the $\log(T)-\log(\rho)$ plane. The interaction of the reduced mass loss with the magnetic field efficiently brakes the surface of the rotating model, producing a strong shear and hence a very efficient mixing that makes the star evolve nearly homogeneously. The core characteristics of our models indicate that solar metallicity models of magnetic VMSs may evolve to PISNe (and pulsation PISNe).

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.

Stellar evolution models of massive stars are very sensitive to the adopted mass-loss scheme. The magnitude and evolution of mass-loss rates significantly affect the main sequence evolution, and the properties of post-main sequence objects, including their rotational velocities. Driven by potential discrepancies between theoretically predicted and observationally derived mass-loss rates in the OB star range, we particularly aim to investigate the response to mass-loss rates that are lower than currently adopted, in parallel with the mass-loss behavior at the "first" bi-stability jump. We perform 1D hydrodynamical model calculations of single $20 - 60 \, M_{\odot}$ Galactic ($Z = 0.014$) stars where the effects of stellar winds are already significant during the main sequence phase. We develop an experimental wind routine to examine the behavior and response of the models under the influence of different mass-loss rates. This observationally guided, simple and flexible wind routine is not a new mass-loss description but a useful tool based on the Wind-momentum Luminosity Relation and other scaling relations, and provides a meaningful base for various tests and comparisons. The main result of this study indicates a dichotomy when accounting for currently debated problems regarding mass-loss rates of hot massive stars. In a fully diffusive approach, and for commonly adopted initial rotational velocities, lower mass-loss rates than theoretically predicted require to invoke an additional source of angular momentum loss (either due to bi-stability braking, or yet unidentified) to brake down surface rotational velocities. [...]