UV wind line variability in OB stars appears to be universal. We review the evidence that the variability is due to large, dense, optically thick structures rooted in or near the photosphere. Using repeated observations and a simple model we translate observed profile variations into optical depth variations and, consequently, variations in measured mass loss rates. Although global rates may be stable, measured rates vary. Consequently, profile variations infer how mass loss rates determined from UV wind lines vary. These variations quantify the intrinsic error inherent in any mass loss rate derived from a single observation. These derived rates can differ by factors of 3 or more. Our results also imply that rates from non-simultaneous observations (such as UV and ground based data) need not agree. Finally, we use our results to examine the nature of the structures responsible for the variability.

The first magnetic field in a star other than the Sun was detected in 1947 in the star 78 Vir. Today, we know that about 10% of these intermediate-mass and high-mass stars have strong, large-scale surface magnetic fields whose origin has remained a mystery till today. It has been suggested that merging of main-sequence and pre-main-sequence stars could produce such strong fields. The massive star τ Sco is a well-known member of the group of magnetic stars and is a blue straggler given its apparently young age compared to that of other members of the Upper Scorpius association. Here, we present 3D magnetohydrodynamic simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain τ Sco’s magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. Such stars are promising progenitors of magnetars, which may give rise to some of the enigmatic fast radio bursts, and their supernova explosions may be affected by the strong magnetic fields.

Spectropolarimetic campaigns have established that large-scale magnetic fields are present at the surfaces of approximately 10% of massive dwarf stars. However, there is a dearth of magnetic field measurements for their deep interiors. Asteroseismology of gravity-mode pulsations combined with rotating magneto-hydrodynamical calculations of the early-B main-sequence star HD 43317 constrain its magnetic field strength to be approximately 5 × 105 G just outside its convective core. This proof-of-concept study for magneto-asteroseismology opens a new window into the observational characterisation of magnetic fields inside massive stars.

While we have growing numbers of massive star observations, it remains unclear how efficient the key physical processes such as mass loss, convection and rotation are, or indeed how they impact each other. We reconcile this with detailed stellar evolution models, yet these models have their own drawbacks with necessary assumptions for 3-dimensional processes like rotation which need to be adapted into 1-dimensional models. The implementation of empirical mass-loss prescriptions in stellar evolution codes can lead to the extrapolation of base rates to unconstrained evolutionary stages leading to a range of uncertain fates. In short, there remain many free parameters and physical processes which need to be calibrated in order to align our theory better with upcoming observations. We have tested various processes such as mass loss and internal mixing, including rotational mixing and convective overshooting, against multiple observational constraints such as using eclipsing binaries, the Humphreys-Davidson limit, and the final masses of Wolf-Rayet stars, across a range of metallicities. In fact, we developed a method of disentangling the effects of mixing and mass loss in the ‘Mass-Luminosity Plane’ allowing direct calibration of these processes. In all cases, it is important to note that a combined appreciation for both stellar winds and internal mixing are important to reproduce observations.

When a star is rapidly rotating, it deviates from spherical symmetry causing non-uniform distributions of the surface gravity and temperature across the surface. These three-dimensional effects lead to an inclination dependence of many observable spectroscopic parameters, however this is often neglected when analyzing rapidly rotating systems. Using spamms, we generate synthetic spectra that account for the 3D geometry of the system and fit them with 1D models to investigate how much the 3D effects can change the derived stellar parameters. We show that these 3D effects can lead to observed temperature differences of thousands of kelvin for the same star viewed at different inclinations, and a systematic underestimation of the helium abundance.

We report recent ESPaDOnS and HARPSpol spectropolarimetric observations from our ongoing magnetic survey of the brightest twenty-five classical Cepheids. Stokes V magnetic signatures are detected in eight of fifteen targets observed to date. The Stokes V profiles show a diversity of morphologies with weak associated longitudinal field measurements of order 1 G. Many of the Stokes V profiles are difficult to interpret in the context of the normal Zeeman effect. They consist of approximately unipolar single or double lobe(s) of positive or negative circular polarization. We hypothesize that these unusual signatures are due to the Zeeman effect modified by atmospheric velocity or magnetic field gradients. In contrast, the Stokes V profiles of Polaris and MY Pup appear qualitatively similar to the complex magnetic signatures of non-pulsating cool supergiants, possibly due to the low pulsation amplitudes of these two stars.

Rotation is one of the important parameters affecting the evolution and final fate of massive stars but the origin of fast rotators remains unclear (imprint of the star formation process, result of binary interactions). In this work, we aim at investigating the binary status, photometric variability, and runaway status of a statistically meaningful sample of Galactic fast-rotating O stars. We perform a comprehensive multi-epoch analysis of new high-quality spectroscopic observations gathered by the IACOB and OWN surveys. Notably, we find that the total percentage of spectroscopic binaries in the investigated sample range between 25 and 40%, in agreement with previous finding for the case of O-type stars with lower projected rotational velocities. On the contrary, the fraction of runaway stars among fast rotators (∼35–50%) is significantly higher than in the case of slow rotators (∼20–30%). By combining all these observational results we will evaluate each scenario about the origin of fast rotators.

The current angular momentum (AM) transport models fail to reproduce asterosieimic observations. One of the best candidates to explain this discrepancy is the magnetic field in radiative zones with its various possible topologies, for instance axisymmetric toroidal magnetic field. If such azimuthal field is strong enough, the Tayler’s instability could occur which induces a magnetic torque that allows a very efficient transport of AM and could trigger dynamo action in radiative layers. If such important field does not emerge at the surface, spectropolarimetry is blind. In this case, the only way to detect and characterise the field is by using magneto-asteroseismology. It consists in searching for the characteristic signatures of magnetic field in the observed frequency spectra of stellar oscillations.

The young open cluster NGC 6231 hosts a rich population of O-type binary stars. We study several of these eccentric short-period massive eclipsing binaries and assess their fundamental parameters. The properties of these systems make them interesting targets to study tidally induced apsidal motion. The analysis of apsidal motion offers a powerful means to obtain information about the internal structure of the stars. Indeed, since the rate of apsidal motion in a binary system is proportional to the internal structure constants of the stars composing it, its value gives direct insight into the internal structure and evolutionary state of these stars. Stellar evolution models are constructed based on the observationally-determined fundamental parameters and a theoretical rate of apsidal motion is inferred. The results are striking: Adopting standard stellar evolution models yields a theoretical rate of apsidal motion much larger than the observational value. This discrepancy results from the standard models predicting too low an efficiency of internal mixing and thus too homogenous stars in terms of density. By enforcing the theoretical rates of apsidal motion to match the observational values, enhanced mixing is required, through a large overshooting parameter and/or additional turbulent/rotational mixing. Our analysis leads to the conclusion that the chemically mixed cores in those massive stars must be more extended than anticipated from standard models.

We study the apsidal motion in close eccentric massive binaries. Measuring the rate of apsidal motion in such a system gives insight into the internal structure and evolutionary state of the stars. We focus on CPD-41° 7742, for which independent studies in the past showed large discrepancies in the longitude of periastron of the orbit, hinting at the presence of apsidal motion. We perform a consistent analysis of all observational data to solve this apparent discrepancy and report the first determination of apsidal motion in this system. This study confirms the need for enhanced mixing in the stellar evolution models of the primary star to reproduce the observational properties. This points towards larger convective cores than usually considered.

During the evolution of massive stars, their properties change significantly. But stellar parameters of massive stars have the biggest uncertainties in stellar astrophysics, specifically in the post-main sequence stages where blue supergiant stars are located. These stars experience mass loss events during their evolution which are supposed to be related to strange-modes instabilities. In this work, we explore the stability of oscillation modes in massive stars for different masses, and a range of mass-loss rates, with the aim to provide clues about the connection between strange modes instabilities and mass-loss events.

The galactic binary star LB-1 contains a recently stripped B-type star. Comparing its properties to detailed binary star models shows that tidal braking and magnetic torques lead to low surface rotational velocities in the stripped donors after Roche-lobe overffiow. Models without magnetic torques cannot reproduce the observed low surface rotation.

The VLT/FLAMES Tarantula Survey (Evans et al. 2011) identified a group of slowly-rotating nitrogen-rich O-type stars that cannot be explained by current evolutionary models. Here we present high-quality VLT/UVES observations of four of these stars that allow a detailed quantitative spectroscopic analysis. We present the analysis of the spectra with a genetic algorithm, and discuss the future steps to be taken to further investigate the cause of the nitrogen enrichment.

The first stars in the Universe have inherited their composition from primordial nucleosynthesis, so they have no metal. These stars, which are also named population III (pop III) stars, began the process of reionization in the Universe and contributed to the metal enrichment with heavy elements. Previous studies showed that they should have been rotating fast due to small or no angular momentum loss, reaching easily the critical velocity since they are massive and have very low stellar winds, thus their mass loss is very low or zero. Our aim is to study how the production of primary nitrogen is affected due to high rotation in the pop III stars. So, we compared grids of pop III stars with zero, average, and high rotation. All these models have been computed using Geneva code (GENEC) in the mass range of 9M⊙ ≤ Mini ≤ 120M⊙. Due to the rotational mixing, the carbon produced in the He-burning core is diffused towards the H-burning shell, triggering the CNO cycle and producing primary nitrogen. In some models the transition of the shell from a pp-chain H-burning to a CNO H-burning induces a strong energy release and a complete change of the stellar structure and the nucleosynthesis. The production of nitrogen is boosted for the high rotation models.

We present results from 3D MHD simulations of the magnetospheres from massive stars with a dipole magnetic axis that has an arbitrary obliquity angle (β) to the stars rotation axis. As an initial direct application, we examine the global structure of co-rotating disks for tilt angles β=0, 45 and 90 degrees using ζ Pup stellar parameters as a prototype. We find that for models with rapid stellar rotation (∼ 0.5 critical rotation), accumulation surfaces closely resemble the form predicted by the analytic Rigidly Rotating Magnetosphere (RRM) model, but with a mass distribution and outer disk termination set by centrifugal breakout processes. However, some significant differences are found including warping resulting from the dynamic nature of the MHD models in contrast to static RRM models. These MHD models can be used to synthesize rotational modulation of photometric absorption and H-alpha emission for a direct comparison with observations.

Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries.

Luminous Blue Variable stars (LBVs) are rare and enigmatic. Often cited as evolutionary stages in the single-star evolution, the idea that binary evolution produces the LBV state was already considered, 30 years ago.

It is now commonly accepted that a significant part of massive stars are born in multiple systems. One aspect that also emerged is that massive stars have on average at least two companions, i.e. they are triples. This immediately implies that a number of LBVs should have evolved as part of multiple systems.

While some LBVs are confirmed as binaries, different methods were used to derive their multiplicity, with different results. We report on a systematic search for multiplicity using spectroscopy, interferometry in a sample of 20 LBVs. Spectroscopy provides us with a bias-corrected binary fraction of $\[62_{ - 24}^{ + 38}\]$%, and a percentage of 50–70% is found from interferometry. This has a high impact on the way that these objects might be formed.

. We present UVIT/Astrosat UV photometry of the RSG population of the Small Cloud galaxy (SMC). As RSGs are extremely faint in the far-UV, these observations directly probe potential companion stars. From a sample of 861 SMC RSGs, we find 88 have detections at far-UV wavelengths: a clear signature of binarity. Stellar parameters are determined for both components, which allows us to study - for the first time - the mass-ratio (q) distribution of RSG binary systems. We find a flat mass-ratio distribution best describes the observations up to MRSG ∼15M⊙. We account for our main observing bias (i.e. the limiting magnitude of the UVIT survey) to determine the intrinsic RSG binary fraction of 18.8 ± 1.5 %, for mass-ratios in the range 0.3.<q<1.0 and orbital periods approximately in the range 3<log P[days]<8.

Since massive stars form preferentially as members of close binary systems, we use dense grids of detailed binary evolution models to explore how binary evolution shapes the main-sequence morphology of young star clusters. We propose that binary mergers might be the origin of the blue main sequence stars in young star clusters. Our results imply that stars may either form by accretion, or through a binary merger, and that both paths lead to distinctly different spins, magnetic fields, and stellar mass distributions.

We present results from a recent study of the spin rate properties of a sample of more than 400 Galactic O-type stars surveyed by the IACOB and OWN projects. We combine vsini, Teff, and logg estimates with information about the spectroscopic binarity status for 285 of the stars in the sample, and provide a renewed overview about how the empirical distribution of projected rotational velocities in the O-star domain depends on mass, evolutionary and binary status. The obtained distributions are then compared with predictions of state-of-the-art population synthesis simulations including binary interaction, and used to provide hints about the initial velocity distribution of stars with masses in the range 15-80 M⊙.