A global disk oscillation implemented in the viscous decretion disk (VDD) model has been used to reproduce most of the observed properties of the well known Be star $\zeta$ Tau. 48 Librae shares several similarities with $\zeta$ Tau -- they are both early-type Be stars, they display shell characteristics in their spectra, and they exhibit cyclic $V/R$ variations -- but has some marked differences as well, such as a much denser and more extended disk, a much longer $V/R$ cycle, and the absence of the so-called triple-peak features. We aim to reproduce the photometric, polarimetric, and spectroscopic observables of 48 Librae with a self-consistent model, and to test the global oscillation scenario for this target. Our calculations are carried out with the three-dimensional NLTE radiative transfer code HDUST. We employ a rotationally deformed, gravity-darkened central star, surrounded by a disk whose unperturbed state is given by the VDD model. A two-dimensional global oscillation code is then used to calculate the disk perturbation, and superimpose it on the unperturbed disk. A very good, self-consistent fit to the time-averaged properties of the disk is obtained with the VDD. The calculated perturbation has a period $P = 12$ yr, which agrees with the observed period, and the behaviour of the $V/R$ cycle is well reproduced by the perturbed model. The perturbed model improves the fit to the photometric data and reproduces some features of the observed spectroscopic data. Some suggestions to improve the synthesized spectroscopy in a future work are given.

Context. Various Be stars exhibit intensity variations of the violet and red emission peaks in their HI lines observed in emission. This so-called $V/R$ phenomenon is usually explained by the precession of a one-armed spiral density perturbation in the circumstellar disk. That global-disk oscillation scenario was confirmed, both observationally and theoretically, in the previous series of two papers analyzing the Be shell star {\zeta} Tauri. The vertically averaged (2D) global-disk oscillation model used at the time was able to reproduce the $V/R$ variations observed in H{\alpha}, as well as the spatially resolved interferometric data from AMBER/VLTI. Unfortunately, that model failed to reproduce the $V/R$ phase of Br15 and the amplitude of the polarization variation, suggesting that the inner disk structure predicted by the model was incorrect. Aims. The first aim of the present paper is to quantify the temporal variations of the shell-line characteristics of {\zeta} Tauri. The second aim is to better understand the physics underlying the $V/R$ phenomenon by modeling the shell-line variations together with the $V/R$ and polarimetric variations. The third aim is to test a new 2.5D disk oscillation model, which solves the set of equations that describe the 3D perturbed disk structure but keeps only the equatorial (i.e., 2D) component of the solution. This approximation was adopted to allow comparisons with the previous 2D model, and as a first step toward a future 3D model. Results. The new 2.5D formalism improves the agreement with the observed $V/R$ variations of H{\alpha} and Br15, under the proviso that a large value of the viscosity parameter, {\alpha} = 0.8, be adopted. Nonetheless, it remains challenging for the models to reproduce consistently the amplitude and the average level of the polarization data, whatever formalism is adopted.

We investigate the surface nitrogen content of the six magnetic O stars known to date as well as of the early B-type star tau Sco. We compare these abundances to predictions of evolutionary models to isolate the effects of magnetic field on the transport of elements in stellar interiors. We conduct a quantitative spectroscopic analysis of the sample stars with state-of-the-art atmosphere models. We rely on high signal-to-noise ratio, high resolution optical spectra obtained with ESPADONS at CFHT and NARVAL at TBL. Atmosphere models and synthetic spectra are computed with the code CMFGEN. Values of N/H together with their uncertainties are determined and compared to predictions of evolutionary models. We find that the magnetic stars can be divided into two groups: one with stars displaying no N enrichment (one object); and one with stars most likely showing extra N enrichment (5 objects). For one star (Theta1 Ori C) no robust conclusion can be drawn due to its young age. The star with no N enrichment is the one with the weakest magnetic field, possibly of dynamo origin. It might be a star having experienced strong magnetic braking under the condition of solid body rotation, but its rotational velocity is still relatively large. The five stars with high N content were probably slow rotators on the zero age main sequence, but they have surface N/H typical of normal O stars, indicating that the presence of a (probably fossil) magnetic field leads to extra enrichment. These stars may have a strong differential rotation inducing shear mixing. Our results should be viewed as a basis on which new theoretical simulations can rely to better understand the effect of magnetism on the evolution of massive stars.

We report the detection of a magnetic field on the Of?p star HD108. Spectropolarimetric observations conducted in 2007, 2008 and 2009 respectively with NARVAL@TBL and ESPaDOnS@CFHT reveal a clear Zeeman signature in the average Stokes V profile, stable on timescales of days to months and slowly increasing in amplitude on timescales of years. We speculate that this timescale is the same as that on which Ha emission is varying and is equal to the rotation period of the star. The corresponding longitudinal magnetic field, measured during each of the three seasons, increases slowly from 100 to 150G, implying that the polar strength of the putatively-dipolar large-scale magnetic field of HD108 is at least 0.5kG and most likely of the order of 1-2 kG. The stellar and wind properties are derived through a quantitative spectroscopic analysis with the code CMFGEN. The effective temperature is difficult to constrain because of the unusually strong HeI4471 and HeI5876 lines. Values in the range 33000-37000 K are preferred. A mass loss rate of about 1e-7 Msun/yr (with a clumping factor f=0.01) and a wind terminal velocity of 2000 km/s are derived. The wind confinement parameter eta_star is larger than 100, implying that the wind of HD108 is magnetically confined. Stochastic short-term variability is observed in the wind-sensitive lines but not in the photospheric lines, excluding the presence of pulsations. Material infall in the confined wind is the most likely origin for lines formed in the inner wind. Wind-clumping also probably causes part of the Ha variability. The projected rotational velocity of HD108 is lower than 50 km/s, consistent with the spectroscopic and photometric variation timescales of a few decades. Overall, HD108 is very similar to the magnetic O star HD191612 except for an even slower rotation.

We have investigated the stellar and wind properties of a sample of five late-type O dwarfs in order to address the weak wind problem. A grid of TLUSTY models was used to obtain the stellar parameters, and the wind parameters were determined by using the CMFGEN code. We found that the spectra have mainly a photospheric origin. A weak wind signature is seen in CIV 1549, from where mass-loss rates consistent with previous CMFGEN results regarding O8-9V stars were obtained. A discrepancy of roughly 2 orders of magnitude is found between these mass-loss rates and the values predicted by theory (Mdot(Vink)), confirming a breakdown or a steepening of the modified wind momentum-luminosity relation at log L/Lsun < 5.2. We have estimated the carbon abundance for the stars of our sample and concluded that its uncertainty cannot cause the weak wind problem. Upper limits on Mdot were established for all objects using lines of different ions, namely, PV 1118,28, CIII 1176, NV 1239,43, Si IV 1394,03, and NIV 1718. All the values obtained are also in disagreement with theoretical predictions, bringing support to the reality of weak winds. Together with CIV 1549, the use of NV 1239,43 results in the lowest mass-loss rates: the upper limits indicate that Mdot must be less than about -1.0 dex Mdot(Vink). Regarding the other transitions, the upper limits still point to low rates: Mdot must be less than about $(-0.5 \pm 0.2)$ dex Mdot(Vink). We have studied the behavior of the Halpha line with different mass-loss rates. We have also explored ways to fit the observed spectra with Mdot(Vink). By using large amounts of X-rays, we verified that few wind emissions take place, as in weak winds. However, unrealistic X-rays luminosities had to be used (log Lx/Lbol > -3.5) (abridged).

We report here the detection of a weak magnetic field of 50 - 100 G on the O9.7 supergiant zeta Ori A, using spectropolarimetric observations obtained with NARVAL at the 2m Telescope Bernard Lyot atop Pic du Midi (France). zeta Ori A is the third O star known to host a magnetic field (along with theta^1 Ori C and HD 191612), and the first detection on a 'normal' rapidly-rotating O star. The magnetic field of zeta Ori A is the weakest magnetic field ever detected on a massive star. The measured field is lower than the thermal equipartition limit (about 100 G). By fitting NLTE model atmospheres to our spectra, we determined that zeta Ori A is a 40 Msun star with a radius of 25 Rsun and an age of about 5 - 6 Myr, showing no surface nitrogen enhancement and losing mass at a rate of about 2x10^(-6) Msol/yr. The magnetic topology of zeta Ori A is apparently more complex than a dipole and involves two main magnetic polarities located on both sides of the same hemisphere; our data also suggest that zeta Ori A rotates in about 7.0 d and is about 40 degrees away from pole-on to an Earth-based observer. Despite its weakness, the detected magnetic field significantly affects the wind structure; the corresponding Alfven radius is however very close to the surface, thus generating a different rotational modulation in wind lines than that reported on the two other known magnetic O stars. The rapid rotation of zeta Ori A with respect to theta^1 Ori C appears as a surprise, both stars having similar unsigned magnetic fluxes (once rescaled to the same radius); it may suggest that the sub-equipartition field detected on zeta Ori A is not a fossil remnant (as opposed to that of theta^1 Ori C and HD 191612), but the result of an exotic dynamo action produced through MHD instabilities.